Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 574,494 A→G 7.6% intergenic (+155/+35) quuD → / ← nmpC DLP12 prophage; putative antitermination protein/DLP12 prophage; truncated outer membrane porin (pseudogene);IS, phage, Tn; Phage or Prophage Related; outer membrane porin protein; locus of qsr prophage

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913574,4940AG 7.6% 412.7 / 18.4 132intergenic (+155/+35)quuD/nmpCDLP12 prophage; putative antitermination protein/DLP12 prophage; truncated outer membrane porin (pseudogene);IS, phage, Tn; Phage or Prophage Related; outer membrane porin protein; locus of qsr prophage
Reads supporting (aligned to +/- strand):  ref base A (60/62);  new base G (2/8);  total (62/70)
Fisher's exact test for biased strand distribution p-value = 1.02e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.69e-01

Only 100 of 139 total aligned reads displayed.

AATTGTTTACTAAAAATTATTAAAAATGGGGCGTTGAGACGCCCCCAAAAATAAAGGGTAATATATAACAGAAGGTTTATATAGTTAGAAGCAAGGTTGTGCTTCTAAAGGAAGTGGCTTGAGGGAGCCACTTATATGTTGGGGAGGCAACGCCTCCCGCAACATATCTTTTTCGTAATCAGATTAGAACTGGTAAACCAGACCTACAGCAACGATGTCATCAGTGCTTACACCGAGTGCTTTAGGGAAGGTGCGAATAAGCGGGGAAATTCTTCTCGGCTGACTCAGTCATTTCAT  >  NC_000913/574346‑574642
                                                                                                                                                     |                                                                                                                                                    
aaTTGTTTACTACAAATTATTAAAAATGGGGCGTTGAGACGCCCCCAAAAATAAAGGGTAATATATAACAGAAGGTTTATATAGTTAGAAGCAAGGTTGTGCTTCTAAAGGAAGTGGCTTGAGGGAGCCACTTATATGTTGGGGAGGCaa                                                                                                                                                     <  1:458300/150‑1 (MQ=255)
aaTTGTTTACTAAAAATTATTAAAAATGGGGCGTTGAGACGCCCCCAAAAATAAAGGGTAATATATAACAGAAGGTTTATATAGTTAGAAGCAAGGTTGTGCTTCTAAAGGAAGTGGCTTGAGGGAGCCACTTATATGTTGGGGAGGCaa                                                                                                                                                     <  1:483043/150‑1 (MQ=255)
aaTTGTTTACTAAAAATTATTAAAAATGGGGCGTTGAGACGCCCCCAAAAATAAAGGGTAATATATAACAGAAGGTTTATATAGTTAGAAGCAAGGTTGTGCTTCTAAAGGAAGTGGCTTGAGGGAGCCACTTATATGTTGGGGAGGCaa                                                                                                                                                     <  1:70459/150‑1 (MQ=255)
aaTTGTTTACTAAAAATTATTAAAAATGGGGCGTTGAGACGCCCCCAAAAATAAAGGGTAATATATAACAGAAGGTTTATATAGTTAGAAGCAAGGTTGTGCTTCTAAAGGAAGTGGCTTGAGGGAGCCACTTATATGTTGGGGAGGCaa                                                                                                                                                     >  1:1830520/1‑150 (MQ=255)
 aTTGTTTACTAAAAATTATTAAAAATGGGGCGTTGAGACGCCCCCAAAAATAAAGGGTAATATATAACAGAAGGTTTATATAGTTAGAAGCAAGGTTGTGCTTCTAAAGGAAGTGGCTTGAGGGAGCCACTTATATGTTGGGGAGGCAAc                                                                                                                                                    <  2:921705/150‑1 (MQ=255)
   tGTTTACTAAAAATTATTAAAAATGGGGCGTTGAGACGCCCCCAAAAATAAAGGGTAATATATAACAGAAGGTTTATATAGTTAGAAGCAAGGTTGTGCTTCTAAAGGAAGTGGCTTGAGGGAGCCACTTATATGTTGGGGAGGCAACGc                                                                                                                                                  <  2:1188675/150‑1 (MQ=255)
       tACTAAAAATTATTAAATATGGGGCGTTGAGACGCCCCCAAAAATAAAGGGTAATATATAACAGAAGGTTTAGATAGTTAGAAGCAAGGTTGTGCTTCTAAAGGAAGTGGCTTGAGGGAGCCACGTATATGTTGGGGAGGCAACGCCTcc                                                                                                                                              >  2:1058480/1‑150 (MQ=255)
       tACTAAAAATTATTAAAAATGGGGCGTTGAGACGCCCCCAAAAATAAGGGGTAATATATAACAGAAGGGTTATATAGTTAGAAGCAAGGTTGTGCTTCTAAAGGACGTGGCTTGAGGGAGCCACTTATATGTTGGGGAGGCAACGCCTcc                                                                                                                                              >  2:13119/1‑150 (MQ=255)
         cTAAAAATTATTAAAAATGGGGCGTTGAGACGCCCCCAAAAATAAAGGGTAATATATAACAGAAGGTTTATATAGTTAGAAGCAAGGTTGTGCTTCTAAAGGAAGTGGCTTGAGGGAGCCACTTATATGTTGGGGAGGCGTTGCCTcccc                                                                                                                                            <  2:2148068/150‑2 (MQ=255)
         cTAAAAATTATTAAAAATGGGGCGTTGAGACGCCCCCAAAAATAAAGGGTAATATATAACAGAAGGTTTATATAGTTAGAAGCAAGGTTGTGCTTCTAAAGGAAGTGGCTTGAGGGAGCCACTTATATGTTGGGGAGGCGTTGCCTcccc                                                                                                                                            <  2:1967552/150‑2 (MQ=255)
         cTAAAAATTATTAAAAATGGGGCGTTGAGACGCCCCCAAAAATAAAGGGTAATATATAACAGAAGGTTTATATAGTTAGAAGCAAGGTTGTGCTTCTAAAGGAAGTGGCTTGAGGGAGCCACTTATATGTTGGGGAGGCGTTGCCTcccc                                                                                                                                            <  2:1962579/150‑2 (MQ=255)
         cTAAAAATTATTAAAAATGGGGCGTTGAGACGCCCCCAAAAATAAAGGGTAATATATAACAGAAGGTTTATATAGTTAGAAGCAAGGTTGTGCTTCTAAAGGAAGTGGCTTGAGGGAGCCACTTATATGTTGGGGAGGCGTTGCCTcccc                                                                                                                                            <  2:671336/150‑2 (MQ=255)
          tAAAAATTATTAAAAATGGGGCGTTGAGACGCCCCCAAAAATAAAGGGTAATATATAACAGAAGGTTTATATAGTTAGAAGCAAGGTTGTGCTTCTAAAGGAAGTGGCTTGAGGGAGCCACTTATATGTTGGGGAGGCAACGCCTCCCGc                                                                                                                                           >  2:1228104/1‑150 (MQ=255)
            aaaaTTATTAAAAATGGGGCGTTGAGACGCCCCCAAAAATAAAGGGTAATATATAACAGAAGGTTTATATAGTTAGAAGCAAGGTTGTGCTTCTAAAGGAAGTGGCTTGAGGGAGCCACTTATATGTTGGGGAGGCAACGCCTCCCGCaa                                                                                                                                         >  1:932572/1‑150 (MQ=255)
             aaaTTATTAAAAATGGGGCGTTGAGACGCCCCCAAAAATAAAGGGTAATATATAACAGAAGGTTTATATAGTTAGAAGCAAGGTTGTGCTTCTAAAGGAAGTGGCTTGAGGGAGCCACTTATATGTTGGGGAGGCAACGCCTCCCGCAAc                                                                                                                                        <  2:697518/150‑1 (MQ=255)
                      aaaaTGGGGCGTTGAGACGCCCCCAAAAATAAAGGGTAATATATAACAGAAGGTTTATATAGTTAGAAGCAAGGTTGTGCTTCTAAAGGAAGTGGCTTGAGGGAGCCACTTATATGTTGCGGGAGGCGTTGCCTCCCCAACatataagtg                                                                                                                               <  1:2324617/150‑6 (MQ=38)
                      aaaaTGGGGCGTTGAGACGCCCCCAAAAATAAAGGGTAATATATAACAGAAGGTTTATATAGTTAGAAGCAAGGTTGTGCTTCTAAAGGAAGTGGCTTGAGGGAGCCACTTATATGTTGGGGAGGCAACGCCTCCCGCAACATATCtttt                                                                                                                               >  2:1543057/1‑150 (MQ=255)
                      aaaaTGGGGCGTTGAGACGCCCCCAAAAATAAAGGGTAATATATAACAGAAGGTTTATATAGTTAGAAGCAAGGTTGTGCTTCTAAAGGAAGTGGCTTGAGGGAGCCACTTATATGTTGGGGAGGCAACGCCTCCCGCAACATATCtttt                                                                                                                               >  2:476069/1‑150 (MQ=255)
                       aaaTGGGGCGTTGAGACGCCCCCAAAAATAAAGGGTAATATATAACAGAAGGTTTATATAGTTAGAAGCAAGGTTGTGCTTCTAAAGGAAGTGGCTTGAGGGAGCCACTTATATGTTGGGGAGGCAACGCCTCCCGCAACATATCttttt                                                                                                                              >  2:2148437/1‑150 (MQ=255)
                       aaaTGGGGCGTTGAGACGCCCCCAAAAATAAAGGGTAATATATAACAGAAGGTTTATATAGTTAGAAGCAAGGTTGTGCTTCTAAAGGAAGTGGCTTGAGGGAGCCACTTATATGTTGGGGAGGCAACGCCTCCCGCAACATATCttttt                                                                                                                              >  2:2221776/1‑150 (MQ=255)
                       aaaTGGGGCGTTGAGACGCCCCCAAAAATAAAGGGTAATATATAACAGAAGGTTTATATAGTTAGAAGCAAGGTTGTGCTTCTAAAGGAAGTGGCTTGAGGGAGCCACTTATATGTTGGGGAGGCAACGCCTCCCGCAACATATCttttt                                                                                                                              >  2:36508/1‑150 (MQ=255)
                       aaaTGGGGCGTTGAGACGCCCCCAAAAATAAAGGGTAATATATAACAGAAGGTTTATATAGTTAGAAGCAAGGTAGTGCTTCTAAAGGAAGTGGCTTGAGGGAGCCACTTATATGTTGGGGAGGCAACGCCTCCCGCAACATATCttatt                                                                                                                              >  2:1670496/1‑150 (MQ=255)
                       aaaTGGGGCGTTGAGACGCCCCCAAAAATAAAGGGCAATATATAACAGAAGGTGTATATAGTTAGAAGCAAGGTTGTGCTTCTAAAGGAAGTGGCTTGAGGGAGCCACTTATATGTTGGGGAGGCAACGCCTCCCGCACCATATCttttt                                                                                                                              >  2:1670495/1‑150 (MQ=255)
                        aaTGGGGCGTTGAGCCGCCCCCAAAAATAAAGGGTAATATATAACAGACGGTTTATATAGTTAGAAGCAAGGTTGTGCTTCTAAAGGAAGTGGCTTGAGGGAGCAACTTATATTTTGGGGAGGCAACGCCTCCCGCAACATATCTTTTTc                                                                                                                             <  1:1894592/150‑1 (MQ=255)
                         aTGGGGCGTTGAGACGCCCCCAAAAATAAAGGGTAATATATAACAGAAGGTTTATATAGTTAGAAGCAAGGTTGTGCTTCTAAAGGAAGTGGCTTGAGGGAGCCACTTATATGTTGGGGAGGCAACGCCTCCCGCAACATATCTTTTTCg                                                                                                                            <  2:1659420/150‑1 (MQ=255)
                         aTGGGGCGTAGAGACGCCCCCAAAAAAAAACCCTAATATATAACACAAGCTTTATATAGTTAGAACCAAGGTTGTGCTTCTAAAGGAAGTGGCTTGAGGGAGCCACTTATATGTTGGGGAGGCAACGCCTCCCGCAACATATCTTTTTCg                                                                                                                            <  2:1075349/150‑1 (MQ=255)
                            gggCGTTGAGACGCCCCCAAAAATAAAGGGTAATATATAACAGAAGGTTTATATAGTTAGAAGCAAGGTTGTGCTTCTAAAGGAAGTGGCTTGAGGGAGCCACTTATATGTTGGGGAGGCAACGCCTCCCGCAACATATCTTTTTCGTaa                                                                                                                         >  1:1374746/1‑150 (MQ=255)
                               cGTTGAGACGCCCCCAAAAATAAAGGGTAATATATAACAGAAGGTTTATATAGTTAGAAGCAAGGTTGTGCTTCTAAAGGAAGTGGCTTGAGGGAGCCACTTATATGTTGGGGAGGCAACGCCTCCCGCAACATATCTTTTTCGTAATCa                                                                                                                      >  1:880912/1‑150 (MQ=255)
                               cGTTGAGACGCCCCCAAAAATAAAGGGTAATATATAACAGAAGGTTTATATAGTTAGAAGCAAGGTTGTGCTTCTAAAGGAAGTGGCTTGAGGGAGCCACTTATATGTTGGGGAGGCAACGCCTCCCGCAACATATCTTTTTCGTAATCa                                                                                                                      >  1:880860/1‑150 (MQ=255)
                               cGTTGAGACGCCCCCAAAAATAAAGGGTAATATATAACAGAAGGTTTATATAGTTAGAAGCAAGGTTGTGCTTCTAAAGGAAGTGGCTTGAGGGAGCCACTTATATGTTGGGGAGGCAACGCCTCCCGCAACATATCTTTTTCGTAATCa                                                                                                                      <  2:1642847/150‑1 (MQ=255)
                               cGTTGAGACGCCCCCAAAAATAAAGGGTAATATATAACAGAAGGTTTATATAGTTAGAAGCAAGGTTGTGCTTCTAAAGGAAGTGGCTTGAGGGAGCCACTTATATGTTGGGGAGGCAACGCCTCCCGCAACATATCTTTTTCGTAATCa                                                                                                                      >  1:1184185/1‑150 (MQ=255)
                               cGTTGAGACGCCCCCAAAAATAAAGGGTAATATATAACAGAAGGTTTATATAGTTAGAAGCAAGGTTGTGCTTCTAAAGGAAGTGGCTTGAGGGAGCCACTTATATGTTGGGGAGGCAACGCCTCCCGCAACATATCTTTTTCGTAATCa                                                                                                                      >  1:1184164/1‑150 (MQ=255)
                                       cGCCCCCAAAAATAAAGGGTAATATATAACAGAAGGTTTATATAGTTAGAAGCAAGGTTGTGCTTCTAAAGGAAGTGGCTTGAGGGAGCCACTTATATGTTGGGGAGGCAACGCCTCCCGCAACATATCTTTTTCGTAATCAGATTAGaa                                                                                                              >  2:375591/1‑150 (MQ=255)
                                             cAAAAATAAAGGGTAATATATAACAGAAGGTTTATATAGTTAGAAGCAAGGTTGTGCTTCTAAAGGAAGTGGCTTGAGGGAGCCACTTATATGTTGGGGAGGCAACGCCTCCCACAACATATCTTTTTCGTAATCAGATTAGAACTGGTa                                                                                                        >  2:1938173/1‑150 (MQ=255)
                                               aaaaTAAAGGGTAATATATAACAGAAGGTTTATATAGTTAGAAGCAAGGTTGTGCTTCTAAAGGAAGTGGCTTGAGGGAGCCACTTATATGTTGGGGAGGCAACGCCTCCCGCAACATATCTTTTTCGTCATCAGATTAGAACTGGTaag                                                                                                      >  2:1294917/1‑149 (MQ=255)
                                               aaaaTAAAGGGTAATATATAACAGAAGGTTTATATAGTTAGAAGCAAGGTTGTGCTTCTAAAGGAAGTGGCTTGAGGGAGCCACTTATATGTTGGGGAGGCAACGCCTCCCGCAACATATCTTTTTCGTAATCAGATTAGAACTGGTaaa                                                                                                      <  1:634781/150‑1 (MQ=255)
                                               aaaaTAAAGCCTAATATATAACACAAGCATTATATAGTTAGAAGCAACGTTCTGCTTCTAAAGGAAGTGGCTTGAGGGAGCCACTTATATGTTGGGGAGGCAACGCCTCCCGCAACATATCTTTTTCGTAATCAGATTAGAACTGGTaaa                                                                                                      <  2:1676011/150‑1 (MQ=255)
                                                 aaTAAAGGGTAATATATAACAGAAGGTTTATATAGTTAGAAGCAAGGTTGTGCTTCTAAAGGAAGTGGCTTGAGGGAGCCACTTATATGTTGGGGAGGCAACGCCTCCCGCAACATATCTTTTTCGTAATCAGATTAGAACTGGTAAAcc                                                                                                    <  2:257256/150‑1 (MQ=255)
                                                    aaaGGGTAATATATAACAGAAGGTTTATATAGTTAGAAGCAAGGTTGTGCTTCTAAAGGAAGTGGCTTGAGGGAGCCACTTATATGTTGGGGAGGCAACGCCTCCCGCAACATATCTTTTTCGTAATCAGATTAGAACTGGTAAACCAGa                                                                                                 <  1:385545/150‑1 (MQ=255)
                                                       gggTAATATATAACAGAAGGTTTATATAGTTAGAAGCAAGGTTGTGCTTCTAAAGGAAGTGGCTTGAGGGAGCCACTTATATGTTGGGGAGGCAACGCCTCCCGCAACATATCTTTTTCGTAATCAGATTAGAACTGGTAAACCAGACCt                                                                                              <  2:18360/150‑1 (MQ=255)
                                                         gTAATATATAACAGACGGTTTATATAGTTAGAAGCAAGGTTGTGCTTCTAAAGGAAGTGGCTTGAGGGAGCCACTTATATGTTGGGGAGGCAACGCCTCCCGCAACATATCTTTTTCGTAATCAGATTAGAACTGGTAAACCAGACCTAc                                                                                            <  1:496681/150‑1 (MQ=255)
                                                         gTAATATATAACAGAAGGTTTATATAGTTAGAAGCAAGGTTGTGCTTCTAAAGGAAGTGGCTTGAGGGAGCCACTTATATGTTGGGGAGGCAACGCCTCCCGCAACATATCTTTTTCGTAATCAGATTAGAACTGGTAAACCAGACCTAc                                                                                            <  2:2159649/150‑1 (MQ=255)
                                                          tAATATATAACAGAAGGTTTATATAGTTAGAAGCAAGGTTGTGCTTCTAAAGGAAGTGGCTTGAGGGAGCCACTTATATGTTGGGGAGGCAACGCCTCCCGCAACATATCTTTTTCGTAATCAGATTAGAACTGGTAAACCAGACCTACa                                                                                           >  2:1445374/1‑150 (MQ=255)
                                                          tAATATATAACAGAAGGTTTATATAGTTAGAAGCAAGGTTGTGCTTCTAAAGGAAGTGGCTTGAGGGAGCCACTTATATGTTGGGGAGGCAACGCCTCCCGCAACATATCTTTTTCGTAATCAGATTAGAACTGGTAAACCAGACCTACa                                                                                           >  2:503874/1‑150 (MQ=255)
                                                          tAATATATAACAGAAGGTTTATATAGTTAGAAGCAAGGTTGTGCTTCTAAAGGAAGTGGCTTGAGGGAGCCACTTATATGTTGGGGAGGCAACGCCTCCCGCAACATATCTTTTTCGTAATCAGATTAGAACTGGTAAACCAGACCTACa                                                                                           <  1:1071024/150‑1 (MQ=255)
                                                              atataACAGAAGGTTTATATAGTTAGAAGCAAGGTTGTGCTTCTAAAGGAAGTGGCTTGAGGGAGCCACTTATATGTTGGGGAGGCAACGCCTCCCGCAACATATCTTTTTCGTAATCAGATTAGAACTGGTAAACCAGACCTACAGCaa                                                                                       <  2:1726504/150‑1 (MQ=255)
                                                                         gttttATATAGTTAGAAGCAAGGTTGTTCTTCTATAGGAAGTGGCTTGAGGGAGCCAGTTATATGTTGGGTAGGCAACGCCTCCCGCAACATATCTTTTTAGTAATCAGATTAGAACTGGTAAACCAGACCTACAGCAACGATGtcatca                                                                            <  1:2297217/148‑1 (MQ=255)
                                                                         ggTTTATATAGTTAGAAGCAAGGTTGTGCTTCTAAAGGAAGTGGCTTGAGGGAGCCACTTATATGTTGGGGAGGCAACGCCTCCCGCAACATATCTTTTTCGTAATCAGATTAGAACTGGTAAACCAGACCTACAGCAACGATGtcatca                                                                            <  2:94325/150‑1 (MQ=255)
                                                                         ggTTTATATAGTTAGAAGCAAGGTTGTGCTTCTAAAGGAAGTGGCTTGAGGGAGCCACTTATATGTTGGGGAGGCAACGCCTCCCGCAACATATCTTTTTCGTAATCAGATTAGAACTGGTAAACCAGACCTACAGCAACGATGtcatca                                                                            >  2:2086467/1‑150 (MQ=255)
                                                                         ggTTTATATAGTTAGAAGCAAGGTTGTGCTTCTAAAGGAAGTGGCTTGAGGGAGCCACTTATATGTTGGGGAGGCAACGCCTCCCGCAACATATCTTTTTCGTAATCAGATTAGAACTGGTAAACCAGACCTACAGCAACGATGtcatca                                                                            >  2:2086604/1‑150 (MQ=255)
                                                                           tttATATAGTTAGAAGCAAGGTTGTGCTTCTAAAGGAAGTGGCTTGAGGGAGCCACTTATATGTTGGGGAGGCAACGCCTCCCGCAACATATCTTTTTCGTAATCAGATTAGAACTGGTAAACCAGACCTACAGCAACGATGTCATCAGt                                                                          >  1:449397/1‑150 (MQ=255)
                                                                                    ttAGAAGCAAGGTTGTGCTTCTAAAGGAAGTGGCTTGAGGGAGCCACTTATATGTTGGGGAGGCAACGCCTCCCGCAACATATCTTTTTCGTAATCAGATTAGAACTGGTAAACCAGACCTACAGCAACGATGTCATCAGTGCTTACAcc                                                                 >  2:1314520/1‑150 (MQ=255)
                                                                                    ttAGAAGCAAGGTTGTGCTTCTAAAGGAAGTGGCTTGAGGGAGCCACTTATATGTTGGGGAGGCAACGCCTCCCGCAACATATCTTTTTCGTAATCAGATTAGAACTGGTAAACCAGACCTACAGCAACGATGTCATCAGTGCTTACAcc                                                                 >  2:550215/1‑150 (MQ=255)
                                                                                       gAAGCAAGGTTGTGCTTCTAAAGGAAGTGGCTTGAGGGAGCCACTTATATGTTGGGGAGGCAACGCCTCCCGCAACATATCTTTTTCGTAATCAGATTAGAACTGGTAAACCAGACCTACAGCAACGATGTCATCAGTGCTTACACCGAg                                                              <  2:313908/150‑1 (MQ=255)
                                                                                          gCAAGGTTGTGCTTCTAAAGGAAGTGGCTTGAGGGAGCCACTTATATGTTGGGGAGGCAACGCCTCCCGCAACATATCTTTTTCGTAATCAGATTAGAACTGGTAAACCAGACCTACAGCAACGATGTCATCAGTGCTTACACCGAGTGc                                                           <  1:822525/150‑1 (MQ=255)
                                                                                          gCAAGGTTGTGCTTCTAAAGGAAGTGGCTTGAGGGAGCCACTTATATGTTGGGGAGGCAACGCCTCCCGCAACATATCTTTTTCGTAATCAGATTAGAACTGGTAAACCAGACCTACAGCAACGATGTCATCAGTGCTTACACCGAGTGc                                                           <  1:1590047/150‑1 (MQ=255)
                                                                                           cAAGGTTGTGCTTCTAAAGGAAGTGGCTTGAGGGAGCCACTTATATGTTGGGGAGGCAACGCCTCCCGCAACATATCTTTTTCGTAATAAGATTAGAACTGGTAAACCAGACCTACAGCAACGATGTCATCAGTGCTTACACCGAGTGCt                                                          >  2:30882/1‑150 (MQ=255)
                                                                                              ggTTGTGCTTCTAAAGGAAGTGGCTTGAGGGAGCCACTTATATGTTGGGGAGGCAACGCCTCCCGCAACATATCTTTTTCGTAATCAGATTAGAACTGGTAAACCAGACCTACAGCAACGATGTCATCAGTGCTTACACCGAGTGCTTTa                                                       <  1:1013767/150‑1 (MQ=255)
                                                                                              ggTTGTGCTTCTAAAGGAAGTGGCTTGAGGGAGCCACTTATATGTTGGGGAGGCAACGCCTCCCGCAACATATCTTTTTCGTAATCAGATTAGAACTGGTAAACCAGACCTACAGCAACGATGTCATCAGTGCTTACACCGAGTGCTTTa                                                       <  2:336923/150‑1 (MQ=255)
                                                                                              ggTTGTGCTTCTAAAGGAAGTGGCTTGAGGGAGCCACTTATATGTTGGGGAGGCAACGCCTCCCGCAACATATCTTTTTCGTAATCAGATTAGAACTGGTAAACCAGACCTACAGCAACGATGTCATCAGTGCTTACACCGAGTGCTTTa                                                       <  2:2178383/150‑1 (MQ=255)
                                                                                              ggTTGTGCTTCTAAAGGAAGTGGCTTGAGGGAGCCACTTATATGTTGGGGAGGCAACGCCTCCCGCAACATATCTTTTTCGTAATCAGATTAGAACTGGTAAACCAGACCTACAGCAACGATGTCATCAGTGCTTACACCGAGTGCTTTa                                                       <  2:2178105/150‑1 (MQ=255)
                                                                                              ggTTGTGCTTCTAAAGGAAGTGGCTTGAGGGAGCCACTTATATGTTGGGGAGGCAACGCCTCCCGCAACATATCTTTTTCGTAATCAGATTAGAACTGGTAAACCAGACCTACAGCAACGATGTCATCAGTGCTTACACCGAGTGCTTTa                                                       <  2:1943440/150‑1 (MQ=255)
                                                                                              ggTTGTGCTTCTAAAGGAAGTGGCTTGAGGGAGCCACTTATATGTTGGGGAGGCAACGCCTCCCGCAACATATCTTTTTCGTAATCAGATTAGAACTGGTAAACCAGACCTACAGCAACGATGTCATCAGTGCTTACACCGAGTGCTTTa                                                       <  2:1681323/150‑1 (MQ=255)
                                                                                               gTTGTGCTTCTAAAGGAAGTGGCTTGAGGGAGCCACTTATATGTTGGGGAGGCAACGCCTCCCGCAACATATCTTTTTCGTAATCAGATTAGAACTGGTAAACCAGACCTACAGCAACGATGTCATCAGTGCTTACACCGAGTGCTTTAg                                                      <  1:927482/150‑1 (MQ=255)
                                                                                                gtgtgCTTCTAAAGGAAGTGGCTTGAGGGAGCCACTTATATGTTGGGGAGGCAACGCCTCCCGCAACATATCTTTTTCGTAATCAGATTAGAACTGGTAAACCAGACCTACAGCAACGATGTCATCAGTGCTTACACCGAGTGCTTTAgg                                                     <  2:1781357/149‑1 (MQ=255)
                                                                                                   tgCTTCTAAAGGAAGTGGCTTGAGGGAGCCACTTATATGTTGGGGAGGCAACGCCTCCCGCAACATATCTTTTTCGTAATCAGATTAGAACTGGTAAACCAGACCTACAGCAACGATGTCATCAGTGCTTACACCGAGTGCTTTAGGGaa                                                  <  1:555128/150‑1 (MQ=255)
                                                                                                   tgCTTCTAAAGGAAGTGGCTTGAGGGAGCCACTTATATGTTGGGGAGGCAACGCCTCCCGCAACATATCTTTTTCGTAATCAGATTAGAACTGGTAAACCAGACCTACAGCAACGATGTCATCAGTGCTTACACCGAGTGCTTTAGGGaa                                                  <  1:555253/150‑1 (MQ=255)
                                                                                                   tgCTTCTAAAGGAAGTGGCTTGAGGGAGCCACTTATATGTTGGGGAGGCAACGCCTCCCGCAACATATCTTTTTCGTAATCAGATTAGAACTGGTAAACCAGACCTACAGCAACGATGTCATCAGTGCTTACACCGAGTGCTTTAGGGaa                                                  <  1:2080565/150‑1 (MQ=255)
                                                                                                   tgCTTCTAAAGGAAGTGGCTTGAGGGAGCCACTTATATGTTGGGGAGGCAACGCCTCCCGCAACATATCTTTTTCGTAATCAGATTAGAACTGGTAAACCAGACCTACAGCAACGATGTCATCAGTGCTTACACCGAGTGCTTTAGGGaa                                                  <  1:2080620/150‑1 (MQ=255)
                                                                                                     cTTCTAAAGGAAGTGGCTTGAGGGAGCCACTTATATGTTGGGGAGGCAACGCCTCCCGCAACATATCTTTTTCGTAATCAGATTAGAACTGGTAAACCAGACCTACAGCAACGATGTCATCAGTGCTTACACCGAGTGCTTTAGGGAAgg                                                <  1:507838/150‑1 (MQ=255)
                                                                                                     cTTCTAAAGGAAGTGGCTTGAGGGAGCCACTTATATGTTGGGGAGGCAACGCCTCCCGCAACATATCTTTTTCGTAATCAGATTAGAACTGGTAAACCAGACCTACAGCAACGATGTCATCAGTGCTTACACCGAGTGCTTTAGGGAAgg                                                <  1:507349/150‑1 (MQ=255)
                                                                                                            aGGAAGTGGCTTGAGGGAGCCACTTATATGTTGGGGAGGCAACGCCTCCCGCAACATATCTTTTTCGTAATCAGATTAGAACTGGTAAACCAGACCTACAGCAACGATGTCATCAGTGCTTACACCGAGTGCTTTAGGGAAGGTGCGAAt                                         >  2:2113/1‑150 (MQ=255)
                                                                                                                aGTGGCTTGAGGGAGCCACTTATATGTTGGGGAGGCAACGCCTCCCGCAACATATCTTTTTCGTAATCAGATTAGAACTGGTAAACCAGACCTACAGCAACGATGTCATCAGTGCTTACACCGAGTGCTTTAGGGAAGGTGCGAATAAGc                                     <  2:579321/150‑1 (MQ=255)
                                                                                                                   ggCTTGAGGGAGCCACTTATATGTTGGGGAGGCAACGCCTCCCGCAACATATCTTTTTCGTAATCAGATTAGAACTGGTAAACCAGACCTACAGCAACGATGTCATCAGTGCTTACACCGAGTGCTTTAGGGAAGGTGCGAATAAGCggg                                  <  1:130674/150‑1 (MQ=255)
                                                                                                                    gCTTGAGGGAGCCACTTATATGTTGGGGAGGCAACGCCTCCCGCAACATATCTTTTTCGTAATCAGATTAGAACTGGTAAACCAGACCTACAGCAACGATGTCATCAGTGCTTACACCGAGTGCTTTAGGGAAGGTGCGAATAAGCgggg                                 >  1:757708/1‑150 (MQ=255)
                                                                                                                      ttGAGGGAGCCACTTATATGTTGGGGAGGCAACGCCTCCCGCAACATATCTTTTTCGTAATCAGATTAGAACTGGTAAACCAGACCTACAGCAACGATGTCATCAGTGCTTACACCGAGTGCTTTAGGGAAGGTGCGAATAAGCGGGGaa                               >  1:2239305/1‑150 (MQ=255)
                                                                                                                      ttGAGGGAGCCACTTATATGTTGGGGAGGCAACGCCTCCCGCAACATATCTTTTTCGTAATCAGATTAGAACTGGTAAACCAGACCTACAGCAACGATGTCATCAGTGCTTACACCGAGTGCTTTAGGGAAGGTGCGAATAAGCGGGGaa                               >  2:401782/1‑150 (MQ=255)
                                                                                                                      ttGAGGGAGCCACTTATATGTTGGGGAGGCAACGCCTCCCGCAACATATCTTTTTCGTAATCAGATTAGAACTGGTAAACCAGACCTACAGCAACGATGTCATCAGTGCTTACACCGAGTGCTTTAGGGAAGGTGCGAATAAGCGGGGaa                               >  2:1164768/1‑150 (MQ=255)
                                                                                                                       tGAGGGAGCCACTTATATGTTGGGGAGGCAACGCCTCCCGCAACATATCTTTTTCGTAATCAGATTAGAACTGGTAAACCAGACCTACAGCAACGATGTCATCAGTGCTTACACCGAGTGCTTTAGGGAAGGTGCGAATAAGCGGGGaaa                              >  1:2140156/1‑150 (MQ=255)
                                                                                                                         aGGGAGCCACTTATATGTTGGGGAGGCAACGCCTCCCGCAACATATCTTTTTCGTAATCAGATTAGAACTGGTAAACCAGACCTACAGCAACGATGTCATCAGTGCTTACACCGAGTGCTTTAGGGAAGGTGCGAATAAGCGGGGAAAtt                            <  1:1230237/150‑1 (MQ=255)
                                                                                                                            gAGCCACTTATATGTTGGGGAGGCAACGCCTCCCGCAACATATCTTTTTCGTAATCAGATTAGAACTGGTAAACCAGACCTACAGCAACGATGTCATCAGTGCTTACACCGAGTGCTTTAGGGAAGGTGCGAATAAGCGGGGAAAttctt                         >  2:1814933/1‑150 (MQ=255)
                                                                                                                               ccACTTATATGTTGGGGAGGCAACGCCTCCCGCAACATATCTTTTTCGTAATCAGATTAGAACTGGTAAACCAGACCTACAGCAACGATGTCATCAGTGCTTACACCGAGTGCTTTAGGGAAGGTGCGAATAAGCGGGGAAATTCTtctc                      >  1:172008/1‑150 (MQ=255)
                                                                                                                                cACTTATATGTTGGGGAGGCAACGCCTCCCGCAACATATCTTTTTCGTAATCAGATTAGAACTGGTAAACCAGACCTACAGCAACGATGTCATCAGTGCTTACACCGAGTGCTTTAGGGAAGGTGCGAATAAGCGGGGAAATTCTTCTCg                     <  2:1754714/150‑1 (MQ=255)
                                                                                                                                cACTTATATGTTGGGGAGGCAACGCCTCCCGCAACATATCTTTTTCGTAATCAGATTAGAACTGGTAAACCAGACCTACAGCAACGATGTCATCAGTGCTTACACCGAGTGCTTTAGGGAAGGTGCGAATAAGCGGGGAAATTCTTCTCg                     <  2:1686117/150‑1 (MQ=255)
                                                                                                                                  cTTATATGTTGGGGAGGCAACGCCTCCCGCAACATATCTTTTTCGTAATCAGATTAGAACTGGTAAACCAGACCTACAGCAACGATGTCATCAGTGCTTACACCGAGTGCTTTAGGGAAGGTGCGAATAAGCGGGGAAATTCTTCTCGGc                   <  1:2221776/150‑1 (MQ=255)
                                                                                                                                  cTTATATGTTGGGGAGGCAACGCCTCCCGCAACATATCTTTTTCGTAATCAGATTAGAACTGGTAAACCAGACCTACAGCAACGATGTCATCAGTGCTTACACCGAGTGCTTTAGGGAAGGTGCGAATAAGCGGGGAAATTCTTCTCGGc                   <  1:36508/150‑1 (MQ=255)
                                                                                                                                   ttATATGTTGGGGAGGCAACGCCTCCCGCAACATATCTTTTTCGTTATCAGATTAGAACTGGTAAACCAGACCTACAGCAACGATGTCATCAGTGCTTACACCGAGTGCTTTAGGGAAGGTGCGAATAAGCGGGGAAATTCTTCTCGGCt                  <  1:1881339/150‑1 (MQ=255)
                                                                                                                                   ttATATGTTGGGGAGGCAACGCCTCCCGCAACATATCTTTTTCGTAATCAGATTAGAACTGGTAAACCAGACCTACAGCAACGATGTCATCAGTGCTTACACCGAGTGCTTTAGGGAAGGTGCGAATAAGCGGGGAAATTCTTCTCGGCt                  <  1:183293/150‑1 (MQ=255)
                                                                                                                                   ttATATGTTGGGGAGGCAACGCCTCCCGCAACATATCTTTTTCGTAATCAGATTAGAACTGGTAAACCAGACCTACAGCAACGATGTCATCAGTGCTTACACCGAGTGCTTTAGGGAAGGTGCGAATAAGCGGGGAAATTCTTCTCGGCt                  <  1:286976/150‑1 (MQ=255)
                                                                                                                                          ttGGGGAGGCAACGCCTCCCGCAACATATCTTTTTCGTAATCAGATTCGAACTGGTAAACCAGACCTACAGCAACGATGTCATCAGTGCTTACACCGAGTGCTTTAGGGAAGGTGCGAATAAGCGGGGAAATTCTTCTCGGCTGACTCAg           >  2:260913/1‑150 (MQ=255)
                                                                                                                                          atGGGGAGGCAACGCCTCCCGCAACATATCTTTTTCGTAATCAGATTAGAACTGGTAAACCAGACCTACAGCAACGATGTCATCAGTGCTTACACCGAGTGCTTTAGGGACGGTGCGAATAAGCGGGGAAATTCTTCTCGGCTGACTCAg           >  2:2327714/2‑150 (MQ=255)
                                                                                                                                           tGGGGAGGCAACGCCTCCCGCAACATATCTTTTTCGTAATCAGATTAGAACTGGTAAACCAGACCTACAGCAACGATGTCATCAGTGCTTACACCGAGTGCTTTAGGGAAGGTGCGAATAAGCGGGGAAATTCTTCTCGGCTGACTCAGt          >  2:1666451/1‑150 (MQ=255)
                                                                                                                                              gggAGGCAACGCCTCCCGCAACATATCTTTTTCGTAATCAGATTAGAACTGGTAAACCAGACCTACAGCAACGATGTCATCAGTGCTTACACCGAGTGCTTTAGGGAAGGTGCGAATAAGCGGGGAAATTCTTCTCGGCTGACTCAGtca        <  2:1992237/150‑1 (MQ=255)
                                                                                                                                               ggAGGCAACGCCTCCCGCAACATATCTTTTTCGTAATCAGATTAGAACTGGTAAACCAGACCTACAGCAACGATGTCATCAGTGCTTACACCGAGTGCTTTAGGGAAGGTGCGAATAAGCGGGGAAATTCTTCTCGGCTGACTCAGtcat       <  1:1732852/150‑1 (MQ=255)
                                                                                                                                                gAGGCAACGCCTCCCGCAACATATCTTTTTCGTAATCAGATTAGAACTGGTAAACCAGACCTACAGCAACGATGTCATCAGTGCTTACACCGAGTGCTTTAGGGAAGGTGCGAATAAGCGGGGAAATTCTTCTCGGCTGACTCAGtcatt      >  2:2239898/1‑150 (MQ=255)
                                                                                                                                                gAGGCAACGCCTCCCGCAACATATCTTTTTCGTAATCAGATTAGAACTGGTAAACCAGACCTACAGCAACGATGTCATCAGTGCTTACACCGAGTGCTTTAGGGAAGGTGCGAATAAGCGGGGAAATTCTTCTCGGCTGACTCAGtcatt      >  2:772611/1‑150 (MQ=255)
                                                                                                                                                gAGGCAACGCCTCCCGCAACAAATCTTTTTCGTAATCAGATTAGAACTGGTAAACCAGACCTACAGCAACGATGTCATCAGTGCTTACACCGAGTGCTTTAGGGAAGGTGCGAATAAGCGGGGAAATTCTTCTCGGCTGACTCAGtcatt      >  1:1015309/1‑150 (MQ=255)
                                                                                                                                                  ggCAACGCCTCCCGCAACATATCTTTTTCGTAATCAGATTAGAACTGGTAAACCAGACCTACAGCAACGATGTCATCAGTGCTTACACCGAGTGCTTTAGGGAAGGTGCGAATAAGCGGGGAAATTCTTCTCGGCTGACTCAGtcatttc    >  1:1918926/1‑150 (MQ=255)
                                                                                                                                                  ggCAACGCCTCCCGCAACATATCTTTTTCGTAATCAGATTAGAACTGGTAAACCAGACCTACAGCAACGATGTCATCAGTGCTTACACCGAGTGCTTTAGGGAAGGTGCGAATAAGCGGGGAAATTCTTCTCGGCTGACTCAGtcatttc    >  1:335143/1‑150 (MQ=255)
                                                                                                                                                    cAACGCCTCCCGCAACATATCTTTTTCGTAATCAGATTAGAACTGGTAAACCAGACCTACAGCAACGATGTCATCAGTGCTTACACCGAGTGCTTTAGGGAAGGTGCGAATAAGCGGGGAAATTCTTCTCGGCTGACTCAGtcatttcat  <  1:408049/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                    
AATTGTTTACTAAAAATTATTAAAAATGGGGCGTTGAGACGCCCCCAAAAATAAAGGGTAATATATAACAGAAGGTTTATATAGTTAGAAGCAAGGTTGTGCTTCTAAAGGAAGTGGCTTGAGGGAGCCACTTATATGTTGGGGAGGCAACGCCTCCCGCAACATATCTTTTTCGTAATCAGATTAGAACTGGTAAACCAGACCTACAGCAACGATGTCATCAGTGCTTACACCGAGTGCTTTAGGGAAGGTGCGAATAAGCGGGGAAATTCTTCTCGGCTGACTCAGTCATTTCAT  >  NC_000913/574346‑574642

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