Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,011,930 T→C 6.6% intergenic (+59/‑71) uup → / → pqiA ATP‑binding protein with possible role in replication/intermembrane transport protein PqiA

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,011,9300TC 6.6% 320.9 / 9.5 106intergenic (+59/‑71)uup/pqiAATP‑binding protein with possible role in replication/intermembrane transport protein PqiA
Reads supporting (aligned to +/- strand):  ref base T (50/49);  new base C (5/2);  total (55/51)
Fisher's exact test for biased strand distribution p-value = 4.40e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.79e-01

Only 100 of 106 total aligned reads displayed.

GTGCTTGCTGATATGGCTGCTGCAGAGCAGGAGCTGGAGCAAGCCTTTGAACGCTGGGAGTATCTTGAAGCGTTAAAAAATGGTGGCTGATCGCAAAATAGTCGATTTAAACAGAGACAGCTGACAACGGTAAAATTGCCTGATGAACTGGGTTTATCAGGCATGAAAACCAACACTTTGTTAGTCATCTGTAAAATAGCGCATCATTAAGGAGTACCAATGTGCGAACATCATCATGCCGCGAAGCACATCCTGTGCTCGCAGTGTGACATGCTGGTGGCGTTACCGCGCCTTGAG  >  NC_000913/1011782‑1012078
                                                                                                                                                    |                                                                                                                                                    
gTGCTTGCTGATATGGCTGCTGCAGAGCAGGAGCTGGAGCAAGCCTTTGAACGCTGGGAGTATCTTGAAGCGTTAAAAAATGGTGGCTGATCGCAAAATAGTCGATTTAAACAGAGACAGCTGACAACGGTAAAATTGCCTGATGAACTg                                                                                                                                                     <  2:533678/150‑1 (MQ=255)
  gcttgctGATATGGCTGCTGCAGAGCAGGAGCTGGAGCAAGCCTTTGAACGCTGGGAGTATCTTGAAGCGTTAAAAAATGGTGGCTGATCGCAAAATAGTCGATTTAAACAGAGACAGCTGACAACGGTAAAATTGCCTGATGAACTggg                                                                                                                                                   >  1:907331/1‑150 (MQ=255)
   cttgctGATATGGCTGCTGCAGAGCAGGAGCTGGAGCAAGCCTTTGAACGCTGGGAGTATCTTGAAGCGTTAAAAAATGGTGGCTGATCGCAAAATAGTCGATTTAAACAGAGACAGCTGACAACGGTAAAATTGCCTGATGAACTGGGt                                                                                                                                                  <  1:1391270/150‑1 (MQ=255)
    ttgctGATATGGCTGCTGCAGAGCAGGAGCTGGAGCAAGCCTTTGAACGCTGGGAGTATCTTGAAGCGTTAAAAAATGGTGGCTGATCGCAAAATAGTCGATTTAAACAGAGACAGCTGACAACGGTAAAATTGCCTGATGAACTGGGtt                                                                                                                                                 <  2:752246/150‑1 (MQ=255)
     tgctGATATGGCTGCTGCAGAGCAGGAGCTGGAGCAAGCCTTTGAACGCTGGGAGTATCTTGAAGCGTTAAAAAATGGTGGCTGATCGCAAAATAGTCGATTTAAACAGAGACAGCTGACAACGGTAAAATTGCCTGATGAACTGGGttt                                                                                                                                                >  1:580963/1‑150 (MQ=255)
     tgctGATATGGCTGCTGCAGAGCAGGAGCTGGAGCAAGCCTTTGAACGCTGGGAGTATCTTGAAGCGTTAAAAAATGGTGGCTGATCGCAAAATAGTCGATTTAAACAGAGACAGCTGACAACGGTAAAATTGCCTGATGAACTGGGttt                                                                                                                                                >  1:1406605/1‑150 (MQ=255)
       ctGATATGGCTGCTGCAGAGCAGGAGCTGGAGCAAGCCTTTGAACGCTGGGAGTATCTTGAAGCGTTAAAAAATGGTGGCTGATCGCAAAATAGTCGATTTAAACAGAGACAGCTGACAACGGTAAAATTGCCTGATGAACTGGGTTTAt                                                                                                                                              >  1:1526368/1‑150 (MQ=255)
        tGATATGGCTGCTGCAGAGCAGGAGCTGGAGCAAGCCTTTGAACGCTGGGAGTATCTTGAAGCGTTAAAAAATGGTGGCTGATCGCAAAATAGTCGATTTAAACAGAGACAGCTGACAACGGTAAAATTGCCTGATGAACTGGGTTTATc                                                                                                                                             <  1:1564242/150‑1 (MQ=255)
         gATATGGCTGCTGCAGAGCAGGAGCTGGAGCAAGCCTTTGAACGCTGGGAGTATCTTGAAGCGTTAAAAAATGGTGGCTGATCGCAAAATAGTCGATTTAAACAGAGACAGCTGACAACGGTAAAATTGCCTGATGAACTGGGTTTATCa                                                                                                                                            >  1:1452034/1‑150 (MQ=255)
          ttatGGCTGCTGCAGAGCAGGAGCTGGAGCAAGCCTTTGAACGCTGGGAGTATCTTGAAGCGTTAAAAAATGGTGGCTGATCGCAAAATAGTCGATTTAAACAGAGACAGCTGACAACGGTAAAATTGCCTGATGAACTGTGTTTATCAg                                                                                                                                           <  2:946864/149‑1 (MQ=255)
           tatGGCTGCTGCAGAGCAGGAGCTGGAGCAAGCCTTTGAACGCTGGGAGTATCTTGAAGCGTTAAAAAATGGTGGCTGATCGCAAAATAGTCGATTTAAACAGAGACAGCTGACAACGGTAAAATTGCCTGATGAACTGGGTTTATCAgg                                                                                                                                          <  1:763340/150‑1 (MQ=255)
            atGGCTGCTGCAGAGCAGGAGCTGGAGCAAGCCTTTTAACGCTGGGAGTATCTTGAAGCGTTAAAAAATGGTGGCTGATCGCAAAATATTCGATTTAAACAGAGACAGCTGACAACGGTAAAATTGCCTGATGAACTGGGTTTATCAGGc                                                                                                                                         <  2:1563333/150‑1 (MQ=255)
             tGGCTGCTGCAGAGCAGGAGCTGGAGCAAGCCTTTGAACGCTGGGAGTATCTTGAAGCGTTAAAAAATGGTGGCTGATCGCAAAATAGTCGATTTAAACAGAGACAGCTGACAACGGTAAAATTGCCTGATAAACCCAGTTCATCAGGCa                                                                                                                                        <  1:1542268/150‑1 (MQ=255)
              ggCTGCTGCAGAGCAGGAGCTGGAGCAAGCCTTTGAACGCTGGGAGTATCTTGAAGCGTTAAAAAATGGTGGCTGATCGCAAAATAGTCGATTTAAACAGAGACAGCTGACAACGGTAAAATTGCCTGATGAACTGGGTTTATCAGGCAt                                                                                                                                       >  1:1025734/1‑150 (MQ=255)
              ggCTGCTGCAGAGCAGGAGCTGGAGCAAGCCTTTGAACGCTGGGAGTATCTTGAAGCGTTAAAAAATGGTGGCTGATCGCAAAATAGTCGATTTAAACAGAGACAGCTGACAACGGTAAAATTGCCTGATAAACCCAGTTCATCAGGCaa                                                                                                                                       >  2:1244004/1‑149 (MQ=255)
                  gctgcAGAGCAGGAGCTGGAGCAAGCCTTTGAACGCTGGGAGTATCTTGAATCGTTAAAAAATGGTGGCTGATCGCAAAATAGTCGATTTAAACAGAGACAGCTGACAACGGTAAAATTGCCTGATGAACTGGGTTTATCAGGCATGaaa                                                                                                                                   <  2:596921/150‑1 (MQ=255)
                  gctgcAGAGCAGGAGCTGGAGCAAGCCTTTGAACGCTGGGAGTATCTTGAAGCGTTAAAAAATGGTGGCTGATCGCAAAATAGTCGATTTAAACAGAGACAGCTGACAACGGTAAAATTGCCTGATGAACTGGGTTTATCAGGCATGaaa                                                                                                                                   >  1:29567/1‑150 (MQ=255)
                  gctgcAGAGCAGGAGCTGGAGCAAGCCTTTGAACGCTGGGAGTATCTTGAAGCGTTAAAAAATGGTGGCTGATCGCAAAATAGTCGATTTAAACAGAGACAGCTGACAACGGTAAAATTGCCTGATAAACCCAGTTCATCAGGCaatttt                                                                                                                                   <  1:1244004/150‑6 (MQ=255)
                   ctgcAGAGCAGGAGCTGGAGCAAGCCTTTGAACGCTGGGAGTATCTTGAAGCGTTAAAAAATGGTGGCTGATCGCAAAATAGTCGATTTAAACAGAGACAGCTGACAACGGTAAAATTGCCTGATGAACTGGGTTTATCAGGCATGaaaa                                                                                                                                  >  1:1012334/1‑150 (MQ=255)
                   ctgcAGAGCAGGAGCTGGAGCAAGCCTTTGAACGCTGGGAGTATCTTGAAGCGTTAAAAAATGGTGGCTGATCGCAAAATAGTCGATTTAAACAGAGACAGCTGACAACGGTAAAATTGCCTGATGAACTGGGTTTATCAGGCATGaaaa                                                                                                                                  <  1:51350/150‑1 (MQ=255)
                    tgcAGAGCAGGAGCTGGAGCAAGCCTTTGAACGCTGGGAGTATCTTGAAGCGTTAAAAAATGGTGGCTGATCGCAAAATAGTCGATTTAAACAGAGACAGCTGACAACGGTAAAATTGCCTGATGAACTGGGTTTATCAGGCATGAAAAc                                                                                                                                 >  2:360318/1‑150 (MQ=255)
                     gcAGAGCAGGAGCTGGAGCAAGCCTTTGAACGCTGGGAGTATCTTGAAGCGTTAAAAAATGGTGGCTGATCGCAAAATAGTCGATTTAAACAGAGACAGCTGACAACGGTAAAATTGCCTGATGAACTGGGTTTATCAGGCATGAAAAcc                                                                                                                                <  2:619743/150‑1 (MQ=255)
                     gcAGAGCAGGAGCTGGAGCAAGCCTTTGAACGCTGGGAGTATCTTGAAGCGTTAAAAAATGGTGGCTGATCGCAAAATAGTCGATTTAAACAGAGACAGCTGACAACGGTAAAATTGCCTGATGAACTGGGTTTATCAGGCATGAAAAcc                                                                                                                                <  1:1586663/150‑1 (MQ=255)
                     gcAGAGCAGGAGCTGGAGCAAGCCTTTGAACGCTGGGAGTATCTTGAAGCGTTAAAAAATGGTGGCTGATCGCAAAATAGTCGATTTAAACAGAGACAGCTGACAACGGTAAAATTGCCTGATGAACTGGGTTTATCAGGCATGAAAAcc                                                                                                                                >  1:287310/1‑150 (MQ=255)
                      cAGAGCAGGAGCTGGAGCAAGCCTTTGAACGCTGGGAGTATCTTGAAGCGTTAAAAAATGGTGGCTGATCGCAAAATAGTCGATTTAAACAGAGACAGCTGACAACGGTAAAATTGCCTGATGAACTGGGTTTATCAGGCATGAAAACCa                                                                                                                               <  1:540768/150‑1 (MQ=255)
                      cAGAGCAGGAGCTGGAGCAAGCCTTTGAACGCTGGGAGTATCTTGAAGCGTTAAAAAATGGTGGCTGATCGCAAAATAGTCGATTTAAACAGAGACAGCTGACAACGGTAAAATTGCCTGATGAACTGGGTTTATCAGGCATGAAAACCa                                                                                                                               >  2:1154003/1‑150 (MQ=255)
                             ggAGCTGGAGCAAGCCTTTGAACGCTGGGAGTATCTTGAAGCGTTAAAAAATGGTGGCTGATCGCAAAATAGTCGATTTAAACAGAGACAGCTGACAACGGTAAAATTGCCTGATGAACTGGGTTTATCAGGCATGAAAACCAACACttt                                                                                                                        >  2:1555581/1‑150 (MQ=255)
                             ggAGCTGGAGCAAGCCTTTGAACGCTGGGAGTATCTTGAAGCGTTAAAAAATGGTGGCTGATCGCAAAATAGTCGATTTAAACAGAGACAGCTGACAACGGTAAAATTGCCTGATGAACTGGGTTTATCAGGCATGAAAACCAACACttt                                                                                                                        >  1:702544/1‑150 (MQ=255)
                             ggAGCTGGAGCAAGCCTTTGAACGCTGGGAGTATCTTGAAGCGTTAAAAAATGGTGGCTGATCGCAAAATAGTCGATTTAAACAGAGACAGCTGACAACGGTAAAATTGCCTGATGAACTGGGTTTATCAGGCATGAAAACCAACACttt                                                                                                                        >  2:1457334/1‑150 (MQ=255)
                              gAGCTGGAGCAAGCCTTTGAACGCTGGGAGTATCTTGAAGCGTTAAAAAATGGTGGCTGATCGCAAAATAGTCGATTTAAACAGAGACAGCTGACAACGGTAAAATTGCCTGATGAACTGGGTTTATCAGGCATGAAAACCAACACTTTg                                                                                                                       <  2:580963/150‑1 (MQ=255)
                                gCTGGAGCAAGCCTTTGAACGCTGGGAGTATCTTGAAGCGTTAAAAAATGGTGGCTGATCGCAAAATAGTCGATTTAAACAGAGACAGCTGACAACGGTAAAATTGCCTGATGAACTGGGTTTATCAGGCATGAAAACCAACACTTTGtt                                                                                                                     >  2:1009147/1‑150 (MQ=255)
                                 cTGGAGCAAGCCTTTGAACGCTGGGAGTATCTTGAAGCGTTAAAAAATGGTGGCTGATCGCAAAATAGTCGATTTAAACAGAGACAGCTGACAACGGTAAAATTGCCTGATGAACTGGGTTTATCAGGCATGAAAACCAACACTTTGTTa                                                                                                                    >  1:1132708/1‑150 (MQ=255)
                                  tGGAGCAAGCCTTTGAACGCTGGGAGTATCTTGAAGCGTTAAAAAATGGTGGCTGATCGCAAAATAGTCGATTTAAACAGAGACAGCTGACAACGGTAAAATTGCCTGATGAACTGGGTTTATCAGGCATGAAAACCAACACTTTGTTAg                                                                                                                   <  2:398150/150‑1 (MQ=255)
                                   ggAGCAAGCCTTTGAACGCTGGGAGTATCTTGAAGCGTTAAAAAATGGTGGCTGATCGCAAAATAGTCGATTTAAACAGAGACAGCTGACAACGGTAAAATTGCCTGATGAACTGGGTTTATCAGGCATGAAAACCAACACTTTGTTAGt                                                                                                                  >  2:210271/1‑150 (MQ=255)
                                   ggAGCAAGCCTTTGAACGCTGGGAGTATCTTGAAGCGTTAAAAAATGGTGGCTGATCGCAAAATAGTCGATTTAAACAGAGACAGCTGACAACGGTAAAATTGCCTGATGAACTGGGTTTATCAGGCATGAAAACCAACACTTTGTTAGt                                                                                                                  >  2:204402/1‑150 (MQ=255)
                                   ggAGCAAGCCTTTGAACGCTGGGAGTATCTTGAAGCGTTAAAAAATGGTGGCTGATCGCAAAATAGTCGATTTAAACAGAGACAGCTGACAACGGTAAAATTGCCTGATGAACTGGGTTTATCAGGCATGAAAACCAACACTTTGTTAGt                                                                                                                  >  2:455723/1‑150 (MQ=255)
                                    gAGCAAGCCTTTGAACGCTGGGAGTATCTTGAAGCGTTAAAAAATGGTGGCTGATCGCAAAATAGTCGATTTAAACAGAGACAGCTGACAACGGTAAAATTGCCTGATGAACTGGGTTTATCAGGCATGAAAACCAACACTTTGTTAGTc                                                                                                                 >  1:1354848/1‑150 (MQ=255)
                                     aGCAAGCCTTTGAACGCTGGGAGTATCTTGAAGCGTTAAAAAATGGTGGCTGATCGCAAAATAGTCGATTTAAACAGAGACAGCTGACAACGGTAAAATTGCCTGATGAACTGGGTTTATCAGGCATGAAAACCAACACTTTGTTAGTCa                                                                                                                >  2:372089/1‑150 (MQ=255)
                                           ccTTTGAACGCTGGGAGTATCTTGAAGCGTTAAAAAATGGTGGCTGATCGCAAAATAGTCGATTTAAACAGAGACAGCTGACAACGGTAAAATTGCCTGATGAACTGGGTTTATCAGGCATGAAAACCAACACTTTGTTAGTCATCTGTa                                                                                                          >  2:1486601/1‑150 (MQ=255)
                                            cTTTGAACGCTGGGAGTATCTTGAAGCGTTAAAAAATGGTGGCTGATCGCAAAATAGTCGATTTAAACAGAGACAGCTGACAACGGTAAAATTGCCTGATGAACTGGGTTTATCAGGCATGAAAACCAACACTTTGTTAGTCATCTGTaa                                                                                                         <  1:483769/150‑1 (MQ=255)
                                                        ggAGTATCTTGAAGCGTTAAAAAATGGTGGCTGATCGCAAAATAGTCGATTTAAACAGAGACAGCTGACAACGGTAAAATTGCCTGATGAACTGGGTTTATCAGGCATGAAAACCAACACTTTGTTAGTCATCTGTAAAATAGCGcatca                                                                                             <  2:690315/150‑1 (MQ=255)
                                                         gAGTATCTTGAAGCGTTAAAAAATGGTGGCTGATCGCAAAATAGTCGATTTAAACAGAGACAGCTGACAACGGTAAAATTGCCTGATGAACTGGGTTTATCAGGCATGAAAACCAACACTTTGTTAGTCATCTGTAAAATAGCGcatcat                                                                                            <  1:1032617/150‑1 (MQ=255)
                                                           gTATCTTGAAGCGTTAAAAAATGGTGGCTGATCGCAAAATAGTCGATTTAAACAGAGACAGCTGACAACGGTAAAATTGCCTGATGAACTGGGTTTATCAGGCATGAAAACCAACACTTTGTTAGTCATCTGTAAAATAGCGCATCATTa                                                                                          <  2:1221573/150‑1 (MQ=255)
                                                                ttGAAGCGTTAAAAAATGGTGGCTGATCGCAAAATAGTCGATTTAAACAGAGACAGCTGACAACGGTAAAATTGCCTGATGAACTGGGTTTATCAGGCATGAAAACCAACACTTTGTTAGTCATCTGTAAAATAGCGCATCATTAAGGAg                                                                                     <  2:268654/150‑1 (MQ=255)
                                                                  gAAGCGTTAAAAAATGGTGGCTGATCGCAAAATAGTCGATTTAAACAGAGACAGCTGACAACGGTAAAATTGCCTGATGAACTGGGTTTATCAGGCATGAAAACCAACACTTTGTTAGTCATCTGTAAAATAGCGCATCATTAAGGAGTa                                                                                   >  2:735588/1‑150 (MQ=255)
                                                                  gAAGCGTTAAAAAATGGTGGCTGATCGCAAAATAGTCGATTTAAACAGAGACAGCTGACAACGGTAAAATTGCCTGATGAACTGGGTTTATCAGGCATGAAAACCAACACTTTGTTAGTCATCTGTAAAATAGCGCATCATTAAGGAGTa                                                                                   <  1:612874/150‑1 (MQ=255)
                                                                    aGCTTTAAAAAATGGTGGCTGATCGCAAAATAGTCGATTTAAACAGAGACAGCTGACAACGGTAAAATTGCCTGATGAACTGGGTTTATCAGGCATGAAAACCAACACTTTGTTAGTCATCTGTAAAATAGCGCATCATTAAGGAGTAcc                                                                                 <  2:310136/150‑1 (MQ=255)
                                                                      cGTTAAAAAATGGTGGCTGATCGCAAAATAGTCGATTTAAACAGAGACAGCTGACAACGGTAAAATTGCCTGATGAACTGGGTTTATCAGGCATGAAAACCAACACTTTGTTAGTCATCTGTAAAATAGCGCATCATTAAGGAGTACCaa                                                                               <  1:360318/150‑1 (MQ=255)
                                                                       gTTAAAAAATGGTGGCTGATCGCAAAATAGTCGATTTAAACAGAGACAGCTGACAACGGTAAAATTGCCTGATGAACTGGGTTTATCAGGCATGAAAACCAACACTTTGTTAGTCATCTGTAAAATAGCGCATCATTAAGGAGTACCAAt                                                                              >  1:621429/1‑150 (MQ=255)
                                                                        ttAAAAAATGGTGGCTGATCGCAAAATAGTCGATTTAAACAGAGACAGCTGACAACGGTAAAATTGCCTGATGAACTGGGTTTATCAGGCATGAAAACCAACACTTTGTTAGTCATCTGTAAAATAGCGCATCATTAAGGAGTACCAAtg                                                                             <  2:1066774/150‑1 (MQ=255)
                                                                          aaaaaaTGGTGGCTGATCGCAAAATAGTCGATTTAAACAGAGACAGCTGACAACGGTAAAATTGCCTGATGAACTGGGTTTATCAGGCATGAAAACCAACACTTTGTTAGTCATCTGTAAAATAgcg                                                                                                  >  2:1625410/1‑127 (MQ=255)
                                                                          aaaaaaTGGTGGCTGATCGCAAAATAGTCGATTTAAACAGAGACAGCTGACAACGGTAAAATTGCCTGATGAACTGGGTTTATCAGGCATGAAAACCAACACTTTGTTAGTCATCTGTAAAATAgcg                                                                                                  <  1:1625410/127‑1 (MQ=255)
                                                                          aaaaaaTGGTGGCTGATCGCAAAATAGTCGATTTAAACAGAGACAGCTGACAACGGTAAAATTGCCTGATGAACTGGGTTTATCAGGCATGAAAACCAACACTTTGTTAGTCATCTGTAAAATAGCGCATCATTAAGGAGTACCAAtgtg                                                                           <  2:907331/150‑1 (MQ=255)
                                                                           aaaaaTGGTGGCTGATCGCAAAATAGTCGATTTAAACAGAGACAGCTGACAACGGTAAAATTGCCTGATGAACTGGGTTTATCAGGCATGAAAACCAACACTTTGTTAGTCATCTGTAAAATAGCGCATCATTAAGGAGTACCAATGTGc                                                                          >  2:648804/1‑150 (MQ=255)
                                                                            aaaaTGGTGGCTGATCGCAAAATAGTCGATTTAAACAGAGACAGCTGACAACGGTAAAATTGCCTGATGAACTGGGTTTATCAGGCATGAAAACCAACACTTTGTTAGTCATCTGTAAAATAGCGCATCATTAAGGAGTACCAATGTGCg                                                                         <  2:1212759/150‑1 (MQ=255)
                                                                               aTGGTGGCTGATCGCAAAATAGTCGATTTAAACAGAGACAGCTGACAACGGTAAAATTGCCTGATGAACTGGGTTTATCAGGCATGAAAACCAACACTTTGTTAGTCATCTGTAAAATAGCGCATCATTAAGGAGTACCAATGTGCGAAc                                                                      <  2:287310/150‑1 (MQ=255)
                                                                                   tggCTGATCGCAAAATAGTCGATTTAAACAGAGACAGCTGACAACGGTAAAATTGCCTGATGAACTGGGTTTATCAGGCATGAAAACCAACACTTTGTTAGTCATCTGTAAAATAGCGCATCATTAAGGAGTACCAATGTGCGAAcatca                                                                  >  1:98511/1‑150 (MQ=255)
                                                                                   tggCTGATCGCAAAATAGTCGATTTAAACAGAGACAGCTGACAACGGTAAAATTGCCTGATGAACTGGGTTTATCAGGCATGAAAACCAACACTTTGTCAGTCATCTGTAAAATAGCGCATCATTAAGGAGTACCAATGTGCGAAcatca                                                                  <  2:1574315/150‑1 (MQ=255)
                                                                                      cTGATCGCAAAATAGTCGATTTAAACAGAGACAGCTGACAACGGTAAAATTGCCTGATGAACTGGGTTTATCAGGCATGAAAACCAACACTTTGTTAGTCATCTGTAAAATAGCGCATCATTAAGGAGTACCAATGTGCGAAcatcatca                                                               >  2:1038452/1‑150 (MQ=255)
                                                                                      cTGATCGCAAAATAGTCGATTTAAACAGAGACAGCTGACAACGGTAAAATTGCCTGATGAACTGGGTTTATCAGGCATGAAAACCAACACTTTGTTAGTCATCTGTAAAATAGCGCATCATTAAGGAGTACCAATGTGAGAAcatcatca                                                               <  2:104143/150‑1 (MQ=255)
                                                                                         aTCGCAAAATAGTCGATTTAAACAGAGACATCTGACAACTGTAAAATTGCCTGATGAACTGGGTTTATCAGGCATGAAAACCAACACTTTGTTAGTCATCTGTAAAATAGCGCATCATTAAGGAGTACCAATGTGCGAACATCATCATGc                                                            <  1:698614/150‑1 (MQ=255)
                                                                                            gCAAAATAGTCGATTTAAACAGAGACAGCTGACAACGGTAAAATTGCCTGATGAACTGGGTTTATCAGGCATGAAAACCAACACTTTGTTAGTCATCTGTAAAATAGCGCATCATTAAGGAGTACCAATGTGCGAACATCATCATGCcgc                                                         >  1:1380389/1‑150 (MQ=255)
                                                                                               aaaTAGTCGATTTAAACAGAGACAGCTGACAACGGTAAAATTGCCTGATGAACTGGGTTTATCAGGCATGAAAACCAACACTTTGTTAGTCATCTGTAAAATAGCGCATCATTAAGGAGTACCAATGTGCGAACATCATCATGCCGCGaa                                                      >  2:42867/1‑150 (MQ=255)
                                                                                                aaTAGTCGATTTAAACAGAGACAGCTGACAACGGTAAAATTGCCTGATGAACTGGGTTTATCAGGCATGAAAACCAACACTTTGTTAGTCATCTGTAAAATAGCGCATCATTAAGGAGTACCAATGTGCGAACATCATCATGCCGCGAAg                                                     <  2:702544/150‑1 (MQ=255)
                                                                                                  tAGTCGATTTAAACAGAGACAGCTGACAACGGTAAAATTGCCTGATGAACTGGGTTTATCAGGCATGAAAACCAACACTTTGTTAGTCATCTGTAAAATAGCGCATCATTAAGGAGTACCAATGTGCGAACATCATCATGCCGCGAAGct                                                   <  2:1526368/150‑2 (MQ=255)
                                                                                                  tAGTCGATTTAAACAGAGACAGCTGACAACGGTAAAATTGCCTGATGAACTGGGTTTATCAGGCATGAAAACCAACACTTTGTTAGTCATCTGTAAAATAGCGCATCATTAAGGAGTACCAATGTGCGAACATCATCATGCCGCGAAGca                                                   >  2:1142369/1‑150 (MQ=255)
                                                                                                    gTCGATTTAAACAGAGACAGCTGACAACGGTAAAATTGCCTGATGAACTGGGTTTATCAGGCATGAAAACCAACACTTTGTTAGTCATCTGTAAAATAGCGCATCATTAAGGAGTACCAATGTGCGAACATCATCATGCCGCGAAGcaca                                                 <  2:1382685/150‑1 (MQ=255)
                                                                                                       gATTTAAACAGAGACAGCTGACAACGGTAAAATTGCCTGATGAACTGGGTTTATCAGGCATGAAAACCAACACTTTGTTAGTCATCTGTAAAATAGCGCATCATTAAGGAGTACCAATGTGCGAACATCATCATGCCGCGAAGCACATcc                                              >  1:1140684/1‑150 (MQ=255)
                                                                                                         tttAAACAGAGACAGCTGACAACGGTAAAATTGCCTGATGAACTGGGTTTATCAGGCATGAAAACCAACACTTTGTTAGTCATCTGTAAAATAGCGCATCATTAAGGAGTACCAATGTGCGAACATCATCATGCCGCGAAGCACATCCtg                                            <  2:1591834/150‑1 (MQ=255)
                                                                                                           tAAACAGAGACAGCTGACAACGGTAAAATTGCCTGATGAACTGGGTTTATCAGGCATGAAAACCAACACTTTGTTAGTCATCTGTAAAATAGCGCATCATTAAGGAGTACCAATGTGCGAACATCATCATGCCGCGAAGCACATCCtgtg                                          >  2:1241546/1‑150 (MQ=255)
                                                                                                           tAAACAGAGACAGCTGACAACGGTAAAATTGCCTGATGAACTGGGTTTATCAGGCATGAAAACCAACACTTTGTTAGTCATCTGTAAAATAGCGCATCATTAAGGAGTACCAATGTGCGAACATCATCATGCCGCGAAGCACATCCtgtg                                          <  2:212290/150‑1 (MQ=255)
                                                                                                            aaaCAGAGACAGCTGACAACGGTAAAATTGCCTGATGAACTGGGTTTATCAGGCATGAAAACCAACACTTTGTTAGTCATCTGTAAAATAGCGCATCATTAAGGAGTACCAATGTGCGAACATCATCATGCCGCGAAGCACATCCTGTGc                                         >  2:743699/1‑150 (MQ=255)
                                                                                                              aCAGAGACAGCTGACAACGGTAAAATTGCCTGATGAACTGGGTTTATCAGGCATGAAAACCAACACTTTGTTAGTCATCTGTAAAATAGCGCATCATTAAGGAGTACCAATGTGCGAACATCATCATGCCGCGAAGCACATCCTGTGCTc                                       <  1:1546891/150‑1 (MQ=255)
                                                                                                                agagACAGCTGACAACGGTAAAATTGCCTGATGAACTGGGTTTATCAGGCATGAAAACCAACACTTTGTTAGTCATCTGTAAAATAGCGCATCATTAAGGAGTACCAATGTGCGAACATCATCATGCCGCGAAGCACATCCTGTGCTCGc                                     <  1:816227/150‑1 (MQ=255)
                                                                                                                 gagaCAGCTGACAACGGTAAAATTGCCTGATGAACTGGGTTTATCAGGCATGAAAACCAACACTTTGTTAGTCATCTGTAAAATAGCGCATCATTAAGGAGTACCAATGTGCGAACATCATCATGCCGCGAAGCACATCCTGTGCTCGCa                                    <  2:1012334/150‑1 (MQ=255)
                                                                                                                  agaCAGCTGACAACGGTAAAATTGCCTGATGAACTGGGTTTATCAGGCATGAAAACCAACACTTTGTTAGTCATCTGTAAAATAGCGCATCATTAAGGAGTACCAATGTGCGAACATCATCATGCCGCGAAGCACATCCTGTGCTCGCAg                                   <  1:1332872/150‑1 (MQ=255)
                                                                                                                  agaCAGCTGACAACGGTAAAATTGCCTGATGAACTGGGTTTATCAGGCATGAAAACCAACACTTTGTTAGTCATCTGTAAAATAGCGCATCATTAAGGAGTACCAATGTGCGAACATCATCATGCCGCGAAGCACATCCTGTGCTCGCAg                                   <  1:743699/150‑1 (MQ=255)
                                                                                                                     cAGCTGACAACGGTAAAATTGCCTGATGAACTGGGTTTATCAGGCATGAAAACCAACACTTTGTTAGTCATCTGTAAAATAGCGCATCATTAAGGAGTACCAATGTGCGAACATCATCATGCCGCGAAGCACATCCTGTGCTCGCAgtgt                                >  2:1590361/1‑150 (MQ=255)
                                                                                                                       gCTGACAACGGTAAAATTGCCTGATGAACTGGGTTTATCAGGCATGAAAACCAACACTTTGTTAGTCATCTGTAAAATAGCGCATCATTAAGGAGTACCAATGTGCGAACATCATCATGCCGCGAAGCACATCCTGTGCTCGCAGTGTGa                              >  2:128950/1‑150 (MQ=255)
                                                                                                                       gCTGACAACGGTAAAATTGCCTGATGAACTGGGTTTATCAGGCATGAAAACCAACACTTTGTTAGTCATCTGTAAAATAGCGCATCATTAAGGAGTACCAATGTGCGAACATCATCATGCCGCGAAGCACATCCTGTGCTCGCAGTGTGa                              >  1:595179/1‑150 (MQ=255)
                                                                                                                         tGACAACGGTAAAATTGCCTGATGAACTGGGTTTATCAGGCATGAAAACCAACACTTTGTTAGTCATCTGTAAAATAGCGCATCATTAAGGAGTACCAATGTGCGAACATCATCATGCCGCGAAGCACATCCTGTGCTCGCAGTGTGACa                            <  1:1555581/150‑1 (MQ=255)
                                                                                                                          gACAACGGTAAAATTGCCTGATGAACTGGGTTTATCAGGCATGAAAACCAACACTTTGTTAGTCATCTGTAAAATAGCGCATCATTAAGGAGTACCAATGTGCGAACATCATCATGCCGCGAAGCACATCCTGTGCTCGCAGTGTGACAt                           >  1:909972/1‑150 (MQ=255)
                                                                                                                            cAACGGTAAAATTGCCTGATGAACTGGGTTTATCAGGCATGAAAACCAACACTTTGTTAGTCATCTGTAAAATAGCGCATCATTAAGGAGTACCAATGTGCGAACATCATCATGCCGCGAAGCACATCCTGTGCTCGCAGTGTGACATGc                         >  1:535157/1‑150 (MQ=255)
                                                                                                                               cGGTAAAATTGCCTGATGAACTGGGTTTATCAGGCATGAAAACCAACACTTTGTTAGTCATCTGTAAAATAGCGCATCATTAAGGAGTACCAATGTGCGAACATCATCATGCCGCGAAGCACATCCTGTGCTCGCAGTGTGACATGCtgg                      >  1:1063761/1‑150 (MQ=255)
                                                                                                                               cGGTAAAATTGCCTGATGAACTGGGTTTATCAGGCATGAAAACCAACACTTTGTTAGTCATCTGTAAAATAGCGCATCATTAAGGAGTACCAATGTGCGAACATCATCATGCCGCGAAGCACATCCTGTGCTCGCAGTGTGACATGCtgg                      >  2:808288/1‑150 (MQ=255)
                                                                                                                               cGGTAAAATTGCCTGATGAACTGGGTTTATCAGGCATGAAAACCAACACTTTGTTAGTCATCTGTAAAATAGCGCATCATTAAGGAGTACCAATGTGCGAACAGCATCATGCCGCGAAGCACATCCTGTGCTCGCAGTGTGACATGCtgg                      >  1:665282/1‑150 (MQ=255)
                                                                                                                                 gttttcatGCCTGATAAACCCAGTTCATCAGGCATGAAAACCAACACTTTGTTAGTCATCTGTAAAATAGCGCATCATTAAGGAGTACCAATGTGCGAACATCATCATGCCGCGAAGCACATCCTGTGCTCGCAGTGTGACATGCtggtg                    >  1:897851/8‑150 (MQ=255)
                                                                                                                                  ttttcatGCCTGATAAACCCAGTTCATCAGGCATGAAAACCAACACTTTGTTAGTCATCTGTAAAATAGCGCATCATTAAGGAGTACCAATGTGCGAACATCATCATGCCGCGAAGCACATCCTGTGCTCGCAGTGTGACATGCtggtgg                   >  2:1578638/7‑150 (MQ=255)
                                                                                                                                  tAAAATTGCCTGATGAACTGGGTTTATCAGGCATGAAAACCAACACTTTGTTAGTCATCTGTAAAATAGCGCATCATTAAGGAGTACCAATGTGCGAACATCATCATGCCGCGAAGCACATCCTGTGCTCGCAGTGTGACATGCtggtgg                   <  1:1154003/150‑1 (MQ=255)
                                                                                                                                   aaaaTTGCCTGATGAACTGGGTTTATCAGGCATGAAAACCAACACTTTGTTAGTCATCTGTAAAATAGCGCATCATTAAGGAGTACCAATGTGCGAACATCATCATGCCGCGAAGCACATCCTGTGCTCGCAGTGTGACATGCTGGTGGc                  <  2:837618/150‑1 (MQ=255)
                                                                                                                                      catGCCTGATAAACCCAGTTCATCAGGCATGAAAACCAACACTTTGTTAGTCATCTGTAAAATAGCGCATCATTAAGGAGTACCAATGTGCGAACATCATCATGCCGCGAAGCACATCCTGTGCTCGCAGTGTGACATGCTGGTGGCGtt               >  2:401368/3‑150 (MQ=255)
                                                                                                                                      aTTGCCTGATGAACTGGGTTTATCAGGCATGAAAACCAACACTTTGTTAGTCATCTGTAAAATAGCGCATCATTAAGGAGTACCAATGTGCGAACATCATCATGCCGCGAAGCACATCCTGTGCTCGCAGTGTGACATGCTGGTGGCGtt               >  1:775418/1‑150 (MQ=255)
                                                                                                                                       ttGCCTGATGAACTGGGTTTATCAGGCATGAAAACCAACACTTTGTTAGTCATCTGTAAAATAGCGCATCATTAAGGAGTACCAATGTGCGAACATCATCATGCCGCGAAGCACATCCTGTGCTCGCAGTGTGACATGCTGGTGGCGTTa              >  1:1058490/1‑150 (MQ=255)
                                                                                                                                        tGCCTGATAAACCCAGTTCATCAGGCATGAAAACCAACACTTTGTTAGTCATCTGTAAAATAGCGCATCATTAAGGAGTACCAATGTGCGAACATCATCATGCCGCGAAGCACATCCTGTGCTCGCAGTGTGACATGCTGGTGGCGTTAc             >  1:174781/1‑150 (MQ=255)
                                                                                                                                           cTGATGAACTGGGTTTATCAGGCATGAAAACCAACACTTTGTTAGTCATCTGTAAAATAGCGCATCATTAAGGAGTACCAATGTGCGAACATCATCATGACGCGAAGCACATCCTGTGCTCGCAGTGTGACATGCTGGTGGCGTTACcgc          >  1:1449190/1‑150 (MQ=255)
                                                                                                                                            tgatgaACTGGGTTTATCATGCAAGAAAACAAAAACTTTGTTAGTCATCTGTAAAATAGAGCATCATTAAGGAGTAACAATGTGAGAACATCATCATGCCGCGAAGCACATCCTGTGCTCGCAGTGTGACATGCTGGTGGCGTTACcgcg         >  2:1343662/1‑150 (MQ=255)
                                                                                                                                            tgatgaACTGGGTTTATCAGGCATGAAAACCAACACTTTGTTAGTCATCTGTAAAATAGCGCATCATTAAGGAGTACCAATGTGCGAACATCATCATGCCGCGAAGCACATCCTGTGCTCGCAGTGTGACATGCTGGTGGCGTTACcgcg         <  2:1354848/150‑1 (MQ=255)
                                                                                                                                             gatgaACTGGGTTTATCAGGCATGAAAACCAACACTTTGTTAGTCAGCGGTATAATAGCGCATCATTAAGGAGTACCAAGGTGCGAACATCATCATGCCGCGACGCACATCCTGTGCTCGCAGTGTGACATGCTGGTGGCGTTGCCgcgc        >  2:947038/1‑150 (MQ=255)
                                                                                                                                              atgaACTGGGTTTATCAGGCATGAAAACCAACACTTTGTTAGTCATCTGTAAAATAGCGCATCATTAAGGAGTACCAATGTGCGAACATCATCATGCCGCGAAGCACATCCTGTGCTCGCAGTGTGACATGCTGGTGGCGTTACCGCGcc       >  1:838218/1‑150 (MQ=255)
                                                                                                                                                   cTGGGTTTATCAGGCATGAAAACCAACACTTTGTTAGTCATCTGTAAAATAGCGCATCATTAAGGAGTACCAATGTGCGAACATCATCATGCCGCGAAGCACATCCTGTGCTCGCAGTGTGACATGCTGGTGGCGTTACCGCGCCTTGAg  <  1:1139467/150‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                    
GTGCTTGCTGATATGGCTGCTGCAGAGCAGGAGCTGGAGCAAGCCTTTGAACGCTGGGAGTATCTTGAAGCGTTAAAAAATGGTGGCTGATCGCAAAATAGTCGATTTAAACAGAGACAGCTGACAACGGTAAAATTGCCTGATGAACTGGGTTTATCAGGCATGAAAACCAACACTTTGTTAGTCATCTGTAAAATAGCGCATCATTAAGGAGTACCAATGTGCGAACATCATCATGCCGCGAAGCACATCCTGTGCTCGCAGTGTGACATGCTGGTGGCGTTACCGCGCCTTGAG  >  NC_000913/1011782‑1012078

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