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)
          gCAACGTATCACTTGTCAATAAATAGAAAGTGTAGATGTATTAGGCGCAGCTATTTAGTCTATATGGCAAGTAGTAATTATTTACTTC‑TTTTTTTCTACg                                                                                   >  1:348703/1‑100 (MQ=255)
          gCAACGTATCACTTGTCAATAAATAGAAAGTGTAGATGTATTAGGCGCAGCTATTTAGTCTATATGGCAAGTAGTAATTATTTACTTC‑TTTTTTTCTACg                                                                                   <  1:2579652/100‑1 (MQ=255)
           cAACGTATCACTTGTCAATAAATAGAAAGTGTAGATGTATTAGGCGCAGCTATTTAGTCTATATGGCAAGTAGTAATTATTTACTTC‑TTTTTTTCTACGa                                                                                  >  2:2194851/1‑100 (MQ=255)
              cttATCACTTGTCAATAAATAGAAAGTGTAGATGTATTAGGCGCAGCTATTTAGTCTATATGGCAAGTAGTAATTATTTACTTC‑TTTTTTTCTACGAACa                                                                               <  1:1368607/98‑1 (MQ=255)
                tATCACTTGTCAATAAATAGAAAGTGTAGATGTATTAGGCGCAGCTATTTAGTCTATATGGCAAGTAGTAATTATTTACTTC‑TTTTTTTCTACGAACaag                                                                             >  2:1368612/1‑98 (MQ=255)
                 aTCACTTGTCAATAAATAGAAAGTGTAGATGTATTAGGCGCAGCTATTTAGTCTATATGGCAAGTAGTAATTATTTACTTC‑TTTTTTTCTACGAACAgcg                                                                            <  2:2035705/100‑1 (MQ=255)
                 aTCACTTGTCAATAAATAGAAAGTGTAGATGTATTAGGCGCAGCTATTTAGTCTATATGGCAAGTAGTAATTATTTACTTC‑TTTTTTTCTACGAACAgcg                                                                            <  2:2222301/100‑1 (MQ=255)
                   cACTTGTCAATAAATAGAAAGTGTAGATGTATTAGGCGCAGCTATTTAGTCTATATGGCAAGTAGTAATTATTTACTTC‑TTTTTTTCTACGAACAGCGCt                                                                          >  2:505666/1‑100 (MQ=255)
                   cACTTGTCAATAAATAGAAAGTGTAGATGTATTAGGCGCAGCTATTTAGTCTATATGGCAAGTAGTAATTATTTACTTC‑TTTTTTTCTACGAACAGCGCt                                                                          >  1:1990500/1‑100 (MQ=255)
                    aCTTGTCAATAAATAGAAAGTGTAGATGTATTAGGCGCAGCTATTTAGTCTATATGGCAAGTAGTAATTATTTACTTC‑TTTTTTTCTACGAACAGCGCTa                                                                         >  1:3061172/1‑100 (MQ=255)
                       tGTCAATAAATAGAAAGTGTAGATGTATTAGGCGCAGCTATTTAGTCTATATGGCAAGTAGTAATTATTTACTTC‑TTTTTTTCTACGAACAGCGCTATCg                                                                      >  2:2733939/1‑100 (MQ=255)
                       tGTCAATAAATAGAAAGTGTAGATGTATTAGGCGCAGCTATTTAGTCTATATGGCAAGTAGTAATTATTTACTTC‑TTTTTTTCTACGAACAGCGCTATCg                                                                      >  1:875769/1‑100 (MQ=255)
                       tGTCAATAAATAGAAAGTGTAGATGTATTAGGCGCAGCTATTTAGTCTATATGGCAAGTAGTAATTATTTACTTC‑TTTTTTTCTACGAACAGCGCTATCg                                                                      >  2:764757/1‑100 (MQ=255)
                         tCAATAAATAGAAAGTGTAGATGTATTAGGCGCAGCTATTTAGTCTATATGGCAAGTAGTAATTATTTACTTC‑TTTTTTTCTACGAACAGCGCTATCGaa                                                                    <  1:1648951/100‑1 (MQ=255)
                          cAATAAATAGAAAGTGTAGATGTATTAGGCGCAGCTATTTAGTCTATATGGCAAGTAGTAATTATTTACTTC‑TTTTTTTCTACGAACAGCGCTATCGaaa                                                                   >  1:1697870/1‑100 (MQ=255)
                          cAATAAATAGAAAGTGTAGATGTATTAGGCGCAGCTATTTAGTCTATATGGCAAGTAGTAATTATTTACTTC‑TTTTTTTCTACGAACAGCGCTATCGaaa                                                                   >  2:21089/1‑100 (MQ=255)
                           aataaataGAAAGTGTAGATGTATTAGGCGCAGCTATTTAGTCTATATGGCAAGTAGTAATTATTTACTTCTTTTTTTTCTACGAACAGCGCTATCGaaa                                                                   <  2:614389/100‑1 (MQ=255)
                           aataaataGAAAGTGTAGATGTATTAGGCGCAGCTATTTAGTCTATATGGCAAGTAGTAATTATTTACTTC‑TTTTTTTCTACGAACAGCGCTATCGAAAg                                                                  >  1:2217482/1‑100 (MQ=255)
                           aataaataGAAAGTGTAGATGTATTAGGCGCAGCTATTTAGTCTATATGGCAAGTAGTAATTATTTACTTC‑TTTTTTTCTACGAACAGCGCTATCAAAAg                                                                  >  1:1480553/1‑100 (MQ=255)
                            ataaataGAAAGTGTAGATGTATTAGGCGCAGCTATTTAGTCTATATGGCAAGTAGTAATTATTTACTTC‑TTTTTTTCTACGAACAGCGCTATCGAAAGc                                                                 <  2:2936524/100‑1 (MQ=255)
                                ataGAAAGTGTAGATGTATTAGGCGCAGCTATTTAGTCTATATGGCAAGTAGTAATTATTTACTTC‑TTTTTTTCTACGAACAGCGCTATCGAAAGCTTGt                                                             >  2:3033516/1‑100 (MQ=255)
                                  ccgatctgtAGATGTATTAGGCGCAGCTATTTAGTCTATATGGCAAGTAGTAATTATTTACTTC‑TTTTTTTCTACGAACAGCGCTATCGAAAGCTTGttt                                                           <  2:2057657/94‑1 (MQ=255)
                                  aGAAAGTGTAGATGTATTAGGCGCAGCTATTTAGTCTATATGGCAAGTAGTAATTATTTACTTC‑TTTTTTTCTACGAACAGCGCTATCGAAAGCTTGttt                                                           <  1:1625487/100‑1 (MQ=255)
                                    aaaGTGTAGATGTATTAGGCGCAGCTATTTAGTCTATATGGCAAGTAGTAATTATTTACTTCTTTTTTTT‑TTTTTATTGCGCTATCGAAAGCTTAtttt                                                          >  1:2552850/1‑99 (MQ=255)
                                     aaGTGTAGATGTATTAGGCGCAGCTATTTAGTCTATATGGCAAGTAGTAATTATTTGCTTCTTTTTTTT‑TTTTTTTAACGCTATCGAAAGCTTTTTTTa                                                         >  1:1960116/1‑99 (MQ=255)
                                      aGTGTAGATGTATTAGGCGCAGCTATTTAGTCTATATGGCAAGTAGTAATTATTTACTTC‑TTTTTTTCTACGAACAGCGCTATCGAAAGCTTGTTTTACa                                                       >  2:1109015/1‑100 (MQ=255)
                                       gtgtAGATGTATTAGGCGCAGCTATTTAGTCTATATGGCAAGTAGTAATTATTTACTTC‑TTTTTTTCTACGAACAGCGCTATCGAAAGCTTGTTTTACAt                                                      >  2:2952213/1‑100 (MQ=255)
                                       gtgtAGATGTATTAGGCGCAGCTATTTAGTCTATATGGCAAGTAGTAATTATTTACTTC‑TTTTTTTCTACGAACAGCGCTATCGAAAGCTTGTTTTACAt                                                      >  2:1839733/1‑100 (MQ=255)
                                         gtAGATGTATTAGGCGCAGCTATTTAGTCTATATGGCAAGTAGTAATTATTTACTTC‑TTTTTTTCTACGAACAGCGCTATCGAAAGCTTGtttagatcgg                                                    >  1:2057649/1‑93 (MQ=255)
                                         gtAGATGTATTAGGCGCAGCTATTTAGTCTATATGGCAAGTAGTAATTATTTACTTC‑TTTTTTTCTACGAACAGCGCTATCGAAAGCTTGTTTTACATCt                                                    <  1:446288/100‑1 (MQ=255)
                                          tAGATGTATTAGGCGCAGCTATTTAGTCTATATGGCAAGTAGTAATTATTTACTTC‑TTTTTTTCTACGAACAGCGCTATCGAAAGCTTGTTTTACATCtt                                                   >  1:1735448/1‑100 (MQ=255)
                                           aGATGTATTAGGCGCAGCTATTTAGTCTATATGGCAAGTAGTAATTATTTACTTC‑TTTTTTTCTACGAACAGCGCTATCGAAAGCTTGTTTTACATCTTc                                                  <  1:1262821/100‑1 (MQ=255)
                                            gATGTATTAGGCGCAGCTATTTAGTCTATATGGCAAGTAGTAATTATTTACTTC‑TTTTTTTCTACGAACAGCGCTATCGAAAGCTTGTTTTACATCTTCa                                                 <  2:395980/100‑1 (MQ=255)
                                            gATGTATTAGGCGCAGCTATTTAGTCTATATGGCAAGTAGTAATTATTTACTTC‑TTTTTTTCTACGAACAGCGCTATCGAAAGCTTGTTTTACATCTTCa                                                 >  1:1227233/1‑100 (MQ=255)
                                                 aTTAGGCGCAGCTATTTAGTCTATATGGCAAGTAGTAATTATTTACTTC‑TTTTTTTCTACGAACAGCGCTATCGAAAGCTTGTTTTACATCTTCAttttt                                            >  1:327895/1‑100 (MQ=255)
                                                 aTTAGGCGCAGCTATTTAGTCTATATGGCAAGTAGTAATTATTTACTTC‑TTTTTTTCTACGAACAGCGCTATCGAAAGCTTGTTTTACATCTTCAttttt                                            <  2:1553266/100‑1 (MQ=255)
                                                 aTTAGGCGCAGCTATTTAGTCTATATGGCAAGTAGTAATTATTTACTTC‑TTTTTTTCTACGAACAGCGCTATCGAAAGCTTGTTTTACATCTTCAttttt                                            <  2:1780508/100‑1 (MQ=255)
                                                 aTTAGGCGCAGCTATTTAGTCTATATGGCAAGTAGTAATTATTTACTTC‑TTTTTTTCTACGAACAGCGCTATCGAAAGCTTGTTTTACATCTTCAttttt                                            <  1:2250347/100‑1 (MQ=255)
                                                  ttAGGCGCAGCTATTTAGTCTATATGGCAAGTAGTAATTATTTACTTC‑TTTTTTTCTACGAACAGCGCTATCGAAAGCTTGTTTTACATCTTCAtttttt                                           <  1:1097358/100‑1 (MQ=255)
                                                    aGGCGCAGCTATTTAGTCTATATGGCAAGTAGTAATTATTTACTTC‑TTTTTTTCTACGAACAGCGCTATCGAAAGCTTGTTTTACATCTTCATTTTTTaa                                         >  2:2579383/1‑100 (MQ=255)
                                                    aGGCGCAGCTATTTAGTCTATATGGCAAGTAGTAATTATTTACTTC‑TTTTTTTCTACGAACAGCGCTATCGAAAGCTTGTTTTACATCTTCATTTTTTaa                                         >  2:2705844/1‑100 (MQ=255)
                                                     ggCGCAGCTATTTAGTCTATATGGCAAGTAGTAATTATTTACTTC‑TTTTTTTCTACGAACAGCGCTATCGAAAGCTTGTTTTACATCTTCATTTTTTaaa                                        >  2:1175173/1‑100 (MQ=255)
                                                     ggCGCAGCTATTTAGTCTATATGGCAAGTAGTAATTATTTACTTC‑TTTTTTTCTACGAACAGCGCTATCGAAAGCTTGTTTTACATCTTCATTTTTTaaa                                        <  2:1051910/100‑1 (MQ=255)
                                                      gcgcAGCTATTTAGTCTATATGGCAAGTAGTAATTATTTACTTC‑TTTTTTTCTACGAACAGCGCTATCGAAAGCTTGTTTTACATCTTCATTTTTTaaa                                        <  1:770092/99‑1 (MQ=255)
                                                      gcgcAGCTATTTAGTCTATATGGCAAGTAGTAATTATTTACTTC‑TTTTTTTCTACGAACAGCGCTATCGAAAGCTTGTTTTACATCTTCATTTTTTAAAg                                       <  2:948688/100‑1 (MQ=255)
                                                       ctcAGCTATTTAGTCTATATGGCAAGTAGTAATTATTTACTTC‑TTTTTTTCTACGAACAGCGCTATCGAAAGCTTGTTTTACATCTTCATTTTTTaa                                         <  1:1497327/95‑1 (MQ=255)
                                                       cgcAGCTATTTAGTCTATATGGCAAGTAGTAATTATTTACTTCTTTTTTTTCTACGAACAGCGCTATCGAAAGCTTGTTTTACATCTTCATTTTTTAAAg                                       >  1:2428201/1‑100 (MQ=255)
                                                       cgcAGCTATTTAGTCTATATGGCAAGTAGTAATTATTTACTTC‑TTTTTTTCTACGAACAGCGCTATCGAAAGCTTGTTTTACATCTTCATTTTTTAAAGa                                      >  1:1110136/1‑100 (MQ=255)
                                                        gcAGCTATTTAGTCTATATGGCAAGTAGTAATTATTTACTTC‑TTTTTTTCTACGAACAGCGCTATCGAAAGCTTGTTTTACATCTTCATTTTTTAAAGAc                                     <  1:2324517/100‑1 (MQ=255)
                                                          aGCTATTTAGTCTATATGGCAAGTAGTAATTATTTACTTC‑TTTTTTTCTACGAACAGCGCTATCGAAAGCTTGTTTTACATCTTCATTTTTTAAAGACCa                                   <  1:1109011/100‑1 (MQ=255)
                                                             tatttattCTATATGGCAAGTAGTAATTATTTACTTC‑TTTTTTTCTACGAACAGCGCTATCGAAAGCTTGTTTTACATCTTCATTTTTTAAAGACCAt                                  <  1:531089/98‑1 (MQ=255)
                                                             tATTTAGTCTATATGGCAAGTAGTAATTATTTACTTC‑TTTTTTTCTACGAACAGCGCTATCGAAAGCTTGTTTTACATCTTCATTTTTTAAAGACCAtt                                 >  1:1080292/1‑99 (MQ=255)
                                                             tATTTAGTCTATATGGCAAGTAGTAATTATTTACTTC‑TTTTTTTCTACGAACAGCGCTATCGAAAGCTTGTTTTACATCTTCATTTTTTAAAGACCATTg                                <  1:1175169/100‑1 (MQ=255)
                                                              aTTTAGTCTATATGGCAAGTAGTAATTATTTACTTC‑TTTTTTTCTACGAACAGCGCTATCGAAAGCTTGTTTTACATCTTCATTTTTTAAAGACCATTGa                               <  2:1886892/100‑1 (MQ=255)
                                                              aTTTAGTCTATATGGCAAGTAGTAATTATTTACTTC‑TTTTTTTCTACGAACAGCGCTATCGAAAGCTTGTTTTACATCTTCATTTTTTAAAGACCATTGa                               <  2:1966992/100‑1 (MQ=255)
                                                                ttAGTCTATATGGCAAGTAGTAATTATTTACTTC‑TTTTTTTCTACGAACAGCGCTATCGAAAGCTTGTTTTACATCTTCATTTTTTAAAGACCATTGAtt                             <  2:63832/100‑1 (MQ=255)
                                                                ttAGTCTATATGGCAAGTAGTAATTATTTACTTC‑TTTTTTTCTACGAACAGCGCTATCGAAAGCTTGTTTTACATCTTCATTTTTTAAAGACCATTGAtt                             <  2:2792445/100‑1 (MQ=255)
                                                                  aGTCTATATGGCAAGTAGTAATTATTTACTCC‑TTTTTTTCTACGAACAGCGCTATCGAAAGCTTGTTTTACATCTTCATTTTTTAAAGACCATTGATTTa                           >  1:1266102/1‑100 (MQ=255)
                                                                   gTCTATATGGCAAGTAGTAATTATTTACTTC‑TTTTTTTCTACGAACAGCGCTATCGAAAGCTTGTTTTACATCTTCATTTTTTAAAGACCATTGATTTaa                          >  2:2680866/1‑100 (MQ=255)
                                                                   gTCTATATGGCAAGTAGTAATTATTTACTTC‑TTTTTTTCTACGAACAGCGCTATCGAAAGCTTGTTTTACATCTTCATTTTTTAAAGACCATTGATTTaa                          <  1:2117028/100‑1 (MQ=255)
                                                                    tCTATATGGCAAGTAGTAATTATTTACTTCTTTTTTTT‑TATGAACAGCGCTATCGAAAGCTTGTTTTACATCTTCATTTTTTAAAGACCATTGATTTaaa                         >  1:2239442/1‑100 (MQ=255)
                                                                     cTATATGGCAAGTAGTAATTATTTACTTCTTTTTTTT‑TTTTTTTTGCGCTATCGAAAGCTTGTTTTACATCTTCATTTTTTAAAGACCATTGATTTaa                          >  1:2458332/1‑98 (MQ=255)
                                                                     cTATATGGCAAGTAGTAATTATTTACTTC‑TTTTTTTCTACGAACAGCGCTATCGAAAGCTTGTTTTACATCTTCATTTTTTAAAGACCATTGATTTaaaa                        >  1:1821354/1‑100 (MQ=255)
                                                                      tataTGGTAAGTAGTAATTATTTACTTC‑TTTTTTTCTACGAACAGCGCTATCGAAAGCTTGTTTTACATCTTCATTTTTTAAAGACCATTGATTTAAAAt                       >  2:1452883/1‑100 (MQ=255)
                                                                      tataTGGCAAGTAGTAATTATTTACTTC‑TTTTTTTCTACGAACAGCGCTATCGAAAGCTTGTTTTACATCTTCATTTTTTAAAGACCATTGATTTAAAAt                       <  1:2579369/100‑1 (MQ=255)
                                                                      tataTGGCAAGTAGTAATTATTTACTTC‑TTTTTTTCTACGAACAGCGCTATCGAAAGCTTGTTTTACATCTTCATTTTTTAAAGACCATTGATTTAAAAt                       <  1:829081/100‑1 (MQ=255)
                                                                       atatGGCAAGTAGTAATTATTTACTTC‑TTTTTTTCTACGAACAGCGCTATCGAAAGCTTGTTTTACATCTTCATTTTTTAAAGACCATTGATTTAAAATa                      <  2:1227237/100‑1 (MQ=255)
                                                                         atGGCAAGTAGTAATTATTTACTTCTTTTTTTT‑TTTTTTTTTCGCTATCGAAAGCTTGTTTTACATCTTCATTTTTTAAAGACCATTGATTTAAAATAGt                    >  1:1413666/1‑100 (MQ=255)
                                                                         atGGCAAGTAGTAATTATTTACTTC‑TTTTTTTCTACGAACAGCGCTATCGAAAGCTTGTTTTACATCTTCATTTTTTAAAGACCATTGATTTAAAATAGt                    >  2:2675191/1‑100 (MQ=255)
                                                                          tGGCAAGTAGTAATTATTTACTTC‑TTTTTTTCTACGAACAGCGCTATCGAAAGCTTGTTTTACATCTTCATTTTTTAAAGACCATTGATTTAAAATAGTa                   >  1:696227/1‑100 (MQ=255)
                                                                            gCAAGTAGTAATTATTTACTTC‑TTTTTTTCTACGAACAGCGCTATCGAAAGCTTGTTTTACATCTTCATTTTTTAAAGACCATTGATTTAAAATAGTa                   <  2:2053087/98‑1 (MQ=255)
                                                                             cAAGTAGTAATTATTTACTTC‑TTTTTTTCTACGAACAGCGCTATCGAAAGCTTGTTTTACATCTTCATTTTTTAAAGACCATTGATTTAAAATAGTATCa                >  1:446746/1‑100 (MQ=255)
                                                                             cAAGTAGTAATTATTTACTTC‑TTTTTTTCTACGAACAGCGCTATCGAAAGCTTGTTTTACATCTTCATTTTTTAAAGACCATTGATTTAAAATAGTATCa                <  2:3060912/100‑1 (MQ=255)
                                                                             cAAGTAGTAATTATTTACTTC‑TTTTTTTCTACGAACAGCGCTATCGAAAGCTTGTTTTACATCTTCATTTTTTAAAGACCATTGATTTAAAATAGTATCa                >  1:2779915/1‑100 (MQ=255)
                                                                              aaGTAGTAATTATTTACTTC‑TTTTTTTCTACGAACAGCGCTATCGAAAGCTTGTTTTACATCTTCATTTTTTAAAGACCATTGATTTAAAATAGTATCaa               <  2:2235001/100‑1 (MQ=255)
                                                                               agtagtAATTATTTACTTC‑TTTTTTTCTACGAACAGCGCTATCGAAAGCTTGTTTTACATCTTCATTTTTTAAAGACCATTGATTTAAAATAGTATCAAc              <  2:261247/100‑1 (MQ=255)
                                                                                gtagtaATTATTTACTTC‑TTTTTTTCTACGAACAGCGCTATCGAAAGCTTGTTTTACATCTTCATTTTTTAAAGACCATTGATTTAAAATAGTATCAAcc             <  1:162490/100‑1 (MQ=255)
                                                                                gtagtaATTATTTACTTC‑TTTTTTTCTACGAACAGCGCTATCGAAAGCTTGTTTTACATCTTCATTTTTTAAAGACCATTGATTTAAAATAGTATCAAcc             <  1:2680850/100‑1 (MQ=255)
                                                                                gtagtaATTATTTACTTC‑TTTTTTTCTACGAACAGCGCTATCGAAAGCTTGTTTTACATCTTCATTTTTTAAAGACCATTGATTTAAAATAGTATCAAcc             <  2:749589/100‑1 (MQ=255)
                                                                                 tagtaATTATTTACTTC‑TTTTTTTCTACGAACAGCGCTATCGAAAGCTTGTTTTACATCTTCATTTTTTAAAGACCATTGATTTAAAATAGTATCAACCa            >  1:530464/1‑100 (MQ=255)
                                                                                    taATTATTTACTTC‑TTTTTTTCTACGAACAGCGCTATCGAAAGCTTGTTTTACATCTTCATTTTTTAAAGACCATTGATTTAAAATAGTATCAACCAAAg         <  1:2666487/100‑1 (MQ=255)
                                                                                    taATTATTTACTTC‑TTTTTTTCTACGAACAGCGCTATCGAAAGCTTGTTTTACATCTTCATTTTTTAAAGACCATTGATTTAAAATAGTATCAACCAAAg         <  2:1682722/100‑1 (MQ=255)
                                                                                     aaTTATTTACTTC‑TTTTTTTCTACGAACAGCGCTATCGAAAGCTTGTTTTACATCTTCATTTTTTAAAGACCATTGATTTAAAATAGTATCAACCAAAGc        >  2:1687099/1‑100 (MQ=255)
                                                                                      attattTACTTC‑TTTTTTTCTACGAACAGCGCTATCGAAAGCTTGTTTTACATCTTCATTTTTTAAAGACCATTGATTTAAAATAGTATCAACCAAAGCa       <  1:2138926/100‑1 (MQ=255)
                                                                                      attattTACTTC‑TTTTTTTCTACGAACAGCGCTATCGAAAGCTTGTTTTACATCTTCATTTTTTAAAGACCATTGATTTAAAATAGTATCAACCAAAGCa       <  2:1888534/100‑1 (MQ=255)
                                                                                      attattTACTTC‑TTTTTTTCTACGAACAGCGCTATCGAAAGCTTGTTTTACATCTTCATTTTTTAAAGACCATTGATTTAAAATAGTATCAACCAAAGCa       <  2:490547/100‑1 (MQ=255)
                                                                                      attattTACTTC‑TTTTTTTCTACGAACAGCGCTATCGAAAGCTTGTTTTACATCTTCATTTTTTAAAGACCATTGATTTAAAATAGTATCAACCAAAGCa       >  1:3013379/1‑100 (MQ=255)
                                                                                           ttACTTC‑TTTTTTTCTACGAACAGCGCTATCGAAAGCTTGTTTTACATCTTCATTTTTTAAAGACCATTGATTTAAAATAGTATCAACCAAAGCATCTTg  >  1:941097/1‑100 (MQ=255)
                                                                                           ttACTTC‑TTTTTTTCTACGAACAGCGCTATCGAAAGCTTGTTTTACATCTTCATTTTTTAAAGACCATTGATTTAAAATAGTATCAACCAAAGCATCTTg  <  2:1990507/100‑1 (MQ=255)
                                                                                           ttACTTC‑TTTTTTTCTACGAACAGCGCTATCGAAAGCTTGTTTTACATCTTCATTTTTTAAAGACCATTGATTTAAAATAGTATCAACCAAAGCATCTTg  <  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: ‑