Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,586,311 T→A 42.8% H59L (CAC→CTC)  ydeP ← putative oxidoreductase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,586,3110TA42.8% 73.5 / 205.0 158H59L (CAC→CTC) ydePputative oxidoreductase
Reads supporting (aligned to +/- strand):  ref base T (47/43);  new base A (29/38);  total (77/81)
Fisher's exact test for biased strand distribution p-value = 3.33e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.76e-01

Only 100 of 161 total aligned reads displayed.

ACCGTATTCTCAGCAAAGAAAGAGGCGTTTACCTGCTTATCCGTGACTTCCCAGGCGATTGCTTTTGCGCCGTTTTCACAAATGTCGAATGACGCACTGTGCTTAGGATCTGGCCATGCACAACCCGGACAGTCAAAGCCCTCTGGCTTATTCATGTCAAACATGGCAATAACATCCTGGCGTATATCCATCTGCTTA  >  NC_000913/1586212‑1586409
                                                                                                   |                                                                                                  
  cGTATTCTCAGCAAAGAAAGAGGCGTTTACCTGCTTATCCGTGACTTCCCAGGCGATTGCTTTTGCGCCGTTTTCACAAATGTCGAATGACGCACTGAGc                                                                                                  <  1:303684/100‑1 (MQ=255)
  cGTATTCTCAGCAAAGAAAGAGGCGTTTACCTGCTTATCCGTGACTTCCCAGGCGATTGCTTTTGCGCCGTTTTCACAAATGTCGAATGACGCACTGAGc                                                                                                  <  1:1913555/100‑1 (MQ=255)
   gTATTCTCAGCAAAGAAAGAGGCGTTTACCTGCTTATCCGTGACTTCCCAGGCGATTGCTTTTGCGCCGTTTTCACAAATGTCGAATGACGCACTGTGCt                                                                                                 >  1:2385680/1‑100 (MQ=255)
   gTATTCTCAGCAAAGAAAGAGGCGTTTACCTGCTTATCCGTGACTTCCCAGGCGATTGCTTTTGCGCCGTTTTCACAAATGTCGAATGACGCACTGAGCt                                                                                                 >  1:285942/1‑100 (MQ=255)
     attaggaGCAAAGAATGATGGTTTTAGCTGGTTACCCGTGTCTGCCCAGGAGATTGCTTGTGCGCAGTTTTCACAAATGTCGCATGACGTACTGAGc                                                                                                  <  1:2718703/91‑1 (MQ=255)
     aTTCTCAGCAAAGAAAGAGGCGTTTACCTGCTTATCAGTGACTTCCCAGGCGATTGCTTTTGCGCCGTTTTCACAAATGTCGAATGACGCACTGTGCTTa                                                                                               <  1:3464436/100‑1 (MQ=255)
      ttCTCAGCAAAGAAAGAGGCGTTTACCTGCTTATCCGTGACTTCCCAGGCGATTGCTTTTGCGCCGTTTTCACAAATGTCGAATGACGCACTGTGCTTAg                                                                                              <  2:538307/100‑1 (MQ=255)
       tctcAGCAAAGAAAGAGGCGTTTACCTGCTTATCCGTGACTTCCCAGGCGATTGCTTTTGCGCCGTTTTCACAAATGTCGAATGACGCACTGTGCTTAgg                                                                                             >  2:1708150/1‑100 (MQ=255)
       tctcAGCAAAGAAAGAGGCGTTTACCTGCTTATCCGTGACTTCCCAGGCGATTGCTTTTGCGCCGTTTTCACAAATGTCGAATGACGCACTGCGCTTAgg                                                                                             >  1:1381257/1‑100 (MQ=255)
        ctcAGCAAAGAAAGAGGCGTCTACCTGCTTATCCGTGACTTCCCAGGCGATTGCTTTTGCGCCGTTTTCACAAATGTCGAATGACGCACTGTGCTTAGGa                                                                                            >  1:2607708/1‑100 (MQ=255)
        ctcAGCAAAGAAAGAGGCGTCTACCTGCTTATCCGTGACTTCCCAGGCGATTGCTTTTGCGCCGTTTTCACAAATGTCGAATGACGCACTGTGCTTAGGa                                                                                            >  1:2607495/1‑100 (MQ=255)
         tcAGCAAAGAAAGAGGCGTTTACCTGCTTATCCGTGACTTCCCAGGCGATTGCTTTTGCGCCGTTTTCACAAATGTCGAATGACGCACTGTGCTTAGGAt                                                                                           <  2:2517579/100‑1 (MQ=255)
          cAGCAAAGAAAGAGGCGTTTACCTGCTTATCCGTGACTTCCCAGGCGATTGCTTTTGCGCCGTTTTCACAAATGTCGAATGACGCACTGTGCTTAGGATc                                                                                          <  2:1084353/100‑1 (MQ=255)
          cAGCAAAGAAAGAGGCGTTTACCTGCTTATCCGTGACTTCCCAGGCGATTGCTTTTGCGCCGTTTTCACAAATGTCGAATGACGCACTGTGCTTAGGATc                                                                                          <  2:3840214/100‑1 (MQ=255)
           aGCAAAGAAAGAGGCGTTTACCTGCTTATCCGTGACTTCCCAGGCGATTGCTTTTGCGCCGTTTTCACAAATGTCGAATGACGCACTGTGCTTAGGATCt                                                                                         >  2:2546989/1‑100 (MQ=255)
           aGCAAAGAAAGAGGCGTTTACCTGCTTATCCGTGACTTCCCAGGCGATTGCTTTTGCGCCGTTTTCACAAATGTCGAATGACGCACTGTGCTTAGGATCt                                                                                         >  2:3910780/1‑100 (MQ=255)
           aGCAAAGAAAGAGGCGTTTACCTGCTTATCCGTGACTTCCCAGGCGATTGCTTTTGCGCCGTTTTCACAAATGTCGAATGACGCACTGAGCtagatcgga                                                                                         >  2:561781/1‑92 (MQ=255)
            gCAAAGAAAGAGGCGTTTACCTGCTTATCCGTGACTTCCCAGGCGATTGCTTTTGCGCCGTTTTCACAAATGTCGAATGACGCACTGAGCTTAGGATCTg                                                                                        >  1:300670/1‑100 (MQ=255)
             cAAAGAAAGAGGCGTTTACCTGCTTATCCGTGACTTCCCAGGCGATTGCTTTTGCGCCGTTTTCACAAATGTCGAATGACGCACTGAGCTTAGGATCTgg                                                                                       <  2:1565995/100‑1 (MQ=255)
             cAAAGAAAGAGGCGTTTACCTGCTTATCCGTGACTTCCCAGGCGATTGCTTTTGCGCCGTTTTCACAAATGTCGAATGACGCACTGAGCTTAGGATCTgg                                                                                       >  2:280462/1‑100 (MQ=255)
              aaagaaagAGGCGTTTACCTGCTTATCCGTGACTTCCCAGGCGATTGCTTTTGCGCCGTTTTCACAAATGTCGAATGACGCACTGTGCTTAGGATCTGGc                                                                                      >  2:2934363/1‑100 (MQ=255)
              aaagaaagAGGCGTTTACCTGCTTATCCGTGACTTCCCAGGCGATTGCTTTTGCGCCGTTTTCACAAATGTCGAATGACGCACTGAGCTTAGGATCTGGc                                                                                      >  1:314808/1‑100 (MQ=255)
               aagaaagaGGCGTTTACCTGCTTATCCGTGACTTCCCAGGCGATTGCTTTTGCGCCGTTTTCACAAATGTCGAATGACGCACTGTGCTTAGGATCTGGcc                                                                                     >  2:3022039/1‑100 (MQ=255)
                agaaagaGGCGTTTACCTGCTTATCCGTGACTTCCCAGGCGATTGCTTTTGCGCCGTTTTCACAAATGTCGAATGACGCACTGTGCTTAGGATCTGGCCa                                                                                    >  1:689436/1‑100 (MQ=255)
                 gaaagaGGCGTTTACCTGCTTATCCGTGACTTCCCAGGCGATTGCTTTTGCGCCGTTTTCACAAATGTCGAATGACGCACTGTGCTTAGGATCTGGCCAt                                                                                   >  1:1430422/1‑100 (MQ=255)
                 gaaagaGGCGTTTACCTGCTTATCCGTGACTTCCCAGGCGATTGCTTTTGCGCCGTTTTCACAAATGTCGAATGACGCACTGAGCTTAGGATCTGGCCAt                                                                                   >  2:1156098/1‑100 (MQ=255)
                  aaagaGGCGTTTACCTGCTTATCCGTGACTTCCCAGGCGATTGCTTTTGCGCCGTTTTCACAAATGTCGAATGACGCACTGAGCTTAGGATCTGGCCATg                                                                                  <  1:673014/100‑1 (MQ=255)
                    agagGCGTTTACCTGCTTATCCGTGACTTCCCAGGCGATTGCTTTTGCGCCGTTTTCACAAATGTCGAATGACGCACTGTGCTTAGGATCTGGCCATGca                                                                                <  2:1494987/100‑1 (MQ=255)
                    agagGCGTTTACCTGCTTATCCGTGACTTCCCAGGCGATTGCTTTTGCGCCGTTTTCACAAATGTCGAATGACGCACTGTGCTTAGGATCTGGCCATGca                                                                                <  2:2406326/100‑1 (MQ=255)
                    agagGCGTTTACCTGCTTATCCGTGACTTCCCAGGCGATTGCTTTTGCGCCGTTTTCACAAATGTCGAATGACGCACTGTGCTTAGGATCTGGCCATGca                                                                                <  2:3981246/100‑1 (MQ=255)
                     gagGCGTTTACCTGCTTATCCGTGACTTCCCAGGCGATTGCTTTTGCGCCGTTTTCACAAATGTCGAATGACGCACTGTGCTTAGGATCTGGCCATGcac                                                                               >  1:3489947/1‑100 (MQ=255)
                     gagGCGTTTACCTGCTTATCCGTGACTTCCCAGGCGATTGCTTTTGCGCCGTTTTCACAAATGTCGAATGACGCACTGTGCTTAGGATCTGGCCATGcac                                                                               >  1:664280/1‑100 (MQ=255)
                     gagGCGTTTACCTGCTTATCCGTGACTTCCCAGGCGATTGCTTTTGCGCCGTTTTCACAAATGTCGAATGACGCACTGAGCTTAGGATCTGGCCATGcac                                                                               <  2:755240/100‑1 (MQ=255)
                      agGCGTTTACCTGCTTATCCGTGACTTCCCAGGCGATTGCTTTTGCGCCGTTTTCACAAATGTCGAATGACGCACTGTGCTTAGGATCTGGCCATGcaca                                                                              >  1:2610206/1‑100 (MQ=255)
                      agGCGTTTACCTGCTTATCCGTGACTTCCCAGGCGATTGCTTTTGCGCCGTTTTCACAAATGTCGAATGACGCACTGAGCTTAGGATCTGGCCATGcaca                                                                              <  2:3832127/100‑1 (MQ=255)
                      agGCGTTTACCTGCTTATCCGTGACTTCCCAGGCGATTGCTTTTGCGCCGTTTTCACAAATGTCGAATGACGCACTGAGCTTAGGATCTGGCCATGcaca                                                                              <  2:2750711/100‑1 (MQ=255)
                       ggCGTTTACCTGCTTATCCGTGACTTCCCAGGCGATTGCTTTTGCGCCGTTTTCACAAATGTAGAATGTCGCACTGTACTTAGGATCTGGCGATGCACaa                                                                             >  1:2886358/1‑100 (MQ=255)
                        gCGTTTACCTGCTTATCCGTGACTTCCCAGGCGATTGCTTTTGCGCCGTTTTCACAAATGTCGAATGACGCACTGTGCTTAGGATCTGGCCATGCACAAc                                                                            <  2:3123306/100‑1 (MQ=255)
                           tttACCTGCTTATCCGTGACTTCCCAGGCGATTGCTTTTGCGCCGTTTTCACAAATGTCGAATGACGCACTGAGCTTAGGATCTGGCCATGCACAACCCg                                                                         >  2:3249365/1‑100 (MQ=255)
                            ttACCTGCTTATCCGTGACTTCCCAGGCGATTGCTTTTGCGCCGTTTTCACAAATGTCGAATGACGCACTGTGCTTAGGATCTGGCCATGCACAACCCgg                                                                        <  2:1526553/100‑1 (MQ=255)
                            ttACCTGCTTATCCGTGACTTCCCAGGCGATTGCTTTTGCGCCGTTTTCACAAATGTCGAATGACGCACTGAGCTTAGGATCTGGCCATGCACAACCCgg                                                                        <  2:2868421/100‑1 (MQ=255)
                             tACCTGCTTATCCGTGACTTCCCAGGCGATTGCTTTTGCGCCGTTTTCACAAATGTCGAATGACGCACTGTGCTTAGGATCTGGCCATGCACAACCCGGa                                                                       >  1:3886286/1‑100 (MQ=255)
                              aCCTGCTTATCCGTGACTTCCCAGGCGATTGCTTTTGCGCCGTTTTCACAAATGTCGAATGACGCACTGTGCTTAGGATCTGGCCATGCACAACCCGGAc                                                                      >  2:2072471/1‑100 (MQ=255)
                               ccTGCTTATCCGTGACTTCCCAGGCGATTGCTTTTGCGCCGTTTTCACAAATGTCGAATGACGCACTGTGCTTAGGATCTGGCCATGCACAACCCGGACa                                                                     >  2:1881570/1‑100 (MQ=255)
                               ccTGCTTATCCGTGACTTCCCAGGCGATTGCTTTTGCGCCGTTTTCACAAATGTCGAATGACGCACTGTGCTTAGGATCTGGCCATGCACAACCCGGACa                                                                     >  2:1883371/1‑100 (MQ=255)
                                cTGCTTATCCGTGACTTCCCAGGCGATTGCTTTTGCGCCGTTTTCACAAATGTCGACTGACGCACTGTGCTTAGGATCTGGCCATGCACAACCCGGACAg                                                                    >  2:3539519/1‑100 (MQ=255)
                                 tGCTTATCCGTGACTTCCCAGGCGATTGCTTTTGCGCCGTTTTCACAAATGTCGAATGACGCACTGTGCTTAGGATCTGGCCATGCACAACCCGGACAGt                                                                   >  1:2819625/1‑100 (MQ=255)
                                  gCTTATCCGTGACTTCCCAGGCGATTGCTTTTGCGCCGTTTTCACAAATGTCGAATGACGCACTGAGCTTAGGATCTGGCCATGCACAACCCGGACAGTc                                                                  <  2:3087578/100‑1 (MQ=255)
                                  gCTTATCAGTGACTTCCCAGGCGATTGCTTTTGCGCCGTTTTCACAAATGTCGAATGACGCACTGAGCTTAGGATCTGGCCATGCACAACCCGGACAGTc                                                                  >  2:2231442/1‑100 (MQ=255)
                                   cTTATCCGTGACTTCCCAGGCGATTGCTTTTGCGCCGTTTTCACAAATGTCGAATGACGCACTGTGCTTAGGATCTGGCCATGCACAACCCGGACAGTCa                                                                 >  1:891255/1‑100 (MQ=255)
                                     tATCCGTGACTTCCCAGGCGATTGCTTTTGCGCCGTTTTCACAAATGTCGAATGACGCACTGTGCTTAGGATCTGGCCATGCACAACCCGGACAGTCaaa                                                               >  2:1360144/1‑100 (MQ=255)
                                     tATCCGTGACTTCCCAGGCGATTGCTTTTGCGCCGTTTTCACAAATGTCGAATGACGCACTGTGCTTAGGATCTGGCCATGCACAACCCGGACAGTCaaa                                                               >  2:2982418/1‑100 (MQ=255)
                                     tATCCGTGACTTCCCAGGCGATTGCTTTTGCGCCGTTTTCACAAATGTCGAATGACGCACTGTGCTTAGGATCTGGCCATGCACAACCCGGACAGTCaaa                                                               >  2:742140/1‑100 (MQ=255)
                                         cGTGACTTCCCAGGCGATTGCTTTTGCGCCGTTTTCACAAATGTCGAATGACGCACTGAGCTTAGGATCTGGCCATGCACAACCCGGACAGTCAAAGccc                                                           >  1:1255758/1‑100 (MQ=255)
                                             aCTTCCCAGGCGATTGCTTTTGCGCCGTTTTCACAAATGTCGAATGACGCACTGTGCTTAGGATCTGGCCATGCACAACCCGGACAGTCAAAGCCCTCTg                                                       <  1:824808/100‑1 (MQ=255)
                                             aCTTCCCAGGCGATTGCTTTTGCGCCGTTTTCACAAATGTCGAATGACGCACTGAGCTTAGGATCTGGCCATGCACAACCCGGACAGTCAAAGCCCTCTg                                                       <  2:1133699/100‑1 (MQ=255)
                                             aCTTCCCAGGCGATTGCTTTTGCGCCGTTTTCACAAATGTCGAATGACGCACTGAGCTTAGGATCTGGCCATGCACAACCCGGACAGTCAAAGCCCTCTg                                                       <  2:1942192/100‑1 (MQ=255)
                                             aCTTCCCAGGCGATTGCTTTTGCGCCGTTTTCACAAATGTCGAATGACGCACTGAGCTTAGGATCTGGCCATGCACAACCCGGACAGTCAAAGCCCTCTg                                                       <  2:4085450/100‑1 (MQ=255)
                                             aCTTCCCAGGCGATTGCTTTTGCGCCGTTTTCACAAATGTCGAATGACGCACTGAGCTTAGGATCTGGCCATGCACAACCCGGACAGTCAAAGCCCTCTg                                                       <  2:33720/100‑1 (MQ=255)
                                              cTTCCCAGGCGATTGCTTTTGCGCCGTTTTCCCAAATGTCGAATGACGCACTGTGCTTAGGATCTGGCCATGCACAACCCGGACAGTCAAAGCCCTCTgg                                                      >  1:1038367/1‑100 (MQ=255)
                                              cTTCCCAGGCGATTGCTTTTGCGCCGTTTTCACAAATGTCGAATGACGCACTGTGCTTAGGATCTGGCCATGCACAACCCGGACAGTCAAAGCCCTCTgg                                                      >  2:267061/1‑100 (MQ=255)
                                               ttCCCAGGCGATTGCTTTTGCGCCGTTTTCACAAATGTCGAATGACGCACTGTGCTTAGGATCTGGCCATGCACAACCCGGACAGTCAAAGCCCTCTGGc                                                     <  2:4026149/100‑1 (MQ=255)
                                               ttCCCAGGCGATTGCTTTTGCGCCGTTTTCACAAATGTCGAATGACGCACTGAGCTTAGGATCTGGCCATGCACAACCCGGACAGTCAAAGCCCTCTGGc                                                     <  1:3247677/100‑1 (MQ=255)
                                                tCCCAGGCGATTGCTTTTGCGCCGTTTTCACAAATGTCGAATGACGCACTGTGCTTAGGATCTGGCCATGCACAACCCGGACAGTCAAAGCCCTCTGGCt                                                    <  1:811149/100‑1 (MQ=255)
                                                 cccAGGCGATTGCTTTTGCGCCGTTTTCACAAATGTCGAATGACGCACTGTGCTTAGGATCTGGCCATGCACAACCCGGACAGTCAAAGCCCTCTGGCtt                                                   >  1:212771/1‑100 (MQ=255)
                                                      gCGATTGCTTTTGCGCCGTTTTCACAAATGTCGAATGACGCACTGTGCTTAGGATCTGGCCATGCACAACCCGGACAGTCAAAGCCCTCTGGCTTATTCa                                              <  1:2315860/100‑1 (MQ=255)
                                                      gCGATTGCTTTTGCGCCGTTTTCACAAATGTCGAATGACGCACTGTGCTTAGGATCTGGCCATGCACAACCCGGACAGTCAAAGCCCTCTGGCTTATTCa                                              <  2:635694/100‑1 (MQ=255)
                                                      gCGATTGCTTTTGCGCCGTTTTCACAAATGTCGAATGACGCACTGTGCTTAGGATCTGGCCATGCACAACCCGGACAGTCAAAGCCCTCTGGCTTATTCa                                              <  2:2739016/100‑1 (MQ=255)
                                                      gCGATTGCTTTTGCGCCGTTTTCACAAATGTCGAATGACGCACTGAGCTTAGGATCTGGCCATGCACAACCCGGACAGTCAAAGCCCTTTGGCTTATTCa                                              <  1:2231436/100‑1 (MQ=255)
                                                      gCGATTGCTTTTGCGCCGTTTCACAAATGTCGAATGACGCACTGTGCTTAGGATCTGGCCATGCACAACCCGGACAGTCAAAGCCCTCTGGCTTATTCAt                                             <  1:181638/100‑1 (MQ=255)
                                                       cGATTGCTTTTGCGCCGTTTTCACAAATGTCGAATGACGCACTGTGCTTAGGATCTGGCCATGCACAACCCGGACAGTCAAAGCCCTCTGGCTTATTCAt                                             >  2:180805/1‑100 (MQ=255)
                                                        gATTGCTTTTGCGCCGTTTTCACAAATGTCGAATGACGCACTGTGCTTAGGATCTGGCCATGCACAACCCGGACAGTCAAAGCCCTCTGGCTTATTCATg                                            <  1:3022027/100‑1 (MQ=255)
                                                           tGCTTTTGCGCCGTTTTCACAAATGTCGAATGACGCACTGAGCTTAGGATCTGGCCATGCACAACCCGGACAGTCAAAGCCCTCTGGCTTATTCATGTCa                                         >  2:3580931/1‑100 (MQ=255)
                                                            gCTTTTGCGCCGTTTTCACAAATGTCGAATGACGCACTGAGCTTAGGATCTGGCCATGCACAACCCGGACAGTCAAAGCCCTCTGGCTTATTCATGTCaa                                        >  1:3748113/1‑100 (MQ=255)
                                                            gCTTTTGCGCCATTTTCACAACTGCCGAATGACGCACTGAGCTTAGGATCTGGCCATGCACAACCCGGACAGTCAAAGCCCTCTGGCTTATTCATGTCaa                                        <  2:2767930/100‑1 (MQ=255)
                                                             cTTTTGCGCCGTTTTCACAAATGTCGAATGACGCACTGAGCTTAGGATCTGGCCATGCACAACCCGGACAGTCAAAGCCCTCTGGCTTATTCATGTCaaa                                       >  2:782399/1‑100 (MQ=255)
                                                                  gcgcCGTTTTCACAAATGTCGAATGACGCACTGTGCTTAGGATATGGCCATGCAAAACCCGGACAGTCAACGCCATCTGGCTTATTCCTGTCAAAc                                      >  1:3896214/1‑96 (MQ=255)
                                                                   cgcCGTTTTCACAAATGTCGAATGACGCACTGTGCTTAGGATCTGGCCATGCACAACCCGGACAGTCAAAGCCCTCTGGCTTATTCATGTCAAACATGGc                                 >  1:3978007/1‑100 (MQ=255)
                                                                   cgcCGTTTTCACAAATGTCGAATGACGCACTGAGCTTAGGATCTGGCCATGCACAACCCGGACAGTCAAAGCCCTCTGGCTTATTCATGTCAAACATGGc                                 <  2:1327985/100‑1 (MQ=255)
                                                                    gcCGTTTTCACAAATGTCGAATGACGCACTGTGCTTAGGATCTGGCCATGCACAACCCGGACAGTCAAAGCCCTCTGGCTTATTCATGTCAAACATGGCa                                >  1:902375/1‑100 (MQ=255)
                                                                    gcCGTTTTCACAAATGTCGAATGACGCACTGAGCTTAGGATCTGGCCATGCACAACCCGGACAGTCAAAGCCCTCTGGCTTATTCATGTCAAACATGGCa                                >  2:1150511/1‑100 (MQ=255)
                                                                    gcCGTTTTCACAAATATCGAATGACGCACTGTGCTTAGGATCTGGCCATGCACAACCCGGACAGTCAAAGCCCTCTGGCTTATTCATGTCAAACATGGCa                                <  2:4202588/100‑1 (MQ=255)
                                                                       gTTTTCACAAATGTCGAATGACGCACTGTGCTTAGGATCTGGCCATGCACAACCCGGACAGTCAAAGCCCTCTGGCTTATTCATGTCAAACATGGCAATa                             <  1:1003573/100‑1 (MQ=255)
                                                                         tttCACAAATGTCGAATGACGCACTGTGCTTAGGATCTGGCCATGCACAACCCGGACAGTCAAAGCCCTCTGGCTTATTCATGTCAAACATGGCAATAAc                           <  1:2988416/100‑1 (MQ=255)
                                                                              caAATGTCGAATGACGCACTGAGCTTAGGATCTGGCCCTGCACAACCCGGACAGTCAAAGCCCTCTTGCTTATTCATGTCAAACATGGCAATAAACTCCt                      >  2:3160387/1‑100 (MQ=255)
                                                                                 aTGTCGAATGACGCACTGAGCTTAGGATCTGGCCATGCACAACCCGGACAGTCAAAGCCCTCTGGCTTATTCATGTCAAACATGGCAATAACATCCTGGc                   >  2:1210960/1‑100 (MQ=255)
                                                                                  tGTCGAATGACGCACTGTGCTTAGGATCTGGCCATGCACAACCCGGACAGTCAAAGCCCTCTGGCTTATTCATGTCAAACATGGCAATAACATCCTGGCg                  <  1:961026/100‑1 (MQ=255)
                                                                                  tGTCGAATGACGCACTGTGCTTAGGATCTGGCCATGCACAACCCGGACAGTCAAAGCCCTCTGGCTTATTCATGTCAAACATGGCAATAACATCCTGGCg                  <  1:508895/100‑1 (MQ=255)
                                                                                  tGTCGAATGACGCACTGAGCTTAGGATCTGGCCATGCACAACCCGGACAGTCAAAGCCCTCTGGCTTATTCATGTCAAACATGGCAATAACATCCTGGCg                  <  2:3400409/100‑1 (MQ=255)
                                                                                   gTCGAATGGCGCACTGAGCTTAGGATCTGGCCATGCACAACCCGGACAGTCAAAGCCCTCTGGCTTATTCATGTCAAACATGGCAATAACATCCTGGCg                  <  1:1169474/99‑1 (MQ=255)
                                                                                    tCGAATGACGCACTGTGCTTAGGATCTGGCCATGCACAACCCGGACAGTCAAAGCCCTCTGGCTTATTCATGTCAAACATGGCAATAACATCCTGGCGta                >  2:3261543/1‑100 (MQ=255)
                                                                                    tCGAATGACGCACTGTGCTTAGGATCTGGCCATGCACAACCCGGACAGTCAAAGCCCTCTGGCTTATTCATGTCAAACATGGCAATAACATCCTGGCGta                <  1:4151189/100‑1 (MQ=255)
                                                                                      gAATGACGCACTGAGCTTAGGATCTGGCCATGCACAACCCGGACAGTCAAAGCCCTCTGGCTTATTCATGTCAAACATGGCAATAACATCCTGGCGtata              <  1:1166725/100‑1 (MQ=255)
                                                                                         tGACGCACTGTGCTTAGGATCTGGCCATGCACAACCCGGACAGTCAAAGCCCTCTGGCTTATTCATGTCAAACATGGCAATAACATCCTGGCGTATATcc           <  1:2982406/100‑1 (MQ=255)
                                                                                         tGACGCACTGAGCTTAGGATCTGGCCATGCACAACCCGGACAGTCAAAGCCCTCTGGCTTATTCATGTCAAACATGGCAATAACATCCTGGCGTATATcc           >  2:4098326/1‑100 (MQ=255)
                                                                                          gACGCACTGAGCTTAGGATCAGGCCATGCACAACCCGGACAGTCAAAGCCCTCTGGCTTATTCATGTCAAACATGGCAATAACATCCTGGCGTATATCCa          >  1:2611841/1‑100 (MQ=255)
                                                                                           aCGCACTGTGCTTAGGATCTGGCCATGCACAACCCGGACAGTCAAAGCCCTCTGGCTTATTCATGTCAAACATGGCAATAACATCCTGGCGTATATCCAt         >  2:240543/1‑100 (MQ=255)
                                                                                               aCTTTGCTTAGGATCTGGCCATGCACAACCCGGACAGTCAAAGCCCTCTGGCTTATTCATGTCAAACATGGCAATAACATCCTGGCGTATATCCATCTGc     <  1:2312288/100‑1 (MQ=255)
                                                                                                  gtgCTTAGGATCTGGCCATGCACAACCCGGACAGTCAAAGCCCTCTGGCTTATTCATGTCAAACATGGCAATAACATCCTGGCGTATATCCATCTGCTTa  <  1:1711923/100‑1 (MQ=255)
                                                                                                  gtgCTTAGGATCTGGCCATGCACAACCCGGACAGTCAAAGCCCTCTGGCTTATTCATGTCAAACATGGCAATAACATCCTGGCGTATATCCATCTGCTTa  >  2:1975786/1‑100 (MQ=255)
                                                                                                   |                                                                                                  
ACCGTATTCTCAGCAAAGAAAGAGGCGTTTACCTGCTTATCCGTGACTTCCCAGGCGATTGCTTTTGCGCCGTTTTCACAAATGTCGAATGACGCACTGTGCTTAGGATCTGGCCATGCACAACCCGGACAGTCAAAGCCCTCTGGCTTATTCATGTCAAACATGGCAATAACATCCTGGCGTATATCCATCTGCTTA  >  NC_000913/1586212‑1586409

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