Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,012,854 G→A 5.7% T54M (ACG→ATG) ‡ fliE ← flagellar basal‑body protein FliE

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,012,8540GA 5.7% 276.3 / 3.1 87T54M (ACG→ATG) ‡fliEflagellar basal‑body protein FliE
Reads supporting (aligned to +/- strand):  ref base G (48/34);  new base A (2/3);  total (50/37)
Fisher's exact test for biased strand distribution p-value = 6.47e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GCATGCTCATCACTTCCTGATACGCCGCCACCAGCTTATTACGCACCTGAATCCCCATTTGCATAGAAACTGAGGCTTTTTGCATATCGGTCATCACATCGTTTAACGCCACGCCGGGTTCACCGAGAGTGAATTTTTCTGCCTGCGTGCGGGCAGCTGTTTGTGTATCACTAATGCGATCGAGCGCGGCGTGCAGCTGCCCGGCAAAACTAATGGTCGGTTGCGGCAGTGATTCCTGCGCACGCGCACTCATCGCCGTAGCCTGTAACTGGCTGATAACCCCTTCAATCCCCTG  >  NC_000913/2012708‑2013002
                                                                                                                                                  |                                                                                                                                                    
gCATGCTCATCACTTCCTGATACGCCGCCACCAGCTTATCACGCACCTGAATCCCCATTTGCATAGAAACTGAGGCTTTTTGCATATCGGTCATCACATCGTTTAACGCCACGCCGGGTTCACCGAGAGTGAATTTTTCTGCCtgcgtgc                                                                                                                                                   >  2:1157790/1‑150 (MQ=255)
 cATGCTCATCACTTCCTGATACGCCGCCACCAGCTTATTACGCACCTGAATCCCCATTTGCATAGAAACTGAGGCTTTTTGCATATCGGTCATCACATCGTTTAACGCCACGCCGGGTTCACCGAGAGTGAATTTTTCTGCCtgcgtgcg                                                                                                                                                  >  1:259201/1‑150 (MQ=255)
       catcaCTTCCTGATACGCCGCCACCAGCTTATTACGCACCTGAATCCCCATTTGCATAGAAACTGAGGCTTTTTGCATATCGGTCATCACATCGTTTAACGCCACGCCGGGTTCACCGAGAGTGAATTTTTCTGCCTGCGTGCGGGCAGc                                                                                                                                            <  2:920385/150‑1 (MQ=255)
        atcaCTTCCTGATACGCCGCCACCAGCTTATTACGCACCTGAATCCCCATTTGCATAGAAACTGAGGCTTTTTGCATATCGGTCATCACATCGTTTAACGCCACGCCGGGTTCACCGAGAGTGAATTTTTCTGCCTGCGTGCGGGCAGCt                                                                                                                                           >  2:1101338/1‑150 (MQ=255)
         tcaCTTCCTGATACGCCGCCACCAGCTTATTACGCACCTGAATCCCCATTTGCATAGAAACTGAGGCTTTTTGCATATCGGTCATCACATCGTTTAACGCCACGCCGGGTTCACCGAGAGTGAATTTTTCTGCCCGCACgcaggcagaaa                                                                                                                                          <  2:1402256/150‑4 (MQ=255)
           aCTTCCTGATACGCCGCCACCAGCTTATTACGCACCTGAATCCCCATTTGCATAGAAACTGAGGCTTTTTGCATATCGGTCATCACATCGTTTAACGCCACGCCGGGTTCACCGAGAGTGAATTTTTCTGCCTGCGTGCGGGCAGCTGtt                                                                                                                                        >  2:1276126/1‑150 (MQ=255)
            cTTCCTGATACGCCGCCACCAGCTTATTACGCACCTGAATCCCCATTTGCATAGAAACTGAGGCTTTTTGCATATCGGTCATCACATCGTTTAACGCCACGCCGGGTTCACCGAGAGTGAATTTTTCTGCCCGCACgcaggcagaaaaat                                                                                                                                       <  1:423290/150‑7 (MQ=255)
             ttCCTGATACGCCGCCACCAGCTTATTACGCACCTGAATCCCCATTTGCATAGAAACTGAGGCTTTTTGCATATGGGTCATCACATCGTTTAACGCCACGCCGTGTTCACCGAGAGTGAATTTTTCTGCCCGCACgcaggcagaaaaatt                                                                                                                                      <  2:261582/150‑8 (MQ=255)
             ttCCTGATACGCCGCCACCAGCTTATTACGCACCTGAATCCCCATTTGCATAGAAACTGAGGCTTTTTGCATATCGGTCATCACATCGTTTAACGCCACGCCGGGTTCACCGAGAGTGAATTTTTCTGCCTGCGTGCGGGCAGCTGTTtg                                                                                                                                      <  1:107484/150‑1 (MQ=255)
             ttCCTGATACGCCGCCACCAGCTTATTACGCACCTGAATCCCCATTTGCATAGAAACTGAGGCTTTTTGCATATCGGTCATCACATCGTTTAACGCCACGCCGGGTTCACCGAGAGTGAATTTTTCTGCCTGCGTGCGGGCAGCTGTTtg                                                                                                                                      <  2:868485/150‑1 (MQ=255)
             ttCCTGATACGCCGCCACCAGCTTATTACGCACCTGAATCCCCATTTGCATAGAAACTGAGGCTTTTTGCATATCGGTCATCACATCGTTTAACGCCACGCCGGGTTCACCGAGAGTGAATTTTTCTGCCTGCGTGCGGGCAGCTGTTtg                                                                                                                                      <  1:540135/150‑1 (MQ=255)
             taccTGATACGCCGCCACCAGCTTTTTACGCACCTGAATCCCCATTTGCATAGAAACTGAGGCTTTTTGCATATCGGTCATCACATCGTTTAACGCCACGCGGGGTTCACGGAGAGTGAATTTTTCTGCCTGCGTGCGGGCAGCTGTTtg                                                                                                                                      <  1:856446/148‑1 (MQ=255)
               ccTGATACGCCGCCACCAGCTTATTACGCACCTGAATCCCCATTTGCATAGAAACTGAGGCTTTTTGCATATCGGTCATCACATCGTTTAACGCCACGCCGGGTTCACCGAGAGTGAATTTTTCTGCCTGCGTGCGGGCAGCTGTTtgtg                                                                                                                                    <  1:737599/150‑1 (MQ=255)
                cTGATACGCCGCCACCAGCTTATTACGCACCTGAATCCCCATTTGCATAGAAACTGAGGCTTTTTGCATATCGGTCATCACATCGTTTAACGCCACGCCGGGTTCACCGAGAGTGAATTTTTCTGCCTGCGTGCGGGCAGCTGTTTgtgt                                                                                                                                   >  2:1168418/1‑150 (MQ=255)
                 tGATACGCCGCCACCAGCTTATTACGCACCTGAATCCCCATTTGCATAGAAACTGAGGCTTTTTGCATATCGGTCATCACATCGTTTAACGCCACGCCGGGTTCACCGAGAGTGAATTTTTCTGCCTGCGTGCGGGCAGCTGTTTGTGTa                                                                                                                                  >  1:1495747/1‑150 (MQ=255)
                 tGATACGCCGCCACCAGCTTATTACGCAACTGAATACCCATTTGCATAGAAACTGAGGCTTTTTGCATATCGGTCATCACATCGTTTAACTCCACGCCGGGTTCACCGAGAGTCAATTTTTCTGCCTGCGTGCGGGCAGCTGTTTGTGTa                                                                                                                                  >  2:241518/1‑150 (MQ=255)
                  gATACGCCGCCACCAGCTTATTACGCACCTGAATCCCCATTTGCATAGAAACTGAGGCTTTTTGCATATCGGTCATCACATCGTTTAACGCCACGCCGGGTTCACCGAGAGTGAATTTTTCTGCCTGCGTGCGGGCAGCTGTTTGTGTAt                                                                                                                                 >  1:702235/1‑150 (MQ=255)
                   aTACGCCGCCACCAGCTTATTACGCACCTGAATCCCCATTTGCATAGAAACTGAGGCTTTTTGCATATCGGTCATCACATCGTTAAACGCCACGCCGGGTTCACCGAGAGTGAATTTTTCTGCCTGCGTGCGGGCAGCTGTTTGTGTATc                                                                                                                                <  1:1182138/150‑1 (MQ=255)
                    tACGCCGCCACCAGCTTATTACGCACCTGAATCCCCATTTGCATAGAAACTGAGGCTTTTTGCATATCGGTCATCACATCGTTTAACGCCACGCCGGGTTCACCGAGAGTGAATTTTTCTGCCTGCGTGCGGGCAGCTGTTTGTGTATCa                                                                                                                               >  1:464924/1‑150 (MQ=255)
                      cgccgcCACCAGCTTATTACGCACCTGAATCCCCATTTGCATAGAAACTGAGGCTTTTTGCATATCGGTCATCACATCGTTTAACGCCACGCCGGGTTCACCGAGAGTGAATTTTTCTGCCTGCGTGCGGGCAGCTGTTTGTGTATCACt                                                                                                                             >  1:747052/1‑150 (MQ=255)
                            caccaGCTTATTACGCACCTGAATCCCCATTTGCATAGAAACTGAGGCTTTTTGCATATCGGTCATCACATCGTTTAACGCCACGCCGGGTTCACCGAGAGTGAATTTTTCTGCCTGCGTGCGGGCAGCTGTTTGTGTATCACTAATGCg                                                                                                                       >  1:923744/1‑150 (MQ=255)
                              ccaGCTTATTACGCACCTGAATCCCCATTTGCATAGAAACTGAGGCTTTTTGCATATCGGTCATCACATCGTTTAACGCCACGCCGGGTTCACCGAGAGTGAATTTTTCTGCCTGCGTGCGGGCAGCTGTTTGTGTATCACTAATGCGAt                                                                                                                     >  2:1508762/1‑150 (MQ=255)
                                    tattaCGCACCTGAATCCCCATTTGCATAGAAACTGAGGCTTTTTGCATATCGGTCATCACATCGTTTAACGCCACGCCGGGTTCACCGAGAGTGAATTTTTCTGCCTGCGTGCGGGCAGCTGTTTGTGTATCACTAATGCGATCGAgcg                                                                                                               >  1:265563/1‑150 (MQ=255)
                                    tattaCGCACCTGAATCCCCATTTGCATAGAAACTGAGGCTTTTTGCATATCGGTCATCACATCGTTTAACGACACGCCGGGTTCACCGAGAGTGAATTTTTCTGCCTGCGTGCGGGCAGCTGTTTGTGTATCACTAATGCGATCGAgcg                                                                                                               >  1:500520/1‑150 (MQ=255)
                                       taCGCACCTGAATCCCCATTTGCATAGAAACTGAGGCTTTTTGCATATCGGTCATCACATCGTTTAACGCCACGCCGGGTTCACCGAGAGTGAATTTTTCTGCCTGCGTGCGGGCAGCTGTTTGTGTATCACTAATGCGATCGAGCgcgg                                                                                                            >  1:1251323/1‑150 (MQ=255)
                                       taCGCACCTGAATCCCCATTTGCATAGAAACTGAGGCTTTTTGCATATCGGTCATCACATCGTTTAACGCCACGCCGGGTTCACCGAGAGTGAATTTTTCTGCCTGCGTGCGGGCAGCTGTTTGTGTATCACTAATGCGATCGAGCgcgg                                                                                                            >  1:668974/1‑150 (MQ=255)
                                            aCCTGAATCCCCATTTGCATAGAAACTGAGGCTTTTTGCATATCGGTCATCACATCGTTTAACGCCACGCCGGGTTCACCGAGAGTGAATTTTTCTGCCTGCGTGCGGGCAGCTGTTTGTGTATCACTAATGCGATCGAGCGCGGCGTGc                                                                                                       <  2:1513719/150‑1 (MQ=255)
                                               tGAATCCCCATTTGCATAGAAACTGAGGCTTTTTGCATATCGGTCATCACATCGTTTAACGCCACGCCGGGTTCACCGAGAGTGAATTTTTCTGCCTGCGTGCGGGCAGCTGTTTGTGTATCACTAATGCGATCGAGCGCGGCGTGCAGc                                                                                                    >  2:1139927/1‑150 (MQ=255)
                                                gAATCCCCATTTGCATAGAAACTGAGGCTTTTTGCATATCGGTCATCACATCGTTTAACGCCACGCCGGGTTCACCGAGAGTGAATTTTTCTGCCTGCGTGCGGGCAGCTGTTTGTGTATCACTAATGCGATCGAGCGCGGCGTGCAGCt                                                                                                   >  1:498033/1‑150 (MQ=255)
                                                    ccccATTTGCATAGAAACTGAGGCTTTTTGCATATCGGTCATCACATCGTTTAACGCCACGCCGGGTTCACCGAGAGTGAATTTTTCTGCCTGCGTGCGGGCAGCTGTTTGTGTATCACTAATGCGATCGAGCGCGGCGTGCAGCTGccc                                                                                               >  1:20957/1‑150 (MQ=255)
                                                     cccATTTGCATAGAAACTGAGGCTTTTTGCATATCGGTCATCACATCGTTTAACGCCACGCCGGGTTCACCGAGAGTGAATTTTTCTGCCTGCGTGCGGGCAGCTGTTTGTGTATCACTAATGCGATCGAGCGCGGCGTGCAGCTGCCCg                                                                                              >  2:1419379/1‑150 (MQ=255)
                                                     cccATTTGCATAGAAACTGAGGCTTTTTGCATATCGGTCATCACATCGTTTAACGCCACGCCGGGTTCACCGAGAGTGAATTTTTCTGCCTGCGTGCGGGCAGCTGTTTGTGTATCACTAATGCGATCGAGCGCGGCGTGCAGCTGCCCg                                                                                              >  2:457627/1‑150 (MQ=255)
                                                       cATTTGCATAGAAACTGAGGCTTTTTGCATATCGGTCATCACATCGTTTAACGCCACGCCGGGTTCACCGAGAGTGAATTTTTCTGCCTGCGTGCGGGCAGCTGTTTGTGTATCACTAATGCGATCGAGCGCGGCGTGCAGCTGCCCGGc                                                                                            <  2:265563/150‑1 (MQ=255)
                                                        aTTTGCATAGAAACTGAGGCTTTTTGCATATCGGTCATCACATCGTTTAACGCCACGCCGGGTTCACCGAGAGTGAATTTTTCTGCCTGCGTGCGGGCAGCTGTTTGTGTATCACTAATGCGATCGAGCGCGGCGTGCAGCTGCCCGGCa                                                                                           >  1:966093/1‑150 (MQ=255)
                                                            gCATAGAAACTGAGGCTTTTTGCATATCGGTCATCACATCGTTTAACGCCACGCCGGGTTCACCGAGAGTGAATTTTTCTGCCTGCGTGCGGGCAGCTGTTTGTGTATCACTAATGCGATCGAGCGCGGCGTGCAGCTGCCCGGCAAAAc                                                                                       >  2:238797/1‑150 (MQ=255)
                                                              aTAGAAACTGAGGCTTTTTGCATATCGGTCATCACATCGTTTAACGCCACGCCGGGTTCACCGAGAGTGAATTTTTCTGCCTGCGTGCGGGCAGCTGTTTGTGTATCACTAATGCGATCGAGCGCGGCGTGCAGCTGCCCGGCAAAACTa                                                                                     >  2:923192/1‑150 (MQ=255)
                                                                aGAAACTGAGGCTTTTTGCATATCGGTCATCACATCGTTTAACGCCACGCCGGGTTCACCGAGAGTGAATTTTTCTGCCTGCGTGCGGGCAGCTGTTTGTGTATCACTAATGCGATCGAGCGCGGCGTGCAGCTGCCCGGCAAAACTAAt                                                                                   <  1:1412780/150‑1 (MQ=255)
                                                                 gAAACTGAGGCTTTTTGCATATCGGTCATCACATCGTTTAACGCCACGCCGGGTTCACCGAGAGTGAATTTTTCTGCCTGCGTGCGGGCAGCTGTTTGTGTATCACTAATGCGATCGCGCGCGGCGTGCAGCTGCCCGGCAAAACTAATg                                                                                  >  2:148497/1‑150 (MQ=255)
                                                                     cTGAGGCTTTTTGCATATCGGTCATCACATCGTTTAACGCCACGCCGGGTTCACCGAGAGTGAATTTTTCTGCCTGCGTGCGGGCAGCTGTTTGTGTATCACTAATGCGATCGAGCGCGGCGTGCAGCTGCCCGGCAAAACTAATGGTCg                                                                              >  1:1324916/1‑150 (MQ=255)
                                                                     atGAGGCTTTTTGCATATCGGTCATCACATCGTTTAACGCCACGCCGGGTTCACCGAGAGTGAATTTTTCTGCCTGCGTGCGGGCAGCTGTTTGTTTATCACTAATGCTATCGAGCGCGGCGTGCAGCTGCCCGGGAAAACTAATGGTCg                                                                              <  2:689078/149‑1 (MQ=255)
                                                                           cTTTTTGCATATCGGTCATCACATCGTTTAACGCCACGCCGGGTTCACCGAGAGTGAATTTTTCTGCCTGCGTGCGGGCAGCTGTTTGTGTATCACTAATGCGATCGAGCGCGGCGTGCAGCTGCCCGGCAAAACTAATGGTCGGTTGCg                                                                        <  2:464924/150‑1 (MQ=255)
                                                                            tttttGCATATCGGTCATCACATCGTTTAACGCCACGCCGGGTTCACCGAGAGTGAATTTTTCTGCCTGCGTGCGGGCAGCTGTTTGTGTATCACTAATGCGATCGAGCGCGGCGTGCAGCTGCCCGGCAAAACTAATGGTCGGTTGCgg                                                                       >  2:1331695/1‑150 (MQ=255)
                                                                               ttGCATATCGGTCATCACATCGTTTAACGCCACGCCGGGTTCACCGAGAGTGAATTTTTCTGCCTGCGTGCGGGCAGCTGTTTGTGTATCACTAATGCGATCGAGCGCGGCGTGCAGCTGCCCGGCAAAACTAATGGTCGGTTGCGGCAg                                                                    <  2:1495747/150‑1 (MQ=255)
                                                                                tGCATATCGGTCATCACATCGTTTAACGCCACGCCGGGTTCACCGAGAGTGAATTTTTCTGCCTGCGTGCGGGCAGCTGTTTGTGTATCACTAATGCGATCGAGCGCGGCGTGCAGCTGCCCGGCAAAACTAATGGTCGGTTGCGGCAGt                                                                   <  1:375052/150‑1 (MQ=255)
                                                                                  cATATCGGTCATCACATCGTTTAACGCCACGCCGGGTTCACCGAGAGTGAATTTTTCTGCCTGCGTGCGGGCAGCTGTTTGTGTATCACTAATGCGATCGAGCGCGGCGTGCAGCTGCCCGGCAAAACTAATGGTCGGTTGCGGCAGTGa                                                                 >  1:212882/1‑150 (MQ=255)
                                                                                   atatCGGTCATCACATCGTTTAACGCCACGCCGGGTTCACCGAGAGTGAATTTTTCTGCCTGCGTGCGGGCAGCTGTTTGTGTATCACTAATGCGATCGAGCGCGGCGTGCAGCTGCCCGGCAAAACTAATGGTCGGTTGCGGCAGTGAt                                                                <  1:715050/150‑1 (MQ=255)
                                                                                   atatCGGTCATCACATCGTTTAACGCCACGCCGGGTTCACCGAGAGTGAATTTTTCTGCCTGCGTGCGGGCAGCTGTTTGTGTATCACTAATGCGATCGAGCGCGGCGTGCAGCTGCCCGGCAAAACTAATGGTCGGTTGCGGCAGTGAt                                                                <  2:593321/150‑1 (MQ=255)
                                                                                    tatCGGTCATCACATCGTTTAACGCCACGCCGGGTTCACCGAGAGTGAATTTTTCTGCCTGCGTGCGGGCAGCTGTTTGTGTATCACTAATGCGATCGAGCGCGGCGTGCAGCTGCCCGGCAAAACTAATGGTCGGTTGCGGCAGTGAtt                                                               >  2:548203/1‑150 (MQ=255)
                                                                                        ggTCATCACATCGTTTAACGCCACGCCGGGTTCACCGAGAGTGAATTTTTCTGCCTGCGTGCGGGCAGCTGTTTGTGTATCACTAATGCGATCGAGCGCGGCGTGCAGCTGCCCGGCAAAACTAATGGTc                                                                               <  2:686457/130‑1 (MQ=255)
                                                                                        ggTCATCACATCGTTTAACGCCACGCCGGGTTCACCGAGAGTGAATTTTTCTGCCTGCGTGCGGGCAGCTGTTTGTGTATCACTAATGCGATCGAGCGCGGCGTGCAGCTGCCCGGCAAAACTAATGGTc                                                                               >  1:686457/1‑130 (MQ=255)
                                                                                        ggTCATCACATCGTTTAACGCCACGCCGGGTTCACCGAGAGTGAATTTTTCTGCCTGCGTGCGGGCAGCTGTTTGTGTATCACTAATGCGATCGAGCGCGGCGTGCAGCTGCCCGGCAAAACTAATGGTCGGTTGCGGCAGTGATTCCTg                                                           <  2:668974/150‑1 (MQ=255)
                                                                                             tcaCATCGTTTAACGCCACGCCGGGTTCACCGAGAGTGAATTTTTCTGCCTGCGTGCGGGCAGCTGTTTGTGTATCACTAATGCGATCGAGCGCGGCGTGCAGCTGCCCGGCAAAACTAATGGTCGGTTGCGGCAGTGATTCCTgcgcac                                                      >  2:916662/1‑150 (MQ=255)
                                                                                               acaTCGTTTAACGCCACGCCGGGTTCACCGAGAGTGAATTTTTCTGCCTGCGTGCGGGCAGCTGTTTGTGTATCACTAATGCGATCGAGCGCGGCGTGCAGCTGCCCGGCAAAACTAATGGTCGGTTGCGGCAGTGATTCCTgcgcacgc                                                    >  2:1075914/1‑150 (MQ=255)
                                                                                               acaTCGTTTAACGCCACGCCGGGTTCACCGAGAGTGAATTTTTCTGCCTGCGTGCGGGCAGCTGTTTGTGTATCACTAATGCGATCGAGCGCGGCGTGCAGCTGCCCGGCAAAACTAATGGTCGGTTGCGGCAGTGATTCCTgcgcacgc                                                    >  2:243885/1‑150 (MQ=255)
                                                                                                caTCGTTTAACGCCACGCCGGGTTCACCGAGAGTGAATTTTTCTGCCTGCGTGCGGGCAGCTGTTTGTGTATCACTAATGCGATCGAGCGCGGCGTGCAGCTGCCCGGCAAAACTAATGGTCGGTTGCGGCAGTGATTCCTgcgcacgcg                                                   >  2:680008/1‑150 (MQ=255)
                                                                                                caTCGTTTAACGCCACGCCGGGTTCACCGAGAGTGAATTTTTCTGCCTGCGTGCGGGCAGCTGTTTGTGTATCACTAATGCGATCGAGCGCGGCGTGCAGCTGCCCGGCAAAACTAATGGTCGGTTGCGGCAGTGATTCCTgcgcacgcg                                                   >  1:719441/1‑150 (MQ=255)
                                                                                                 aTCGTTTAACGCCACGCCGGGTTCACCGAGAGTGAATTTTTCTGCCTGCGTGCGGGCAGCTGTTTGTGTATCACTAATGCGATCGAGCGCGGCGTGCAGCTTCCCGGCAAAACTAATGGTCGGTTGCGGCAGTGATTCCTgcgcacgcgc                                                  >  1:293122/1‑150 (MQ=255)
                                                                                                   cGTTTAACGCCACGCCGGGTTCACCGAGAGTGAATTTTTCTGCCTGCGTGCGGGCAGCTGTTTGTGTATCACTAATGCGATCGAGCGCGGCGTGCAGCTGCCCGGCAAAACTAATGGTCGGTTGCGGCAGTGATTCCTgcgcacgcgcac                                                >  1:535217/1‑150 (MQ=255)
                                                                                                    gTTTAACGCCACGCCGGGTTCACCGAGAGTGAATTTTTCTGCCTGCGTGCGGGCAGCTGTTTGTGTATCACTAATGCGATCGAGCGCGGCGTGCAGCTGCCCGGCAAAACTAATGGTCGGTTGCGGCAGTGATTCCTGCGCACGCGCACt                                               >  1:1292273/1‑150 (MQ=255)
                                                                                                       tAACGCCACGCCGGGTTCACCGAGAGTGAATTTTTCTGCCTGCGTGCGGGCAGCTGTTTGTGTATCACTAATGCGATCGAGCGCGGCGTGCAGCTGCCCGGCAAAACTAATGGTCGGTTGCGGCAGTGATTCCTGCGCACGCGCACTCAt                                            <  1:1445441/150‑1 (MQ=255)
                                                                                                       tAACGCCACGCCGGGTTCACCGAGAGTGAATTTTTCTGCCTGCGTGCGGGCAGCTGTTTGTGTATCACTAATGCGATCGAGCGCGGCGTGCAGCTGCCCGGCAAAACTAATGGTCGGTTGCGGCAGTGATTCCTGCGCACGCGCACTCAt                                            <  1:1445444/150‑1 (MQ=255)
                                                                                                           gccacgccGGGTTCACCGAGAGTGAATTTTTCTGCCTGCGTGCGGGCAGCTGTTTGTGTATCACTAATGCGATCGAGCGCGGCGTGCAGCTGCCCGGCAAAACGAATGGTCGGTTGCGGCAGTGATTCCTGCGCACGCGCACTCATCGcc                                        >  2:306079/1‑150 (MQ=255)
                                                                                                              acgccGGGTTCACCGAGAGTGAATTTTTCTGCCTGCGTGCGGGCAGCTGTTTGTGTATCACTAATGCGATCGAGCGCGGCGTGCAGCTGCCCGGCAAAACTAATGGTCGGTTGCGGCAGTGATTCCTGCGCACGCGCACTCATCGCCGTa                                     >  1:1159146/1‑150 (MQ=255)
                                                                                                               cgccGGGTTCACCGAGAGTGAATTTTTCTGCCTGCGTGCGGGCAGCTGTTTGTGTATCACTAATGCGATCGAGCGCGGCGTGCAGCTGCCCGGCAAAACTAATGGTCGGTTGCGGCAGTGATTCCTGCGCACGCGCACTCATCGCCGTAg                                    <  1:1599300/150‑1 (MQ=255)
                                                                                                                gccGGGTTCACCGAGAGTGAATTTTTCTGCCTGCGTGCGGGCAGCTGTTTGTGTATCACTAATGCGATCGAGCGCGGCGTGCAGCTGCCCGGCAAAACTAATGGTCGGTTGCGGCAGTGATTCCTGCGCACGCGCACTCATCGCCGTAGc                                   <  1:457627/150‑1 (MQ=255)
                                                                                                                 ccGGGTTCACCGAGAGTGAATTTTTCTGCCTGCGTGCGGGCAGCTGTTTGTGTATCACTAATGCGATCGAGCGCGGCGTGCAGCTGCCCGGCAAAACTAATGGTCGGTTGCGGCAGTGATTCCTGCGCACGCGCACTCATCGCCGTAGcc                                  <  2:212882/150‑1 (MQ=255)
                                                                                                                 ccGGGTTCACCGAGAGTGAATTTTTCTGCCTGCGTGCGGGCAGCTGTTTGTGTATCACTAATGCGATCGAGCGCGGCGTGCAGCTGCCCGGCAAAACTAATGGTCGGTTGCGGCAGTGATTCCTGCGCACGCGCACTCATCGCCGTAGcc                                  <  1:1276126/150‑1 (MQ=255)
                                                                                                                 ccGGGTTCACCGAGAGTGAATTTTTCTGCCTGCGTGCGGGCAGCTGTTTGTGTATCACTAATGCGATCGAGCGCGGCGTGCAGCTGCCCGGCAAAACTAATGGTCGGTTGCGGCAGTGATTCCTGCGCACGCGCACTCATCGCCGTAGcc                                  >  1:1520608/1‑150 (MQ=255)
                                                                                                                   gggTTCACCGAGAGTGAATTTTTCTGCCTGCGTGCGGGCAGCTGTTTGTGTATCACTAATGCGATCGAGCGCGGCGTGCAGCTGCCCGGCAAAACTAATGGTCGGTTGCGGCAGTGATTCCTGCGCACGCGCACTCATCGCCGTAGCCTg                                <  2:966093/150‑1 (MQ=255)
                                                                                                                        cACCGAGAGTGAATTTTTCTGCCTGCGTGCGGGCAGCTGTTTGTGTATCACTAATGCGATCGAGCGCGGCGTGCAGCTGCCCGGCAAAACTAATGGTCGGTTGCGGCAGTGATTCCTGCGCACGCGCACTCATCGCCGTAGCCTGTAACt                           >  1:1344634/1‑150 (MQ=255)
                                                                                                                        cACCGAGAGTGAATTTTTCTGCCTGCGTGCGGGCAGCTGTTTGTGTATCACTAATGCGATCGAGCGCGGCGTGCAGCTGCCCGGCAAAACTAATGGTCGGTTGCGGCAGTGATTCCTGCGCACGCGCACTCATCGCCGTAGCCTGTAACt                           >  1:1376925/1‑150 (MQ=255)
                                                                                                                        cACCGAGAGTGAATTTTTCTGCCTGCGTGCGGGCAGCTGTTTGTGTATCACTAATGCGATCGAGCGCGGCGTGCAGCTGCCCGGCAAAACTAATGGTCGGTTGCGGCAGTGATTCCTGCGCACGCGCACTCATCGCCGTAGCCTGTAACt                           >  1:490637/1‑150 (MQ=255)
                                                                                                                        cACCGAGAGTGAATTTTTCTGCCTGCGTGCGGGCAGCTGTTTGTGTATCACTAATGCGATCGAGCGCGGCGTGCAGCTGCCCGGCAAAACTAATGGTCGGTTGCGGCAGTGATTCCTGCGCACGCGCACTCATCGCCGTAGCCTGTAACt                           <  1:640563/150‑1 (MQ=255)
                                                                                                                          ccGAGAGTGAATTTTTCTGCCTGCGTGCGGGCAGCTGTTTGTGTATCACTAATGCGATCGAGCGCGGCGTGCAGCTGCCCGGCAAAACTAATGGTCGGTTGCGGCAGTGATTCCTGCGCACGCGCACTCATCGCCGTAGCCTGTAACTgg                         >  1:1506797/1‑150 (MQ=255)
                                                                                                                           cGAGAGTGAATTTTTCTGCCTGCGTGCGGGCAGCTGTTTGTGTATCACTAATGCGATCGAGCGCTGCGTGCAGCTGCCCGGCAAAACTAATGGTCGGTTGCGGCAGTGATTCCTGCGCACGCGCACTCATCGCCGTAGCCTGTAACTGGc                        >  1:623763/1‑150 (MQ=255)
                                                                                                                           cGAGAGTGAATTTTTCTGCCTGCGTGCGGGCAGCTGTTTGTGTATCACTAATGCGATCGAGCGCGGCGTGCAGCTGCCCGGCAAAACTAATGGTCGGTTGCGGCAGTGATTCCTGCGCACGCGCACTCATCGCCGTAGCCTGTAACTGGc                        <  1:1139927/150‑1 (MQ=255)
                                                                                                                             agagTGAATTTTTCTGCCTGCGTGCGGGCAGCTGTTTGTGTATCACTAATGCGATCGAGCGCGGCGTGCAGCTGCCCGGCAAAACTAATGGTCGGTTGCGGCAGTGATTCCTGCGCACGCGCACTCATCGCCGTAGCCTGTAACTGGCTg                      <  1:1508762/150‑1 (MQ=255)
                                                                                                                              gagTGAATTTTTCTGCCTGCGTGCGGGCAGCTGTTTGTGTATCACTAATGCGATCGAGCGCGGCGTGCAGCTGCCCGGCAAAACTAATGGTCGGTTGCGGCAGTGATTCCTGCGCACGCGCACTCATCGCCGTAGCCTGTAACTGGCTGa                     >  2:1538993/1‑150 (MQ=255)
                                                                                                                              gagTGAATTTTTCTGCCTGCGTGCGGGCAGCTGTTTGTGTATCACTAATGCGATCGAGCGCGGCGTGCAGCTGCCCGGCAAAACTAATGGTCGGTTGCGGCAGTGATTCCTGCGCACGCGCACTCATCGCCGTAGCCTGTAACTGGCTGa                     <  1:1101338/150‑1 (MQ=255)
                                                                                                                                  gAATTTTTCTGCCTGCGTGCGGGCAGCTGTTTGTGTATCACTAATGCGATCGAGCGCGGCGTGCAGCTGCCCGGCAAAACTAATGGTCGGTTGCGGCAGTGATTCCTGCGCACGCGCACTCATCGCCGTAGCCTGTAAc                            >  1:775498/1‑139 (MQ=255)
                                                                                                                                  gAATTTTTCTGCCTGCGTGCGGGCAGCTGTTTGTGTATCACTAATGCGATCGAGCGCGGCGTGCAGCTGCCCGGCAAAACTAATGGTCGGTTGCGGCAGTGATTCCTGCGCACGCGCACTCATCGCCGTAGCCTGTAAc                            <  2:775498/139‑1 (MQ=255)
                                                                                                                                  gAATTTTTCTGCCTGCGTGCGGGCAGCTGTTTGTGTATCACTAATGCGATCGAGCGCGGCGTGCAGCTGCCCGGCAAAACTAATGGTCGGTTGCGGCAGTGATTCCTGCGCACGCGCACTCATCGCCGTAGCCTGTAACTGGCTGATAAc                 >  2:226234/1‑150 (MQ=255)
                                                                                                                                   caaacagctgcccgcACGCAGGCAGCTGTTTGTGTATCACTAATGCGATCGAGCGCGGCGTGCAGCTGCCCGGCAAAACTAATGGTCGGTTGCGGCAGTGATTCCTGCGCACGCGCACTCATCGCCGTAGCCTGTAACTGGCTGATAAcc                >  1:1082995/8‑150 (MQ=255)
                                                                                                                                   aaTTTTTCTGCCTGCGTGCGGGCAGCTGTTTGTGTATCACTAATGCGATCGAGCGCGGCGTGCAGCTTCCCGGCAAAACTAATGGTCGGTTGCGGCAGTGATTCCTGCGCACGCGCACTCATCGCCGTAGCCTGTAACTGGCTGATAAcc                <  2:293122/150‑1 (MQ=255)
                                                                                                                                     aacagctgcccgcACGCAGGCAGCTGTTTGTGTATCACTAATGCGATCGAGCGCGGCGTGCAGCTGCCCGGCAAAACTAATGGTCGGTTGCGGCAGTGATTCCTGCGCACGCGCACTCATCGCCGTAGCCTGTAACTGGCTGATAAcccc              >  2:1487920/6‑150 (MQ=255)
                                                                                                                                      ttttCTGCCTGCGTGCGGGCAGCTGTTTGTGTATCACTAATGCGATCGAGCGCGGCGTGCAGCTGCCCGGCAAAACTAATGGTCGGTTGCGGCAGTGATTCCTGCGCACGCGCACTCATCGCCGTAGCCTGTAACTGGCTGATAACCCCt             <  2:1324916/150‑1 (MQ=255)
                                                                                                                                        ttCTGCCTGCGTGCGGGCAGCTGTTTGTGTATCACTAATGCGATCGAGCGCGGCGTGCAGCTGCCCGGCAAAACTAATGGTCGGTTGCGGCAGTGATTCCTGCGCACGCGCACTCATCGCCGTAGCCTGTAACTGGCTGATAACCCCTTc           <  2:685438/150‑1 (MQ=255)
                                                                                                                                            gcctgcGTGCGGGCAGCTGTTTGTGTATCACTAATGCGATCGAGCGCGGCGTGCAGCTGCCCGGCAAAACTAATGGTCGGTTGCGGCAGTTATTCCTGCGCACGCGCAGTCATCGCCGTAGCCTGTAACTGGCTGATAACCCCTTCAATc       <  2:500520/150‑1 (MQ=255)
                                                                                                                                            gcctgcGTGCGGGCAGCTGTTTGTGTATCACTAATGCGATCGAGCGCGGCGTGCAGCTGCCCGGCAAAACTAATGGTCGGTTGCGGCAGTGATTCCTGCGCACGCGCACTCATCGCCGTAGCCTGTAACTGGCTGATAACCCCTTCAATc       >  2:1155858/1‑150 (MQ=255)
                                                                                                                                               tgcgtgcgGGCAGCTGTTTGTGTATCACTAATGCGATCGAGCGCTGCGTGCAGCTGCCCGGCAAAACTAATGGTCGGTTGCGGCAGTGATTCCTGCGCACGCGCACTCATCGCCGTAGCCTGTAACTGGCTGATAACCCCTTCAATcccc    <  2:623763/150‑1 (MQ=255)
                                                                                                                                                 cgtgcgGGCAGCTGTTTGTGTATCACTAATGCGATCGAGCGCGGCGTGCAGCTGCCCGGCAAAACTAATGGTCGGTTGCGGCAGTGATTCCTGCGCACGCGCACTCATCGCCGTAGCCTGTAACTGGCTGATACCCCCTTCAATCCCCTg  >  1:840749/1‑150 (MQ=255)
                                                                                                                                                  |                                                                                                                                                    
GCATGCTCATCACTTCCTGATACGCCGCCACCAGCTTATTACGCACCTGAATCCCCATTTGCATAGAAACTGAGGCTTTTTGCATATCGGTCATCACATCGTTTAACGCCACGCCGGGTTCACCGAGAGTGAATTTTTCTGCCTGCGTGCGGGCAGCTGTTTGTGTATCACTAATGCGATCGAGCGCGGCGTGCAGCTGCCCGGCAAAACTAATGGTCGGTTGCGGCAGTGATTCCTGCGCACGCGCACTCATCGCCGTAGCCTGTAACTGGCTGATAACCCCTTCAATCCCCTG  >  NC_000913/2012708‑2013002

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: