Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 217,037 A→T 15.9% intergenic (‑34/+20) 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,0370AT15.9% 289.7 / 41.2 113intergenic (‑34/+20)yaeF/proSputative lipoprotein/prolyl‑tRNA synthetase
Reads supporting (aligned to +/- strand):  ref base A (55/40);  new base T (9/9);  total (64/49)
Fisher's exact test for biased strand distribution p-value = 6.08e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.23e-01

Only 100 of 119 total aligned reads displayed.

GGTACAAGCCGAAAGTAAAAGAAAGCTGGGGAGAAACAGGCGGCAGTACGCCTTTGGTTTATCCATTTTATACAATCCATGTAAAAAAAGGGCCCTGAAATTCAGGACCCTTTCTGGCATCAGCCTTTAATCTGTTTCACCAGATATTCGACGATGTCACCAGTCTTAATTAACTGTTTCTCGCCGTTGCGACGATATT  >  NC_000913/216938‑217136
                                                                                                   |                                                                                                   
ggTACAAGCCGAAAGTAAAAGAAAGCTGGGGAGAAACAGGCGGCAGTACGCCTTTGGTTTATCCATTTTATACAATCCATGTAAAAAAAGGGCCCTGaaa                                                                                                     <  2:1964518/100‑1 (MQ=255)
ggTACAAGCCGAAAGTAAAAGAAAGCTGGGGAGAAACAGGCGGCAGTACGCCTTTGGTTTATCCATTTTATACAATCCATGTAAAAAAAGGGCCCTGaaa                                                                                                     <  2:1230661/100‑1 (MQ=255)
  tACAAGCCGAAAGTAAAAGAAAGCTGGGGAGAAACAGGCGGCAGTACGCCTTTGGTTTATCCATTTTATACAATCCATGTAAAAAAAGGGCCCTGAAAtt                                                                                                   >  2:2529945/1‑100 (MQ=255)
    cAAGCCGAAAGTAAAAGAAAGCTGGGGAGAAACAGGCGGCAGTACGCCTTTGGTTTATCCATTTTATACAATCCATGTAAAAAAAGGGCCCTGAAATTCa                                                                                                 <  2:1404959/100‑1 (MQ=255)
     aaGCCGAAAGTAAAAGAAAGCTGGACAGAAACAGGCGGCAGTACGCCTTTGGTTTATCCATTTTATACAATCCATGTAAAAAAAGGGCCCTGAAATTCAg                                                                                                <  1:1911745/100‑1 (MQ=255)
      aGCCGAAAGTAAAAGAAAGCTGGGGAGAAACAGGCGGCAGTACGCCTTTGGTTTATCCATTTTATACAATCCATGTAAAAAAAGGGCCCTGAAATTCAgg                                                                                               <  2:612204/100‑1 (MQ=255)
       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)
           aaaGTAAAAGAAAGCTGGGGAGAAACAGGCGGCAGTACGCCTTTGGTTTATCCATTTTATACAATCCATGTAAAAAAAGGGCCCTGAAATTCAGGACCCt                                                                                          >  2:804585/1‑100 (MQ=255)
            aaGTAAAAGAAAGCTGGGGAGAAACAGGCGGCAGTACGCCTTTGGTTTATCCATTTTATACAATCCATGTAAAAAAAGGGTCCTGAATTTCAGGGCCCtt                                                                                         <  1:342799/100‑1 (MQ=39)
               tAAAAGAAAGCTGGGGAGAAACAGGCGGCAGTACGCCTTTGGTTTATCCATTTTATACAATCCATGTAAAAAAAGGGCCCTGAAATTCAGGACCCTTTCt                                                                                      >  1:1962199/1‑100 (MQ=255)
               naaaaGAAAGCTGGGGAGAAACAGGCGGCAGTACGCCTTTGGTTTATCCATTTTATACAATCCATGTAAAAAAAGGGCCCTGAAATTCAGGACCCTTTCt                                                                                      >  1:2042504/2‑100 (MQ=255)
                aaaaGAAAGCTGGGGAGAAACAGGCGGCAGTACGCCTTTGGTTTATCCATTTTATACAATCCATGTAAAAAAAGGGTCCTGAATTTCAGGGCCCtttttt                                                                                     <  1:947580/100‑4 (MQ=38)
                aaaaGAAAGCTGGGGAGAAACAGGCGGCAGTACGCCTTTGGTTTATCCATTTTATACAATCCATGTAAAAAAAGGGCCCTGAAATTCAGGACCCTTTCTg                                                                                     <  2:136832/100‑1 (MQ=255)
                aaaaGAAAGCTGGGGAGAAACAGGCGGCAGTACGCCTTTGGTTTATCCATTTTATACAATCCATGTAAAAAAAGGGCCCTGAAATTCAGGACCCTTTCTg                                                                                     >  2:2330734/1‑100 (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)
                  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)
                     aaagCTGGGGAGAAACAGGCGGCAGTACGCCTTTGGTTTATCCATTTTATACAATCCATGTAAAAAAAGGGCCCTGAAATTCAGGACCCTTTCTGGCATc                                                                                >  2:1625280/1‑100 (MQ=255)
                       agCTGGGGAGAAACAGGCGGCAGTACGCCTTTGGTTTATCCATTTTATACAATCCATGTAAAAAAAGGGTCCTGAATTTCAGGGCCCtttttttacatgg                                                                              <  2:2050431/100‑11 (MQ=25)
                         cTGGGGAGAAACAGGCGGCAGTACGCCTTTGGTTTATCCATTTTATACAATCCATGTAAAAAAAGGGCCCTGAAATTCAGGACCCTTTCTGGCATCAGcc                                                                            <  2:860908/100‑1 (MQ=255)
                         cTGGGGAGAAACAGGCGGCAGTACGCCTTTGGTTTATCCATTTTATACAATCCATGTAAAAAAAGGGCCCTGAAATTCAGGACCCTTTCTGGCATCAGcc                                                                            <  2:2256741/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)
                             ggAGAAACAGGCGGCAGTACGCCTTTGGTTTATCCATTTTATACAATCCATGTAAAAAAAGGACCCTGAAATTCAGGACCCTTTCTGGCATCAGCCTTTa                                                                        >  1:2481697/1‑100 (MQ=255)
                             ggAGAAACAGGCGGCAGTACACCTTTGGTTTATCCATTTTATACAATCCATGTAAAAAAAGGGCCCTGAAATTCAGGACCCTTTCTGGCATCAGCCTTTa                                                                        >  2:2437580/1‑100 (MQ=255)
                                gataCAGGCGGCAGTACTCCGTTGGTTTATCCATTTTATACAAGCCATTTAAAAAAAGGGCACTGAAATTCATGACCCTTTCTGGCATCAGCCTTTAAt                                                                      <  1:1386076/96‑1 (MQ=255)
                                  aaCAGGCGGCAGTACGCCTTTGGTTTATCCATTTTATACAATCCATGTAAAAAAAGGGCCCTGAAATTCAGGACCCTTTCTGGCATCAGCCTTTAATCTg                                                                   >  2:2128267/1‑100 (MQ=255)
                                  aaCAGGCGGCAGTACGCCTTTGGTTTATCCATTTTATACAATCCATGTAAAAAAAGGGCCCTGAAATTCAGGACCCTTTCTGGCATCAGCCTTTAATCTg                                                                   >  1:1757089/1‑100 (MQ=255)
                                    cAGGCGGCAGTACGCCTTTGGTTTATCCATTTTATACAATCCATGTAAAAAAAGGGCCCTGAAATTCAGGACCCTTTCTGGCATCAGCCTTTAATCTGtt                                                                 >  2:1520369/1‑100 (MQ=255)
                                    cAGGCGGCAGTACGCCTTTGGTTTATCCATTTTATACAATCCATGTAAAAAAAGGGCCCTGAAATTCAGGACCCTTTCTGGCATCAGCCTTTAATCTGtt                                                                 >  2:499818/1‑100 (MQ=255)
                                        cggcAGTACGCCTTTGGTTTATCCATTTTATACAATCCATGTAAAAAAAGGGCCCTGAAATTCAGGACCCTTTCTGGCATCAGCCTTTAATCTGTTTCAc                                                             <  2:2465617/100‑1 (MQ=255)
                                         ggcAGTACGCCTTTGGTTTATCCATTTTATACAATCCATGTAAAAAAAGGGCCCTGAAATTCAGGACCCTTTCTGGCATCAGCCTTTAATCTGTTTCAcc                                                            <  2:1819393/100‑1 (MQ=255)
                                         ggcAGTACGCCTTTGGTTTATCCATTTTATACAATCCATGTAAAAAAAGGGCCCTGAAATTCAGGACCCTTTCTGGCATCAGCCTTTAATCTGTTTCAcc                                                            >  1:217324/1‑100 (MQ=255)
                                          gcAGTACGCCTTTGGTTTATCCATTTTATACAATCCATGTAAAAAAAGGGCCCTGAAATTCAGGACCCTTTCTGGCATCAGCCTTTAATCTGTTTCACCa                                                           >  2:1329212/1‑100 (MQ=255)
                                          gcAGTACGCCATTGGTTTATCCATTTTATACAATCCATGTAAAAAAAGGGCCCTGAAATTCAGGACCCTTTCTGGCATCAGCCTTTAATCTGTTTCAcc                                                            <  1:2039649/99‑1 (MQ=255)
                                            aGTACGCCTTTGGTTTATCCATTTTATACAATCCATGTAAAAAAAGGGCCCTGAAATTCAGGACCCTTTCTCGCATCAGCCTTTAATCTGTTTCACCAGa                                                         >  2:2025202/1‑100 (MQ=255)
                                             gTACGCCTTTGGTTTATCCATTTTATACAATCCATGTAAAAAAAGGGCCCTGAAATTCAGGACCCTTTCTGGCATCAGCCTTTAATCTGTTTCACCAGat                                                        <  1:878141/100‑1 (MQ=255)
                                             gTACGCCTTTGGTTTATCCATTTTATACAATCCATGTAAAAAAAGGGCCCTGAAATTCAGGACCCTTTCTGACATCAGCCTTTAATCTGTTTCACCAGat                                                        >  1:1214087/1‑100 (MQ=255)
                                                cGCCTTTGGTTTATCCATTTTATACAATCCATGTAAAAAAAGGGCCCTGAAATTCAGGACCCTTTCTGGCATCAGCCTTTAATCTGTTTCACCAGATAtt                                                     >  1:2178497/1‑100 (MQ=255)
                                                cGCCTTTGGTTTATCCATTTTATACAATCCATGTAAAAAAAGGGCCCTGAAATTCAGGACCCTTTCTGGCATCAGCCTTTAATCTGTTTCACCAGATAtt                                                     >  2:1559273/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)
                                                             tCCATTTTATACAATCCATGTAACAAAAGGCCCCTGAAATTCAGGACCCTTTCTGGCATCAGCCTTTAATCTGTTTGACCAGATATTCGACGATGTctc                                         >  1:1330960/1‑97 (MQ=255)
                                                             tCCATTTTATACAATCCATGTAAAAAAAGGGCCCTGAAATTCAGGACCCTTTCTGGCATCAGCCTTTAATCTGTTTCACCAGATATTCGACGATGTCAcc                                        <  1:2562920/100‑1 (MQ=255)
                                                               cATTTTATACAATCCATGTAAAAAAAGGGCCCTGAAATTCAGGACCCTTTCTGGCATCAGCCTTTAATCTGTTTCACCAGATATTCGACGATGTCACCAg                                      >  1:719434/1‑100 (MQ=255)
                                                               cATTTTATACAATCCATGTAAAAAAAGGGCCCTGAAATTCAGGACCCTTTCTGGCATCAGCCTTTAATCTGTTTCACCAGATATTCGACGATGTCACCAg                                      >  1:1774177/1‑100 (MQ=255)
                                                                aTTTTATACAATCCATGTAAAAAAAGGGCCCTGAAATTCAGGACCCTTTCTGGCATCAGCCTTTAATCTGTTTCACCAGATATTCGACGATGTCACCAGt                                     >  2:1591074/1‑100 (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)
                                                                                tgccagaaaGGGTCCTGAATTTCAGGGCCCTTTCTGGCATCAGCCTTTAATCTGTTTCACCAGATATTCGACGATGTCACCAGTCTTAATTAACTGttt                      >  1:2267946/7‑99 (MQ=37)
                                                                                 tAAAAAAAGGGCCCTGAAATTCAGGACCCTTTCTGGCATCAGCCTTTAATCTGTTTCACCAGATATTCGACGATGTCACCAGTCTTAATTAACTGTTtct                    >  1:36626/1‑100 (MQ=255)
                                                                                 tAAAAAAAGGGCCCTGAAATTCAGGACCCTTTCTGGCATCAGCCTTTAATCTGTTTCACCAGATATTCGACGATGTCACCAGTCTTAATTAACTGTTtct                    >  2:2586383/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)
                                                                                        aGGGTCCTGAATTTCAGGGCCCTTTCTGGCATCAGCCTTTAATCTGTTTCACCAGATATTCGACGATGTCACCAGTCTTAATTAACTGTTTCTCGCCGtt             >  1:2173799/1‑100 (MQ=255)
                                                                                        aGGGTCCTGAATATCAGGGCCCTTTCTGGCATCAGCCTTTAATCTGTTTCACCAGATATTCGACGATGTCACCAGTCTTAATTAACTGTTTCTCGCCGtt             >  2:1697927/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)
                                                                                           gtccTGAATTTCAGGGCCCTTTCTGGCATCAGCCTTTAATCTGTTTCACCAGATATTCGACGATGTCACCAGTCTTAATTAACTGTTTCTCGCCGTTGcg          >  1:41779/3‑100 (MQ=255)
                                                                                           gCCCTGAAATTCAGGACCCTTTCTGGCATCAGCCTTTAATCTGTTTCACCAGATATTCGACGATGTCACCAGTCTTAATTAACTGTTTCTCGCCGTTGc           >  2:67885/1‑99 (MQ=255)
                                                                                            cccTGAAATTCAGGACCCTTTCTGGCATCAGCCTTTAATCTGTTTCACCAGATATTCGACGATGTCACCAGTCTTAATTAACTGTTTCTCGCCGTTGcga         >  1:2220313/1‑100 (MQ=255)
                                                                                            cccTGAAATTCAGGACCCTTTCTGGCATCAGCCTTTAATCTGTTTCACCAGATATTCGACGATGTCACCAGTCTTAATTAACTGTTTCTCGCCGTTGcga         >  2:2046107/1‑100 (MQ=255)
                                                                                            cccTGAAATTCAGGACCCTTTCTGGCATCAGCCTTTAATCTGTTTCACCAGATATTCGACGATGTCACCAGTCTTAATTAACTGTTTCTCGCCGTTGcga         >  1:1874524/1‑100 (MQ=255)
                                                                                             ccTGAAATTCAGGACCCTTTCTGGCATCAGCCTTTAATCTGTTTCACCAGATATTCGACGATGTCACCAGTCTTAATTAACTGTTTCTCGCCGTTGcgac        <  2:847248/100‑1 (MQ=255)
                                                                                             ccTGAAATTCAGGACCCTTTCTGGCATCAGCCTTTAATCTGTTTCACCAGATATTCGACGATGTCACCAGTCTTAATTAACTGTTTCTCGCCGTTGcgac        <  2:217326/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)
                                                                                                   |                                                                                                   
GGTACAAGCCGAAAGTAAAAGAAAGCTGGGGAGAAACAGGCGGCAGTACGCCTTTGGTTTATCCATTTTATACAATCCATGTAAAAAAAGGGCCCTGAAATTCAGGACCCTTTCTGGCATCAGCCTTTAATCTGTTTCACCAGATATTCGACGATGTCACCAGTCTTAATTAACTGTTTCTCGCCGTTGCGACGATATT  >  NC_000913/216938‑217136

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: