Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,432,985 A→G 5.9% intergenic (‑43/+27) folC ← / ← accD bifunctional folylpolyglutamate synthase/ dihydrofolate synthase/acetyl‑CoA carboxylase, beta (carboxyltransferase) subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,432,9850AG 5.9% 162.6 / 3.0 68intergenic (‑43/+27)folC/accDbifunctional folylpolyglutamate synthase/ dihydrofolate synthase/acetyl‑CoA carboxylase, beta (carboxyltransferase) subunit
Reads supporting (aligned to +/- strand):  ref base A (37/26);  new base G (2/2);  total (39/29)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.36e-01

GCGCTCAAGGCCGAGATCGATAGTTTTACTGTGCAGGTTTTCCAGATAAGAAAGCCACGAAGCCAGAGGCGACGCGGCTTGAGGAGTGCGTTTGATAATCATGGTATCCGCTGATTCGTTACGGTGAGAATAGCAAAAGGGCAGAGCCAGTGGCCCTGCCCTTATCAGTTATCAGGCCTCAGGTTCCTGATCCGGTACCGGGGGTACCACTACGCCTTCACGCGGCGCTTCAGGATTCGGCGCTGGCAGATTCATCAACTTCGCCAGAATGCTCGCCAGTTTCAGGCG  >  NC_000913/2432841‑2433128
                                                                                                                                                |                                                                                                                                               
