Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 574,494 A→G 10.7% intergenic (+155/+35) ybcQ → / ← nmpC DLP12 prophage; putative antitermination protein/DLP12 prophage; putative outer membrane porin NmpC

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913574,4940AG10.7% 301.4 / 23.0 112intergenic (+155/+35)ybcQ/nmpCDLP12 prophage; putative antitermination protein/DLP12 prophage; putative outer membrane porin NmpC
Reads supporting (aligned to +/- strand):  ref base A (42/58);  new base G (3/9);  total (45/67)
Fisher's exact test for biased strand distribution p-value = 3.55e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.99e-01

Only 100 of 117 total aligned reads displayed.

CAATTGTTTACTAAAAATTATTAAAAATGGGGCGTTGAGACGCCCCCAAAAATAAAGGGTAATATATAACAGAAGGTTTATATAGTTAGAAGCAAGGTTGTGCTTCTAAAGGAAGTGGCTTGAGGGAGCCACTTATATGTTGGGGAGGCAACGCCTCCCGCAACATATCTTTTTCGTAATCAGATTAGAACTGGTAAACCAGACCTACAGCAACGATGTCATCAGTGCTTACACCGAGTGCTTTAGGGAAGGTGCGAATAAGCGGGGAAATTCTTCTCGGCTGACTCAGTCATTTCAT  >  NC_000913/574345‑574642
                                                                                                                                                      |                                                                                                                                                    
cAATTGTTTACTAAAAATTATTAAAAATGGGGCGTTGAGACGCCCCCAAAAATAAAGGGTAATATATAACAGAAGGTTTATATAGTTAGAAGCAAGGTTGTGCTTCTAAAGGAAGTGGCTTGAGGGAGCCCCTTATATGTTGGGGAGGCa                                                                                                                                                      >  2:1323051/1‑150 (MQ=255)
cAATTGTTTACTAAAAATTATTAAAAATGGGGCGTTGAGACGCCCCCAAAAATAAAGGGTAATATATAACAGAAGGTTTATATAGTTAGAAGCAAGGTTGTGCTTCTAAAGGAAGTGGCTTGAGGGAGCCACTTATATGTTGGGGAGGCa                                                                                                                                                      >  2:1387600/1‑150 (MQ=255)
cAATTGTTTACTAAAAATTATTAAAAATGGGGCGTTGAGACGCCCCCAAAAATAAAGGGTAATATATAACAGAAGGTTTATATAGTTAGAAGCAAGGTTGTGCTTCTAAAGGAAGTGGCTTGAGGGAGCCACTTATATGTTGGGGAGGCa                                                                                                                                                      >  2:46856/1‑150 (MQ=255)
cAATTGTTTACTAAAAATTATTAAAAATGGGGCGTTGAGACGCCCCCAAAAATAAAGGGTAATATATAACAGAAGGTTTATATAGTTAGAAGCAAGGTTGTGCTTCTAAAGGAAGTGGCTTGAGGGAGCCACTTATATGTTGGGGAGGCa                                                                                                                                                      >  1:1213155/1‑150 (MQ=255)
     gTTTACTAAAAATTATTAAAAATGGGGCGTTGAGACGCCCCCAAAAATAAAGGGTAATATATAACAGAAGGTTTATATAGTTAGAAGCAAGGTTGTGCTTCTAAAGGAAGTGGCTTGAGGGAGCCACTTATATGTTGGGGAGGCAACGcc                                                                                                                                                 >  2:1210421/1‑150 (MQ=255)
      tttACTAAAAATTATTAAAAATGGGGCGTTGAGACGCCCCCAAAAATAAAGGGTAATATATAACAGAAGGTTTATATAGTTAGAAGCAAGGTTGTGCTTCTAAAGGAAGTGGCTTGAGGGAGCCACTTATATGTTGGGGAGGCAACGCCt                                                                                                                                                <  1:782983/150‑1 (MQ=255)
      tttACTAAAAATTATTAAAAATGGGGCGTTGAGACGCCCCCAAAAATAAAGGGTAATATATAACAGAAGGTTTATATAGTTAGAAGCAAGGTTGTGCTTCTAAAGGAAGTGGCTTGAGGGAGCCACTTATATGTTGGGGAGGCAACGCCt                                                                                                                                                >  2:780356/1‑150 (MQ=255)
         aCTAAAAATTATTAAAAATGGGGCGTTGAGACGCCCCCAAAAATAAAGGGTAATATATAACAGAAGGTTTATATAGTTAGAAGCAAGGTTGTGCTTCTAAAGGAAGTGGCTTGAGGGAGCCACTTATATGTTGGGGAGGCAACGCCTccc                                                                                                                                             <  2:262259/150‑1 (MQ=255)
          cTAAAAATTATTAAAAATGGGGCGTTGAGACGCCCCCAAAAATAAAGGGTAATATATAACAGAAGGTTTATATAGTTAGAAGCAAGGTTGTGCTTCTAAAGGAAGTGGCTTGAGGGAGCCACTTATATGTTGGGGAGGCAACGCCTCCCg                                                                                                                                            <  1:430290/150‑1 (MQ=255)
            aaaaaTTATTAAAAATGGGGCGTTGAGACGCCCCCAAAAATAAAGGGTAATATATAACAGAAGGTTTATATAGTTAGAAGCAAGGTTGTGCTTCTAAAGGAAGTGGCTTGAGGGAGCCACTTATATGTTGGGGAGGCAACGCCTCCCGCa                                                                                                                                          <  2:148385/150‑1 (MQ=255)
            aaaaaTTATTAAAAATGGGGCGTTGAGACGCCCCCAAAAATAAAGGGTAATATATAACAGAAGGTTTATATAGTTAGAAGCAAGGTTGTGCTTCTAAAGGAAGTGGCTTGAGGGAGCCACTTATATGTTGGGGAGGCAACGCCTCCCGCa                                                                                                                                          <  1:767180/150‑1 (MQ=255)
            aaaaaTTATTAAAAATGGGGCGTTGAGACGCCCCCAAAAATAAAGGGTAATATATAACAGAAGGTTTATAAAGTTAGAAGCAAGGTTGTGCTTCTAAAGGAAGTGGCTTGAGGTAGCCACTTATATGTTGGGGAGGCAACGCCTCCCGCa                                                                                                                                          >  1:1492031/1‑150 (MQ=255)
              aaaTTATTAAAACTGGGGCGGTGAGACGCCCCCAAAAATAAAGGGTAATATATAACAGAAGGTTTATATAGTTAGAAGCACGGTTGTGCTTCTAAAGGAAGTGGCTTGAGGGAGCCACTTATATGGTGCGGGAGGCGTTGCCTccccaac                                                                                                                                         <  2:865807/150‑5 (MQ=255)
              aaaTTATTAAAAATGGGGCGTTGAGACGCCTCCAAAAATAAAGGGTAATATATAACAGAAGGTTTATATAGTTAGAAGCAAGGTTGTGCTTCTAAAGGAAGTGGCTTGAGGGAGCCACTTATATGTTGGGGAGGCAACGCCTCCCGCAAc                                                                                                                                        >  1:1257759/1‑150 (MQ=255)
              aaaTTATTAAAAATGGGGCGTTGAGACGCCCCCAAAAATAAAGGGTAATATATAACAGAAGGTTTATATAGTTAGAAGCAAGGTTGTGCTTCTAAAGGAAGTGGCTTGAGGGAGCCACTTATATGTTGGGGAGGCAACGCCTCCCGCAAc                                                                                                                                        <  1:1515225/150‑1 (MQ=255)
                attattAAAAATGGGGCGTTGAGACGCCCCCAAAAATAAAGGGTAATATATAACAGAAGGTTTATATAGTTAGAAGCAAGGTTGTGCTTCTAAAGGAAGTGGCTTGAGGGAGCCACTTATATGTTGCGGGAGGCGTTGCCTCCCCAACat                                                                                                                                      <  2:1295205/150‑1 (MQ=255)
                 ttattaAAAATGGGGCGTTGAGACGCCCCCAAAAATAAAGGGTAATATATAACAGAAGGTTTATATAGTTAGAAGCAAGGTTGTGCTTCTAAAGGAAGTGGCTTGAGGGAGCCACTTATATGTTGGGGAGGCAACGCCTCCCGCAACata                                                                                                                                     <  1:1569096/150‑1 (MQ=255)
                    ttaAAAATGGGGCGTTGAGACGCCCCCAAAAATAAAGGGTAATATATAACAGAAGGTTTATATAGTTAGAAGCAAGGTTGTGCTTCTAAAGGAAGTGGCTTGAGGGAGCCACTTATATGTTGGGGAGGCAACGCCTCCCGCAACATATCt                                                                                                                                  >  1:1241736/1‑150 (MQ=255)
                     taAAAATGGGGCGTTGAGACGCCCCCAAAAATAAAGGGTAATATATAACAGAAGGTTTATATAGTTAGAAGCAAGGTTGTGCTTCTAAAGGAAGTGGCTTGAGGGAGCCACTTATATGTTGCGGGAGGCGTTGCCTCCCCAACatataag                                                                                                                                 <  2:51082/150‑4 (MQ=38)
                     taAAAATGGGGCGTTGAGACGCCCCCAAAAATAAAGGGTAATATATAACAGAAGGTTTATATAGTTAGAAGCAAGGTTGTGCTTCTAAAGGAAGTGGCTTGAGGGAGCCACTTATATGTTGGGGAGGCAACGCCTCCCGCAACATATCtt                                                                                                                                 >  2:1268214/1‑150 (MQ=255)
                      aaaaaTGGGGCGTTGAGACGCCCCCAAAAATAAAGGGTAATATATAACAGAAGGTTTATATAGTTAGAAGCAAGGTTGTGCTTCTAAAGGAAGTGGCTTGAGGGAGCCACTTATATGTTGGGGAGGCAACGCCTCCCGCAACATATCttt                                                                                                                                >  2:1228494/1‑150 (MQ=255)
                       aaaaTGGGGCGTTGAGACGCCCCCAAAAATAAAGGGTAATATATAACAGAAGGTTTATATAGTTAGAAGCAAGGTTGTGCTTCTAAAGGAAGTGGCTTGAAGGAGCCACTTATATGTTGCGGGAGGCGTTGCCTCCCCAACatataagtg                                                                                                                               <  2:724589/150‑6 (MQ=38)
                        aaaTGGGGCGTTGAGACGCCCCCAAAAATAAAGGGTAATATATAACAGAAGGTTTATATAGTTAGAAGCAAGGTTGTGCTTCTAAAGGAAGTGGCTTGAGGGAGCCACTTATATGTTGCGGGAGGCGTTGCCTCCCCAACatataagtgg                                                                                                                              <  1:805832/150‑7 (MQ=38)
                          aTGGGGCGTTGAGACGCCCCCAAAAATAAAGTGTAATATATAACAGAAGGTTTATATAGTTAGAAGCAAGGTTGTGCTTCTAAAGGAAGTGGCTTGAGGGAGCCACTTATATGTTGGGGAGGCAACGCCTCCCGCAACATATCTTTTTCg                                                                                                                            <  2:853659/150‑1 (MQ=255)
                           tGGGGCGTTGAGACGCCCCCAAAAATAAAGGGTAATATATAACAGAAGGTTTATATAGTTAGAAGCAAGGTTGTGCTTCTAAAGGAAGTGGCTTGAGGGAGCCACTTATATGTTGGGGAGGCAACGCCTCCCGCAACATATCTTTTTCGt                                                                                                                           <  1:254558/150‑1 (MQ=255)
                                cGTTGAGACGCCCCCAAAAATAAAGGGTAATATATAACAGAAGGTTTATATAGTTAGAAGCAAGGTTGTGCTTCTAAAGGAAGTGGCTTGAGGGAGCCACTTATATGTTGGGGAGGCAACGCCTCCCGCAACATATCTTTTTCGTAATCa                                                                                                                      <  1:1387600/150‑1 (MQ=255)
                                 gTTGAGACGCCCCCAAAAATAAAGGGTAATATATAACAGAAGGTTTATATAGTTAGAAGCAAGGTTGTGCTTCTAAAGGAAGTGGCTTGAGGGAGCCACTTATATGTTGGGGAGGCAACGCCTCCCGCAACATATCTTTTTCGTAATCAg                                                                                                                     <  1:179560/150‑1 (MQ=255)
                                         gCCCCCAAAAATAAAGGGTAATATATAACAGAAGGTTTATATAGTTAGAAGCAAGGTTGTGCTTCTAAAGGAAGTGGCTTGAGGGAGCCACTTATATGTTGGGGAGGCAACGCCTCCCGCAACATATCTTTTTCGTAATCAGATTAGAAc                                                                                                             <  2:1105822/150‑1 (MQ=255)
                                          cccccAAAAATAAAGGGTAATATATAACAGAAGGTTTATATAGTTAGAAGCAAGGTTGTGCTTCTAAAGGAAGTGGCTTGAGGGAGCCACTTATATGTTGGGGAGGCAACGCCTCCCGCAACATATCTTTTTCGTAATCAGATTAGAACt                                                                                                            <  1:881620/150‑1 (MQ=255)
                                          cccccAAAAATAAAGGGTAATATATAACAGAAGGTTTATATAGTTAGAAGCAAGGTTGTGCTTCTAAAGGAAGTGGCTTGAGGGAGCCACTTATATGTTGGGGAGGCAACGCCTCCCGCAACATATCTTTTTCGTAATCAGATTAGAACt                                                                                                            <  2:144960/150‑1 (MQ=255)
                                            cccAAAAATAAAGGGTAATATATAACAGAAGGTTTATATAGTTAGAAGCAAGGTTGTGCTTCTAAAGGAAGTGGCTTGAGGGAGCCACTTATATGTTGGGGAGGCAACGCCTCCCGCAACATATCTTTTTCGTAATCAGATTAGAACTgg                                                                                                          >  2:446378/1‑150 (MQ=255)
                                             ccAAAAATAAAGGGTAATATATAACAGAAGGTTTATATAGTTAGAAGCAAGGTTGTGCTTCTAAAGGAAGTGGCTTGAGGGAGCCACTTATATGTTGGGGAGGCAACGCCTCCCGCAACATATCTTTTTCGTAATCAGATTAGAACTGGt                                                                                                         <  2:289592/150‑1 (MQ=255)
                                             acAAAAATAAAGGGAAATATATAAAAGAAGGTTTATATAGTTATAAGCAAGGTTGTGGTTCTAAAGGAAGTGGATTGAGGGAGCCACTTATATGTTGGGGAGGCAACGCCTAACGCAACATATCTTTTTAGTAATCAGATTTGAAATGGt                                                                                                         >  2:1325460/2‑150 (MQ=255)
                                                aaaaTAAAGGGTAATATATAACAGAAGGTTTATATAGTTAGAAGCAAGGTTGTGCTTCTAAAGGAAGTGGCTTGAGGGAGCCACTTATATGTTGGGGCGGCAACGCCTCCCGCAACATATCTTTTTCGTAATCAGATTAGAACTGGTaaa                                                                                                      >  2:927579/1‑150 (MQ=255)
                                                     aaaGGGTAATATATAACAGAAGGTTTATATAGTTAGAAGCAAGGTTGTGCTTCTAAAGGAAGTGGCTTGAGGGAGCCACTTATATGTTGGGGAGGCAACGCCTCCCGCAACATATCTTTTTCGTAATCAGATTAGAACTGGTAAACCAGa                                                                                                 <  1:424187/150‑1 (MQ=255)
                                                     aaaGGGTAATATATAACAGAAGGTTTATATAGTTAGAAGCAAGGTTGTGCTTCTAAAGGAAGTGGCTTGAGGGAGCCACTTATATGTTGGGGAGGCAACGCCTCCCGCAACATATCTTTTTCGTAATCAGATTAGAACTGGTAAACCAGa                                                                                                 <  1:1103917/150‑1 (MQ=255)
                                                      aaGTGTATTATATAACATAATTTTTATATAGTTAGAAGCTAGGTTTTTCTTCTAAAGTAAGTGGCTTGAGGTAGACACTTATATTTTGGGGAGGCATCGCCTCCCGCAACATATCTTTTTCGTAATCAGATTTGAAGTTGTATACCTGAc                                                                                                <  2:1492031/150‑1 (MQ=255)
                                                      aaGGGTAATATATAACAGAAGGTTTATATAGTTAGAAGCAAGGTTGTGCTTCTAAAGGAAGTGGCTTGAGGGAGCCACTTATATGTTGGGGAGGCAACGCCTCCCGCAACATATCTTTTTCGTAATCAGATTAGAACTGGTAAACCAGAc                                                                                                <  1:331282/150‑1 (MQ=255)
                                                       aGGGTAATATATAACAGAAGGTTTATATAGTTAGAAGCAAGGTTGTGCTTCTAAAGGAAGTGGCTTGAGGGAGCCACTTATATGTTGGGGAGGCAACGCCTCCCGCAACATATCTTTTTCGTAATCAGATTAGAACTGGTAAACCAGAcc                                                                                               <  2:277255/150‑1 (MQ=255)
                                                       aGGGTAATATATAACAGAAGGTTTATATAGTTAGAAGCAAGGTTGTGCTTCTAAAGGAAGTGGCTTGAGGGAGCCACTTATATGTTGGGGAGGCAACGCCTCCCGCAACATATCTTTTTCGTAATCAGATTAGAACTGGTAAACCAGAcc                                                                                               >  1:1151649/1‑150 (MQ=255)
                                                        gggTAATATATAACAGAAGGTTTATATAGTTAGAAGCAAGGTTGTGCTTCTAAAGGAAGTGGCTTGAGGGAGCCACTTATATGTTGGGGAGGCAACGCCTCCCGCAACATATCTTTTTCGTAATCAGATTAGAACTGGTAAACCAGACCt                                                                                              >  2:484454/1‑150 (MQ=255)
                                                         ggTAATATATAACAGAAGGTTTATATAGTTAGAAGCAAGGTTGTGCTTCTAAAGGAAGTGGCTTGAGGGAGCCACTTATATGTTGGGGAGGCAACGCCTCCCGCAACATATCTTTTTCGTAATCAGATTAGAACTGGTAAACCAGACCTa                                                                                             >  2:546193/1‑150 (MQ=255)
                                                          gTAATATATAACAGAAGGTTTATATAGTTAGAAGCAAGGTTGTGCTTCTAAAGGAAGTGGCTTGAGGGAGCCACTTATATGTTGGGGAGGCAACGCCTCCCGCAACATATCTTTTTCGTAATCAGATTAGAACTGGTAAACCAGACCTAc                                                                                            <  2:822508/150‑1 (MQ=255)
                                                           tAATATATAACAGAAGGTTTATATAGTTAGAAGCAAGGTTGTGCTTCTAAAGGAAGTGGCTTGAGGGAGCCACTTATATGTTGGGGAGGCAACGCCTCCCGCAACATATCTTTTTCGTAATCAGATTAGAACTGGTAAACCAGACCTACa                                                                                           <  1:835029/150‑1 (MQ=255)
                                                             atatatAACAGAAGGTTTATATAGTTAGAAGCAAGGTTGTGCTTCTAAAGGAAGTGGCTTGAGGGAGCCACTTATATGTTGGGGAGGCAACGCCTCCCGCAACATATCTTTTTCGTAATCAGATTAGAAc                                                                                                             >  1:31968/1‑130 (MQ=255)
                                                             atatatAACAGAAGGTTTATATAGTTAGAAGCAAGGTTGTGCTTCTAAAGGAAGTGGCTTGAGGGAGCCACTTATATGTTGGGGAGGCAACGCCTCCCGCAACATATCTTTTTCGTAATCAGATTAGAAc                                                                                                             <  2:31968/130‑1 (MQ=255)
                                                                tataACATAAGTTTTATATATTTAGAAGCAATGTTGTGCTTCTAAAGTAAGTGGATTGAGGTAGACACTTATATGTTGGGGATGCAAAGCATCACGCAACATATCTTTTTAGTAATAATATTAGAACTGGTAAACCAGACCTACAGAAAc                                                                                      >  2:512164/1‑150 (MQ=255)
                                                                tataACAGAAGGTTTATATAGTTAGAAGCAAGGGTGTGCTTCTAAAGGAAGTGGCTTGAGGGCGCCACTTATATGTTGGGGAGGCAACGCCTCCCGCCACATATCTATTTCGTAATCAGATAAGAACTGGTAACCCAGCCCTACAGCAAc                                                                                      >  2:1670606/1‑150 (MQ=255)
                                                                 ataACAGAAGGTTTATATAGTTAGAAGCAAGGTTGTGCTTCTAAAGGAAGTGGCTTGAGGGAGCCACTTATATGTTGGGGAGGCAACGCCTCCCGCAACATATCTTTTTCGTAATCAGATTAGAACTGGTAAACCAGACCTACAGCAACg                                                                                     <  1:46856/150‑1 (MQ=255)
                                                                  taACAGAAGGTTTATATAGTTAGAAGCAAGGTTGTGCTTCTAAAGGAAGTGGCTTGAGGGAGCCACTTATATGTTGGGGAGGCAACGCCTCCCGCAACATATCTTTTTCGTAATCAGATTAGAACTGGTAAACCAGACCTACAGCAACGa                                                                                    <  2:416229/150‑1 (MQ=255)
                                                                  taACAGAAGGTTTATATAGTTAGAAGCAAGGTTGTGCTTCTAAAGGAAGTGGCTTGAGGGAGCCACTTATATGTTGGGGAGGCAACGCCTCCCGCAACATATCTTTTTCGTAATCAGATTAGAACTGGTAAACCAGACCTACAGCAACGa                                                                                    >  1:593025/1‑150 (MQ=255)
                                                                    aCAGAAGGTTTATATAGTTAGAAGCAAGGTTGTGCTTCTAAAGGAAGTGGCTTGAGGGAGCCACTTATATGTTGGGGAGGCAACGCCTCCCGCAACATATCTTTTTCGTAATCAGATTAGAACTGGTAAACCAGACCTACAGCAACGATg                                                                                  <  1:1100334/150‑1 (MQ=255)
                                                                    aCAGAAGGTTTATATAGTTAGAAGCAAGGTTGTGCTTCTAAAGGAAGTGGCTTGAGGGAGCCACTTATATGTTGGGGAGGCAACGCCTCCCGCAACATATCTTTTTCGTAATCAGATTAGAACTGGTAAACCAGACCTACAGCAACGATg                                                                                  >  1:1353484/1‑150 (MQ=255)
                                                                     cAGAAGGTTTATATAGTTAGAAGCAAGGTTGTGCTTCTAAAGGAAGTGGCTTGAGGGAGCCACTTATATGTTGGGGAGGCAACGCCTCCCGCAACATATCTTTTTCGTAATCAGATTAGAACTGGTAAACCAGACCTACAGCAACGATGt                                                                                 <  2:1257759/150‑1 (MQ=255)
                                                                           gTTTATATAGTTAGAAGCAAGGTTGTGCTTCTAAAGGAAGTGGCTTGAGGGAGCCACTTATATGTTGGGGAGGCAACGCCTCCCGCAACATATCTTTTTCGTAATCAGATTAGAACTGGTAAACCAGACCTACAGCAACGATGTCATCAg                                                                           <  2:1213155/150‑1 (MQ=255)
                                                                             ttATATAGTTAGAAGCAAGGTTGTGCTTCTAAAGGAAGTGGCTTGAGGGAGCCACTTATATGTTGGGGAGGCAACGCCTCCCGCAACATATCTTTTTCGTAATCAGATTAGAACTGGTAAACCAGACCTACAGCAACGATGTCATCAGTg                                                                         <  2:1534037/150‑1 (MQ=255)
                                                                              tatataGTTAGAAGCAAGGTTGTGCTTCTAAAGGAAGTGGCTTGAGGGAGCCACTTATATGTTGGGGAGGCAACGCCTCCCGCAACATATCTTTTTCGTAATCAGATTAGAACTGGTAAACCAGACCTACAGCAACGATGTCATCAGTGc                                                                        >  2:1123929/1‑150 (MQ=255)
                                                                              tatataGTTAGAAGCAAGGTTGTGCTTCTAAAGGAAGTGGCTTGAGGGAGCCACTTATATGTTGGGGAGGCAACGCCTCCCGCAACATATCTTTTTCGTAATCAGATTAGAACTGGTAAACCAGACCTACAGCAACGATGTCATCAGTGc                                                                        <  1:1006655/150‑1 (MQ=255)
                                                                                   aGTTAGAAGCAAGGTTGTGCTTCTAAAGGAAGTGGCTTGAGGGAGCCACTTATATGTTGGGGAGGCAACGCCTCCCGCAACATATCTTTTTCGTAATCAGATTAGAACTGGTAAACCAGACCTACAGCAACGATGTCATCAGTGCTTaca                                                                   <  1:1449394/150‑1 (MQ=255)
                                                                                    gTTAGAAGCAAGGTTGTGCTTCTAAAGGAAGGGGCTTGAGGGAGCCACTTATATGTTGGGGAGGCAACGCCTCCCGCAACATATCTTTTTCGTAATCAGATTAGAACTGGTAAACCAGACCTACAGCACTGATGACATCGTTGCtgtagg                                                                  <  2:1049874/150‑6 (MQ=255)
                                                                                      tAGAAGCAAGGTTGTGCTTCTAAAGGAAGTGGCTTGAGGGAGCCACTTATATGTTGGGGAGGCAACGCCTCCCGCAACATATCTTTTTCGTAATCAGATTAGAACTGGTAAACCAGACCTACAGCAACGATGTCATCAGTGCTTACACCg                                                                <  1:1210421/150‑1 (MQ=255)
                                                                                         aaGCAAGGTTGTGCTTCTAAAGGAAGTGGCTTGAGGGAGCCACTTATATGTTGGGGAGGCAACGCCTCCCGCAACATATCTTTTTCGTAATCAGATTAGAACTGGTAAACCAGACCTACAGCAACGATGTCATCAGTGCTTACACCGAGt                                                             <  1:780356/150‑1 (MQ=255)
                                                                                          aGCAAGGTTGTGCTTCTAAAGGAAGTGGCTTGAGGGAGCCACTTATATGTTGGGGAGGCAACGCCTCCCGCAACATATCTTTTTCGTAATCAGATTAGAACTGGTAAACCAGACCTACAGCAACGATGTCATCAGTGCTTACACCGAGTg                                                            <  1:1630115/150‑1 (MQ=255)
                                                                                          aGCAAGGTTGTGCTTCTAAAGGAAGTGGCTTGAGGGAGCCACTTATATGTTGGGGAGGCAACGCCTCCCGCAACATATCTTTTTCGTAATCAGATTAGAACTGGTAAACCAGACCTACAGCAACGATGTCATCAGTGCTTACACCGAGTg                                                            >  2:551665/1‑150 (MQ=255)
                                                                                               ggTTGTGCTTCTAAAGGAAGTGGCTTGAGGGAGCCACTTATATGTTGGGGAGGCAACGCCTCCCGCAACATATCTTTTTCGTAATCAGATTAGAACTGGTAAACCAGACCTACAGCAACGATGTCATCAGTGCTTACACCGAGTGCTTTa                                                       <  2:1241736/150‑1 (MQ=255)
                                                                                                gTTGTGCTTCTAAAGGAAGTGGCTTGAGGGAGCCACTTATATGTTGGGGAGGCAACGCCTCCCGCAACATATCTTTTTCGTAATCAGATTAGAACTGGTAAACCAGACCTACAGCAACGATGTCATCAGTGCTTACACCGAGTGCTTTAg                                                      <  1:888596/150‑1 (MQ=255)
                                                                                                   gtgCTTCTAAAGGAAGTGGCTTGAGGGAGCCACTTATATGTTGGGGAGGCAACGCCTCCCGCAACATATCTTTTTCGTAATCAGATTAGAACTGGTAAACCAGACCTACAGCAACGATGTCATCAGTGCTTACACCGAGTGCTTTAGGGa                                                   >  2:18881/1‑150 (MQ=255)
                                                                                                    tgCTTCTAAAGGAAGTGGCTTGAGGGACCCACTTATATGTTGGGGAGGCAACCCCTCCCGCAACATATATTTTTCGTAATCAGATTAGAACTGGTAAACCAGCCCTACAGCAACGATGTCATCAGTGCTTACACCGAGTGCTTTAGGGaa                                                  >  2:906473/1‑150 (MQ=255)
                                                                                                     gCTTCTAAAGGAAGTGGCTTGAGGGAGCCACTTATATGTTGGGGAGGCAACGCCTCCCGCAACATATCTTTTTCGTAATCAGATTAGAACTGGTAAACCAGACCTACAGCAACGATGTCATCAGTGCTTACACCGAGTGCTTTAGGGAAg                                                 <  2:359094/150‑1 (MQ=255)
                                                                                                      cTTCTAAAGGAAGTGGCTTGAGGGAGCCACTTATATGTTGGGGAGGCAACGCCTCCCGCAACATATCTTTTTCGTAATCAGATTAGAACTGGTAAACCAGACCTACAGCAACGATGTCATCAGTGCTTACACCGAGTGCTTTAGGGAAgg                                                >  1:4837/1‑150 (MQ=255)
                                                                                                        tCTAAAGGAAGTGGCTTGAGGGAGCCACTTATATGTTGGGGAGGCAACGCCTCCCGCAACATATCTTTTTCGTAATCAGATTAGAACTGGTAAACCAGACCTACAGCAACGATGTCATCAGTGCTTACACCGAGTGCTTTAGGGAAGGTg                                              >  1:1622045/1‑150 (MQ=255)
                                                                                                        tCTAAAGGAAGTGGCTTGAGGGAGCCACTTATATGTTGGGGAGGCAACGCCTCCCGCAACATATCTTTTTCGTAATCAGATTAGAACTGGTAAACCAGACCTACAGCAACGATGTCATCAGTGCTTACACCGAGTGCTTTAGGGAAGGTg                                              <  1:497252/150‑1 (MQ=255)
                                                                                                        tCTAAAGGAAGTGGCTTGAGGGAGCCACTTATATGTTGGGGAGGCAACGCCTCCCGCAACATATCTTTTTCGTAATCAGATTAGAACTGGTAAACCAGACCTACAGCAACGATGTCATCAGTGCTTACACCGAGTGCTTTAGGGAAGGTg                                              <  1:546946/150‑1 (MQ=255)
                                                                                                           aaaGGAAGTGGCTTGAGGGAGCCACTTATATGTTGGGGAGGCAACGCCTCCCGCAACATATCTTTTTCGTAATCAGATTAGAACTGGTAAACCAGACCTACAGCAACGATGTCATCAGTGCTTACACCGAGTGCTTTAgg                                                     <  1:1584930/140‑1 (MQ=255)
                                                                                                           aaaGGAAGTGGCTTGAGGGAGCCACTTATATGTTGGGGAGGCAACGCCTCCCGCAACATATCTTTTTCGTAATCAGATTAGAACTGGTAAACCAGACCTACAGCAACGATGTCATCAGTGCTTACACCGAGTGCTTTAgg                                                     >  2:1584930/1‑140 (MQ=255)
                                                                                                             aGGAAGTGGCTTGAGGGAGCCACTTATATGTTGGGGAGGCAACGCCTCCCGCAACATATCTTTTTCGTAATCAGATTAGAACTGGTAAACCAGACCTACAGCAACGATGTCATCAGTGCTTACACCGAGTGCTTTAGGGAAGGTGCGAAt                                         <  2:1328709/150‑1 (MQ=255)
                                                                                                                   tGGCTTGAGGGAGCCACTTATATGTTGGGGAGGCAACGCCTCCCGCAACATATCTTTTTCGTAATCAGATTAGAACTGGTAAACCAGACCTACAGCAACGATGTCATCAGTGCTTACACCGAGTGCTTTAGGGAAGGTGCGAATAAGCgg                                   <  2:1353484/150‑1 (MQ=255)
                                                                                                                    ggCTTGAGGGAGCCACTTATATGTTGGGGAGGCAACGCCTCCCGCAACATATCTTTTTCGTAATCAGATTAGAACTGGTAAACCAGACCTACAGCAACGATGTCATCAGTGCTTACACCGAGTGCTTTAGGGAAGGTGCGAATAAGCggg                                  <  1:85819/150‑1 (MQ=255)
                                                                                                                      cTTGAGGGAGCCACTTATATGTTGGGGAGGCAACGCCTCCCGCAACATATCTTTTTCGTAATCAGATTAGAACTGGTAAACCAGACCTACAGCAACGATGTCATCAGTGCTTACACCGAGTGCTTTAGGGAAGGTGCGAATAAGCGGGGa                                <  2:1151649/150‑1 (MQ=255)
                                                                                                                         gAGGGAGCCACTTATATGTTGGGGAGGCAACGCCTCCCGCAACATATCTTTTTCGTAATCAGATTAGAACTGGTAAACCAGACCTACAGCAACGATGTCATCAGTGCTTACACCGAGTGCTTTAGGGAAGGTGCGAATAAGCGGGGAAAt                             >  1:1268979/1‑150 (MQ=255)
                                                                                                                         gAGGGAGCCACTTATATGTTGGGGAGGCAACGCCTCCCGCAACATATCTTTTTCGTAATCAGATTAGAACTGGTAAACCAGACCTACAGCAACGATGTCATCAGTGCTTACACCGAGTGCCTTAGGGAAGGTGCGAATAAGCGGGGAAAt                             >  1:1361063/1‑150 (MQ=255)
                                                                                                                          aGGGAGCCACTTATATGTTGGGGAGGCAACGCCTCCCGCAACATATCTTTTTCGTAATCAGATTAGAACTGGTAAACCAGACCTACAGCAACGATGTCATCAGTGCTTACACCGAGTGCTTTAGGGAAGGTGCGAATAAGCGGGGAAAtt                            >  2:663019/1‑150 (MQ=255)
                                                                                                                              aGCCACTTATATGTTGGGGAGGCAACGCCTCCCGCAACATATCTTTTTCGTAATCAGATTAGAACTGGTAAACCAGACCTACAGCAACGATGTCATCAGTGCTTACACCGAGTGCTTTAGGGAAGGTGCGAATAAGCGGGGAAAttcttc                        >  1:1609883/1‑150 (MQ=255)
                                                                                                                                   cTTATATGTTGGGGAGGCAACGCCTCCCGCAACATATCTTTTTCGTAATCAGATTAGAACTGGTAAACCAGACCTACAGCAACGATGTCATCAGTGCTTACACCGAGTGCTTTAGGGAAGGTGCGAATAAGCGGGGAAAttct                          <  1:1051426/143‑1 (MQ=255)
                                                                                                                                   cTTATATGTTGGGGAGGCAACGCCTCCCGCAACATATCTTTTTCGTAATCAGATTAGAACTGGTAAACCAGACCTACAGCAACGATGTCATCAGTGCTTACACCGAGTGCTTTAGGGAAGGTGCGAATAAGCGGGGAAATTCTTCTCGGc                   <  2:207548/150‑1 (MQ=255)
                                                                                                                                   cTTATATGTTGGGGAGGCAACGCCTCCCGCAACATATCTTTTTCGTAATCAGATTAGAACTGGTAAACCAGACCTACAGCAACGATGTCATCAGTGCTTACACCGAGTGCTTTAGGGAAGGTGCGAATAAGCGGGGAAATTCTTCTCGGc                   >  1:84462/1‑150 (MQ=255)
                                                                                                                                   cTTATATGTTGGGGAGGCAACGCCTCCCGCAACATATCTTTTTCGTAATCAGATTAGAACTGGTAAACCAGACCTACAGCAACGATGTCATCAGTGCTTACACCGAGTGCTTTAGGGAAGGTGCGAATAAGCGGGGAAATTCTTCTCGGc                   >  1:219887/1‑150 (MQ=255)
                                                                                                                                   cTTATATGTTGGGGAGGCAACGCCTCCCGCAACATATCTTTTTCGTAATCAGATTAGAACTGGTAAACCAGACCTACAGCAACGATGTCATAAGTGCTTACACCGAGTGCTTTAGGGAAGGTGCGAATAAGCGGGGAAAttct                          >  2:1051426/1‑143 (MQ=255)
                                                                                                                                   aagatatGTTGCGGGAGGCGTTGCCTCCCGCAACATATCTTTTTCGTAATCAGATTAGAACTGGTAAACCAGACCTACAGCAACGATGTCATCAGTGCTTACACCGAGTGCTTTAGGGAAGGTGCGAATAAGCGGGGAAATTCTTCTCgg                    >  2:1605711/4‑150 (MQ=39)
                                                                                                                                     tataTGTTGGGGAGGCAACGCCTCCCGCAACATATCTTTTTCGTAATCAGATTAGAACTGGTAAACCAGACCTACAGCAACGATGTCATCAGTGCTTACACCGAGTGCTTTAGGGAAGGTGCGAATAAGCGGGGAAATTCTTCTCGGCTg                 <  1:18881/150‑1 (MQ=255)
                                                                                                                                        atGTTGGGGAGGCAACGCCTCCCGCAACATATCTTTTTCGTAATCAGATTAGAACTGGTAAACCAGACCTACAGCAACGATGTCATCAGTGCTTACACCGAGTGCTTTAGGGAAGGTGCGAATAAGCGGGGAAATTCTTCTCGGCTGACt              >  1:1173424/1‑150 (MQ=255)
                                                                                                                                         tGTTGGGGAGGCAACGCCTCCCGCAACATATCTTTTTCGTAATCAGATTAGAACTGGTAAACCAGACCTACAGCAACGATGTCATCAGTGCTTACACCGAGTGCTTTAGGGAAGGTGCGAATAAGCGGGGAAATTCTTCTCGGCTGACTc             <  1:1123929/150‑1 (MQ=255)
                                                                                                                                          gTTGGGGAGGCAACGCCTCCCGCAACATATCTTTTTCGTAATCAGATTAGAACTGGTAAACCAGACCTACAGCAACGATGTCATCAGTGCTTACACCGAGTGCTTTAGGGAAGGTGCGAATAAGCGGGGAAATTCTTCTCGGCTGACTCa            >  2:988642/1‑150 (MQ=255)
                                                                                                                                            tGGGGAGGCAACGCCTCCCGCAACATATCTTTTTCGTAATCAGATTAGAACTGGTAAACCAGACCTACAGCAACGATGTCATCAGTGCTTACACCGAGTGCTTTAGGGAAGGTGCGAATAAGCGGGGAAATTCTTCTCGGCTGACTCAGt          <  1:1575265/150‑1 (MQ=255)
                                                                                                                                             gcgggAGGCGTTGCCTCCCGCAACATATCTTTTTCGTAATCAGATTAGAACTGGTAAACCAGACCTACAGCAACGATGTCATCAGTGCTTACACCGAGTGCTTTAGGGAAGGTGCGAATAAGCGGGGAAATTCTTCTCgg                    <  1:893256/138‑1 (MQ=255)
                                                                                                                                             gcgggAGGCGTTGCCTCCCGCAACATATCTTTTTCGTAATCAGATTAGAACTGGTAAACCAGACCTACAGCAACGATGTCATCAGTGCTTACACCGAGTGCTTTAGGGAAGGTGCGAATAAGCGGGGAAATTCTTCTCgg                    >  2:893256/3‑140 (MQ=255)
                                                                                                                                               gggAGGCAACGCCTCCCGCAACATATCTTTTTCGTAATCAGATTAGAACTGGTAAACCAGACCTACAGCAACGATGTCATCAGTGCTTACACCGAGTGCTTTAGGGAAGGTGCGAATAAGCGGGGAAATTCTTCTCGGCTGACTCAGtca        <  1:484454/150‑1 (MQ=255)
                                                                                                                                                ggAGGCAACGCCTCCCGCAACATATCTTTTTCGTAATCAGATTAGAACTGGTAAACCAGACCTACAGCAACGATGTCATCAGTGCTTACACCGAGTGCTTTAGGGAAGGTGCGAATAAGCGGGGAAATTCTTCTCGGCTGACTCAGtcat       >  2:5860/1‑150 (MQ=255)
                                                                                                                                                     cAACGCCTCCCGCAACATATCTTTTTCGTAATCAGATTAGAACTGGTAAACCAGACCTACAGCAACGATGTCATCAGTGCTTACACCGAGTGCTTTAGGGAAGGTGCGAATAAGCGGGGAAATTCTTCTCGGCTGACTCAGtcatttcat  >  2:931382/1‑150 (MQ=255)
                                                                                                                                                     cAACGCCTCCCGCAACATATCTTTTTCGTAATCAGATTAGAACTGGTAAACCAGACCTACAGCAACGATGTCATCAGTGCTTACACCGAGTGCTTTAGGGAAGGTGCGAATAAGCGGGGAAATTCTTCTCGGCTGACTCAGtcatttcat  >  2:972166/1‑150 (MQ=255)
                                                                                                                                                      |                                                                                                                                                    
CAATTGTTTACTAAAAATTATTAAAAATGGGGCGTTGAGACGCCCCCAAAAATAAAGGGTAATATATAACAGAAGGTTTATATAGTTAGAAGCAAGGTTGTGCTTCTAAAGGAAGTGGCTTGAGGGAGCCACTTATATGTTGGGGAGGCAACGCCTCCCGCAACATATCTTTTTCGTAATCAGATTAGAACTGGTAAACCAGACCTACAGCAACGATGTCATCAGTGCTTACACCGAGTGCTTTAGGGAAGGTGCGAATAAGCGGGGAAATTCTTCTCGGCTGACTCAGTCATTTCAT  >  NC_000913/574345‑574642

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: