Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,932,094 C→G 5.9% intergenic (+35/+21) purT → / ← eda phosphoribosylglycinamide formyltransferase 2/KHG/KDPG aldolase; 2‑dehydro‑3‑deoxy‑phosphogluconate/4‑hydroxy‑2‑ oxoglutarate aldolase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,932,0940CG 5.9% 345.0 / 10.4 116intergenic (+35/+21)purT/edaphosphoribosylglycinamide formyltransferase 2/KHG/KDPG aldolase; 2‑dehydro‑3‑deoxy‑phosphogluconate/4‑hydroxy‑2‑ oxoglutarate aldolase
Reads supporting (aligned to +/- strand):  ref base C (63/46);  new base G (5/2);  total (68/48)
Fisher's exact test for biased strand distribution p-value = 6.98e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.12e-01

Only 100 of 121 total aligned reads displayed.

GCCGGAAATTGATGGCAGCCGTCGTCTGGGGGTGGCACTGGCTACTGCAGAGAGTGTTGTTGACGCCATTGAACGCGCGAAGCACGCCGCCGGACAGGTAAAAGTACAGGGTTAAACCCGGCAAAAAAACGCTACAAAAATGCCCGATCCTCGATCGGGCATTTTGACTTTTACAGCTTAGCGCCTTCTACAGCTTCACGCGCCAGCTTAGTAATGCGGTCGTAATCGCCCGCTTCCAGCGCATCTGCCGGAACCAGCCAGGAACCACCGATGCACAGCACGCTTTTCAGCGCCAGGTA  >  NC_000913/1931945‑1932243
                                                                                                                                                     |                                                                                                                                                     