gcgcTCAAGGCCGAGATCGATAGTTTTACTGTGCAGGTTTTCCAGATAAGAAAGCCACGAAGCAAGAGGCGCCGCGGCTTGAGGAGTGCGTTTGATAATCATGGTATCCGCTGATTCGTTACGGTGAGAATAGCAAAAGGGCAGAGCCa                                                                                                                                             >  1:1016560/1‑149 (MQ=255)
 cgcTCACGGCCGAGATCGATAGTTTTACTGTGCAGGTTTTCCAGATAAGAAAGCCACGAAGCCAGAGGCGACGCGGCTTGAGGAGTGCGTTTGATAATCATGGTATCCGCTGATTCGTTACGGTGAGAATAGCAAAAGGGCAGAGCCAg                                                                                                                                            <  1:895539/149‑1 (MQ=255)
    tCAAGGCCGAGATCGATAGTTTTACTGTGCAGGTTTTCCAGATAAGAAAGCCACGAAGCCAGAGGCGACGCGGCTTGAGGAGTGCGTTTGATAATCATGGTATCCGCTGATTCGTTACGGTGAGAATAGCAAAAGGGCAGAGCCAGTgg                                                                                                                                         >  1:790967/1‑149 (MQ=255)
       aGGCCGAGATCGATAGTTTTACTGTGCAGGTTTTCCAGATAAGAAAGCCACGAAGCCAGAGGCGACGCGGCTTGAGGAGTGCGTTTGATAATCATGGTATCCGCTGATTCGTTACGGTGAGAATAGCAAAAGGGCAGAGCCAGTGgccc                                                                                                                                      >  2:629195/1‑149 (MQ=255)
       aGGCCGAGATCGATAGTTTTACTGTGCAGGTTTTCCAGATAAGAAAGCCACGAAGCCAGAGGCGACGCGGCTTGAGGAGTGCGTTTGATAATCATGGTATCCGCTGATTCGGTACGGTGAGAATAGCAAAAGGGCAGAGCCAGTGgccc                                                                                                                                      >  1:833631/1‑149 (MQ=255)
        ggCCGAGATCGATAGTTTTACTGTGCAGGTTTTCCAGATAAGAAAGCCACGAAGCCAGAGGCGACGCGGCTTGAGGAGTGCGTTTGATAATCATGGTATCCGCTGATTCGTTACGGTGAGAATAGCAAAAGGGCAGAGCCAGTGgccct                                                                                                                                     >  2:489382/1‑149 (MQ=255)
         gCCGAGATCGATAGTTTTACTGTGCAGGTTTTCCAGATAAGAAAGCCACGAAGCCAGAGGCGACGCGGCTTGAGGAGTGCGTTTGATAATCATGGTATCCGCTGATTCGTTACGGTGAGAATAGCAAAAGGGCAGAGCCAGTGgccctg                                                                                                                                    <  2:22482/149‑1 (MQ=255)
          ccGAGATCGATAGTTTTACTGTGCAGGTTTTCCAGATAAGAAAGCCACGAAGCCAGAGGCGACGCGGCTTGAGGAGTTCGTTTGATAATCATGGTATCCGCTGATTCGTTACGGTGAGAATAGCAAAAGGGCAGGGCCACTGgctctgc                                                                                                                                   <  1:570508/149‑1 (MQ=255)
           cGAGATCGATAGTTTTAATGTGCAGGTTTTCCAGATAAGAAAGCCACGAAGCCAGAGGCGACGCGGCTTGAGGAGTGCGTTTGATAATCATGGTATCCGCTGATTCGTTACGGTGAGAATAGCAAAAGGGCAGAGCCAGTGgccctgcc                                                                                                                                  <  1:572904/149‑1 (MQ=255)
               aTCGATAGTTTTCCTGTGCAGGTGTTCCAGATAAGAAAGCCACTAAGCCAGAGGCGACGCGGCTTGAGGAGTGCGTTTGATAATCATGCAATCCGCTGATTCGTTACGGTGAGAATAGCAAAAGGGCAGAGCCAGTGGCCCTGCCCtta                                                                                                                              <  2:147434/149‑1 (MQ=255)
                tCGATAGTTTTACTGTGCAGGTTTTCCAGATAAGAACGCCACGAAGCCAGAGGCGACGCGGCTTGAGGAGTGCGTTTGATAATCATGGTATCCGCTGATTCGTTACGGTTAGAATAGCAAAAGGGCAGGGCCACTGGCTCTGCCCtttt                                                                                                                             <  1:541042/149‑3 (MQ=255)
                 cGATAGTTTTACTGTGCAGGTTTTCCAGATAAGAAAGCCACGAAGCCAGAGGCGACGCGGCTTGAGGAGTGCGTTTGATAATCATGGTATCCGCTGATTCGTTACGGAGAGAATAGCAAAAGGGCAGAGCCAGTGGCCCTGCCCttatc                                                                                                                            >  2:821297/1‑149 (MQ=255)
                  gATAGTTTTACTGTGCAGGTTTTCCAGATAAGAAAGCCACGAAGCCAGAGGCGACGCGGCTTGAGGAGTGCGTTTGATAATCATGGTATCCTCTGATTCGTTACGGTGAGAATAGCAAAAGGGCAGAGCCAGTGGCCCTGCCCttatca                                                                                                                           >  1:910076/1‑149 (MQ=255)
                  gATAGTTTTACTGTGCAGGTTTTCCAGATAAGAAAGCCACGAAGCCAGAGGCGACGCGGCTTGAGGAGTGCGTTTGATAATCATGGTATCCGCTGATTCGTTACGGTGAGAATAGCAAAAGGGCAGAGCCAGTGGCCCTGCCCttatca                                                                                                                           >  1:617059/1‑149 (MQ=255)
                   aTAGTTTTACTGTGCAGGTTTTCCAGATAAGAAAGCCACGAAGCCAGAGTCGACGCGGCTTGAGGAGTGCGTTTGATAATCATGGTATCCGCTGATTCGTTACGGTGAGAATAGCAAAAGGGCAGAGCCAGTGTCCCTGCCGttatcag                                                                                                                          >  1:691912/1‑149 (MQ=255)
                   aTAGTTTTACTGTGCAGGTTTTCCAGATAAGAAAGCCACGAAGCCAGAGGCGACGCGGCTTGAGGAGTGCGTTTGATAATCATGGTATCCGCTGATTCGTTACGGTGAGAATAGCAAAAGGGCAGAGCCAGTGGCCCTGCCCttatcat                                                                                                                          >  1:216120/1‑148 (MQ=255)
                   aTAGTTTTACTGTGCAGGTTTTCCAGATAAGAAAGCCACGAAGCCAGAGGCGACGCGGCTTGAGGAGTGCGTTTGATAATCATGGTATCCGCTGATTCGTTACGGTGAGAATAGCAAAAGGGCAGAGCCAGTGGCCCTGCCCttatcag                                                                                                                          >  1:393299/1‑149 (MQ=255)
                        tttACTGTGCAGGTTTTCCAGATAAGAAAGCCACGAAGCCAGAGGCGACGCGGCTTGAGGAGTGCGTTTGATAATCATGGTATCCGCTGATTCGTTACGGTGAGAATAGCAAAAGGGCAGAGCCAGTGGCCCTGCCCttatcagttatc                                                                                                                     <  2:617059/149‑1 (MQ=255)
                         ttAATGAGCAGGTTTTCCACTTAAGAAAGCCACCAAGGCAGAGGCGAAGCGGCTTGAGGAGTGCGTTTGATAATCATGGTATCCGCTGATTCGATACGGTGATAATAGCAAAAGGGCAGAGCCAGTGGCCCTGCCCttatcagttatca                                                                                                                    <  1:584424/149‑1 (MQ=255)
                             tgtgCATGTTTTCCAGATAAGAAAGCCACGAAGCCAGATGCGACGCGCCTAGAGGAGTGCGTTTGATAATCATGGTATCCGCTGATTATTTACGGTGAGAATAGCAAAAGGGCAGAGCCAGTGGCCCTGCCCTTATCAGTTATCAGGcc                                                                                                                >  2:254325/1‑149 (MQ=255)
                                     ttttCCAGATATGAAAGCCACGAAGCCAGAGGCGACGCGGCTTGAGGAGGGCGTTTGATAATCATGGTAGCCGCTGATTCGTTACGGTGAGAATAGCAAAAGCGCAGAGCCAGTGGCCCTGCCCTTATCAGTTATCAGGCCTCCGGTTc                                                                                                        <  1:691276/149‑1 (MQ=255)
                                       ttCCAGATAAGAAAGCCACGAAGCCAGAGGCGCCGCGGCTCGAGGAGTGCGTTTGATAATCATGGTATCCGCTGATTCGTTACGGTGAGAATAGCAAAAGGGCAGAGCCAGTGGCCCTGCCCTTATCAGTTATCAGGCCTCAGGTTCCt                                                                                                      >  1:590930/1‑149 (MQ=255)
                                            gATAAGAAAGCCACGAAGCCAGAGGCGACGCGGCTTGAGGAGTGCGTTTGATAATCATGGTATCCGCTGATTCGTTACGGTGAGAATAGCAAAAGGGCAGAGCCAGTGTCCCTGCCCTTATCAGTTATCAGGCCTCAGGTTCCTGATcc                                                                                                 >  2:185846/1‑149 (MQ=255)
                                            gATAAGAAAGCCACGAAGCCAGAGGCGACGCGGCTTGAGGAGTGCGTTTGATAATCATGGTATCCGCTGATTCGTTACGGTGAGAATAGCAAAAGGGCAGAGCCAGTGGCCCTGCCCTTATCAGTTATCAGGCCTCAGGTTCCTGATcc                                                                                                 >  2:182490/1‑149 (MQ=255)
                                               aaGAAAGCCACGAAGCCAGAGGCGACGCGGCTTGAGGAGTGCGTTTGATAATCATGGTATCCGCTGATTCGTTACGGTGAGAATAGCAAAAGGGCAGAGCCAGTGGCCCTGCCCTTATCAGTTATCATGCCTCAGGTTCCTGATCCGGt                                                                                              >  2:298742/1‑149 (MQ=255)
                                                 gAAAGCCACGAAGCCAGAGGCGACGCGGCTTGAGGAGTGCGTTTGATAATCCTGGTATCCGCTGATTCGTTACGGTGAGAATAGCAAAAGGGCAGAGCCAGTGGCCCTGCCCAGATCAGATAGCCGGCCTCAGGTTCCTGATCCGGTAc                                                                                            >  1:557900/1‑149 (MQ=255)
                                                 gAAAGCCACGAAGCCAGAGGCGACGCGGCTTGAGGAGTGCGTTTGATAATCATGGTATCCGCTGATTCGTTACGGTGAGAATAGCAAAAGGGCAGAGCCAGTGGCCCTGCCCTTATCAGTTATCAGGCCTCAGGTTCCTGATCCGGTAc                                                                                            >  1:109988/1‑149 (MQ=255)
                                                      ccACGACGCCAGAGGCGACGCGGCTTGAGGAGTGCGTTTGATAATCATGGTATCCGCTGATTCGTTACGGTGAGAATAGCAAAAGGGCAGAGCCAGTGGCCCTGCCCTTATCAGTTATCAGGCCTCAGGTTCCTGATCCGGTACCgggg                                                                                       >  1:135367/1‑149 (MQ=255)
                                                        ccGAAGCCAGAGGCGACGCGGCTTGAGGAGTGCGTTTGATAATCATGGTATCCGCTGATTCGTTCCGGAGAGAATAGCAACAGGGCAGAGCCAGTGGCCCTGCCCTTATCAGTTATCAGGCCTCAGGTTCCTGATCCGGTACCGGGGGt                                                                                     <  1:568333/148‑1 (MQ=255)
                                                          gAAGCCAGAGGCGACGCGGCTTGAGGAGTGCGTTTGATAATCATGGTATCCGCTGATTCGTTACGGTGAGAATAGCAAAAGGGCAGAGCCAGTGGCCCTGCCCTTATCAGTTATCAGGCCTCAGGTTCCTGATCCGGTACCGGGGGTAc                                                                                   <  1:555181/149‑1 (MQ=255)
                                                            aGCCAGAGGCGACGCGGCGTGAGGAGTGCGTTAGATAATCATGGTATCCGCTGATTCCTTACGGTGAGAATAGCAAAAGGGCAGCGCCAGTGGCCCTGCCCTTATCAGTTATCAGGCCTCAGGATCCTGATCCGGTACCGGGGGTACCa                                                                                 <  2:286359/149‑1 (MQ=255)
                                                             gCCAGAGGCGACGCGGCTTGAGGAGTGCGTTTGATAATCATGGTATCCGCTGATTCGTTACGGTGAGAATAGCAAAAGGGCAGAGCCAGTGGCCCTGCCCTTATCAGTTATCAGGCCTCAGGTTCCTGATCCGGTACCGGGGGTACCAc                                                                                >  1:460774/1‑149 (MQ=255)
                                                                   ggCGACGCGGCTTGAGGAGTGCGTTTGATAATCATGGTATCCGCTGATTCGTTACGGTGAGAATAGCAAAAGGGCAGAGCCAGTGGCCCTGCCCTTATCAGTTATCAGGCCTCAGGTTCCTGATCCGGTACCGGGGGTACCACTACGcc                                                                          >  2:123514/1‑149 (MQ=255)
                                                                    gCGACGCGGCTTGAGTAGTGCGTTTGATAATCATGGTATCCGAAGATTAGTTACGGTGAGAATAGCAAAAGGGCAGAGACAGTGGCCCTGACCTTAACAGTTATCAGGCCTCATGTTCCTGATCCGGTACCGGGGGTACCACTACGCCt                                                                         >  1:110515/1‑149 (MQ=255)
                                                                    gCGACGCGGCTTGAGGAGTGCGTTTGATAATCATGGTATCCGCTGATTCGTTACGGTGAGAATAGCAAAAGGGCAGAGCCAGTGGCCCTGCCCTTATCAGTTATCAGGCCTCAGGTTCCTGATCCGGTACCGCGGGTACCACTACGCCt                                                                         >  2:782827/1‑149 (MQ=255)
                                                                    ccGACGCGGCTTGAGGAGAGCGTTTGATAATCAAGGTATCCGCTGATTCGTTACGGAGAGAATAGCAAAAGGGCAGAGCCAGTGGCCCTGCCCTTATCAGTTATCAGGCCTCAGGTTCCTGATCCGGTACCGGGGGTACCACTACGCCt                                                                         <  1:432171/148‑1 (MQ=255)
                                                                                 aggagTGCGTTTGATAATCATGGTTTCAGCTGATTCGTTACGGTGAGAAAAGCAAAAGGGCAGAGCCAGTGGCCCTGCCCTTATCAGTTATCAGGCCTCAGGTTCCTGATCCGGTACCGGGGGTACCACTACGCCTTCACGCGGCGCtt                                                            <  1:900163/149‑1 (MQ=255)
                                                                                   gagTGCGTTTGATAATCATGGTATCCGCTGATTCGTTACGGTGAGAATAGCAAAAGGGCAGAGCCAGTGGCCCTGCCCTTATCAGTTATCAGGCCTCAGGTTCCTGATCCGGTACCGGGGGTACCACTACGCCTTCACGCGGCGCTTCa                                                          >  2:596889/1‑149 (MQ=255)
                                                                                     gTGCGTTTGATAATCATGGTATCCGCTGATTCGTTACGGTGAGAATAGCAAAAGGGCAGAGCCAGGGGCCCTGCCCTTATCAGTTATCAGGCCTCAGGTTCCTGATCCGGTACCGGGGGTACCACTACGCCTTCACGCGGCGCTTCAgg                                                        <  1:424395/149‑1 (MQ=255)
                                                                                      tGCGTTTGATAATCATGGTATCCGCTGATTCGTTACGGTGAGAATAGCAAAAGGGCAGAGCCAGTGGCCCTGCCCTTATCAGTTATCAGGCCTCAGGTTCCTGATCCGGTACCGGGGGTACCACTACGCCTTCACGCGGCGCTTCAGGa                                                       >  2:382731/1‑149 (MQ=255)
                                                                                                aaTCCTGGTATCCGCTGATTCGTTACGGGGAGAATAGCAAACGGGCCGAGCCAGTGGCCCTGCCCTTATCAGTTATCAGGCCTCAGGTGCCTGATCCGGTACCGGGGGTACCACTACGCCTTCACGCGGCGCTTCAGGATTCGGCGCTg                                             <  1:1002757/149‑1 (MQ=255)
                                                                                                 aTCATGGTATCCGCTGATTCGTTACGGTGCGAATAGCAAAAGGGCAGAGCCAGTGGCCCTGCCCTTATCAGTTATCAGGCCTCCGGTTCCTGTTCCGGTACCGGGGGTACCACTACGCCTTCACGCGGCGCTTCAGGATTCGGCGCTgg                                            <  1:489382/149‑1 (MQ=255)
                                                                                                  tCATGGGAGCCGCTGATTCGTTACGGTGAGAATAGCAAAAGGGCAGAGCCAGTGGCCCTGCCCTTATCAGTTATCAGGCCTCAGGTTCCTGATCCGGTACCGGGGGTACCACTACGCCTTCACGCGGCGCTTCAGGATTCGGCGCTGGc                                           <  1:852260/149‑1 (MQ=255)
                                                                                                   cATGGTATCTGCTGATTCGTTACGGTGAGAATAGCAAAAGGGCAGAGCCAGTGGCCCTGCCCTTATCAGTTATCAGGCCTCAGGTTCCTGATCCGGTACCGGGGGTACCACTACGCCTTCACGCGGCGCTTCAGGATTCGGCGCTGGCa                                          <  2:370847/149‑1 (MQ=255)
                                                                                                    aTGGTATCCGCTGATTCGTTACGGTGAGAATAGCAAAAGGGCAGAGCCAGTGGCCCTGCCCTTATCAGTTATCAGGCCTCAGGTTCCTGATCCGGTACCGGGGGTACCACTACGCCTTCACGCGGCGCTTCAGGATTCGGCGCTGGCAg                                         >  2:771941/1‑149 (MQ=255)
                                                                                                     tGGTATCCGCTGATTCGTTACGGTGAGAATAGCAAAAGGGCAGAGCCAGTGGCCCTGCCCTTATCAGTTATCAGGCCTCAGGTTCCTGATCCGGTACCGGGGGTACCACTACGCCTTCACGCGGCGCTTCAGGATTCGGCGCTGGCAGa                                        >  2:801921/1‑149 (MQ=255)
                                                                                                      ggTATCCGCTGATTCGTTACGGTGAGAATAGCAAAAGGGCAGAGCCAGTGGCCCTGCCCTTATCAGTTATCAGGCCGCAGGTTCCTGATCCGGGACCGGGGGTACCACTACGCCTTCACGCGGCGCTTCAGGATTCGGCGCTGGCAGAt                                       <  2:1019138/149‑1 (MQ=255)
                                                                                                       gTATCCGCTGATTCGTTACGGTGAGAATAGCAAAAGGGCAGAGCCAGTGGCCCTGCCCTTATCAGTTATCAGGCCTCAGGTTCCTGATCCGGTACCGGGGGTACCACTACGCCTTCACGCGGCGCTTCAGGATTCGGCGCTGGCAGAtt                                      >  2:915385/1‑149 (MQ=255)
                                                                                                            cGCTGATTCGTTACGGTGAGAATAGCAAAAGGGCAGAGCCAGTGGCCCTGCCCTTATCAGTTATCAGGCCTCAGGTTCCTGATCCGGTACCGGGGGTACCACTACGCCTTCACGCGGCGCTTCAGGATTCGGCGCTGGCAGATtcatca                                 >  1:961695/1‑149 (MQ=255)
                                                                                                                     gTTACGGTGAGAATAGCAAAAGGGCAGAGCCAGTGGCCCTGCCCTTATCAGTTATCAGGCCTCAGGTTCCTGATCCGGTACCGGGGGTACCACTACGCCTTCACGCGGCGCTTCAGGATTCGGCGCTGGCAGATTCATCAACTACGCCa                        >  2:69825/1‑149 (MQ=255)
                                                                                                                         cGGTGAGAATAGCAAAAGGGCAGAGCCAGTGGCCCTGCCCTTATCAGTTATCAGGCCTCAGGTTCCTGCTCCGGTACCGGGGGTACCACTACGCCTTCACGCCGCGCTTCAGGATTCGGCGCTGGCAGATTCATCAACTTCGCCAGaaa                    >  1:936379/1‑148 (MQ=255)
                                                                                                                              agaATAGCAAAAGGGCAGAGCCAGTGGCCCTGCCCTTATCAGTTATCAGGCCTCAGGTTCCTGATCCGGTACCGGGGGTTCAACTACGCCTTCACGCGGCGCTTCAGGATTCGGCGCTGGCAGCTTCATCCACTTCGCCAGCATGCTCg               >  2:922476/1‑149 (MQ=255)
                                                                                                                              agaATAGCAAAAGGGCAGAGCCAGTGGCCCTGCCCTTATCAGTTATCAGGCCTCAGGTTCCTGATCCGGTACCGGGGGTACCACTACGCCTTCACGCGGCGCTTCAGGATTCGGCGCTGGCAGATTCATCAACTTCGCCAGAATGCTCg               >  1:293170/1‑149 (MQ=255)
                                                                                                                               taaTAGCAAAAGGGCAGAGCCAGTGGCCGTGCACTTATCAGTTATCAGGCCTCAGGGTCCTGATCCGGTACCGGGGGTACCACTACGCCTTCACGCGGCGCTTCAGGATTCGGCGCTGGCAGATTCATCAACTTCGCCAGAATGCTCGc              <  2:990193/148‑1 (MQ=255)
                                                                                                                                aaTAGCAAAAGGGCAGAGCCAGTGGCCCTGCCCTTATCAGTTATCAGGCCTCAGGTTCCTGATCCGGTACCGGGGGTACCACTACGCCTTCACGCGGCGCTTCAGGATTCGGCGCTGGCAGATTCATCAACTTCGCCAGAATGCTCGcc             >  1:303275/1‑149 (MQ=255)
                                                                                                                                  tAGCAAAAGGGCAGAGCCAGTGGCCCTGCCCTTATCAGTTATCAGGCCTCAGGTTCCTGATCCGGTACCGGGGGTACCACTACGCCTTCACGCGGCGCTTCAGGATTCGGCGCTGGCAGATTCATCAACTTCGCCAGAATGCTCGCCAg           <  2:109988/149‑1 (MQ=255)
                                                                                                                                    gCAAAAGGGCAGAGCCAGTGGCCCTGCCCTTATCAGTTATCAGGCCTCAGGTTCCTGATCCGGTACCGGGGGTACCACTACGCCTTCACGCGGGGCTTCAGGATTCGGCGCTGGCCGCTTCATCAACTTCGCCAGACTGCTCGCCAGtt         >  1:123988/1‑149 (MQ=255)
                                                                                                                                      aaaaGGTCAGATCAATTGGCCCTGCCCTTATCAGTTATCAGGCCTCAGTTTCCTGATCCGGTACCGGGGGTACCACTACGCCTTCACTCTGCGATTCAGGATTCGTCGCTGGCATATTCATCAACTTCGCCAGAATGCTCGCCAGTTTc       <  2:188097/149‑1 (MQ=255)
                                                                                                                                      aaaaGGGCAGAGCCAGTGGCCCTGCCCTTATCAGTTATCAGGCCTCAGGTTCCTGATCCGGTACCGGGGGTACCACTACGCCTTCACGCGGCGCTTCAGGATTCGGCGCTGGCAGATTCATCAACTTCGCCAGAATGCTCGCCAGTTTc       <  1:407877/149‑1 (MQ=255)
                                                                                                                                      aaaaGGGCAGAGCCAGTGGCCCTGCCCTTATCAGTTATCAGGCCTCAGGTTCCTGATCCGGTACCGGGGGTACCACTACGCCTTCACGCGGCGCTTCAGGATTCGGCGCTGGCAGATTCATCAACTTCGCCAGAATGCTCGCCAGTTTc       <  1:483373/149‑1 (MQ=255)
                                                                                                                                        aaGGGCTGAGCCAGTGGCCCTGCCATTATCAGTTATCCGGCCTCAGGTTCCTGATCCGGTACCGGGGGTACCACTACGCCTTCACGCGGCGCTTCAGGATTCGGCGCTGGCAGATTCATCAACTTCGCCAGAATGCTCGCCAGTTTCAg     <  2:961695/149‑1 (MQ=255)
                                                                                                                                        aaGGGCAGAGCCAGTGGCCCTGCCCTTATCAGTTATCAGGCCTCAGGTTCCTGATCCGGTACCGGGGGTACCACTACGCCTTCACGCGGCGCTTCAGGATTCGGCGCTGGCAGATTCATCAACTTCGCCAGAATGCTCGCCAGTTTCAg     <  1:262182/149‑1 (MQ=255)
                                                                                                                                        aaGGGCAGAGCCAGTGGCCCTGCCCTTATCAGTTATCAGGCCTCAGGTTCCTGATCCGGTACCGGGGGTACCACTACGCCTTCACGCGGCGCTTCAGGATTCGGCGCTGGCAGATTCATCAACTTCGCCAGAATGCTCGCCAGTTTCAg     <  1:802171/149‑1 (MQ=255)
                                                                                                                                         agggcagggcCACTGGCTCTGCCCTTATCAGTTATCAGGCCTCAGGTTGCTGAGCCGGTACCGGGGGTTCCACTACGCCTTCAAGCGGCGCGGAAGGATTCGGCGCTGGCAGATTCATCAACATCGCCCGAATGATCACAAGATTCAgg    >  1:422588/1‑149 (MQ=255)
                                                                                                                                          gggCAGGGCCACTGGCTCTGCCCTGATAAGTTATCAGGCCTAAGGTTCCTGATCCGGTACCGGGGGTACCACAAAGCCTTCACGCGCCGCTTCAGCATTCGGCGCTGGCTGCTTCATCCACTTCGCCCGAAAGCTCGCCAGTTTCCGgc   >  2:127326/1‑149 (MQ=255)
                                                                                                                                          gggCAGAGCCAGTGCCCTGCCCTTATCAGTTATCAGGCCTCAGGTTCCTGATCCGGTACCGGGGGTACCACTACGCCTTCACGCGGCGCTTCAGGATTCGGCGCTGGCAGATTCATCAACTTCGCCAGAATGCTCGCCAGTTTCAGgcg  >  1:316470/1‑149 (MQ=255)
                                                                                                                                           ggCAGGGCCACTGGCTCTGCCCTTTTCAGTTATCAGGCCTCAGGTTCCTGATCCGGTACCGGGGGTACCACTACGCCTTCACGCGGCGCTTCAGGATTCGGCGCTGGCAGATTCATCAACTTCGCCAGAATGCTCGCCAGTTTCAGgcg  >  1:545123/1‑149 (MQ=255)
                                                                                                                                            gCAGAGCCAGTGGCCCTGCCCTTATCAGTTATCAGGCCTCAGGTTCCTGATCCGGTACCGGGGGTACCACTACGCCTTCACGCGGCGCTTCAGGATTCGGCGCTGGCAGATTCATCAACTTCGCCAGAATGCTCGCCAGTTTCAg     >  2:250850/1‑145 (MQ=255)
                                                                                                                                            gCAGAGCCAGTGGCCCTGCCCTTATCAGTTATCAGGCCTCAGGTTCCTGATCCGGTACCGGGGGTACCACTACGCCTTCACGCGGCGCTTCAGGATTCGGCGCTGGCAGATTCATCAACTTCGCCAGAATGCTCGCCAGTTTCAg     <  1:250850/145‑1 (MQ=255)
                                                                                                                                                |                                                                                                                                               
GCGCTCAAGGCCGAGATCGATAGTTTTACTGTGCAGGTTTTCCAGATAAGAAAGCCACGAAGCCAGAGGCGACGCGGCTTGAGGAGTGCGTTTGATAATCATGGTATCCGCTGATTCGTTACGGTGAGAATAGCAAAAGGGCAGAGCCAGTGGCCCTGCCCTTATCAGTTATCAGGCCTCAGGTTCCTGATCCGGTACCGGGGGTACCACTACGCCTTCACGCGGCGCTTCAGGATTCGGCGCTGGCAGATTCATCAACTTCGCCAGAATGCTCGCCAGTTTCAGGCG  >  NC_000913/2432841‑2433128

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 15 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: