Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 217,044 A→G 15.7% intergenic (‑41/+13) yaeF ← / ← proS putative lipoprotein/prolyl‑tRNA synthetase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913217,0440AG15.7% 256.0 / 35.5 113intergenic (‑41/+13)yaeF/proSputative lipoprotein/prolyl‑tRNA synthetase
Reads supporting (aligned to +/- strand):  ref base A (54/41);  new base G (9/9);  total (63/50)
Fisher's exact test for biased strand distribution p-value = 6.14e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.80e-01

Only 100 of 117 total aligned reads displayed.

GCCGAAAGTAAAAGAAAGCTGGGGAGAAACAGGCGGCAGTACGCCTTTGGTTTATCCATTTTATACAATCCATGTAAAAAAAGGGCCCTGAAATTCAGGACCCTTTCTGGCATCAGCCTTTAATCTGTTTCACCAGATATTCGACGATGTCACCAGTCTTAATTAACTGTTTCTCGCCGTTGCGACGATATTTATAT  >  NC_000913/216945‑217141
                                                                                                   |                                                                                                 
gCCGAAAGTAAAAGAAAGCTGGGGAGAAACAGGCGGCAGTACGCCTTTGGTTTATCCATTTTATACAATCCATGTAAAAACAGGGCCCTGAAATTCAGGa                                                                                                   >  2:1895185/1‑100 (MQ=255)
gCCGAAAGTAAAAGAAAGCTGGGGAGAAACAGGCGGCAGTACGCCTTTGGTTTATCCATTTTATACAATCCATGTAAAAAAAGGGCCCTGAAATTCAGGa                                                                                                   >  1:1148560/1‑100 (MQ=255)
gCCGAAAGTAAAAGAAAGCTGGGGAGAAACAGGCGGCAGTACGCCTTTGGTTTATCCATTTTATACAATCCATGTAAAAAAAGGGCCCTGAAATTCAGGa                                                                                                   >  1:595537/1‑100 (MQ=255)
gCCGAAAGTAAAAGAAAGCTGGGGAGAAACAGGCGGCAGTACGCCTTTGGTTTATCCATTTTATACAATCCATGTAAAAAAAGGGCCCTGAAATTCAGGa                                                                                                   >  1:1836116/1‑100 (MQ=255)
 ccGAAAGTAAAAGAAAGCTGGGGAGAAACAGGCGGCAGTACGCCTTTGGTTTATCCATTTTATACAATCCATGTAAAAAAAGGGCCCTGAAATTCAGGAc                                                                                                  <  1:608582/100‑1 (MQ=255)
  cGAAAGTAAAAGAAAGCTGGGGAGAAACAGGCGGCAGTACGCCTTTGGTTTATCCATTTTATACAATCCATGTAAAAAAAGGGTCCTGAATTTCAGGGcc                                                                                                 >  1:1458667/1‑100 (MQ=255)
     aaGTAAAAGAAAGCTGGGGAGAAACAGGCGGCAGTACGCCTTTGGTTTATCCATTTTATACAATCCATGTAAAAAAAGGGTCCTGAATTTCAGGGCCCtt                                                                                              <  1:342799/100‑1 (MQ=39)
        naaaaGAAAGCTGGGGAGAAACAGGCGGCAGTACGCCTTTGGTTTATCCATTTTATACAATCCATGTAAAAAAAGGGCCCTGAAATTCAGGACCCTTTCt                                                                                           >  1:2042504/2‑100 (MQ=255)
         aaaaGAAAGCTGGGGAGAAACAGGCGGCAGTACGCCTTTGGTTTATCCATTTTATACAATCCATGTAAAAAAAGGGCCCTGAAATTCAGGACCCTTTCTg                                                                                          >  2:2330734/1‑100 (MQ=255)
         aaaaGAAAGCTGGGGAGAAACAGGCGGCAGTACGCCTTTGGTTTATCCATTTTATACAATCCATGTAAAAAAAGGGCCCTGAAATTCAGGACCCTTTCTg                                                                                          <  2:136832/100‑1 (MQ=255)
         aaaaGAAAGCTGGGGAGAAACAGGCGGCAGTACCCCTTGCGTTTATCCATTTTATACATTCCATGTAAAAAAAGGGTCCTGAATTTCAGGGCCCtttttt                                                                                          <  1:1411246/100‑4 (MQ=37)
          ctagacagCTGGGGAGTAACAGGCCGCACTACTCCTCCAGTTTATCCATTTTATGCAATCCAGGTCAGAAAAGGGCCCTGAAATTCAGGACCCCTTCt                                                                                           >  1:2353516/3‑98 (MQ=255)
          aaagaaagCTGGGGAGAAACAGGCGGCAGTACGCCTTTGGTTTATCCATTTTATACAATCCATGTAAAAAAAGGGTCCTGAATTTCAGGGCCCttttttt                                                                                         <  2:1845317/100‑5 (MQ=38)
          aaagaaagCTGGGGAGAAACAGGCGGCAGTACGCCTTTGGTTTATCCATTTTATACAATCCATGTAAAAAAAGGGCCCTGAAATTCAGGACCCTTTCTgg                                                                                         <  2:279531/100‑1 (MQ=255)
          aaagaaagCTGGGGAGAAACAGGCGGCAGTACGCCTTTGGTTTATCCATTTTATACAATCCATGTAAAAAAAGGGCCCTGAAATTCAGGACCCTTTCTgg                                                                                         >  2:1848507/1‑100 (MQ=255)
          aaagaaagCTGGGGAGAAACAGGCGGCAGTACGCCTTTGGTTTATCCATTTTATACAATCCAGGTAAATAAAGGGCCCTGAAATTCAGGACCCTTTCTgg                                                                                         >  1:1749517/1‑100 (MQ=255)
           aagaaagCTGGGGAGAAACAGGCGGCAGTACGCCTTTGGTTTATCCATTTTATACAATCCATGTAAAAAAAGGGCCCTGAAATTCAGGACCCTTTCTGGc                                                                                        >  1:1160161/1‑100 (MQ=255)
            agaaagCTGGGGAGAAACAGGCGGCAGTACGCCTTTGGTTTATCCATTTTATACAATCCATGTAAAAAAAGGGTCCTGAATTTCAGGGCCCtttttttac                                                                                       <  1:1462356/100‑7 (MQ=37)
             gaaagCTGGGGAGAAACAGGCGGCAGTACGCCTTTGGTTTATCCATTTTATACAATCCATGTAAAAAAAGGGTCCTGAATTTCAGGGCCCtttttttaca                                                                                      <  2:1458675/100‑8 (MQ=37)
              aaagCTGGGGAGAAACAGGCGGCAGTACGCCTTTGGTTTATCCATTTTATACAATCCATGTAAAAAAAGGGTCCTGAATTTCAGGGCCCtttttttacat                                                                                     <  2:1736358/100‑9 (MQ=37)
              aaagCTGGGGAGAAACAGGCGGCAGTACGCCTTTGGTTTATCCATTTTATACAATCCATGTAAAAAAAGGGTCCTGAATTTCAGGGCCCtttttttacat                                                                                     <  2:311565/100‑9 (MQ=37)
                agCTGGGGAGAAACAGGCGGCAGTACGCCTTTGGTTTATCCATTTTATACAATCCATGTAAAAAAAGGGTCCTGAATTTCAGGGCCCtttttttacatgg                                                                                   <  2:2050431/100‑11 (MQ=25)
                  cTGGGGAGAAACAGGCGGCAGTACGCCTTTGGTTTATCCATTTTATACAATCCATGTAAAAAAAGGGCCCTGAAATTCAGGACCCTTTCTGGCATCAGcc                                                                                 <  2:2256741/100‑1 (MQ=255)
                  cTGGGGAGAAACAGGCGGCAGTACGCCTTTGGTTTATCCATTTTATACAATCCATGTAAAAAAAGGGCCCTGAAATTCAGGACCCTTTCTGGCATCAGcc                                                                                 <  2:860908/100‑1 (MQ=255)
                   tGGGGAGAAACAGGCGGCAGTACGCCTTTGGTTTATCCATTTTATACAATCCATGTAAAAAAAGGGCCCTGAAATTCAGGACCCTTTCTGGCATCAGcc                                                                                 <  1:114829/99‑1 (MQ=255)
                   tGGGGAGAAACAGGCGGCAGTACGCCTTTGGTTTATCCATTTTATACAATCCATGTAAAAAAAGGGCCCTGAAATTCAGGACCCTTTCTGGCATCAGCCt                                                                                >  1:2527676/1‑100 (MQ=255)
                    ggggAGAAACAGGCGGCAGTACGCCTTTGGTTTATCCATTTTATACAATCCATGTAAAAAAAGGGCCCTGAAATTCAGGACCCTTTCTGGCATCAGCCtt                                                                               >  1:2022712/1‑100 (MQ=255)
                      ggAGAAACAGGCGGCAGTACACCTTTGGTTTATCCATTTTATACAATCCATGTAAAAAAAGGGCCCTGAAATTCAGGACCCTTTCTGGCATCAGCCTTTa                                                                             >  2:2437580/1‑100 (MQ=255)
                       gagaAACAGGCGGCAGTACGCCTTTGGTTTATCCATTTTATACAATCCATGTAAAAAAAGGGCCCTGAAATTCAGGACCCTTTCTGGCATCAGCCTTTaa                                                                            <  2:535018/100‑1 (MQ=255)
                       gagaAACAGGCGGCAGTACGCCTTTGGTTTATCCATTTTATACAATCCATGTAAAAAAAGGGCCCTGAAATTCAGGACCCTTTCTGGCATCAGCCTTTaa                                                                            <  2:2342995/100‑1 (MQ=255)
                        agaAACAGGCGGCAGTACGCCTTTGGTTTATCCATTTTATACAATCCATGTAAAAAAAGGGCCCTGAAATTCAGGACCCTTTCTGGCATCAGCCTTTAAt                                                                           <  2:1612203/100‑1 (MQ=255)
                         gataCAGGCGGCAGTACTCCGTTGGTTTATCCATTTTATACAAGCCATTTAAAAAAAGGGCACTGAAATTCATGACCCTTTCTGGCATCAGCCTTTAAt                                                                           <  1:1386076/96‑1 (MQ=255)
                         gaAACAGGCGGCAGTACGCCTTTGGTTTATCCATTTTATACAATCCATGTAAAAAAAGGGCCCTGAAATTCAGGACCCTTTCTGGCATCAGCCTTTAATc                                                                          <  1:2121572/100‑1 (MQ=255)
                           aaCAGGCGGCAGTACGCCTTTGGTTTATCCATTTTATACAATCCATGTAAAAAAAGGGCCCTGAAATTCAGGACCCTTTCTGGCATCAGCCTTTAATCTg                                                                        >  2:2128267/1‑100 (MQ=255)
                             cAGGCGGCAGTACGCCTTTGGTTTATCCATTTTATACAATCCATGTAAAAAAAGGGCCCTGAAATTCAGGACCCTTTCTGGCATCAGCCTTTAATCTGtt                                                                      >  2:1520369/1‑100 (MQ=255)
                             cAGGCGGCAGTACGCCTTTGGTTTATCCATTTTATACAATCCATGTAAAAAAAGGGCCCTGAAATTCAGGACCCTTTCTGGCATCAGCCTTTAATCTGtt                                                                      <  1:1885214/100‑1 (MQ=255)
                                  ggcAGTACGCCTTTGGTTTATCCATTTTATACAATCCATGTAAAAAAAGGGCCCTGAAATTCAGGACCCTTTCTGGCATCAGCCTTTAATCTGTTTCAcc                                                                 >  1:217324/1‑100 (MQ=255)
                                  ggcAGTACGCCTTTGGTTTATCCATTTTATACAATCCATGTAAAAAAAGGGCCCTGAAATTCAGGACCCTTTCTGGCATCAGCCTTTAATCTGTTTCAcc                                                                 <  2:1819393/100‑1 (MQ=255)
                                   gcAGTACGCCTTTGGTTTATCCATTTTATACAATCCATGTAAAAAAAGGGCCCTGAAATTCAGGACCCTTTCTGGCATCAGCCTTTAATCTGTTTCACCa                                                                >  2:1329212/1‑100 (MQ=255)
                                   gcAGTACGCCATTGGTTTATCCATTTTATACAATCCATGTAAAAAAAGGGCCCTGAAATTCAGGACCCTTTCTGGCATCAGCCTTTAATCTGTTTCAcc                                                                 <  1:2039649/99‑1 (MQ=255)
                                      gTACGCCTTTGGTTTATCCATTTTATACAATCCATGTAAAAAAAGGGCCCTGAAATTCAGGACCCTTTCTGGCATCAGCCTTTAATCTGTTTCACCAGat                                                             <  1:878141/100‑1 (MQ=255)
                                      gTACGCCTTTGGTTTATCCATTTTATACAATCCATGTAAAAAAAGGGCCCTGAAATTCAGGACCCTTTCTGACATCAGCCTTTAATCTGTTTCACCAGat                                                             >  1:1214087/1‑100 (MQ=255)
                                         cGCCTTTGGTTTATCCATTTTATACAATCCATGTAAAAAAAGGGCCCTGAAATTCAGGACCCTTTCTGGCATCAGCCTTTAATCTGTTTCACCAGATAtt                                                          >  2:1559273/1‑100 (MQ=255)
                                         cGCCTTTGGTTTATCCATTTTATACAATCCATGTAAAAAAAGGGCCCTGAAATTCAGGACCCTTTCTGGCATCAGCCTTTAATCTGTTTCACCAGATAtt                                                          >  1:2178497/1‑100 (MQ=255)
                                           ccTTTGGTTTATCCATTTTATACAATCCATGTAAAAAAAGGGCCCTGAAATTCAGGACCCTTTCTTGCATCAGCCTTTAATCTGTTTCACCAGATATTcg                                                        <  1:609028/100‑1 (MQ=255)
                                             tttGGTTTATCCATTTTATACAATCCATGTAAAAAAAGGGCCCTGAAATTCATGACCCTTTCTGGCATCAGCCTTTAATCTGTTTCACCAGATATTcgac                                                      <  1:855582/100‑1 (MQ=255)
                                             tttGGTTTATCCATTTTATACAATCCATGTAAAAAAAGGGCCCTGAAATTCAGGACCCTTTCTGGCATCAGCCTTTAATCTGTTTCACCAGATATTcgac                                                      <  2:569158/100‑1 (MQ=255)
                                               tGGTTTATCCATTTTATACAATCCATGTAAAAAAAGGGCCCTGAAATTCAGGACCCTTTCTGGCATCAGCCTTTAATCTGTTTCACCAGATATTcgacga                                                    <  2:1576689/100‑1 (MQ=255)
                                                ggTTTATCCATTTTATACAATCCATGTAAAAAAAGGGCCCTGAAATTCAGGACCCTTTCTGGCATCAGCCTTTAATCTGTTTCACCAGATATTCGACGAt                                                   >  1:1166293/1‑100 (MQ=255)
                                                ggTTTATCCATTTTATACAATCCATGTAAAAAAAGGGCCCTGAAATTCAGGACCCTTTCTGGCATCAGCCTTTAATCTGTTTCACCAGATATTCGACGAt                                                   >  2:60760/1‑100 (MQ=255)
                                                gatttATCCATTTAATACAATCCATGTAAAAAAAGGGCCCTGAAATTCAGGACCCTTTCTGGCATCAGCCTTTAATCTGTTTCACCAGATATTCGACGAt                                                   >  2:1745296/3‑100 (MQ=255)
                                                 gTTTATCCATTTTATACAATCCATGTAAAAAAAGGGCCCTGAAATTCAGGACCCTTTCTGGCATCAGCCTTTAATCTGTTTCACCAGATATTCGACGATg                                                  >  2:423534/1‑100 (MQ=255)
                                                    tATCCATTTTATACAATCCATGTAAAAAAAGGGCCCTGAAATTCAGGACCCTTTCTGGCATCAGCCTTTAATCTGTTTCACCAGATATTCGACGATGTCa                                               <  1:84713/100‑1 (MQ=255)
                                                     aTCCATTTTATACAATCCATGTAAAAAAAGGGCCCTGAAATTCAGGACCCTTTCTGGCATCAGCCTTTAATCTGTTTCACCAGATATTCGACGATGTCa                                               <  1:872487/99‑1 (MQ=255)
                                                     aTCCATTTTATACAATCCATGTAAAAAAAGGGCCCTGAAATTCAGGACCCTTTCTGGCATCAGCCTTTAATCTGTTTCACCAGATATTCGACGATGTCAc                                              <  2:787325/100‑1 (MQ=255)
                                                      tCCATTTTATACAATCCATGTAACAAAAGGCCCCTGAAATTCAGGACCCTTTCTGGCATCAGCCTTTAATCTGTTTGACCAGATATTCGACGATGTctc                                              >  1:1330960/1‑97 (MQ=255)
                                                      tCCATTTTATACAATCCATGTAAAAAAAGGGCCCTGAAATTCAGGACCCTTTCTGGCATCAGCCTTTAATCTGTTTCACCAGATATTCGACGATGTCAcc                                             <  1:2562920/100‑1 (MQ=255)
                                                        cATTTTATACAATCCATGTAAAAAAAGGGCCCTGAAATTCAGGACCCTTTCTGGCATCAGCCTTTAATCTGTTTCACCAGATATTCGACGATGTCACCAg                                           >  1:1774177/1‑100 (MQ=255)
                                                         aTTTTATACAATCCATGTAAAAAAAGGGCCCTGAAATTCAGGACCCTTTCTGGCATCAGCCTTTAATCTGTTTCACCAGATATTCGACGATGTCACCAGt                                          >  2:1591074/1‑100 (MQ=255)
                                                          ttttATACAATCCATGTAAAAAAAGGGCCCTGAAATTCAGGACCCTTTCTGGCATCAGCCTTTAATCTGTTTCACCAGATATTCGACGATGTCACCAGTc                                         <  2:1962209/100‑1 (MQ=255)
                                                             tataCAATCCATGTAAAAAAAGGGCCCTGAAATTCAGGACCCTTTCTGGCATCAGCCTTTAATCTGTTTCACCAGATATTCGACGATGTCACCAGTCtta                                      >  2:210543/1‑100 (MQ=255)
                                                                aCAATCCATGTAAAAAAAGGGCCCTGAAATTCAGGACCCTTTCTGGCATCAGCCTTTAATCTGTTTCACCAGATATTCGACGATGTCACCAGTCttaatt                                   <  1:2124024/100‑1 (MQ=255)
                                                                aCAAACCAAGTAAAAAAAAGGCCCTGAGATTCAGGACCCTTTCTGGCATCAGCCTTTAATCAGTTTTTCCAGATATTCGACGTTGTCACCAGTCttact                                    >  1:938486/1‑97 (MQ=255)
                                                                  aaTCCATGTAAAAAAAGGGCCCTGAAATTCAGGACCCTTTCTGGCATCAGCCTTTAATCTGTTTCACCAGATATTCGACGATGTCACCAgatcggaagag                                 >  2:1801393/1‑90 (MQ=255)
                                                                   aTCCATGTAAAAAAAGGGCCCTGAAATTCAGGACCCTTTCTGGCATCAGCCTTTAATCTGTTTCACCAGATATTCGACGATGTCACCAGTCTTAATTAAc                                >  2:1684375/1‑100 (MQ=255)
                                                                    tacATGTAAAAAAAGGGCTCTGAAATTCAGGACCCTTTTTGGCATCAGCCTTTAATCTGTTTCTCCAGATATTCGACGATGTCACTAGTCttaatt                                   <  1:545618/94‑1 (MQ=255)
                                                                     ccATGTAAAAAAAGGGCCCTGAAATTCAGGACCCTTTCTGGCATCAGCCTTTAATCTGTTTCACCAGATATTCGACGATGTCACCAGTCTTAATTAACTg                              <  2:262652/100‑1 (MQ=255)
                                                                         tgccagaaaGGGTCCTGAATTTCAGGGCCCTTTCTGGCATCAGCCTTTAATCTGTTTCACCAGATATTCGACGATGTCACCAGTCTTAATTAACTGttt                           >  1:2267946/7‑99 (MQ=37)
                                                                         tgccagaaaGGGTCCTGAATTTCAGGGCCCTTTCTGGCATCAGCCTTTAATCTGTTTCACCAGATATTCGACGATGTCACCAGTCTTAATTAACTGTTtc                          >  1:223174/7‑100 (MQ=37)
                                                                         gTAAAAAAAGGGCCCTGAAATTCAGGACCCTTTCTGGCATCAGCCTTTAATCTGTTTCACCAGATATTCGACGATGTCACCAGTCTTAATTAACTGTTtc                          >  1:852196/1‑100 (MQ=255)
                                                                          tAAAAAAAGGGCCCTGAAATTCAGGACCCTTTCTGGCATCAGCCTTTAATCTGTTTCACCAGATATTCGACGATGTCACCAGTCTTAATTAACTGTTtct                         >  2:2586383/1‑100 (MQ=255)
                                                                          tAAAAAAAGGGCCCTGAAATTCAGGACCCTTTCTGGCATCAGCCTTTAATCTGTTTCACCAGATATTCGACGATGTCACCAGTCTTAATTAACTGTTtct                         >  1:36626/1‑100 (MQ=255)
                                                                           aaaaaaaGGGCCCTGAAATTCAGGACCCTTTCTGGCATCAGCCTTTAATCTGTTTCACCAGATATTCGACGATGTCACCAGTCTTAATTAACTGTTtctc                        >  2:712344/1‑100 (MQ=255)
                                                                           aaaaaaaGGGCCCTGAAATTCAGGACCCTTTCTGGCATCAGCCTTTAATCTGTTTCACCAGATATTCGACGATGTCACCAGTCTTAATTAACTGTTtctc                        >  2:2358597/1‑100 (MQ=255)
                                                                            cagaaaGGGTCCTGAATTTCAGGGCCCTTTCTGGCATCAGCCTTTAATCTGTTTCACCAGATATTCGACGATGTCACCAGTCTTAATTAACTGTTTCTCg                       >  1:1142697/4‑100 (MQ=38)
                                                                            aaaaaaGGGCCCTGAAATTCAGGACCCTTTCTGGCATCAGCCTTTAATCTGTTTCACCAGATATTCGACGATGTCACCAGTCTTAATTAACTGTTTCTCg                       >  1:2287992/1‑100 (MQ=255)
                                                                             aaaaaGGGCCCTGAAATTCAGGACCCTTTCTGGCATCAGCCTTTAATCTGTTTCACCAGATATTCGACGATGTCACCAGTCTTAATTAACTGTTTCTCGc                      <  2:875253/100‑1 (MQ=255)
                                                                              aaaaGGGCCCTGAAATTCAGGACCCTTTCTGGCATCAGCCTTTAATCTGTTTCACCAGATATTCGACGATGTCACCAGTCTTAATTAACTGTTTCTCGcc                     >  2:1770699/1‑100 (MQ=255)
                                                                               aaaGGGCCCTGAAATTCAGGACCCTTTCTGGCATCAGCCTTTAATCTGTTTCACCAGATATTCGACGATGTCACCAGTCTTAATTAACTGTTTCTCGCCg                    <  1:1892133/100‑1 (MQ=255)
                                                                                 aGGGTCCTGAATATCAGGGCCCTTTCTGGCATCAGCCTTTAATCTGTTTCACCAGATATTCGACGATGTCACCAGTCTTAATTAACTGTTTCTCGCCGtt                  >  2:1697927/1‑100 (MQ=255)
                                                                                  gggTCCTGAATTTCAGGGCCCTTTCTGGCATCAGCCTTTAATCTGTTTCACCAGATATTCGACGATGTCACCAGTCTTAATTAACTGTTTCTCGCCGTTg                 >  2:882635/1‑100 (MQ=255)
                                                                                  gggTCCTGAATTTCAGGGCCCTTTCTGGCATCAGCCTTTAATCTGTTTCACCAGATATTCGACGATGTCACCAGTCTTAATTAACTGTTTCTCGCCGTTg                 >  2:1088483/1‑100 (MQ=255)
                                                                                  gggCCCTGAAATTCAGGACCCTTTCTGGCATCAGCCTTTAATCTGTTTCACCAGATATTCGACGATGTCACCAGTCTTAATTAACTGTTTCTCGCCGTTg                 <  1:1767129/100‑1 (MQ=255)
                                                                                  gggCCCTGAAATTCAGGACCCTTTCTGGCATCAGCCTTTAATCTGTTTCACCAGATATTCGACGATGTCACCAGTCTTAATTAACTGTTTCTCGCCGTTg                 >  1:415348/1‑100 (MQ=255)
                                                                                   ggCCCTGAAATTCAGGACCCTTTCTGGCATCAGCCTTTAATCTGTTTCACCAGATATTCGACGATGTCACCAGTCTTAATTAACTGTTTCTCGCCGTTGc                <  2:2357970/100‑1 (MQ=255)
                                                                                    gCCCTGAAATTCAGGACCCTTTCTGGCATCAGCCTTTAATCTGTTTCACCAGATATTCGACGATGTCACCAGTCTTAATTAACTGTTTCTCGCCGTTGc                >  2:67885/1‑99 (MQ=255)
                                                                                     cccTGAAATTCAGGACCCTTTCTGGCATCAGCCTTTAATCTGTTTCACCAGATATTCGACGATGTCACCAGTCTTAATTAACTGTTTCTCGCCGTTGcga              >  2:2046107/1‑100 (MQ=255)
                                                                                     cccTGAAATTCAGGACCCTTTCTGGCATCAGCCTTTAATCTGTTTCACCAGATATTCGACGATGTCACCAGTCTTAATTAACTGTTTCTCGCCGTTGcga              >  1:1874524/1‑100 (MQ=255)
                                                                                     cccTGAAATTCAGGACCCTTTCTGGCATCAGCCTTTAATCTGTTTCACCAGATATTCGACGATGTCACCAGTCTTAATTAACTGTTTCTCGCCGTTGcga              >  1:2220313/1‑100 (MQ=255)
                                                                                      ccTGAAATTCAGGACCCTTTCTGGCATCAGCCTTTAATCTGTTTCACCAGATATTCGACGATGTCACCAGTCTTAATTAACTGTTTCTCGCCGTTGcgac             <  2:2127829/100‑1 (MQ=255)
                                                                                      ccTGAAATTCAGGACCCTTTCTGGCATCAGCCTTTAATCTGTTTCACCAGATATTCGACGATGTCACCAGTCTTAATTAACTGTTTCTCGCCGTTGcgac             <  2:847248/100‑1 (MQ=255)
                                                                                       cTGAAATTCAGGACCCTTTCTGGCATCAGCCTTTAATCTGTTTCACCAGATATTCGACGATGTCACCAGTCTTAATTAACTGTTTCTCGCCGTTGcgacg            <  2:804731/100‑1 (MQ=255)
                                                                                          aaaTTCAGGACCCTTTCTGGCATCAGCCTTTAATCTGTTTCACCAGATATTCGACGATGTCACCAGTCTTAATTAACTGTTTCTCGCCGTTGcgac             <  1:548952/96‑1 (MQ=255)
                                                                                          aaaTTCAGGACCCTTTCTGGCATCAGCCTTTAATCTGTTTCACCAGATATTCGACGATGTCACCAGTCTTAATTAACTGTTTCTCGCCGTTGCGACGata         <  2:2132338/100‑1 (MQ=255)
                                                                                          aaaTTCACGACCCTTTCTGGCATCAGCCTTTAATCTGTTTCACCAGATATTCGACGATGTCACCAGTCTTAATTAACTGTTTCTCGCCGTTGCGACGata         <  2:37193/100‑1 (MQ=255)
                                                                                            aTTCAGGACCCTTTCTGGCATCAGCCTTTAATCTGTTTCACCAGATATTCGACGATGTCACCAGTCTTAATTAACTGTTTCTCGCCGTTGCGACGATAtt       >  2:1035242/1‑100 (MQ=255)
                                                                                              tCAGGACCCTTTCTGGCATCAGCCTTTAATCTGTTTCACCAGATATTCGACGATGTCACCAGTCTTAATTAACTGTTTCTCGCCGTTGCGACGATATTta     <  1:293796/100‑1 (MQ=255)
                                                                                               cAGGACCCTTTCTGGCATCAGCCTTTAATCTGTTTCACCAGATATTCGACGATGTCACCAGTCTTAATTAACTGTTTCTCGCCGTTGCGACGATATTtat    >  2:1656324/1‑100 (MQ=255)
                                                                                                 ggACCCTTTCTGGCATCAGCCTTTAATCTGTTTCACCAGATATTCGACGATGTCACCAGTCTTAATTAACTGTTTCTCGCCGTTGCGACGATATTTatat  >  2:1076107/1‑100 (MQ=255)
                                                                                                 ggACCCTTTCTGCCATCAGCCTTTAATCTGTTTCACCAGATATTCGACGATGTCACCAGTCTTAATTAACTGTTTCTCGCCGTTGCGACGATATTTatat  <  1:444719/100‑1 (MQ=255)
                                                                                                   |                                                                                                 
GCCGAAAGTAAAAGAAAGCTGGGGAGAAACAGGCGGCAGTACGCCTTTGGTTTATCCATTTTATACAATCCATGTAAAAAAAGGGCCCTGAAATTCAGGACCCTTTCTGGCATCAGCCTTTAATCTGTTTCACCAGATATTCGACGATGTCACCAGTCTTAATTAACTGTTTCTCGCCGTTGCGACGATATTTATAT  >  NC_000913/216945‑217141

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