gCCGGAAATTGATGGCAGCCGTCGTCTGGGGGTGGCACTGGCTACTGCAGAGAGTGTTGTTGACGCCATTGAACGCGCGAAGCACGCCGCCGGACAGGTAAAAGTACAGGGTTAAACCCGGCAAAAAAACGCTACAAAAATGCCCGATcc                                                                                                                                                       >  2:1867955/1‑150 (MQ=255)
  cGGAAATTGATGGCAGCCGTCGTCTGGGGGTGGCACTGGCTACTGCAGAGAGTGTTGTTGACGCCATTGAACGCGCGAAGCACGCCGCCGGACAGGTAAAAGTACAGGGTTAAACCCGGCAAAAAAACGCTACAAAAATGCCCGATCCTc                                                                                                                                                     >  1:1204188/1‑150 (MQ=255)
  cGGAAATTGATGGCAGCCGTCGTCTGGGGGTGGCACTGGCTACTGCAGAGAGTGTTGTTGACGCCATTGAACGCGCGAAGCACGCCGCCGGACAGGTAAAAGTACAGGGTTAAACCCGGCAAAAAAACGCTACAAAAATGCCCGATCCTc                                                                                                                                                     >  1:382957/1‑150 (MQ=255)
           aTGGCAGCCGTCGTCTGGGGGTGGCACTGGCTACTGCAGAGAGTGTTGTTGACGCCATTGAACGCGCGAAGCACGCCGCCGGACAGGTAAAAGTACAGGGTTAAACCCGGCAAAAAAACGCTACAAAAATGCCCGATCCTCGATCGGGCa                                                                                                                                            <  2:12664/150‑1 (MQ=255)
           aTGGCAGCCGTCGTCTGGGGGTGGCACTGGCTACTGCAGAGAGTGTTGTTGACGCCATTGAACGCGCGAAGCACGCCGCCGGACAGGTAAAAGTACAGGGTTAAACCCGGCAAAAAAACGCTACAAAAATGCCCGATCCTCGATCGGGCa                                                                                                                                            >  1:897030/1‑150 (MQ=255)
           aTGGCAGCCGTCGTCTGGGGGTGGCACAGGGTACTGCAGAGAGAGTTGTTGACTCCATTGAACGCGCGAAGGACGCCTCCGGACAGGTTAAAGTACAGGGTAAAACCCGGCAAAAAACCGCTACAAAAATGCCAGATCCTCGATCGGGCa                                                                                                                                            >  1:114389/1‑150 (MQ=255)
             ggCAGCCGGCGTCTGGGGGTGGCACTGGCTACTGCAGAGAGTGTTGTTGACGCCATTGAACGCGCGAAGCACGCCGCCGGACAGGTAAAAGTACAGGGTTAAACCCGGCAAAAAAACGCTACAAAAATGCCCGATCCTCGATCGGGCAtt                                                                                                                                          <  2:906484/150‑1 (MQ=255)
              gCAGCCGTCGTCTGGGGGTGGCACTGGCTACTGCAGAGAGTGTTGTTGACGCCATTGAACGCGCGAAGCACGCCGCCGGACAGGTAAAAGTACAGGGTTAAACCCGGCAAAAAAACGCTACAAAAATGCCCGATCCTCGATCGGGCAttt                                                                                                                                         >  2:1675612/1‑150 (MQ=255)
                 gCCGTCGTCTGGGGGTGGCACTGGCTACTGCAGAGAGTGTTGTTGACGCCATTGAACGCGCGAAGCACGCCGCCGGACAGGTAAAAGTACAGGGTTAAACCCGGCAAAAAAACGCTACAAAAATGCCCGATCCTCGATCGGGCATTTTGa                                                                                                                                      <  1:729163/150‑1 (MQ=255)
                   cgtcgtCTGGGGGTGGCACTGGCTACTGCAGAGAGTGTTGTTGACGCCATTGAACGCGCGAAGCACGCCGCCGGACAGGTAAAAGTACAGGGTTAAACCCGGCAAAAAAACGCTACAAAAATGCCCGATCCTCGATCGGGCATTGTGACt                                                                                                                                    >  1:842507/1‑150 (MQ=255)
                      cgtcTGGGGGTGGCACTGGCTACTGCAGAGAGTGTTGTTGACGCCATTGAACGCGCGAAGCACGCCGCCGGACAGGTAAAAGTACAGGGTTAAACCCGGCAAAAAAACGCTACAAAAATGCCCGATCCTCGATCGGGCATTTTGACtttt                                                                                                                                 <  1:79545/150‑1 (MQ=255)
                       gtcTGGGGGTGGCACTGGCTACTGCAGAGAGTGTTGTTGACGCCATTGAACGCGCGAAGCACGCCGCCGGACAGGTAAAAGTACAGGGTTAAACCCGGCAAAAAAACGCTACAAAAATGCCCGATCCTCGATCGGGCATTTTGACTTTTa                                                                                                                                >  1:1031994/1‑150 (MQ=255)
                        tcTGGGGGTGGCACTGGCTACTGCAGAGAGTGTTGTTGACGCCATTGAACGCGCGAAGCACGCCGCCGGACAGGTAAAAGTACAGGGTTAAACCCGGCAAAAAAACGCTACAAAAATGCCCGATCCTCGATCGGGCATTTTGACTATTAc                                                                                                                               >  1:96365/1‑150 (MQ=255)
                                 ggCACTGGCTACTGCAGAGAGTGTTGTTGACGCCATTGAACGCGCGAAGCACGCCGCCGGACAGGTAAAAGTACAGGGTTAAACCCGGCAAAAAAACGCTACAAAAATGCCCGCTCCTCGATCGGGCATTTTGACTTTTACAGCTTAgcg                                                                                                                      <  1:338160/150‑1 (MQ=255)
                                 ggCACTGGCTACTGCAGAGAGTGTTGTTGACGCCATTGAACGCGCGAAGCACGCCGCCGGACAGGTAAAAGTACAGGGTTAAACCCGGCAAAAAAACGCTACAAAAATGCCCGATCCTCGATCGGGCATTTTGACTTTTACAGCTTAgcg                                                                                                                      <  2:425494/150‑1 (MQ=255)
                                   cACTGGCTACTGCAGAGAGTGTTGTTGACGCCATTGAACGCGCGAAGCACGCCGCCGGACAGGTAAAAGTACAGGGTTAACCCCGGCAAAAAAACGCTACAAAAATGCCCGATCCTCGCTCGGGCATTTTGACTTTTACAGCTTAGCGcc                                                                                                                    >  1:36984/1‑150 (MQ=255)
                                    aCTGGCTACTGCAGAGAGTGTTGTTGACGCCATTGAACGCGCGAAGCACGCCGCCGGACAGGTAAAAGTACAGGGTTAAACCCGGCAAAAAAACGCTACAAAAATGCCCGATCCTCGATCGGGCATTTTGACTTTTACAGCTTAGCGCCt                                                                                                                   <  2:1836636/150‑1 (MQ=255)
                                         cTACTGCAGAGAGTGTTGTTGACGCCATTGAACGCGCGAAGCACGCCGCCGGACAGGTAAAAGTACAGGGTTAAACCCGGCAAAAAAACGCTACAAAAATGCCCGATCCTCGATCGGGCATTTTGACTTTTACAGCTTAGCGCCTTCTAc                                                                                                              >  1:825284/1‑150 (MQ=255)
                                          tACTGCAGAGAGTGTTGTTGACGCCATTGAACGCGCGAAGCACGCCGCCGGACAGGTAAAAGTACAGGGTTAAACCCGGCAAAAAAACGCTACAAAAATGCCCGATCCTCGATCGGGCATTTTGACTTTTACAGCTTAGCGCCTTCTACa                                                                                                             <  2:865332/150‑1 (MQ=255)
                                          tACTGCAGAGAGTGTTGTTGACGCCATTGAACGCGCGAAGCACGCCGCCGGACAGGTAAAAGTACAGGGTTAAACCCGGCAAAAAAACGCTACAAAAATGCCCGATCCTCGATCGGGCATTTTGACTTTTACAGCTTAGCGCCTTCTACa                                                                                                             <  1:508538/150‑1 (MQ=255)
                                          tACTGCAGAGAGTGTTGTTGACGCCATTGAACGCGCGAAGCACGACGTCGGACAGGTAAAAGTACAGGGTTAAACCCGGCAAAAAAACGCTACAAAAATGCCCGATCCTCGATCGGGCATTTTGACTTTTACAGCTTAGCGCCTTCTACa                                                                                                             <  2:611981/150‑1 (MQ=255)
                                             tGCAGAGAGTGTTGTTGACGCCATTGAACGCGCGAAGCACGCCGCCGGACAGGTAAAAGTACAGGGTTAAACCCGGCAAAAAAACGCTACAAAAATGCCCGCTCCTCGATCGGGCATTTTGACTTTTACAGCTTAGCGCCTTCTACAGCt                                                                                                          >  2:121864/1‑150 (MQ=255)
                                              gaagagagTGTGGTTGACGCCATTGAACGCGCGAAGCACGCCGCCGGACAGGTAAAAGTACAGGGTTAAACCCGGCAAAAAAACGCTACAAAAATGCCCGATCCTCGATCGGGCATTTTGACTTTTACAGCTTAGCGCCTTCTACAGCtt                                                                                                         <  1:1243841/148‑1 (MQ=255)
                                              gCAGAGAGTGTTGTTGACGCCATTGAACGCGCGAAGCACGCCGCCGGACAGGTAAAAGTACAGGGTTAAACCCGGCAAAAAAACGCTACAAAAATGCCCGATCCTCGATCGGGCATTTTGCCTTTTACAGCTTAGCGCCTGCTACAGCtg                                                                                                         >  2:1262027/1‑149 (MQ=255)
                                              gCAGAGAGTGTTGTTGACGCCATTGAACGCGCGAAGCACGCCGCCGGACAGGTAAAAGTACAGGGTTAAACCCGGCAAAAAAACGCTACAAAAATGCCCGATCCTCGATCGGGCATTTTGACTTTTACAGCTTAGCGCCTTCTACAGCtt                                                                                                         >  1:1008176/1‑150 (MQ=255)
                                               cAGAGAGTGTTGTTGCCGCCATTGAACGCGCGAAGCACGCCGCCGGACAGGTAAAAGTACAGGGTTAAACCCGGCAAAAAAACGCTACAAAAATGCCCGATCCTCGATCGGGGATTTTGACTTTTACAGCTTAGCGCCTTCTACAGCTTc                                                                                                        >  2:327488/1‑150 (MQ=255)
                                                 gagagTGTTGTTGACGCCATTGAACGCGCGAAGCACGCCGCCGGACAGGTAAAAGTACAGGGTTAAACCCGGCAAAAAAACGCTACAAAAATGCCCGATCCTCGATCGGGCATTTTGACTTTTACAGCTTAGCGCCTTCTACAGCTTCAc                                                                                                      >  1:1310187/1‑150 (MQ=255)
                                                           ttgACGCCATTGAACGCGCGAAGCACGCCGCCGGACAGGTAAAAGTACAGGGTTAAACCCGGCAAAAAAACGCTACAAAAATGCCCGATCCTCGATCGGGCATTTTGACTTTTACAGCTTAGCGCCTTCTACAGCTTCACGCGCCAGCtt                                                                                            <  2:1204188/150‑1 (MQ=255)
                                                            tgACGCCATTGAACGCGCGAAGCACGCCGCCGGACAGGTAACAGTACAGGGTTAAACCCGGCAAAAAAACGCTACAAAAATGCCCGATCCTCGAGCGGGCATTTTGACTTTTACAGCTTAGCGCCTTCTACAGCTTCACGAGCCAGCTTa                                                                                           >  2:1190639/1‑150 (MQ=255)
                                                                  cATTGAACGCGCGAAGCACGCCGCCGGACAGGTAAAAGTACAGGGTTAAACCCGGCAAAAAAACGCTACAAAAATGCCCGATCCTCGATCGGGCATTTTGACTTTTACAGCTTAGCGCCTTCTACAGCTTCACGCGCCAGCTTAGTAATg                                                                                     >  1:1956070/1‑150 (MQ=255)
                                                                    ttGAACGCGCGAAGCACGCCGCCGGACAGGTAAAAGTACAGGGTTAAACCCGGCAAAAAAACGCTACAAAAATGCCCGATCCTCGTTCGGGCATTTTGACTTTTACAGCTTAGCGCCTTCTACAGCTTCACGCGCCAGCTTAGTAATGCg                                                                                   <  2:1163464/150‑1 (MQ=255)
                                                                        ccgcgcgCAGCACGCCGCCGGACAGGTAAAAGGACAGGGGTAAAGCCGGCAAAAAAACGCTACAAAAATGCCCGAGCCTCGCGCGGGCATTTTGACTTTTACCGCTTAGCGCCTTCTACAGCTTCACGCGCCAGCGTAGTAATGCGGTCg                                                                               <  2:976366/149‑1 (MQ=255)
                                                                        aCGCGCGAAGCACGCCGCCGGACAGGTAAAAGTACAGGGTTAAACCCGGCAAAAAAACGCTACAAAAATGCCCGATCCTCGATCGGGCATTTTGACTTTTACAGCTTAGCGCCTTCTACAGCTTCACGCGCCAGCTTAGTAATGCGGTCg                                                                               <  2:1381087/150‑1 (MQ=255)
                                                                        aCGCGCGAAGCACGCCGCCGGACAGGTAAAAGTACAGGGTTAAACCCGGCAAAAAAACGCTACAAAAATGCCCGATCCTCGATCGGGCATTTTGACTTTTACAGCTTAGCGCCTTCTACAGCTTCACGCGCCAGCTTAGTAATGCGGTCg                                                                               <  2:1468924/150‑1 (MQ=255)
                                                                        aCGCGCGAAGCACGCCGCCGGACAGGTAAAAGTACAGGGTTAAACCCGGCAAAAAAACGCTACAAAAATGCCCGATCCTCGATCGGGCATTTTGACTTTTACAGCTTAGCGCCTTCTACAGCTTCACGCGCCAGCTTAGTAATGCGGTCg                                                                               <  2:1891810/150‑1 (MQ=255)
                                                                         cgcgcgAAGCACGCCGCCGGACAGGTAAAAGTACAGGGTTAAACCCGGCAAAAAAACGCTACAAAAATGCCCGATCCTCGATCGGGCATTTTGACTTTTACAGCTTAGCGCCTTCTACAGCTTCACGCGCCAGCTTAGTAATGCGGTCGt                                                                              >  1:249577/1‑150 (MQ=255)
                                                                         cgcgcgAAGCACGCCGCCGGACAGGTAAAAGTACAGGGTTAAACCCGGCAAAAAAACGCTACAAAAATGCCCGAGGGTCGATCGGGCATTTTGACTTTTACCGCTTAGCGCCTTCTACAGCTTCACGCGCCAGCTTAGTAATGCGGTCgg                                                                              >  1:998051/1‑149 (MQ=255)
                                                                           cgcgAAGCACGCCGCCGGACAGGTAAAAGTACAGGGTTAAACCCGGCAAAAAAACGCTACAAAAATGCCCGATCCTCGATCGGGCATTTTGACTTTTACAGCTTAGCGCCTTCTACAGCTTCACGCGCCAGCTTAGTAATGCGGTCGTaa                                                                            >  2:504874/1‑150 (MQ=255)
                                                                                  cACGCCGCCGGACAGGTAAAAGTACAGGGTTAACCCCGGCAAAAAAACGCTACAAAAATGCCCGATCCTCGATCGGGCATTTTGACTTTTACAGCTTAGCGCCTTCTACAGCTTCACGCGCCAGCTTAGTAATGCGGTCGTAATCGCCCg                                                                     >  1:1550069/1‑150 (MQ=255)
                                                                                   aCGCCGCCGGACAGGTAAAAGTACAGGGTTAAACCCGGCAAAAAAACGCTACAAAAATGCCCGATCCTCGATCGGGCATTTTGACTTTTACAGCTTAGCGACTTCTACAGCTACACGCGCCAGCTTAGAAATGCGGTCGAGATCGccccc                                                                    >  1:294682/1‑148 (MQ=255)
                                                                                      ccgccgGACAGGTAAAAGTACAGGGTTAAACCCGGCAAAAAAACGCTACAAAAATGCCCGATCCTCGATCGGGCATTTTGACTTTTACAGCTTAGCGCCTTCTACAGCTTCACGCGCCAGCTTAGTAATGCGGTCGTAATCGCCCGCTTc                                                                 >  2:976432/1‑150 (MQ=255)
                                                                                      ccgccgGACAGGTAAAAGTACAGGGTTAAACCCGGCAAAAAAACGCTACAAAAATGCCCGATCCTCGATCGGGCATTTTGACTTTTACAGCTTAGCGCCTTCTACAGCTTCACGCGCCAGCTTAGTAATGCGGTCGTAATCGCCCGCTTc                                                                 >  1:192873/1‑150 (MQ=255)
                                                                                      ccgccgGACAGGTAAAAGTACAGGGTTAAACCCGGCAAAAAAACGCTACAAAAATGCCCGATCCTCGATCGGGCATTTTGACTTTTACAGCTTAGCGCCTTCTACAGCTTCACGCGCCAGCTTAGTAATGCGGTCGTAATCGCCCGCTTc                                                                 >  1:1748190/1‑150 (MQ=255)
                                                                                       cgccgGACAGGTAAAAGTACAGGGTTAAACCCGGCAAAAAACCGCTACAAAAATGCCCGCTCCTCGATCGGGCATTTTGACTTTTACAGCTTAGCGCCTTCTACAGCTTCTCGCGCCAGCTTAGTAATGCGGTCGTAATCGCCCGCTTcc                                                                >  1:294947/1‑150 (MQ=255)
                                                                                       cgccgGACAGGTAAAAGTACAGGGTTAAACCCGGCAAAAAAACGCTACAAAAATGCCCGATCCTCGATCGGGCATTTTGACTTTTACAGCTTAGCGCCTTCTACAGCTTCACGCGCCAGCTTAGTAATGCGGTCGTAATCGCCCGCTTcc                                                                >  1:1589583/1‑150 (MQ=255)
                                                                                        gccgGACAGGTAAAAGTACAGGGTTAAACCCGTCAAAAAAACGCTACAAAAATGCCCGATCCTCGATCGGGCATTTTGACTTTTACAGCTTAGCGCCTTCTACAGCTTCACGCGCCAGCTTAGTAATGCGGTCGTAATCGCCCGCTTCCa                                                               >  1:718155/1‑150 (MQ=255)
                                                                                                ggTAAAAGTACAGGGTTAAACCCGGCAAAAAAACGCTACAAAAATGCCCGATCCTCGATCGGGCATTTTGCCTTTTACAGCTTAGCGCCTTCTACAGCTTCACGCGCCAGCTTAGTAATGCGGTCGTAATCGCCCGCTTCCAGCGCATCt                                                       >  2:599552/1‑150 (MQ=255)
                                                                                                ggTAAAAGTACAGGGTTAAACCCGGCAAAAAAACGCTACAAAAATGCCCGATCCTCGATCGGGCATTTTGACTTTTACAGCTTAGCGCCTTCTACAGCTTCACGCGCCAGCTTAGTAATGCGGTCGTAATCGCCCGCTTCCAGCGCATCt                                                       >  1:663603/1‑150 (MQ=255)
                                                                                                ggTAAAAGTACAGGGTTAAACCCGGCAAAAAAACGCTACAAAAATGCCCGATCCTCGATCGGGCATTTTGACTTTTACAGCTTAGCGCCTTCTACAGCTTCACGCGCCAGCTTAGTAATGCGGTCGTAATCGCCCGCTTCCAGCGCATCt                                                       >  1:731718/1‑150 (MQ=255)
                                                                                                ggTAAAAGTACAGGGTTAAACCCGGCAAAAAAACGCTACAAAAATGCCCGATCCTCGATCGGGCATTTTGACTTTTACAGCTTAGCGCCTTCTACAGCTTCACGCGCCAGCTTAGTAATGCGGTCGTAATCGCCCGCTTCCAGCGCATCt                                                       >  2:1618941/1‑150 (MQ=255)
                                                                                                ggTAAAAGTACAGGGTTAAACCCGGCAAAAAAACGCTACAAAAATGCCCGATCCTCGATCGGGCATTTTGACATTTACAGCTTAGCGACTTCTACAGCTTCACGCGCCAGCTTAGTAATTCGGTCGTAATCGCCCGCTTCCAGCGCATCt                                                       >  2:627149/1‑150 (MQ=255)
                                                                                                ggTAAAAGTACAGGGTTAAACCCGGCAAAAAAACGCTACAAAAATGCCCGATCCACGATCGGGCATTTTGACTTTTACAGCTTAGCGCCTTCTACAGCTTCACGCGCCAGCTTAGTAATGCGGTCGTAATCGCCCGCTTCCAGCGCATCt                                                       >  2:309769/1‑150 (MQ=255)
                                                                                                 gTAAAAGTACAGGGTTAAACCCGGCAAAAAAACGCTACAAAAATGCCCGATCCTCGATCGGGCATTTTGACTTTTACAGCTTAGCGCCTTCTACAGCTTCACGCGCCAGCTTAGTAATGCGGTCGTAATCGCCCGCTTCCAGCGCATCTg                                                      <  1:1867955/150‑1 (MQ=255)
                                                                                                  tAAAAGTACAGGGTTAAACCCGGCAAAAAAACGCTACAAAAATGCCCGATCCTCGATCGGGCATTTTGACTTTTACAGCTTAGCGCCTTCTACAGCTTCACGCGCCAGCTTAGTAATGCGGTCGTAATCGCCCGCTTCCCGCGCATCTGc                                                     >  2:976601/1‑150 (MQ=255)
                                                                                                  tAAAAGTACAGGGTTAAACCCGGCAAAAAAACGCTACAAAAATGCCAGATCCTCGAACGGGAATTTTGACTGTTACAGCTTAGCGCCTTCTACAGCTTCACGCGCCAGCTTAGTAATGCGGACGTAATCGCCCGCTTCCAGCGCAGCTGc                                                     >  1:1386622/1‑150 (MQ=255)
                                                                                                   aaaaGTAGAGGGGTAAACCGGGCAAAAAAACGCTACAAAAATGCCCGATCCTCGATCGGGCATTTTGACTTTTACTGCTTAGCGCCTTCTACAGCATCACGCGCCAGCTTAGTAATGCGGTCGTAATCGCCCGCTTCCAGCGCATCTGcc                                                    <  2:395040/150‑1 (MQ=255)
                                                                                                   aaaaGTACAGGGTTAAACCCGGCAAAAAAACGCTACAAAAATGCCCGATCCTCGATCGGGCATTTTGACTTTTACAGCTTAGCGCCTTCTACAGCTTCACGCGCCAGCTTAGTAATGCGGTCGTAATCGCCCGCTTCCAGCGCATCTAcc                                                    >  1:629796/1‑150 (MQ=255)
                                                                                                   aaaaGTACAGGGTTAAACCCGGCAAAAAAAAGCTACAAAAATGCCCGATCCTCGATCGGGCATTTTGACTTTTACAGCTTAGCGCCTTCTACAGCTTCACGCGCCAGCTTAGTAATGCTGTCGTAATCGCCCGCTTCCAGCGCATCTGcc                                                    <  1:1413930/150‑1 (MQ=255)
                                                                                                     aaGTACAGGGTTAAACCCGGCAAAAAAACGCTACAAAAATGCCCGATCCTCGATCGGGCATTTTGACTTTTACAGCTTAGCGCCTTCTACAGCTTCACGCGCCAGCTTAGTAATGCGGTCGTAATCGCCCGCTTCCAGCGCATCTGCCgg                                                  <  2:1273848/150‑1 (MQ=255)
                                                                                                           aGGGTTAAACCCGGCAAAAAAACGCTACAAAAATGCCCGATCCTCGATCGGGCATTTTGACTTTTACAGCTTAGCGCCTTCTACAGCTTCACGCGCCAGCTTAGTAATGCGGTCGTAATCGCCCGCTTCCAGCGCATCTGCCGGCAccag                                            >  2:1282011/1‑150 (MQ=255)
                                                                                                            gggTTAAACCCGGCAAAAAAACGCTACAAAAATGCCCGATCCTCGATCGGGCATTTTGACTTTTACAGCTTAGCGCCTTCTACAGCTTCACGCGCCAGCTTAGTAATGCGGTCGTAATCGCCCGCTTCCAGCGCATCTGCCGGAAccagc                                           <  2:1191378/150‑1 (MQ=255)
                                                                                                            gggTTAAACCCGGCAAAAAAACGCTACAAAAATGCCCGATCCTCGATCGGGCATTTTGACTTTTACAGCTTAGCGCCTTCTACAGCTTCACGCGCCAGCTTAGTAATGCGGTCGTAATCGCCCGCTTCCAGCGCATCTGCCGGAAccagc                                           >  1:406020/1‑150 (MQ=255)
                                                                                                             ggTTAAACCCGGCAAAAAAACGCTACAAAAATGCCCGATCCTCGATCGGGCATTTTGACTTTTACAGCTTAGCGCCTTCTACAGCTTCACGCGCCAGCTTAGTAATGCGGTCGTAATCGCCCGCTTCCAGCGCATCTGCCGGAAccagcc                                          <  2:249577/150‑1 (MQ=255)
                                                                                                             ggTTAAACCCGGCAAAAAAACGCTACAAAAATGCCCGATCCTCGATCGGGCATTTTGACTTTTACAGCTTAGCGCCTTCTACAGCTTCACGCGCCAGCTTAGTAATGCGGTCGTAATCGCCCGCTTCCAGCGCATCTGCCGGAAccagcc                                          >  2:1572880/1‑150 (MQ=255)
                                                                                                              gTTAAACCCGGCAAAAAAACGCTACAAAAATGCCCGATCCTCGATCGGGCATTTTGACTTTTACAGCTTAGCGCCTTCTACAGCTTCACGCGCCAGCTTAGTAATGCGGTCGTAATCGCCCGCTTCCAGCGCATCTGCCGGAAccagcca                                         >  2:1731535/1‑150 (MQ=255)
                                                                                                               ttAAACCCGGCAAAAAAACGCTACAAAAATGCCCGATCCTCGATCGGGCATTTTGACTTTTACAGCTTAGCGCCTTCTACAGCTTCACGCGCCAGCTTAGTAATGCGGTCGTAATCGCCCGCTTCCAGCGCATCTGCCGGAAccagccag                                        <  1:1352267/150‑1 (MQ=255)
                                                                                                                 aaaCCCGGCAAAAAAACGCTACAAAAATGCCCGATCCTCGCTCGGGCATTTTGACTTTTACAGCTTAGCGCCTTCTACAGCGTCACGCGCCAGCTTAGTAATGCGGTCGTAATCGCCCGCTTCCAGCGCATCTGCCGGAACCAGCCAGGa                                      >  1:115942/1‑150 (MQ=255)
                                                                                                                 aaaCCCGGCAAAAAAACGCTACAAAAATGCCCGATCCTCGATCGGGCATTTTGACTTTTACAGCTTAGCGCCTTCTACAGCTTCACGCGCCAGCTTAGTAATGCGGTCGTAATCGCCCGCTTCCAGCGCATCTGCCGGAACCAGCCAGGa                                      <  2:1358244/150‑1 (MQ=255)
                                                                                                                  aaCCCGGCAAAAAAACGCTACAAAAATGCCCGATCCTCGATCGGGCATTTTGACTTTTACAGCTTAGCGCCTTCTACAGCTTCACGCGCCAGCTTATTAATGCGGTCGTAATCGCCCGCATCCAGCGCCTCTGCCGGAACCAGACAGGaa                                     >  2:1270447/1‑150 (MQ=255)
                                                                                                                     ccGGCAAAAAAACGCTACAAAAATGCCCGATCCTCGATCGGGCATGTTGACTTTTACAGCTTAGCGCCTTCTACAGCTTCACGCGCCAGCTTAGTAATGCGGTCGTAATCGCCCGCTTCCAGCGCATCTGCCGGAACCAGCCAGGAacca                                  >  1:1258426/1‑150 (MQ=255)
                                                                                                                     ccGGCAAAAAAACGCTACAAAAATGCCCGAGCCTCGATCGGGCATTTTGACTTTTACAGCTTAGCGCCTTCTACAGCTTCACGCGCCAGCTTAGTAATGCGGTCGTAATCGCCCGCTTCCAGCGCATCTGCCGGAACCAGCCAGGAacca                                  <  1:1341316/150‑1 (MQ=255)
                                                                                                                      cGGCAAAAAAACGCTACAAAAATGCCCGATCCTCGATCGGGCATTTTGACTTTTACAGCTTAGCGCCTTCTACAGCTTCACGCGCCAGCTTAGTAATGCGGTCGTAATCGCCCGCTTCCAGCGCATCTGCCGGAACCAGCCAGGAaccac                                 >  1:556805/1‑150 (MQ=255)
                                                                                                                      cGGCAAAAAAACGCTACAAAAATGCCCGATCCTCGATCGGGCATTTTGACTTTTACAGCTTAGCGCCGTCTACAGCTTCACGCGCCAGCTTAGTAATGCGGTCGTAATCGCCCGCTTCCAGCGCATCGGCCGGAACCAGCCAGGAcccac                                 >  1:1812725/1‑150 (MQ=255)
                                                                                                                        ccAAAAAAACGCTACAAAAATGCCCGATCCTCGATCGGGCATTTTGACTTTTACAGCTTAGCGCCTTCTACAGCTTCACGCGCCAGCTTAGTAATGCGGTCGTAATCGCCCGCTTCCAGCGCATCTGCCGGAACCAGCCAGGAACCACCg                               <  2:890660/149‑1 (MQ=255)
                                                                                                                         cAAAAAAACGCTACAAAAATGCCCGATCCTCGATCGGGCATTTTGACTTTTACAGCTTAGCGCCTTCTACAGCTTCACGCGCCAGCTTAGTAATGCGGTCGTAATCGCCCGCTTCCAGCGCATCTGCCGGAACCAGCCAGGACCCACCGa                              >  1:870987/1‑150 (MQ=255)
                                                                                                                         cAAAAAAACGCTACAAAAATGCCCGATCCTCGATCGGGCATTTTGACTTTTACAGCTTAGCGCCTTCTACAGCTTCACGCGCCAGCTTAGTAATGCGGTCGTAATCGCCCGCTTCCAGCGCATCTGCCGGAACCAGCCAGGAACCACCGa                              >  2:356575/1‑150 (MQ=255)
                                                                                                                           aaaaaaCGCTACAAAAATGCCCGATCCTCGATCGGGCATTTTGACTTTTACAGCTTAGCGCCTTCTACAGCTTCACGCGCCAGCTTAGTAATGCGGTCGTAATCGCCCGCTTCCAGCGCATCTGCCGGAACCAGCCAGGAACCACCGATg                            <  2:1550069/150‑1 (MQ=255)
                                                                                                                           aaaaaaCGCTACAAAAATGCCCGATCCTCGATCGGGCATTTTGACTTTTACAGCTTAGCGCCTTCTACAGCTTCACGCGCCAGCTTAGTAATGCGGTCGTAATCGCCCGCTTCCAGCGCATCTGCCGGAACCAGCCAGGAACCACCGATg                            >  1:1060364/1‑150 (MQ=255)
                                                                                                                              aaaCGCTACAAAAATGCCCGATCCTCGATCGGGCATTTTGACTTTTACAGCTTAGCGCCTTCTACAGCTTCACGCGCCAGCTTAGTAATGCGGTCGTAATCGCCCGCTTCCAGCGCATCTGCCGGAACCAGCCAGGAACCACCGATGcac                         <  2:525705/150‑1 (MQ=255)
                                                                                                                              aaaCGCTACAAAAATGCCCGATCCTCGATCGGGCATTTTGACTTTTACAGCTTAGCGCCTTCTACAGCTTCACGCGCCAGCTTAGTAATGCGGTCGTAATCGCCCGCTTCCAGCGCATCTGCCGGAACCAGCCAGGAACCACCGATGcac                         <  1:1551850/150‑1 (MQ=255)
                                                                                                                              aaaCGATACAAAAATGCCCGATCCTCGATCGGGCATTTTGACTTTTACAGCTTAGCGCCTTCTACAGCTTCACGCGCCAGCTTAGTAATGCGGTCGTAATCGCCCGCTTCCAGCGCATCTGCCGGAACCAGCCAGGAACCACCGATGcac                         <  2:1483080/150‑1 (MQ=255)
                                                                                                                                aCGCTACAAAAATGCCCGATCCTCGATCGGGCATTTTGACTTTTACAGCTTAGCGCCTTCTACAGCTTCACGCGCCAGCTTAGTAATGCGGTCGTAATCGCCCGCTTCCAGCGCATCTGCCGGAACCAGCCAGGAACCACCGATGCACAg                       <  1:1569644/150‑1 (MQ=255)
                                                                                                                                aCGCTACAAAAATGCCCGATCCTCGATCGGGCATTTTGACTTTTACAGCTTAGCGCCTTCTACAGCTTCACGCGCCAGCTTAGTAATGCGGTCGTAATCGCCCGCTTCCAGCGCATCTGCCGGAACCAGCCAGGAACCACCGATGCACAg                       <  1:356575/150‑1 (MQ=255)
                                                                                                                                   cTACAAAAATGCCCGATCCTCGATCGGGCATTTTGACTTTTACAGCTTAGCGCCTTCTACAGCTTCACGCGCCAGCTTAGTAATGCGGTCGTAATCGCCCGCTTCCAGCGCATCTGCCGGAACCAGCCAGGAACCACCGATGCACAGCAc                    >  2:1640913/1‑150 (MQ=255)
                                                                                                                                   aagtcaaaaTGCCCGATCGAGGATCGGGCATTTTGACTTTTACAGCTTAGCGCCTTCTACAGCTTCACGCGCCAGCTTAGTAATGCGGTCGTAATCGCCCGCTTCCAGCGCATCTGCCGGAACCAGCCAGGAACCACCGATGCACAGCAc                    >  1:1897382/6‑150 (MQ=255)
                                                                                                                                     aCAAAAATGCCCGATCCTCGATCGGGCATTTTGACTTTTACAGCTTAGCGCCTTCTACAGCTTCACGCGCCAGCTTAGTAATGCGGTCGTAATCGCCCGCTTCCAGCGCATCTGCCGGAACCAGCCAGGAACCACCGATGCACAGCACGc                  >  1:1707169/1‑150 (MQ=255)
                                                                                                                                      cAAAAATGCCCGATCCTCGATCGGGCATTTTGACTTTTACAGCTTAGCGCCTTCTACAGCTTCACGCGCCAGCTTAGTAATGCGGTCGTAATCGCCCGCTTCCAGCGCATCTGCCGGGACCAGCCAGGAACCACCGATGCACAGCCCGcg                 >  1:1870845/1‑149 (MQ=255)
                                                                                                                                       caaaaTGCCCGATCGAGGATCGGGCATTTTGACTTTTACAGCTTAGCGCCTTCTACAGCTTCACGCGCCAGCTTAGTAATGCGGTCGTAATCGCCCGCTTCCAGCGCATCTGCCGGAACCAGCCAGGAACCACCGATGCACAGCACGCtt                >  1:1534199/2‑150 (MQ=255)
                                                                                                                                          aaTGCCCGATCCTCGATCGGGCATTTTGACTTTTACAGCTTAGCGCCTTCTACAGCTTCACGCGCCAGCTTAGTAATGCGGTCGTAATCGCCCGCTTCCAGCGCATCTGCCGGAACCAGCCAGGAACCACCGATGCACAGCACGCtttta             >  2:2002041/1‑149 (MQ=255)
                                                                                                                                          aaTGCCCGATCCTCGATCGGGCATTTTGACTTTTACAGCTTAGCGCCTTCTACAGCTTCACGCGCCAGCTTAGTAATGCGGTCGTAATCGCCCGCTTCCAGCGCATCTGCCGGAACCAGCCAGGAACCACCGATGCACAGCACGCTTTTc             <  1:1262027/150‑1 (MQ=255)
                                                                                                                                           aTGCCCGATCCTCGATCGGGCATTTTGACTTTTACAGCTTAGCGCCTTCTACAGCTTCACGCGCCAGCTTAGTAATGCGGTCGTAATCGCCCGCTTCCAGCGCATCTGCCGGAACCAGCCAGGAACCACCGATGCACAGCACGCTTTTCa            <  1:976601/150‑1 (MQ=255)
                                                                                                                                            tGCCCGATCGAGGATCGGGCATTTTGACTTTTACAGCTTAGCGCCTTCTACAGCTTCACGCGCCAGCTTAGTAATGCGGTCGTAATCGCCCGCTTCCAGCGCATCTGCCGGAACCAGCCAGGAACCACCGATGCACCGCACGCTTTTCAg           >  1:225832/1‑150 (MQ=255)
                                                                                                                                            tGCCCGATCCTCGATCGGGCATTTTGACTTTTACAGCTTAGCGCCTTCTACAGCTTCACGCGCCAGCTTAGTAATGCGGTCGTAATCGCCCGCTTCCAGCGCATCTGCCGGAACCAGCCAGGAACCACCGATGCACAGCACGCTTTTCAg           <  1:1393175/150‑1 (MQ=255)
                                                                                                                                              cccGATCCTCGATCGGGCATTTTGACTTTTACAGCTTAGCGCCTTCTACAGCTTCACGCGCCAGCTTAGTAATGCGGGCGTAATCGCCCGCTTCCAGCGCATCTGCCGGAACCAGCCAGGAACCACCGATGCACAGCACGCTTTTCAgcg         <  2:1069387/150‑1 (MQ=255)
                                                                                                                                                cGATCCTCGATCGGGCATTTTGACTTTTACAGCTTAGCGCCTTCTACAGCTTCACGCGCCAGCTTAGTAATGCGGTCGTAATCGCCCGCTTCCAGCGCATCTGCCGGAACCAGCCAGGAACCACCGATGCACAGCCCGCTTTTCCGCGcc       >  1:361359/1‑150 (MQ=255)
                                                                                                                                                 gATCCTCGATCGGGCATTTTGACTTTTACAGCTTAGCGCCTTCTACAGCTTCACGCGCCAGCTTAGTAATGCGGTCGTAATCGCCCGCTTCCAGCGCATCTGCCGGAACCAGCCAGGAACCACCGATGCACAGCACGCTTTTCAGCGCCa      >  1:859723/1‑150 (MQ=255)
                                                                                                                                                   tCCTCGATCGGGCATTTTGACTTTTACATCTTAGCGCCTTCTACTGCTTCACGCGCGAGCTTTGTAATGCGGTCGTAATCGCCCGCTTCCAGCGCATCTGCCGGAACCAGCCAGGAACCACCGATGCACAGCACGCTTTTCAGCGCCAgg    <  1:327488/150‑1 (MQ=255)
                                                                                                                                                    ccTCGATCGTGCATTTTGACTTTTACAGCTAAGCGCAGTCTTCAGCTTCACGGGCCAGCTTAGTAATGCGGTCGTAATCGCCCGCTTCCAGCGCATCTGCCGGAACCAGCCAGGAACCACCGATGCACAGCACGCTTTTCAGCTCCAGGt   <  1:366619/150‑1 (MQ=255)
                                                                                                                                                    ccTCGATCGGGCATTTTGACTTTTACAGCTTAGCGCCTTCTACAGCTTCACGCGCCAGCTTAGTAATGCGGTCGTAATCGCCCGCTTCCAGCGCATCTGCCGGAACCAGCCAGGAACCACCGATGCACAGCACGCTTTTCAGCGCCAGGt   <  1:1675612/150‑1 (MQ=255)
                                                                                                                                                     cTCGATCGGGCATTTTGACTTTTACAGCTTAGCGCCTTCTACAGCTTCACGCGCCAGCTTAGTAATGCGGTCGTAATCGCCCGCTTCCAGCGCATCTGCCGGAACCAGCCAGGAACCACCGATGCACAGCACGCTTTTCAGCGCCAGGTa  >  2:541678/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
GCCGGAAATTGATGGCAGCCGTCGTCTGGGGGTGGCACTGGCTACTGCAGAGAGTGTTGTTGACGCCATTGAACGCGCGAAGCACGCCGCCGGACAGGTAAAAGTACAGGGTTAAACCCGGCAAAAAAACGCTACAAAAATGCCCGATCCTCGATCGGGCATTTTGACTTTTACAGCTTAGCGCCTTCTACAGCTTCACGCGCCAGCTTAGTAATGCGGTCGTAATCGCCCGCTTCCAGCGCATCTGCCGGAACCAGCCAGGAACCACCGATGCACAGCACGCTTTTCAGCGCCAGGTA  >  NC_000913/1931945‑1932243

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: ‑