Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 2,757,662 (T)8→7 coding (939/957 nt) yfjH ← CP4‑57 prophage; uncharacterized protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,757,6550T.95.0% 585.7 / 24.4 141coding (946/957 nt)yfjHCP4‑57 prophage; uncharacterized protein
Reads supporting (aligned to +/- strand):  ref base T (6/1);  new base . (67/67);  total (73/68)
Fisher's exact test for biased strand distribution p-value = 1.17e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.42e-01

Only 100 of 157 total aligned reads displayed.

TGCTCCTTTTGCAACGTATCACTTGTCAATAAATAGAAAGTGTAGATGTATTAGGCGCAGCTATTTAGTCTATATGGCAAGTAGTAATTATTTACTTCTTTTTTTTCTACGAACAGCGCTATCGAAAGCTTGTTTTACATCTTCATTTTTTAAAGACCATTGATTTAAAATAGTATCAACCAAAGCATCTTG  >  NC_000913/2757557‑2757748
                                                                                                  |                                                                                             
  cTCCTTTTGCAACGTATCACTTGTCAATAAATAGAAAGTGTAGATGTATTAGGCGCAGCTATTTAGTCTATATGGCAAGTAGTAATTATTTACTTCtttt                                                                                            <  1:969924/100‑1 (MQ=255)
   tCCTTTTGCAACGTATCACTTGTCAATAAATAGAAAGTGTAGATGTATTAGGCGCAGCTATTTAGTCTATATGGCAAGTAGTAATTATTTACTTCttttt                                                                                           <  1:859716/100‑1 (MQ=255)
    ccTTTTGCAACGTATCACTTGTCAATAAATAGAAAGTGTAGATGTATTAGGCGCAGCTATTTAGTCTATATGGCAAGTAGTAATTATTTACTTCtttttt                                                                                          >  1:2202956/1‑100 (MQ=255)
      ttttGCAACGTATCACTTGTAAATAAATAGAAAGTGTAGATGTATTAGGCGCAGCTATTTAGTCTATATGGCAAGTAGTAATTATTTACTTCtttttttc                                                                                        <  2:2087364/100‑2 (MQ=255)
       tttGCAACGTATCACTTGTCAATAAATAGAAAGTGTAGATGTATTAGGCGCAGCTATTTAGTCTATATGGCAAGTAGTAATTATTTACTTCtttttttct                                                                                       <  2:1145888/100‑3 (MQ=255)
       tttGCAACGTATCACTTGTCAATAAATAGAAAGTGTAGATGTATTAGGCGCAGCTATTTAGTCTATATGGCAAGTAGTAATTATTTACTTCtttttttct                                                                                       >  2:137730/1‑98 (MQ=255)
        ttGCAACGTATCACTTGTCAATAAATAGAAAGTGTAGATGTATTAGGCGCAGCTATTTAGTCTATATGGCAAGTAGTAATTATTTACTTCtttttttcta                                                                                      <  1:2522066/100‑4 (MQ=255)
        ttGCAACGTATCACTTGTCAATAAATAGAAAGTGTAGATGTATTAGGCGCAGCTATTTAGTCCATATGGCAAGTAGTAATTATTTACTTCtttttttcta                                                                                      <  1:1636688/100‑4 (MQ=255)
         tGCAACGTATCACTTGTCAATAAATAGAAAGTGTAGATGTATTAGGCGCAGCTATTTAGTCTATATGGCAAGTAGTAATTATTTACTTCtttttttctac                                                                                     <  2:1488023/100‑5 (MQ=255)
         tGCAACGTATCACTTGTCAATAAATAGAAAGTGTAGATGTATTAGGCGCAGCTATTTAGTCTATATGGCAAGTAGTAATTATTTACTTCtttttttctac                                                                                     <  1:1408431/100‑5 (MQ=255)
          gCAACGTATCACTTGTCAATAAATAGAAAGTGTAGATGTATTAGGCGCAGCTATTTAGTCTATATGGCAAGTAGTAATTATTTACTTCTTTTTTTCTACg                                                                                   >  1:348703/1‑100 (MQ=255)
          gCAACGTATCACTTGTCAATAAATAGAAAGTGTAGATGTATTAGGCGCAGCTATTTAGTCTATATGGCAAGTAGTAATTATTTACTTCTTTTTTTCTACg                                                                                   <  1:2579652/100‑1 (MQ=255)
           cAACGTATCACTTGTCAATAAATAGAAAGTGTAGATGTATTAGGCGCAGCTATTTAGTCTATATGGCAAGTAGTAATTATTTACTTCTTTTTTTCTACGa                                                                                  >  2:2194851/1‑100 (MQ=255)
              cttATCACTTGTCAATAAATAGAAAGTGTAGATGTATTAGGCGCAGCTATTTAGTCTATATGGCAAGTAGTAATTATTTACTTCTTTTTTTCTACGAACa                                                                               <  1:1368607/98‑1 (MQ=255)
                tATCACTTGTCAATAAATAGAAAGTGTAGATGTATTAGGCGCAGCTATTTAGTCTATATGGCAAGTAGTAATTATTTACTTCTTTTTTTCTACGAACaag                                                                             >  2:1368612/1‑98 (MQ=255)
                 aTCACTTGTCAATAAATAGAAAGTGTAGATGTATTAGGCGCAGCTATTTAGTCTATATGGCAAGTAGTAATTATTTACTTCTTTTTTTCTACGAACAgcg                                                                            <  2:2035705/100‑1 (MQ=255)
                 aTCACTTGTCAATAAATAGAAAGTGTAGATGTATTAGGCGCAGCTATTTAGTCTATATGGCAAGTAGTAATTATTTACTTCTTTTTTTCTACGAACAgcg                                                                            <  2:2222301/100‑1 (MQ=255)
                   cACTTGTCAATAAATAGAAAGTGTAGATGTATTAGGCGCAGCTATTTAGTCTATATGGCAAGTAGTAATTATTTACTTCTTTTTTTCTACGAACAGCGCt                                                                          >  2:505666/1‑100 (MQ=255)
                   cACTTGTCAATAAATAGAAAGTGTAGATGTATTAGGCGCAGCTATTTAGTCTATATGGCAAGTAGTAATTATTTACTTCTTTTTTTCTACGAACAGCGCt                                                                          >  1:1990500/1‑100 (MQ=255)
                    aCTTGTCAATAAATAGAAAGTGTAGATGTATTAGGCGCAGCTATTTAGTCTATATGGCAAGTAGTAATTATTTACTTCTTTTTTTCTACGAACAGCGCTa                                                                         >  1:3061172/1‑100 (MQ=255)
                       tGTCAATAAATAGAAAGTGTAGATGTATTAGGCGCAGCTATTTAGTCTATATGGCAAGTAGTAATTATTTACTTCTTTTTTTCTACGAACAGCGCTATCg                                                                      >  2:2733939/1‑100 (MQ=255)
                       tGTCAATAAATAGAAAGTGTAGATGTATTAGGCGCAGCTATTTAGTCTATATGGCAAGTAGTAATTATTTACTTCTTTTTTTCTACGAACAGCGCTATCg                                                                      >  1:875769/1‑100 (MQ=255)
                       tGTCAATAAATAGAAAGTGTAGATGTATTAGGCGCAGCTATTTAGTCTATATGGCAAGTAGTAATTATTTACTTCTTTTTTTCTACGAACAGCGCTATCg                                                                      >  2:764757/1‑100 (MQ=255)
                         tCAATAAATAGAAAGTGTAGATGTATTAGGCGCAGCTATTTAGTCTATATGGCAAGTAGTAATTATTTACTTCTTTTTTTCTACGAACAGCGCTATCGaa                                                                    <  1:1648951/100‑1 (MQ=255)
                          cAATAAATAGAAAGTGTAGATGTATTAGGCGCAGCTATTTAGTCTATATGGCAAGTAGTAATTATTTACTTCTTTTTTTCTACGAACAGCGCTATCGaaa                                                                   >  1:1697870/1‑100 (MQ=255)
                          cAATAAATAGAAAGTGTAGATGTATTAGGCGCAGCTATTTAGTCTATATGGCAAGTAGTAATTATTTACTTCTTTTTTTCTACGAACAGCGCTATCGaaa                                                                   >  2:21089/1‑100 (MQ=255)
                           aataaataGAAAGTGTAGATGTATTAGGCGCAGCTATTTAGTCTATATGGCAAGTAGTAATTATTTACTTCTTTTTTTTCTACGAACAGCGCTATCGaaa                                                                   <  2:614389/100‑1 (MQ=255)
                           aataaataGAAAGTGTAGATGTATTAGGCGCAGCTATTTAGTCTATATGGCAAGTAGTAATTATTTACTTCTTTTTTTCTACGAACAGCGCTATCGAAAg                                                                  >  1:2217482/1‑100 (MQ=255)
                           aataaataGAAAGTGTAGATGTATTAGGCGCAGCTATTTAGTCTATATGGCAAGTAGTAATTATTTACTTCTTTTTTTCTACGAACAGCGCTATCAAAAg                                                                  >  1:1480553/1‑100 (MQ=255)
                            ataaataGAAAGTGTAGATGTATTAGGCGCAGCTATTTAGTCTATATGGCAAGTAGTAATTATTTACTTCTTTTTTTCTACGAACAGCGCTATCGAAAGc                                                                 <  2:2936524/100‑1 (MQ=255)
                                ataGAAAGTGTAGATGTATTAGGCGCAGCTATTTAGTCTATATGGCAAGTAGTAATTATTTACTTCTTTTTTTCTACGAACAGCGCTATCGAAAGCTTGt                                                             >  2:3033516/1‑100 (MQ=255)
                                  ccgatctgtAGATGTATTAGGCGCAGCTATTTAGTCTATATGGCAAGTAGTAATTATTTACTTCTTTTTTTCTACGAACAGCGCTATCGAAAGCTTGttt                                                           <  2:2057657/94‑1 (MQ=255)
                                  aGAAAGTGTAGATGTATTAGGCGCAGCTATTTAGTCTATATGGCAAGTAGTAATTATTTACTTCTTTTTTTCTACGAACAGCGCTATCGAAAGCTTGttt                                                           <  1:1625487/100‑1 (MQ=255)
                                    aaaGTGTAGATGTATTAGGCGCAGCTATTTAGTCTATATGGCAAGTAGTAATTATTTACTTCTTTTTTTTTTTTTATTGCGCTATCGAAAGCTTAtttt                                                          >  1:2552850/1‑99 (MQ=255)
                                     aaGTGTAGATGTATTAGGCGCAGCTATTTAGTCTATATGGCAAGTAGTAATTATTTGCTTCTTTTTTTTTTTTTTTAACGCTATCGAAAGCTTTTTTTa                                                         >  1:1960116/1‑99 (MQ=255)
                                      aGTGTAGATGTATTAGGCGCAGCTATTTAGTCTATATGGCAAGTAGTAATTATTTACTTCTTTTTTTCTACGAACAGCGCTATCGAAAGCTTGTTTTACa                                                       >  2:1109015/1‑100 (MQ=255)
                                       gtgtAGATGTATTAGGCGCAGCTATTTAGTCTATATGGCAAGTAGTAATTATTTACTTCTTTTTTTCTACGAACAGCGCTATCGAAAGCTTGTTTTACAt                                                      >  2:2952213/1‑100 (MQ=255)
                                       gtgtAGATGTATTAGGCGCAGCTATTTAGTCTATATGGCAAGTAGTAATTATTTACTTCTTTTTTTCTACGAACAGCGCTATCGAAAGCTTGTTTTACAt                                                      >  2:1839733/1‑100 (MQ=255)
                                         gtAGATGTATTAGGCGCAGCTATTTAGTCTATATGGCAAGTAGTAATTATTTACTTCTTTTTTTCTACGAACAGCGCTATCGAAAGCTTGtttagatcgg                                                    >  1:2057649/1‑93 (MQ=255)
                                         gtAGATGTATTAGGCGCAGCTATTTAGTCTATATGGCAAGTAGTAATTATTTACTTCTTTTTTTCTACGAACAGCGCTATCGAAAGCTTGTTTTACATCt                                                    <  1:446288/100‑1 (MQ=255)
                                          tAGATGTATTAGGCGCAGCTATTTAGTCTATATGGCAAGTAGTAATTATTTACTTCTTTTTTTCTACGAACAGCGCTATCGAAAGCTTGTTTTACATCtt                                                   >  1:1735448/1‑100 (MQ=255)
                                           aGATGTATTAGGCGCAGCTATTTAGTCTATATGGCAAGTAGTAATTATTTACTTCTTTTTTTCTACGAACAGCGCTATCGAAAGCTTGTTTTACATCTTc                                                  <  1:1262821/100‑1 (MQ=255)
                                            gATGTATTAGGCGCAGCTATTTAGTCTATATGGCAAGTAGTAATTATTTACTTCTTTTTTTCTACGAACAGCGCTATCGAAAGCTTGTTTTACATCTTCa                                                 <  2:395980/100‑1 (MQ=255)
                                            gATGTATTAGGCGCAGCTATTTAGTCTATATGGCAAGTAGTAATTATTTACTTCTTTTTTTCTACGAACAGCGCTATCGAAAGCTTGTTTTACATCTTCa                                                 >  1:1227233/1‑100 (MQ=255)
                                                 aTTAGGCGCAGCTATTTAGTCTATATGGCAAGTAGTAATTATTTACTTCTTTTTTTCTACGAACAGCGCTATCGAAAGCTTGTTTTACATCTTCAttttt                                            >  1:327895/1‑100 (MQ=255)
                                                 aTTAGGCGCAGCTATTTAGTCTATATGGCAAGTAGTAATTATTTACTTCTTTTTTTCTACGAACAGCGCTATCGAAAGCTTGTTTTACATCTTCAttttt                                            <  2:1553266/100‑1 (MQ=255)
                                                 aTTAGGCGCAGCTATTTAGTCTATATGGCAAGTAGTAATTATTTACTTCTTTTTTTCTACGAACAGCGCTATCGAAAGCTTGTTTTACATCTTCAttttt                                            <  2:1780508/100‑1 (MQ=255)
                                                 aTTAGGCGCAGCTATTTAGTCTATATGGCAAGTAGTAATTATTTACTTCTTTTTTTCTACGAACAGCGCTATCGAAAGCTTGTTTTACATCTTCAttttt                                            <  1:2250347/100‑1 (MQ=255)
                                                  ttAGGCGCAGCTATTTAGTCTATATGGCAAGTAGTAATTATTTACTTCTTTTTTTCTACGAACAGCGCTATCGAAAGCTTGTTTTACATCTTCAtttttt                                           <  1:1097358/100‑1 (MQ=255)
                                                    aGGCGCAGCTATTTAGTCTATATGGCAAGTAGTAATTATTTACTTCTTTTTTTCTACGAACAGCGCTATCGAAAGCTTGTTTTACATCTTCATTTTTTaa                                         >  2:2579383/1‑100 (MQ=255)
                                                    aGGCGCAGCTATTTAGTCTATATGGCAAGTAGTAATTATTTACTTCTTTTTTTCTACGAACAGCGCTATCGAAAGCTTGTTTTACATCTTCATTTTTTaa                                         >  2:2705844/1‑100 (MQ=255)
                                                     ggCGCAGCTATTTAGTCTATATGGCAAGTAGTAATTATTTACTTCTTTTTTTCTACGAACAGCGCTATCGAAAGCTTGTTTTACATCTTCATTTTTTaaa                                        >  2:1175173/1‑100 (MQ=255)
                                                     ggCGCAGCTATTTAGTCTATATGGCAAGTAGTAATTATTTACTTCTTTTTTTCTACGAACAGCGCTATCGAAAGCTTGTTTTACATCTTCATTTTTTaaa                                        <  2:1051910/100‑1 (MQ=255)
                                                      gcgcAGCTATTTAGTCTATATGGCAAGTAGTAATTATTTACTTCTTTTTTTCTACGAACAGCGCTATCGAAAGCTTGTTTTACATCTTCATTTTTTaaa                                        <  1:770092/99‑1 (MQ=255)
                                                      gcgcAGCTATTTAGTCTATATGGCAAGTAGTAATTATTTACTTCTTTTTTTCTACGAACAGCGCTATCGAAAGCTTGTTTTACATCTTCATTTTTTAAAg                                       <  2:948688/100‑1 (MQ=255)
                                                       ctcAGCTATTTAGTCTATATGGCAAGTAGTAATTATTTACTTCTTTTTTTCTACGAACAGCGCTATCGAAAGCTTGTTTTACATCTTCATTTTTTaa                                         <  1:1497327/95‑1 (MQ=255)
                                                       cgcAGCTATTTAGTCTATATGGCAAGTAGTAATTATTTACTTCTTTTTTTTCTACGAACAGCGCTATCGAAAGCTTGTTTTACATCTTCATTTTTTAAAg                                       >  1:2428201/1‑100 (MQ=255)
                                                       cgcAGCTATTTAGTCTATATGGCAAGTAGTAATTATTTACTTCTTTTTTTCTACGAACAGCGCTATCGAAAGCTTGTTTTACATCTTCATTTTTTAAAGa                                      >  1:1110136/1‑100 (MQ=255)
                                                        gcAGCTATTTAGTCTATATGGCAAGTAGTAATTATTTACTTCTTTTTTTCTACGAACAGCGCTATCGAAAGCTTGTTTTACATCTTCATTTTTTAAAGAc                                     <  1:2324517/100‑1 (MQ=255)
                                                          aGCTATTTAGTCTATATGGCAAGTAGTAATTATTTACTTCTTTTTTTCTACGAACAGCGCTATCGAAAGCTTGTTTTACATCTTCATTTTTTAAAGACCa                                   <  1:1109011/100‑1 (MQ=255)
                                                             tatttattCTATATGGCAAGTAGTAATTATTTACTTCTTTTTTTCTACGAACAGCGCTATCGAAAGCTTGTTTTACATCTTCATTTTTTAAAGACCAt                                  <  1:531089/98‑1 (MQ=255)
                                                             tATTTAGTCTATATGGCAAGTAGTAATTATTTACTTCTTTTTTTCTACGAACAGCGCTATCGAAAGCTTGTTTTACATCTTCATTTTTTAAAGACCAtt                                 >  1:1080292/1‑99 (MQ=255)
                                                             tATTTAGTCTATATGGCAAGTAGTAATTATTTACTTCTTTTTTTCTACGAACAGCGCTATCGAAAGCTTGTTTTACATCTTCATTTTTTAAAGACCATTg                                <  1:1175169/100‑1 (MQ=255)
                                                              aTTTAGTCTATATGGCAAGTAGTAATTATTTACTTCTTTTTTTCTACGAACAGCGCTATCGAAAGCTTGTTTTACATCTTCATTTTTTAAAGACCATTGa                               <  2:1886892/100‑1 (MQ=255)
                                                              aTTTAGTCTATATGGCAAGTAGTAATTATTTACTTCTTTTTTTCTACGAACAGCGCTATCGAAAGCTTGTTTTACATCTTCATTTTTTAAAGACCATTGa                               <  2:1966992/100‑1 (MQ=255)
                                                                ttAGTCTATATGGCAAGTAGTAATTATTTACTTCTTTTTTTCTACGAACAGCGCTATCGAAAGCTTGTTTTACATCTTCATTTTTTAAAGACCATTGAtt                             <  2:63832/100‑1 (MQ=255)
                                                                ttAGTCTATATGGCAAGTAGTAATTATTTACTTCTTTTTTTCTACGAACAGCGCTATCGAAAGCTTGTTTTACATCTTCATTTTTTAAAGACCATTGAtt                             <  2:2792445/100‑1 (MQ=255)
                                                                  aGTCTATATGGCAAGTAGTAATTATTTACTCCTTTTTTTCTACGAACAGCGCTATCGAAAGCTTGTTTTACATCTTCATTTTTTAAAGACCATTGATTTa                           >  1:1266102/1‑100 (MQ=255)
                                                                   gTCTATATGGCAAGTAGTAATTATTTACTTCTTTTTTTCTACGAACAGCGCTATCGAAAGCTTGTTTTACATCTTCATTTTTTAAAGACCATTGATTTaa                          >  2:2680866/1‑100 (MQ=255)
                                                                   gTCTATATGGCAAGTAGTAATTATTTACTTCTTTTTTTCTACGAACAGCGCTATCGAAAGCTTGTTTTACATCTTCATTTTTTAAAGACCATTGATTTaa                          <  1:2117028/100‑1 (MQ=255)
                                                                    tCTATATGGCAAGTAGTAATTATTTACTTCTTTTTTTTTATGAACAGCGCTATCGAAAGCTTGTTTTACATCTTCATTTTTTAAAGACCATTGATTTaaa                         >  1:2239442/1‑100 (MQ=255)
                                                                     cTATATGGCAAGTAGTAATTATTTACTTCTTTTTTTTTTTTTTTTGCGCTATCGAAAGCTTGTTTTACATCTTCATTTTTTAAAGACCATTGATTTaa                          >  1:2458332/1‑98 (MQ=255)
                                                                     cTATATGGCAAGTAGTAATTATTTACTTCTTTTTTTCTACGAACAGCGCTATCGAAAGCTTGTTTTACATCTTCATTTTTTAAAGACCATTGATTTaaaa                        >  1:1821354/1‑100 (MQ=255)
                                                                      tataTGGTAAGTAGTAATTATTTACTTCTTTTTTTCTACGAACAGCGCTATCGAAAGCTTGTTTTACATCTTCATTTTTTAAAGACCATTGATTTAAAAt                       >  2:1452883/1‑100 (MQ=255)
                                                                      tataTGGCAAGTAGTAATTATTTACTTCTTTTTTTCTACGAACAGCGCTATCGAAAGCTTGTTTTACATCTTCATTTTTTAAAGACCATTGATTTAAAAt                       <  1:2579369/100‑1 (MQ=255)
                                                                      tataTGGCAAGTAGTAATTATTTACTTCTTTTTTTCTACGAACAGCGCTATCGAAAGCTTGTTTTACATCTTCATTTTTTAAAGACCATTGATTTAAAAt                       <  1:829081/100‑1 (MQ=255)
                                                                       atatGGCAAGTAGTAATTATTTACTTCTTTTTTTCTACGAACAGCGCTATCGAAAGCTTGTTTTACATCTTCATTTTTTAAAGACCATTGATTTAAAATa                      <  2:1227237/100‑1 (MQ=255)
                                                                         atGGCAAGTAGTAATTATTTACTTCTTTTTTTTTTTTTTTTTCGCTATCGAAAGCTTGTTTTACATCTTCATTTTTTAAAGACCATTGATTTAAAATAGt                    >  1:1413666/1‑100 (MQ=255)
                                                                         atGGCAAGTAGTAATTATTTACTTCTTTTTTTCTACGAACAGCGCTATCGAAAGCTTGTTTTACATCTTCATTTTTTAAAGACCATTGATTTAAAATAGt                    >  2:2675191/1‑100 (MQ=255)
                                                                          tGGCAAGTAGTAATTATTTACTTCTTTTTTTCTACGAACAGCGCTATCGAAAGCTTGTTTTACATCTTCATTTTTTAAAGACCATTGATTTAAAATAGTa                   >  1:696227/1‑100 (MQ=255)
                                                                            gCAAGTAGTAATTATTTACTTCTTTTTTTCTACGAACAGCGCTATCGAAAGCTTGTTTTACATCTTCATTTTTTAAAGACCATTGATTTAAAATAGTa                   <  2:2053087/98‑1 (MQ=255)
                                                                             cAAGTAGTAATTATTTACTTCTTTTTTTCTACGAACAGCGCTATCGAAAGCTTGTTTTACATCTTCATTTTTTAAAGACCATTGATTTAAAATAGTATCa                >  1:446746/1‑100 (MQ=255)
                                                                             cAAGTAGTAATTATTTACTTCTTTTTTTCTACGAACAGCGCTATCGAAAGCTTGTTTTACATCTTCATTTTTTAAAGACCATTGATTTAAAATAGTATCa                <  2:3060912/100‑1 (MQ=255)
                                                                             cAAGTAGTAATTATTTACTTCTTTTTTTCTACGAACAGCGCTATCGAAAGCTTGTTTTACATCTTCATTTTTTAAAGACCATTGATTTAAAATAGTATCa                >  1:2779915/1‑100 (MQ=255)
                                                                              aaGTAGTAATTATTTACTTCTTTTTTTCTACGAACAGCGCTATCGAAAGCTTGTTTTACATCTTCATTTTTTAAAGACCATTGATTTAAAATAGTATCaa               <  2:2235001/100‑1 (MQ=255)
                                                                               agtagtAATTATTTACTTCTTTTTTTCTACGAACAGCGCTATCGAAAGCTTGTTTTACATCTTCATTTTTTAAAGACCATTGATTTAAAATAGTATCAAc              <  2:261247/100‑1 (MQ=255)
                                                                                gtagtaATTATTTACTTCTTTTTTTCTACGAACAGCGCTATCGAAAGCTTGTTTTACATCTTCATTTTTTAAAGACCATTGATTTAAAATAGTATCAAcc             <  1:162490/100‑1 (MQ=255)
                                                                                gtagtaATTATTTACTTCTTTTTTTCTACGAACAGCGCTATCGAAAGCTTGTTTTACATCTTCATTTTTTAAAGACCATTGATTTAAAATAGTATCAAcc             <  1:2680850/100‑1 (MQ=255)
                                                                                gtagtaATTATTTACTTCTTTTTTTCTACGAACAGCGCTATCGAAAGCTTGTTTTACATCTTCATTTTTTAAAGACCATTGATTTAAAATAGTATCAAcc             <  2:749589/100‑1 (MQ=255)
                                                                                 tagtaATTATTTACTTCTTTTTTTCTACGAACAGCGCTATCGAAAGCTTGTTTTACATCTTCATTTTTTAAAGACCATTGATTTAAAATAGTATCAACCa            >  1:530464/1‑100 (MQ=255)
                                                                                    taATTATTTACTTCTTTTTTTCTACGAACAGCGCTATCGAAAGCTTGTTTTACATCTTCATTTTTTAAAGACCATTGATTTAAAATAGTATCAACCAAAg         <  1:2666487/100‑1 (MQ=255)
                                                                                    taATTATTTACTTCTTTTTTTCTACGAACAGCGCTATCGAAAGCTTGTTTTACATCTTCATTTTTTAAAGACCATTGATTTAAAATAGTATCAACCAAAg         <  2:1682722/100‑1 (MQ=255)
                                                                                     aaTTATTTACTTCTTTTTTTCTACGAACAGCGCTATCGAAAGCTTGTTTTACATCTTCATTTTTTAAAGACCATTGATTTAAAATAGTATCAACCAAAGc        >  2:1687099/1‑100 (MQ=255)
                                                                                      attattTACTTCTTTTTTTCTACGAACAGCGCTATCGAAAGCTTGTTTTACATCTTCATTTTTTAAAGACCATTGATTTAAAATAGTATCAACCAAAGCa       <  1:2138926/100‑1 (MQ=255)
                                                                                      attattTACTTCTTTTTTTCTACGAACAGCGCTATCGAAAGCTTGTTTTACATCTTCATTTTTTAAAGACCATTGATTTAAAATAGTATCAACCAAAGCa       <  2:1888534/100‑1 (MQ=255)
                                                                                      attattTACTTCTTTTTTTCTACGAACAGCGCTATCGAAAGCTTGTTTTACATCTTCATTTTTTAAAGACCATTGATTTAAAATAGTATCAACCAAAGCa       <  2:490547/100‑1 (MQ=255)
                                                                                      attattTACTTCTTTTTTTCTACGAACAGCGCTATCGAAAGCTTGTTTTACATCTTCATTTTTTAAAGACCATTGATTTAAAATAGTATCAACCAAAGCa       >  1:3013379/1‑100 (MQ=255)
                                                                                           ttACTTCTTTTTTTCTACGAACAGCGCTATCGAAAGCTTGTTTTACATCTTCATTTTTTAAAGACCATTGATTTAAAATAGTATCAACCAAAGCATCTTg  >  1:941097/1‑100 (MQ=255)
                                                                                           ttACTTCTTTTTTTCTACGAACAGCGCTATCGAAAGCTTGTTTTACATCTTCATTTTTTAAAGACCATTGATTTAAAATAGTATCAACCAAAGCATCTTg  <  2:1990507/100‑1 (MQ=255)
                                                                                           ttACTTCTTTTTTTCTACGAACAGCGCTATCGAAAGCTTGTTTTACATCTTCATTTTTTAAAGACCATTGATTTAAAATAGTATCAACCAAAGCATCTTg  <  1:3033497/100‑1 (MQ=255)
                                                                                                  |                                                                                             
TGCTCCTTTTGCAACGTATCACTTGTCAATAAATAGAAAGTGTAGATGTATTAGGCGCAGCTATTTAGTCTATATGGCAAGTAGTAATTATTTACTTCTTTTTTTTCTACGAACAGCGCTATCGAAAGCTTGTTTTACATCTTCATTTTTTAAAGACCATTGATTTAAAATAGTATCAACCAAAGCATCTTG  >  NC_000913/2757557‑2757748

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