Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,332,822 G→A 19.4% intergenic (‑87/+40) yhbE ← / ← rpmA EamA family inner membrane putative transporter/50S ribosomal subunit protein L27

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,332,8220GA19.4% 258.0 / 51.2 122intergenic (‑87/+40)yhbE/rpmAEamA family inner membrane putative transporter/50S ribosomal subunit protein L27
Reads supporting (aligned to +/- strand):  ref base G (59/39);  new base A (13/11);  total (72/50)
Fisher's exact test for biased strand distribution p-value = 6.47e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.09e-01

Only 100 of 127 total aligned reads displayed.

CCTGCTGCTTCATACCGTATTCCCTGCAAAGCGCATTTGCGGATTAACGTCGCACCATGCGCATTTTTACCGGTTATCGAATGTAAAAAGCCCCGCAACGTGTTGCGGGGCTTTCATCCGTTACCGGGACGCGAAAAACTTATTCAGCTTCGATGCTGATAAATTTACGGTTTTTCGGGCCTTTAACTTCGAATTTCAC  >  NC_000913/3332723‑3332921
                                                                                                   |                                                                                                   
cctgcTGCTTCATACCGTATTCCCTGCAAAGCGCATTTGCGGATTAACGTCGCACCATGCGCATTTTTACCGGTTATCGAATGTAAAAAGCCCCGCAACg                                                                                                     <  1:442513/100‑1 (MQ=255)
cctgcTGCTTCATACCGTATTCCCTGCAAAGCGCATTTGCGGATTAACGTCGCACCATGCGCATTTTTACCGGTTATCGAATGTAAAAAGCCCCGCAACg                                                                                                     >  1:2091828/1‑100 (MQ=255)
    ctgctTCATACCGTATTCCCTGCAAAGCGCATTTGCGGATTAACGTCGCACCCTGCGCATTTTTACCGGTTCTCGAATGTAAAAAGCCCCGCAACGtttt                                                                                                 >  2:932477/1‑100 (MQ=255)
    ctgctTCATACCGTATTCCCTGCAAAGCGCATTTGCGGATTAACGTCGCACCATGCGCATTTTTACCGGTTATCGAATGTAAAAAGCCCCGCAACGTGtt                                                                                                 >  2:1157603/1‑100 (MQ=255)
     tgctTCATACCGTATTCCCTGCAAAGCGCATTTGCGGATTAACGTCGCACCATGCGCATTTTTACCGGTAACGGATG‑‑‑AAAGCCCCGCAACACGTTGCggg                                                                                            <  1:1331263/100‑1 (MQ=21)
      gctTCATTCCGTATGCCCTGCAAAGCGGATTTGCGGATTAACGTCGCACCATGCGCATTTTTTCCGGTTATTGAATGTAAAAAGCCCCGCAACGTGTTGc                                                                                               <  1:2451838/100‑1 (MQ=255)
      gctTCATACCGTATTCCCTGCAAAGCGCATTTGCGGATTAACGTCGCACCATGCGCATTTTTACCGGTTATCGAATGTAAAAAGCCCCGCAACGTGTTGc                                                                                               >  2:1492473/1‑100 (MQ=255)
      gctTCATACCGTATTCCCTGCAAAGCGCATTTGCGGATTAACGTCGCACCATGCGCATTTTTACCGGTTATCGAATGTAAAAAGCCCCGCAACGTGTTGc                                                                                               >  1:150465/1‑100 (MQ=255)
      gctTCATACCGTATTCCCTGCAAAGCGCATTTGCGGATTAACGTCGCACCATGCGCATTTTTACCGGTTATCGAATGTAAAAAGCCCCGCAACGTGTTGc                                                                                               >  1:1843069/1‑100 (MQ=255)
         tCATACCGTATTCCCTGCAAAGCGCATTTGCGGATTAACGTCGCACCATGCGCATTTTTACCGGTTATCGAATGTAAAAAGCCCCGCAACGTGTTGCggg                                                                                            >  1:1066089/1‑100 (MQ=255)
         tCATACCGTATTCCCTGCAAAGCGCATTTGCGGATTAACGTCGCACCATGCGCATTTTTACCGGTTATCGAATGTAAAAAGCCCCGCAACACGTTGCggg                                                                                            <  1:2283281/100‑1 (MQ=255)
          cATACCGTATTCCCTGCAAAGCGCATTTGCGGATTAACGTCGCACCATGCGCATTTTTACCGGTTATCGAATGTAAAAAGCCCCGCAACGTGTTGCgggg                                                                                           <  1:1757086/100‑1 (MQ=255)
          cATACCGTATTCCCTGCAAAGCGCATTTGCGGATTAACGTCGCACCATGCGCATTTTTACCGGTTATCGAATGTAAAAAGCCCCGCAACACGTTGCgggg                                                                                           <  1:172745/100‑1 (MQ=255)
           aTACCGTATTCCCTGCAAAGCGCATTTGCGGATTAACGTCGCACCATGCGCATTTTTACCGGTTATCGAATGTAGAAAGCCCCGCAACACGTTGCGGGGc                                                                                          <  2:1255889/100‑1 (MQ=39)
           aTACCGTATTCCCTGCAAAGCGCATTTGCGGATTAACGTCGCACCATGCGCATTTTTACCGGTTATCGAATGTAAAAAGCCCCGCAACGTGTTGCGGGGc                                                                                          <  1:1012725/100‑1 (MQ=255)
           aTACCGTATTCCCTGCAAAGCGCATTTGCGGATTAACGTCGCACCATGCGCATTTTTACCGGTTATCGAATGTAAAAAGCCCCGCAACACGTTGCGGGGc                                                                                          <  1:909060/100‑1 (MQ=255)
           aTACCGTATTCCCTGCAAAGCGCATTTGCGGATTAACGTCGCACCATGCGCATTTTTACCGGTTATCGAATGTAAAAAGCCCCGCAACACGTTGCGGGGc                                                                                          <  2:2053084/100‑1 (MQ=255)
           aTACCGTATTCCCTGCAAAGCGCATTTGCGGATTAACGTCGCACCATGCGCATTTTTACCGGTTATCGAATGTAAAAAGCCCCGCAACACGTTGCGGGGc                                                                                          <  2:590739/100‑1 (MQ=255)
            tACCGTATTCCCTGCAAAGCGCATTTGCGGATTAACGTCGCACCATGCGCATTTTTACCGGTTATCGAATGTAAAAAGCCCCGCAACGTGTTGCGGGGc                                                                                          <  2:1133107/99‑1 (MQ=255)
            tACCGTATTCCCTGCAAAGCGCATTTGCGGATTAACGTCGCACCATGCGCATTTTTACCGGTTATCGAATGTAAAAAGCCCCGCAACGTGTTGCGGGGCt                                                                                         <  2:2119521/100‑1 (MQ=255)
            tACCGTATTCCCTGCAAAGCGCATTTGCGGATTAACGTCGCACCATGCGCATTTTTACCGGTTATCGAATGTAAAAAGCCCCGCAACGTGTTGCGGGGCt                                                                                         <  2:841402/100‑1 (MQ=255)
             aCCGTATTCCCTGCAAAGCGCATTTGCGGATTAACGTCGCACCATGCGCATTTTTACCGGTTATCGAATGTAAAAAGCCCCGCAACGTGTTGCGGGGCtt                                                                                        <  1:2522320/100‑1 (MQ=255)
             aCCGTATTCCCTGCAAAGCGCATTTGCGGATTAACGTCGCACCATGCGCATTTTTACCGGTTATCGAATGTAAAAAGCCCCGCAACGTGTTGCGGGGCtt                                                                                        <  2:2289727/100‑1 (MQ=255)
             aCCGTATTCCCTGCAAAGCGCATTTGCGGATTAACGTCGCACCATGCGCATTTTTACCGGTTATCGAATGTAAAAAGCCCCGCAACGTGTTGCGGGGCtt                                                                                        >  1:1368827/1‑100 (MQ=255)
                gTATTCCCTGCAAAGCGCATTTGCGGATTAACGTCGCACCATGCGCATTTTTACCGGTTATCGAATGTAAGAAGCCCCGCAACACGTTGCGGGGCttttt                                                                                     <  2:1812234/100‑3 (MQ=38)
                gTATTCCCTGCAAAGCGCATTTGCGGATTAACGTCGCACCATGCGCATTTTTACCGGTTATCGAATGTAAAAAGCCCCGCAACGTGTTGCGGGGCTTTCa                                                                                     >  2:936110/1‑100 (MQ=255)
                gTATTCCCTGCAAAGCGCATTTGCGGATTAACGTCGCACCATGCGCATTTTTACCGGTTATCGAATGTAAAAAGCCCCGCAACACGTTGCGGGGCttttt                                                                                     <  2:724901/100‑3 (MQ=255)
                 tATTCCCTGCAAAGCGCATTTGCGGATTAACGTCGCACCATGCGCATTTTTACCGGTTATCGAATGTAAAAAGCCCCGCAACACGTTGCGGGGCttttta                                                                                    <  2:2576620/100‑4 (MQ=38)
                 naTTCCCTGCAAAGCGCATTTGCGGATTAACGTCGCACCATGCGCATTTTTACCGGTTATCGAATGTAAAAAGCCCCGCAACGTGTTGCGGGGCTTTCAt                                                                                    >  2:1289893/2‑100 (MQ=255)
                  aTTCCCTGCAAAGCGCATTTGCGGATTAACGTCGCACCATGCGCATTTTTACCGGTTATCGAATGTAAAAAGCCCCGCAACGTGTTGCGGGGCTTTCATc                                                                                   >  1:2273577/1‑100 (MQ=255)
                    tCCCTGCAAAGCGCATTTGCGGATTAACGTCGCACCATGCGCATTTTTACCGGTTATCGAATGTAAAAAGCCCCGCAACGTGTTGCGGGGCTTTCATCCg                                                                                 <  1:2325573/100‑1 (MQ=255)
                     cccTGCAAAGCGCATTTGCGGATTAACGTCGCACCATGCGCATTTTTACCGGTTATCGAATGTAAAAAGCCCCGCAACGTGTTGCGGGGCTTTCATCCGt                                                                                >  2:1403962/1‑100 (MQ=255)
                      ccTGCAAAGCGCATTTGCGGATTAACGTCGCACCATGCGCATTTTTACCGGTTATCGAATGTAAAAAGCCCCGCAACGTGTTGCGGGGCTTTCATCCGtt                                                                               >  2:515055/1‑100 (MQ=255)
                      ccTGCAAAGCGCATTTGCGGATTAACGTCGCACCATGCGCATTTTTACCGGTTATCGAATGTAAAAAGCCCCGCAACGTGTTGCGGGGCTTTCATCCGtt                                                                               >  1:239216/1‑100 (MQ=255)
                      ccTGCAAAGCGCATTTGCGGATTAACGTCGCACCATGCGCATTTTTACCGGTTATCGAATGTAAAAAGCCCCGCAACGTGTTGCGGGGCTTTCATCCGtt                                                                               >  1:2369542/1‑100 (MQ=255)
                      ccTGCAAAGCGCATTTGCGGATTAACGTCGCACCATGCGCATTTTTACCGGTTATCGAATGTAAAAAGCCCCGCAACGTGTTGCGGGGCTTTCATCCGtt                                                                               >  1:2017714/1‑100 (MQ=255)
                       cTGCAAAGCGCATTTGCGGATTAACGTCGCACCATGCGCATTTTTACCGGTTATCGAATGTAAAAAGCCCCGCAACGTGTTGCGGGGCTTTCATCCGTTa                                                                              >  1:773612/1‑100 (MQ=255)
                         gCAAAGCGCATTTGCGGATTAACGTCGCACCATGCGCATTTTTACCGGTTATCGAATGTAAAAAGCCCCGCAACGTGTTGCGGGGCTTTCATCCGTTAcc                                                                            >  2:2297032/1‑100 (MQ=255)
                         gCAAAGCGCATTTGCGGATTAACGTCGCACCATGCGCATTTTTACCGGTTATCGAATGTAAAAAGCCCCGCAACGTGTTGCGGGGCTTTCATCCGTTAcc                                                                            >  1:1623538/1‑100 (MQ=255)
                          cAAAGCGCATTTGCGGATTAACGTCGCACCATGCGCATTTTTACCGGTTATCGAATGTAAAAAGCCCCGCAACGTGTTGCGGGGCTTTCATCCGTTACCg                                                                           <  1:936107/100‑1 (MQ=255)
                           aaaGCGCATTTGCGGATTAACGTCGCACCATGCGCATTTTTACCGGTTATCGAATGTAAAAAGCCCCGCAACGTGTTGCGGGGCTTTCATCCGTTACCgg                                                                          <  1:580477/100‑1 (MQ=255)
                               cgcATTTGCGGATTAACGTCGCACCATGCGCATTTTTACCGGTTATCGAATGTAAAAAGCCCCGCAACGTGTTGCGGGGCTTTCATCCGTTACCGGGAcg                                                                      <  2:136705/100‑1 (MQ=255)
                               cgcATTTGCGGATTAACGTCGCACCATGCGCATTTTTACCGGTTATCGAATGTAAAAAGCCCCGCAACGTGTTGCGGGGCTTTCATCCGTTACCGGGAcg                                                                      <  2:398211/100‑1 (MQ=255)
                                gcATTTGCGGATTAACGTCGCACCATGCGCATTTTTTCCGGTTATCGCATGTAAAAAGGCCCGCAACGTGTTGCGGGGCTTTCATCCGTTACCgg                                                                          >  1:916514/1‑95 (MQ=255)
                                gcATTTGCGGATTAACGTCGCACCATGCGCATTTTTACCGGTTATCGAATGTAAAAAGCCCCGCAACGTGTTGCGGGGCTTTCATCCGTTACCGGGAcgc                                                                     <  1:1606716/100‑1 (MQ=255)
                                  aTTTGCGGATTAACGTCGCACCATGCGCATTTTTACCGGTTATCGAATGTAAAAAGCCCCGCAACGTGTTGCGGGGCTTTCATCCGTTACCGGGACGCGa                                                                   >  2:1916565/1‑100 (MQ=255)
                                    ttGCGGATTAACGTCGCACCATGCGCATTTTTACCGGTTATCGAATGTAAAAAGCCCCGCAACGTGTTGCGGGGCTTTCATCCGTTACCGGGACGCGaaa                                                                 >  2:254797/1‑100 (MQ=255)
                                     tGCGGATTAACGTCGCACCATGCGCATTTTTACCGGTTATCGAATGTAAAAAGCCCCGCAACGTGTTGCGGGGCTTTCATCCGTTACCGGGACGCGaaaa                                                                >  1:2129737/1‑100 (MQ=255)
                                     tGCGGATTAACGTCGCACCATGCGCATTTTTACCGGTTATCGAATGTAAAAAGCCCCGCAACGTGTTGCGGGGCTTTCATCCGTTACCGGGACGCGaaaa                                                                <  1:1024456/100‑1 (MQ=255)
                                     tGCGGATTAACGTCGCACCATGCGCATTTTTACCGGTTATCGAATGTAAAAAGCCCCGCAACGTGTTGCGGGGCTTTCATCCGTTACCGGGACGCGaaa                                                                 >  1:2070181/1‑99 (MQ=255)
                                       cGGATTAACGTCGCACCATGCGCATTTTTACCGGTTATCGAATGTAAAAAGCCCCGCAACGTGTTGCGGGGCTTTCATCCGTTACCGGGACGCGAAAAAc                                                              >  2:534083/1‑100 (MQ=255)
                                         gATTAACGTCGCACCATGCGCATTTTTACCGGTTATCGAATGTAAAAAGCCCCGCAACGTGTTGCGGGGCTTTCATCCGTTACCGGGACGCGAAAAACtt                                                            >  2:181936/1‑100 (MQ=255)
                                          aTTAACGTCGCACCATGCGCATTTTTACAGGTTATCGAATGTAAAAAGCCCCGCAACGTGTTGCGGGGCTTTCATCCGTTACCGGGACGCGAAAAAc                                                              <  1:62818/97‑1 (MQ=255)
                                           ttAACGTCGCACCATGCGCATTTTTACCGGTTATCGAATGTAAAAAGCCCCGCAACGTGTTGCGGGGCTTTCATCCGTTACCGGGACGCGAAAAACTTAt                                                          <  2:89640/100‑1 (MQ=255)
                                             aaCGTCGCACCATGCGCATTTTTACCGGTTATCGAATGTAAAAAGCCCCGCAACGTGTTGCGGGGCTTTCATCCGTTACCGGGACGCGAAAAACTTATTc                                                        >  2:2079406/1‑100 (MQ=255)
                                              aCGTCGCACCATGCGCATTTTTACCGGTTATCGAATGTAAAAAGCCCCGCAACGTGTTGCGGGGCTTTCATCCGTTACCGGGACGCGAAAAACTTATTCa                                                       <  2:1468799/100‑1 (MQ=255)
                                              aCGTCGCACCATGCGCATTTTTACCGGTTATCGAATGTAAAAAGCCCCGCAACGTGTTGCGGGGCTTTCATCCGTTACCGGGACGCGAAAAACTTATTCa                                                       <  2:1384473/100‑1 (MQ=255)
                                              aCGTCGCACCATGCGCATTTTTACCGGTTATCGAATGTAAAAAGCCCCGCAACGTGTTGCGGGGCTTTCATCCGTTACCGGGACGCGAAAAACTTATTCa                                                       <  1:1262982/100‑1 (MQ=255)
                                                   gCACCATGCGCATTTTTACTTTTTATCGAATGTAAAAAGCCCCGCAACGTGTTGCGGGGCTTTCATCCGTTACCGGGACGCGAAACACTTATTCAGCTTc                                                  >  1:1457165/1‑100 (MQ=255)
                                                   gCACCATGCGCATTTTTACCGGTTATCGAATGTAAAAAGCCCCGCAACGTGTTGCGGGGCTTTCATCCGTTACCGGGACGCGAAAAACTTATTCAGCTTc                                                  >  2:1490836/1‑100 (MQ=255)
                                                     aCCATGCGCATTTTTACCGGTTATCGAATGTAAAAAGCCCCGCAACGTGTTGCGGGGCTTTCATCCGTTACCGGGACGCGAAAAACTTATTCAGCTTCGa                                                >  1:2468725/1‑100 (MQ=255)
                                                            gcATTTTTACCGGTTATCGAATGTAAAAAGCCCCGCAACGTGTTGCGGGGCTTTCATCCGTTACCGGGACGCGAAAAACTTATTCAGCTTCGATGCTGAt                                         <  2:1454759/100‑1 (MQ=255)
                                                             cATTTTTACCGGTTATCGAATGTAAAAAGCCCCGCAACGTGTTGCGGGGCTTTCATCCGTTACCGGGACGCGAAAAACTTATTCAGCTTCGATGCTGATa                                        >  2:1717119/1‑100 (MQ=255)
                                                                      cGGTTATCGAATGTAAAAAGCCCCGCAACGTGTTGCGGGGCTTTCATCCGTTACCGGGACGCGAAAAACTTATTCAGCTTCGATGCTGATAAATTTACgg                               >  1:453988/1‑100 (MQ=255)
                                                                       ggTTATCGAATGTAAAAAGCCCCGCAACGTGTTGCGGAGCTTTCATCCGTTACCGGGACGCGAAAAACTTATTCAGCTTCGATGCTGATAAATTTACGGt                              >  1:1378330/1‑100 (MQ=255)
                                                                          tATCGAATGTAAAAAGCCCCGCAACGTGTTGCGGGGCTTTCATCCGTTACCGGGACGCGAAAAACTTATTCAGCTTCGATGCTGATAAATTTa                                  >  1:349749/1‑93 (MQ=255)
                                                                            tCGAATGTAAAAAGCCCCGCAACGTGTTGCGGGGCTTTCATCCGTTACCGGGACGCGAAAAACTTATTCAGCTTCGATGCTGATAAATTTACGGTTTTTc                         >  2:2432718/1‑100 (MQ=255)
                                                                             cGAATGTAAAAAGCCCCGCAACGTGTTGCGGGGCTTTCATCCGTTACCGGGACGCGAAAACTTATTCAGCTTCGATGCTGATAAATTTACGGTTTTTCgg                       >  1:1963802/1‑100 (MQ=255)
                                                                              gAATGTAAAAAGCCCCGCAACGTGTTGCGGGGCTTTCATCCGTTACCGGGACGCGAAAAACTTATTCAGCTTCGATGCTGATAAATTTACGGTTTTTCgg                       >  2:642056/1‑100 (MQ=255)
                                                                              gAATGTAAAAAGCCCCGCAACGTGTTGCGGGGCTTTCATCCGTTACCGGGACGCGAAAAACTTATTCAGCTTCGATGCTGATAAATTTACGGTTTTTCgg                       <  2:643739/100‑1 (MQ=255)
                                                                              gAATGTAAAAAGCCCCGCAACGTGTTGCGGGGCTTTCATCCGTTACCGGGACGCGAAAAACTTATTCAGCTTCGATGCTGATAAATTTACGGTTTTTCgg                       >  2:2216034/1‑100 (MQ=255)
                                                                              gAATGTAAAAAGCCCCGCAACGTGTTGCGGGGCTTTCATCCGTTACCGGGACGCGAAAAACTTATTCAGCTTCGATGCTGATAAATTTACGGTTTTTCgg                       <  2:1680250/100‑1 (MQ=255)
                                                                               taTGTAAAAAGCCCCGCAACGTGTTGCGGGGCTTTCATCCGTTACCGGGACGCGAAAAACTTATTCAGCTTCGATGCTGATAAATTTACGGTTTTTCggg                      <  2:649578/99‑1 (MQ=255)
                                                                               aaTGTAAAAAGCCCCGCAACGTGTTGCGGGGCTTTCATCCGTTACCGGGACGCGAAAAACTTATTCAGCTTCGATGCTGATAAATTTACGGTTTTTCggg                      <  2:258136/100‑1 (MQ=255)
                                                                                 tGTAAAAAGCCCCGCAACGTGTTGCGGGGCTTTCATCCGTTACCGGGACGCGAAAAACTTATTCAGCTTCGATGCTGATAAATTTACGGTTTTTCGGGcc                    <  1:812268/100‑1 (MQ=255)
                                                                                 tGTAAAAAGCCCCGCAACGTGTTGCGGGGCTTTCATCCGTTACCGGGACGCGAAAAACTTATTCAGCTTCGATGATGATAAATTTACGGTTTTTCGGGcc                    >  1:2286065/1‑100 (MQ=255)
                                                                                 ggatgaaaGCCCCGCAACACGTTGCGGGGCTTTCATCCGTTACCGGGACGCGAAAAACTTATTCAGCTTCGATGCTGATAAATTTACGGTTTTTCTGGcc                    >  1:1837906/6‑100 (MQ=38)
                                                                                   atgaaaGCCCCGCAACACGTTGCGGGGCTTTCATCCGTTACCGGGACGCGAAAAACTTATTCAGCTTCGATGCTGATAAATTTACGGTTTTTCGGGCCtt                  >  2:2221270/4‑100 (MQ=38)
                                                                                    aaaaaGCCCCGCAACGTGTTGCGGGGCTTTCATCCGTTACCGGGACGCGAAAAACTTATTCAGCTTCGATGCTGATAAATTTACGGTTTTTCGGGCCttt                 >  1:100311/1‑100 (MQ=255)
                                                                                     aaaaGCCCCGCAACGTGTTGCGGGGCTTTCATCCGTTACCGGGACGCGAAAAACTTATTCAGCTTCGATGCTGATAAATTTACGGTTTTTCGGGCCTTTa                <  1:453496/100‑1 (MQ=255)
                                                                                      aaaGCCCCGCAACACGTTGCGGGGCTTTCATCCGTTACCGGGACGCGAAAAACTTATTCAGCTTCGATGCTGATAAATTTACGGTTTTTCGGGCCTTTaa               >  1:1906866/1‑100 (MQ=255)
                                                                                       aaGCCCCGCAACGTGTTGCGGGGCTTTCATCCGTTACCGGGACGCGAAAAACTTATTCAGCTTCGATGCTGATAAATTTACGGTTTTTCGGGCCTTTAAc              >  2:1228465/1‑100 (MQ=255)
                                                                                        aGCCCCGCAACGTGTTGCGGGGCTTTCATCCGTTACCGGGACGCGAAAAACTTATTCAGCTTCGATGCTGATAAATTTACGGTTTTTCGGGCCTTTAACt             >  2:468559/1‑100 (MQ=255)
                                                                                         gCCCCGCAACGTGTTGCGGGGCTTTCATCCGTTACCGGGACGCGAAAAACTTATTCAGCTTCGATGCTGATAAATTTACGGTTTTTCTGGCCTTTAACtt            >  2:2329712/1‑100 (MQ=255)
                                                                                         gCCCCGCAACACGTTGCGGGGCTTTCATCCGTTACCGGGACGCGAAAAACTTATTCAGCTTCGATGCTGATAAATTTACGGTTTTTCGGGCCTTTAACtt            >  1:1689051/1‑100 (MQ=255)
                                                                                         gCCCCGCAACACGTTGCGGGGCTTTCATCCGTTACCGGGACGCGAAAAACTTATTCAGCTTCGATGCTGATAAATTTACGGTTTTTCGGGCCTTTAACtt            >  2:753403/1‑100 (MQ=255)
                                                                                         gCCCCGCAACACGTTGCGGGGCTTTCATCCGTTACCGGGACGCGAAAAACTTATTCAGCTTCGATGCTGATAAATTTACGGTTTTTCGGGCCTTTAACtt            >  2:1352230/1‑100 (MQ=255)
                                                                                         gCCCCGCAACACGTTGCGGGGCTTTCATCCGTTACCGGGACGCGAAAAACTTATTCAGCTTCGATGCTGATAAATTTACGGTTTTTCGGGCCTTTAACtt            >  2:262503/1‑100 (MQ=255)
                                                                                         gCCCCGCAACACGTTGAGGGGCTTTCATCCGATACCAGGACGCGAAAACCTAATTCAGGTTCGATTCAGATAAATTTACGGTTTTTTGGGCCTTTaa               >  1:528572/1‑97 (MQ=255)
                                                                                         gCCCAGCAACATGTTACGGGGCTTTCCTCCGTTACCGGGACGCGAAAAACTTATTCAGCTTCGATTCTGATAAATGTACGGTTTTTCGGGCCTTTATCtt            >  1:212720/1‑100 (MQ=255)
                                                                                          ccccGCAACACGTTGCGGGGCTTTCATCCGTTACCGGGACGCGAAAAACTTATTCAGCTTCGATGCTGATAAATTTACGGTTTTTCGGGCCTTTAACTTc           >  2:45670/1‑100 (MQ=255)
                                                                                          ccccGCAACACGTTGCGGGGCTTTCATCCGTTACCGGGACGCGAAAAACTTATTCAGCTTCGATGCTGATAAATTTACGGTTTTTCGGGCCTTTAACTTc           >  2:1617981/1‑100 (MQ=255)
                                                                                          ccccGCAACACGTTGCGGGGCTTTCATCCGTTACCGGGACGCGAAAAACTTATTCAGCTTCGACGCTGATAAATTTACGGTTTTTCGGGCCTTTAACTTc           >  2:1099516/1‑100 (MQ=255)
                                                                                           cccGCAACGTGTTGCGGGGCTTTCATCCGTTACCGGGACGCGAAAAACTTATTCAGCTTCGATGCTGATAAATTTACGGTTTTTCGGGCCTTTAACTTCg          >  2:2093250/1‑100 (MQ=255)
                                                                                           cccGCAACGTGTTGCGGGGCTTTCATCCGTTACCGGGACGCGAAAAACTTATTCAGCTTCGATGCTGATAAATTTACGGTTTTTCGGGCCTTTAACTTCg          <  2:379927/100‑1 (MQ=255)
                                                                                           cccGCAACGTGTTGCGGGGCTTTCATCCGTTACCGGGACGCGAAAAACTTATTCAGCTTCGATGCTGATAAATTTACGGTTTTTCGGGCCTTTAACTTCg          >  2:2336697/1‑100 (MQ=255)
                                                                                           cacGCAACGTGTTGCGGGGCTTTCATCCGTTACCGGGACGCGAAAAACTTATTCAGCTTCGATGCTGATAAATTTACGGTTTTTCGGGCCTTTAACTTc           >  2:1751843/3‑99 (MQ=255)
                                                                                            ccGCAACGTGTTGCGGGGCTTTCATCCGTTACCGGGACGCGAAAAACTTATTCAGCTTCGATGCTGATAAATTTACGGTTTTTCGGGCCTTTAACTTCGa         >  2:519276/1‑100 (MQ=255)
                                                                                            ccGCAACGTGTTGCGGGGCTTTCATCCGTTACCGGGACGCGAAAAACTTATTCAGCTTCGATGCTGATAAATTTACGGTTTTTCGGGCCTTTAACTTCGa         <  1:1910066/100‑1 (MQ=255)
                                                                                                   gtgtTGCGGGGCTTTCATCCGTTACCGGGACGCGAAAAACTTATTCAGCTTCGATGCTGATAAATTTACGGTTTTTCGGGCCTTTAACTTCGAATTTCAc  <  2:804617/100‑1 (MQ=255)
                                                                                                   |                                                                                                   
CCTGCTGCTTCATACCGTATTCCCTGCAAAGCGCATTTGCGGATTAACGTCGCACCATGCGCATTTTTACCGGTTATCGAATGTAAAAAGCCCCGCAACGTGTTGCGGGGCTTTCATCCGTTACCGGGACGCGAAAAACTTATTCAGCTTCGATGCTGATAAATTTACGGTTTTTCGGGCCTTTAACTTCGAATTTCAC  >  NC_000913/3332723‑3332921

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: