Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,332,823 T→C 15.1% intergenic (‑88/+39) 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,8230TC15.1% 513.3 / 78.2 220intergenic (‑88/+39)yhbE/rpmAEamA family inner membrane putative transporter/50S ribosomal subunit protein L27
Reads supporting (aligned to +/- strand):  ref base T (98/88);  new base C (10/23);  total (109/111)
Fisher's exact test for biased strand distribution p-value = 2.28e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.75e-01

Only 100 of 226 total aligned reads displayed.

TGCTGCTTCATACCGTATTCCCTGCAAAGCGCATTTGCGGATTAACGTCGCACCATGCGCATTTTTACCGGTTATCGAATGTAAAAAGCCCCGCAACGTGTTGCGGGGCTTTCATCCGTTACCGGGACGCGAAAAACTTATTCAGCTTCGATGCTGATAAATTTACGGTTTTTCGGGCCTTTAACTTCGAATTTCACT  >  NC_000913/3332725‑3332922
                                                                                                  |                                                                                                   
tgctgcTTCATACCGTATTCCCTGCCATGCGCATTTGCGGATTAACGTCGCACCACGTGCATTGTTACCTGTTATGGAATTTAAACAGCCCCGCAACgtg                                                                                                    >  2:1121036/1‑100 (MQ=255)
    gctTCATACCGTATTCCCTGCAAAGCGCATTTGCGGATTAACGTCGCACCATGCGCATTTTTACCGGTTATCGAATGTAAAAAGCCCCGCAACGTGTTGc                                                                                                <  1:1277072/100‑1 (MQ=255)
    gctTCATACCGTATTCCCTGCAAAGCGCATTTGCGGATTAACGTCGCACCATGCGCATTTTTACCGGTTATCGAATGTAAAAAGCCCCGCAACGTGTTGc                                                                                                >  1:1948209/1‑100 (MQ=255)
    gctTCATACCGTATTCCCTGCAAAGCGCATTTGCGGATTAACGTCGCACCATGCGCATTTTTACCGGTTATCGAATGTAAAAAGCCCCGCAACGTGTTGc                                                                                                <  1:2655456/100‑1 (MQ=255)
    gattCATACCGTATTCCCTGCAAAGCGCATTTGCGGATTAACGTCGCACCATGCGCATTTTTACCGGTTATCGAATGTAAAAAGCCCCGCAACGTGTTg                                                                                                 <  1:2468508/97‑1 (MQ=255)
      ttCATACCGTATTCCCTGCAAAGCGCATTTGCGGATTACCGTCGCACCATGCGCATTTTTACCGGTTATCGCATGTAAAAAGCCCCGCACCGTGTTGCgg                                                                                              <  2:3343009/100‑1 (MQ=255)
       tCATACCGTATTCCCTGCAAAGCGCATTTGCGGATTAACGTCGCACCATGCGCATTTTTACCGGTTATCGAATGTAAAAAGCCCCGCAACGTGTTGCggg                                                                                             >  1:1601183/1‑100 (MQ=255)
       tCATACCGTATTCCCTGCAAAGCGCATTTGCGGATTAACGTCGCACCATGCGCATTTTTACCGGTTATCGAATGTAAAAAGCCCCGCAACGTGTTGCggg                                                                                             >  1:2172489/1‑100 (MQ=255)
       tCATACCGTATTCCCTGCAAAGCGCATTTGCGGATTAACGTCGCACCATGCGCATTTTTACCGGTTATCGAATGTAAAAAGCCCCGCAACACGTTGCggg                                                                                             <  2:3806786/100‑1 (MQ=255)
        cATACCGTATTCCCTGCAAAGCGCATTTGCGGATTAACGTCGCACCATGCGCATTTTTACCGGTTATCGAATGTAAAAAGCCCCGCAACACGTTGCgggg                                                                                            <  1:634252/100‑1 (MQ=255)
        cATACCGTATTCCCTGCAAAGCGCATTTGCGGATTAACGTCCCACCATGCGCATTTTTACCGGTTATCGAATGTAAAAAGCCCCGCAACACGTTGCgggg                                                                                            <  2:4123213/100‑1 (MQ=255)
         aTACCGTATTCCCTGCAAAGCGCATTTGCGGATTAACGTCGCACCATGCGCATTTTTACCGGTTATCGAATGTAAAAAGCCCCGCAACGTGTTGCGGGGc                                                                                           <  2:1845761/100‑1 (MQ=255)
         aTACCGTATTCCCTGCAAAGCGCATTTGCGGATTAACGTCGCACCATGCGCATTTTTACCGGTTATCGAATGTAAAAAGCCCCGCAACACGTTGCGGGGc                                                                                           <  1:203006/100‑1 (MQ=255)
         aTACCGTATTCCCTGCAAAGCGCATTTGCGGATTAACGTCGCACCATGCGCATTTTTACCGGTTATCGAATGTAAAAAGCCCCGCAACACGTTGCGGGGc                                                                                           <  1:2124287/100‑1 (MQ=255)
         aTACCGTATTCCCTGCAAAGCGCATTTGCGGATTAACGTCGCACCATGCGCATTTTTACCGGTTATCGAATGTAAAAAGCCCCGCAACACGTTGCGGGGc                                                                                           <  1:1859838/100‑1 (MQ=255)
         aTACCGTATTCCCTGCAAAGCGCATTTGCGGATTAACGTCGCACCATGCGCATTTTTACCGGTTATCGAATGTAAAAAGCCCCGCAACACGTTGCGGGGc                                                                                           <  1:232942/100‑1 (MQ=255)
          tACCGTATTCCCTGCAAAGCGCATTTGCGGATTAACGTCGCACCATGCGCATTTTTACCGGTTATCGAATGTAAAAAGCCCCGCACCGTATTGCGGGGc                                                                                           >  2:2065880/1‑99 (MQ=255)
          tACCGTATTCCCTGCAAAGCGCATTTGCGGATTAACGTCGCACCATGCGCATTTTTACCGGTTATCGAATGTAAAAAGCCCCGCAACGTGTTGCGGGGCt                                                                                          >  2:1626574/1‑100 (MQ=255)
          tACCGTATTCCCTGCAAAGCGCATTTGCGGATTAACGTCGCACCATGCGCATTTTTACCGGTTATCGAATGTAAAAAGCCCCGCAACACGTTGCGGGGCt                                                                                          <  2:1703231/100‑1 (MQ=255)
          tACCGTACTCCCTGCAAAACACATTTTCGCATTAACGTTGCAACATTCGCATTTTTACCGGTTATCGAATGTAAAAAGCCCCGCAACACTTTGCgggg                                                                                            <  2:3890186/98‑1 (MQ=255)
            ccGTATTCCCTGCAAAGCGCATTTGCGGATTAACGTCGCACCATGCGCATTTTTACCGGTTATCGAATGTAAAAAGCCCCGCAACGTGTTGCGGGGCttt                                                                                        >  1:3854658/1‑100 (MQ=255)
            ccGTATTCCCTGCAAAGCGCATTTGCGGATTAACGTCGCACCATGCGCATTTTTACCGGTTATCGAATGTAAAAAGCCCCGCAACGTGTTGCGGGGCttt                                                                                        <  1:1656850/100‑1 (MQ=255)
             cGTATTCCCTGCAAAGCGCATTTGCGGATTAACGTCGCACCATGCGCATTTTTACCGGTTATCGAATGTAAAAAGCCCCGCAACACGTTGCGGGGCtttt                                                                                       <  1:3693414/100‑2 (MQ=255)
              gTATTCCCTGCAAAGCGCATTTGCGGATTAACGTCGCACCATGCGCATTTTTACCGGTTATCGAATGTAAAAAGCCCCGCAACGTGTTGCGGGGCTTTCa                                                                                      <  1:2086107/100‑1 (MQ=255)
                aTTCCCTGCAAAGCGCATTTGCGGATTAACGTCGCACCATGCGCATTTTTACCGGTTATCGAATGTAAAAAGCCCCGCAACGTGTTGCGGGGCTTTCATc                                                                                    >  2:3703021/1‑100 (MQ=255)
                 ttCCCTGCAAAGCGCATTTGCGGATTAACGTCGCACCATGCGCATTTTTACCGGTTATCGAATGTAAAAAGCCCCGCAACGTGTTGCGGGGCTTTCATcc                                                                                   <  1:641076/100‑1 (MQ=255)
                   cgcagggagtcgcACTTGCTGGTCAACGTCGCATTATAAAAATTTTTACCGGTTATCGAACGTAAAAACACCCGCAACGTGTTGCGGGGCTTTCATCCGt                                                                                 >  2:767600/11‑100 (MQ=255)
                    ccTGCAAAGCGCATTTGCGGATTAACGTCGCACCATGCGCATTTTTACCGGTTATCGAATGTAAAAAGCCCCGCAACGTGTTGCGGGGCTTTCATCCGtt                                                                                >  1:1613246/1‑100 (MQ=255)
                     cTGCAAAGCGCATTTGCGGATTAACGTCGCACCATGCGCATTTTTACCGGTTATCGAATGTAAAAAGCCCCGCAACGTGTTGCGGGGCTTTCATCCGTTa                                                                               >  1:3548804/1‑100 (MQ=255)
                      tGCAAAGCGCTTTTGCGGATTAACGTCGCACCATGCGCATTTTTACCGGTTATCGAATGTAAAAAGTCTCGCAACGTGTTGCGGGGCTTTCATCCg                                                                                  <  2:3485344/96‑1 (MQ=255)
                      tGCAAAGCGCATTTGCGGATTAACGTCGCACCATGCGCATTTTTACCGGTTATCGAATGTAAAAAGCCCCGCAACGTGTTGCGGGGCTTTCATCCGTTAc                                                                              >  1:3107883/1‑100 (MQ=255)
                       gCAAAGCGCATTTGCGGATTAACGTCGCACCATGCGCATTTTTACCGGTTATCGAATGTAAAAAGCCCCGCAACGTGTTGCGGGGCTTTCATCCGTTAcc                                                                             <  1:2015747/100‑1 (MQ=255)
                       gCAAAGCGCATTTGCGGATTAACGTCGCACCATGCGCATTTTTACCGGTTATCGAATGTAAAAAGCCCCGCAACGTGTTGCGGGGCTTTCATCCGTTAcc                                                                             <  1:3568848/100‑1 (MQ=255)
                         aaaGCGCATTTGCGGATTAACGTCGCACCATGCGCATTTTTACCGGTTATCGAATGTAAAAAGCCCCGCAACGTGTTGCGGGGCTTTCATCCGTTACCgg                                                                           <  2:3107870/100‑1 (MQ=255)
                          aaGCGCATTTGCGGATTAACGTCGCACCATGCGCATTTTTACCGGTTATCGAATGTAAAAAGCCCCGCAACGTGTTGCGGGGCTTTCATCCGTTACCggg                                                                          >  2:385830/1‑100 (MQ=255)
                           aGCGCATTTGCGGATTAACGTCGCACCATGCGCATTTTTACCGGTTATCGAATGTAAAAAGCCCCGCAACGTGTTGCGGGGCTTTCATCCGTTACCGGGa                                                                         <  1:280012/100‑1 (MQ=255)
                           aGCGCATTTGCGGATTAACGTAGCACAATGCGCGTTTTTACCGGTTATCGAATGTAAAAAACCCCGCAACGTGTTGCGGGGCTTTCATACGTTACCgg                                                                           >  2:2953490/1‑98 (MQ=255)
                                aTTTGCGGATTAACGTCGCACCATGCGCATTTTTACCGGTTATCGAATGTAAAAAGCCCCGCAACGTGTTGCGGGGCTTTCATCCGTTACCGGGACGCGa                                                                    <  1:4195913/100‑1 (MQ=255)
                                  ttGCGGATTAACGTCGCACCATGCGCATTTTTACCGGTTATCGAATGTAAAAAGCCCCGCAACGTGTTGCGGGGCTTTCATCCGTTACCGGGACGCGaaa                                                                  >  2:3157328/1‑100 (MQ=255)
                                    gCGGATTAACGTCGCACCATGCGCATTTTTACCGGTTATCGAATGTAAAAAGCCCCGCAACGTGTTGCGGGGCTTTCATCCGTTACCGGGACGCGaaaaa                                                                >  2:4053804/1‑100 (MQ=255)
                                     cGGATTAACGTCGCACCATGCGCATTTTTACCGGTTATCGAATGTAAAAAGCCCCGCAACGTGTTGCGGGGCTTTCATCCGTTACCGGGACGCGAAAAAc                                                               >  1:1195014/1‑100 (MQ=255)
                                      ggATTAACGTCGCACCATGCGCATTTTTACCGGTTATCGAATGTAAAAAGCCCCGCAACGTGTTGCGGGCCTTTCATCCGTTACCGGGCCGCGAAAAACt                                                              >  2:519387/1‑100 (MQ=255)
                                          tAACGTCGCACCATGCGCATTTTTACCGGTTATCGAATGTAAAAAGCCCCGCAACGTGTTGCGGGGCTTTCATCCGTTACCGGGACGCGAAAAACTTAtt                                                          <  1:3912449/100‑1 (MQ=255)
                                           aaCGTCGCACCATGCGGATTTTTACCGGTGATCGAATGTAAAAAGCCCCGCAACGTGTTGCGGGGCTTTCATCCGTTAAAAGGAAAAGAAAAACTTAtta                                                         <  2:2813225/100‑2 (MQ=255)
                                           aaCGTCGCACCATGCGCATTTTTACCGGTTATCGAATGTAAAAAGCCCCGCAACGTGTTGCGGGGCTTTCATCCGTTACCGGGACGCGAAAAACTTATTc                                                         >  1:2988482/1‑100 (MQ=255)
                                           aaCGTCGCACCATGCGCATTTTTACCGGTTATCGAATGTAAAAAGCCCCGCAACGTGTTGCGGGGCTTTCATCCGTTACCGGGACGCGAAAAACTTATTc                                                         >  1:1416918/1‑100 (MQ=255)
                                              gTCGCACCATGCGCATTTTTACCGGTTATCGAATGTAAAAAGCCCCGCAACGTGTTGCGGGGCTTTCATCCGTTACCGGGACGCGAAAAACTTATTCAGc                                                      <  1:673294/100‑1 (MQ=255)
                                                  cACCATGCGCATTTTTACCGGTTATCGAATGTAAAAAGCCCCGCAACGTGTTGCGGGGCTTTCATCCGTTACCGGGACGCGAAAAACTTATTCAGCTTCg                                                  <  1:2860295/100‑1 (MQ=255)
                                                  cACCATGCGCATTTTTACCGGTTATCGAATGTAAAAAGCCCCGCAACGTGTTGCGGGGCTTTCATCCGTTACCGGGACGCGAAAAACTTATTCAGCTTCg                                                  <  2:3917375/100‑1 (MQ=255)
                                                        gcgcATTTTTACCGGTTATCGAATGTAAAAAGCCCCGCAACGTGTTGCGGGGCTTTCATCCGTTACCGGGACGCGAAAAACTTATTCAGCTTCGATGCTg                                            >  1:1509638/1‑100 (MQ=255)
                                                        gcgcATTTTTACCGGTTATCGAATGTAAAAAGCCCCGCAACGTGTTGCGGGGCTTTCATCCGTTACCGGGACGCGAAAAACTTATTCAGCTTCGATGCTg                                            >  2:1750078/1‑100 (MQ=255)
                                                         cgcATTTTTACCGGTTATCGAATGTAAAAAGCCCCGCAACGTGTTGCGGGGCTTTCATCCGTTACCGGGACGCGAAAAACTTATTCAGCTTCGATGCTGa                                           <  2:347583/100‑1 (MQ=255)
                                                           cATTTTTACCGGTTATCGAATGTAAAAAGCCCCGCAACGTGTTGCGGGGCTTTCATCCGTTACCGGGACGCGAAAAACTTATTCAGCTTCGATGCTGATa                                         >  1:2713959/1‑100 (MQ=255)
                                                             tttttACCGGTTATCGAATGTAAAAAGCCCCGCAACGTGTTGCGGGGCTTTCATCCGTTACCGGGACGCGAAAAACTTATTCAGCTTCGATGCTGATaaa                                       >  2:83762/1‑100 (MQ=255)
                                                               tttACCGGTTATCGAATGTAAAAAGCCCCGCAACGTGTTGCGGGGCTTTCATCCGTTACCGGGACGCGAAAAACTTATTCAGCATCGAAGCTGaataagt                                     <  1:1810886/100‑7 (MQ=255)
                                                                   ccGGTTATCGAATGTAAAAAGCCCCGCAACGTGTTGCGGGGCTTTCATCCGTTACCGGGACGCGAAAAACTTATTCAGCTTCGATGCTGATAAATTTACg                                 <  1:3222042/100‑1 (MQ=255)
                                                                   ccGGTTATCGAATGTAAAAAGCCCCGCAACGTGTTGCGGGGCTTTCATCCGTTACCGGGACGCGAAAAACTTATTCAGCTTCGATGCTGATAAATTTACg                                 >  1:27500/1‑100 (MQ=255)
                                                                     ngTTATCGAATGTAAAAAGCCCCGCAACGTGTTGCGGGGCTTTCATCCGTTACCGGGACGCGAAAAACTTATTCAGCTTCGATGCTGATAAATTTACGGt                               >  1:2557298/2‑100 (MQ=255)
                                                                     ggTTATCGAATGTAAAAAGCCCCGCAACGTGTTGCGGGGCTTTCATCCGTTACCGGGACGCGAAAAACTTATTCAGCTTCGATGCTGATAAATTTACgat                               >  2:2537765/1‑98 (MQ=255)
                                                                        tagcGAATGTAAAAAGCCCCGCAACGTGTCGCGGGGCTTTCATCGGTTGCCGGGACGTGACAACTTTATCCAGCTGCTCTTCTGATTAATGTACggat                              >  2:3586797/4‑96 (MQ=255)
                                                                        tATCGAATGTAAAAAGCCCCGCAACGTGTTGCGGGGCTTTCATCCGTTACCGGGACGCGAAAGACTTATTCAGCTTCGATGCTGATAAATTTACGGtttt                            >  2:2574007/1‑100 (MQ=255)
                                                                         aTCGAATGTAAAAAGCCCCGCAACGTGTTGCGGGGCTTTCATCCGTTACCGGGACGCGAAAAACTTATTCAGCTTCGATGCTGATAAATTTACGGttttt                           >  1:811758/1‑100 (MQ=255)
                                                                          tCGAATGTAAAAAGCCCCGCAACGTGTTGCGGGGCTTTCATCCGTTACCGGGACGCGAAAAACTTATTCAGCTTCGATGCTGATAAATTTACGGTTTTTc                          >  1:1065048/1‑100 (MQ=255)
                                                                           cGAATGTAAAAAGCCCCGCAACGTGTTGCGGGGCTTTCATCCGTTACCGGGACGCGAAAAACTTATTCAGCTTCGATGCTGATAAATTTACGGTTTTTCg                         <  1:50982/100‑1 (MQ=255)
                                                                             aaTGTAAAAAGCCCCGCAACGTTTTGCGTGGTTTTCATCCGTTACAGGGACGCGAAAAACTTATTCAGCTTCGATGCTGATAAATTTACGGTTTTTCg                         <  2:2656812/98‑1 (MQ=255)
                                                                              aTGTAAAAAGCCCCGCAACGTGTTGCGGGGCTTTCATCCGTTACCGGGACGCGAAAAACTTATTCAGCTTCGATGCTGATAAATTTACGGTTTTTCGGGc                      <  1:915483/100‑1 (MQ=255)
                                                                               tGTAAAAAGCCCCGCAACGTGTTGCGGGGCTTTCATCCGTTACCGGGACGCGAAAAACTTATTCAGCTTCGATGCTGATAAATTTACGGTTTTTCGGGcc                     <  1:2228150/100‑1 (MQ=255)
                                                                               tGTAAAAAGCCCCGCAACGTGTTGCGGGGCTTTCATCCGTTACCGGGACGCGAAAAACTTATTCAGCTTCGATGCTGATAAATTTACGGTTTTTCGGGcc                     >  1:4056155/1‑100 (MQ=255)
                                                                                 tAAAAAGCCCCGCAACGTGTTGCGGGGCTTTCATCCGTTACCGGGACGCGAAAACCTTATTCAGCTTCGATGCTGATAGATTTACGGTTTTTCGGGCCtt                   >  2:3550694/1‑100 (MQ=255)
                                                                                 tAAAAAGCCCCGCAACGTGTTGCGGGGCTTTCATCCGTTACCGGGACGCGAAAAACTTATTCAGCTTCGATGCTGATAAATTTACGGTTTTTCGGGcc                     >  1:2866939/1‑98 (MQ=255)
                                                                                 tAAAAAGCCCCGCAACGTGTTGCGGGGCTTTCATCCGTTACCGGGACGCGAAAAACTTATTCAGCTTCGATGCTGATAAATTTACGGTTTTTCGGGCCtt                   <  2:181713/100‑1 (MQ=255)
                                                                                 tAAAAAGCCCCGCAACGTGTTGCGGGGCTTTCATCCGTTACCGGGACGCGAAAAACTTATTCAGCTTCGATGCTGATAAATTTACGGTTTTTCGGGCCtt                   >  2:3378282/1‑100 (MQ=255)
                                                                                 tAAAAAGCCCCGCAACGTGTTGCGGGGCTTTCATCCGTTACCGGGACGCGAAAAACTTATTCAGCTTCGATGCTGATAAATTTACGGTTTTTCGGGCCtt                   >  1:608938/1‑100 (MQ=255)
                                                                                 tAAAAAGCCCCGCAACGTGTTGCGGGGCTTTCATCCGTTACCGGGACGCGAAAAACTTATTCAGCTTCGATGCTGATAAATTTACGGTTTTTCGGGCCtt                   >  2:1240685/1‑100 (MQ=255)
                                                                                  aaaaaGCCCCGCAACGTGTTGCGGGGCTTTCATCCGTTACCGGGACGCGAAAAACTTATTCAGCTTCGATGCTGATAAATTTACGGTTTTTCGGGCCttt                  <  1:1767314/100‑1 (MQ=255)
                                                                                  aaaaaGCCCCGCAACGTGTTGCGGGGCTTTCATCCGTTACCGGGACGCGAAAAACTTATTCAGCTTCGATGCTGATAAATTTACGGTTTTTCGGGCCttt                  >  2:3687885/1‑100 (MQ=255)
                                                                                  aaaaaGCCCCGCAACGTGTTGCGGGGCTTTCATCCGTTACCGGGACGCGAAAAACTTATTCAGCTTCGATGCTGATAAATTTACGGTTTTTCGGGCCttt                  >  1:259995/1‑100 (MQ=255)
                                                                                   aaaaGCCCCGCAACGTGTTGCGGGGCTTTCATCCGTTACCGGGACGCGAAAAACTTATTCAGCTTCGATGCTGATAAATTTACGGTTTTTCGGGCCTTTa                 >  1:3979666/1‑100 (MQ=255)
                                                                                    aaaGCCCCGCAACGTGTTGCGGGGCTTTCATCCGTTACCGGGACGCGAAAAACTTATTCAGCTTCGATGCTGATAAATTTACGGTTTTTCGGGCCTTTa                 >  2:2308755/1‑99 (MQ=255)
                                                                                     aaGCCCCGCAACACGTTGCGGGGCTTTTATCCGTTACCGGGACGCGAAAAACTTATTCAGCTTCGATGCTGATAAATTTACGGTTTTTCGGGCAtttcc                >  2:1376590/1‑97 (MQ=38)
                                                                                       gCCCTGCAACACGTTGCGGGGCTTTCATCCGTTACCGGGACGCGAAAAACTTATTCAGCTTCGATGCTGATAAATTTACGGTTTTTCGGGCCTTTAACtt             >  2:4086379/1‑100 (MQ=255)
                                                                                       gCCCCGCAACGTGTTGCGGGGCTTTCATCCGTTACCGGGACGCGAAAAACTTATTCAGCTTCGATGCTGATAAATTTACGGTTTTTCGGGCCTTTAACtt             >  1:1343466/1‑100 (MQ=255)
                                                                                       gCCCCGCAACGTGTTGCGGGGCTTTCATCCGTTACCGGGACGCGAAAAACTTATTCAGCTTCGATGCTGATAAATTTACGGTTTTTCGGGCCTTTAACtt             >  1:1694571/1‑100 (MQ=255)
                                                                                       gCCCCGCAACACGTTGCGGGGCTTTCATCCGTTACCGGGACGCGAAAAACTTATTCAGCTTCGATGCTGATAAATTTACGGTTTTTCGGGCCTTTAACtt             >  1:2875244/1‑100 (MQ=255)
                                                                                       gCCCCGCAACAAGTTGCGGGGCTTTCATCCGTTACCGGGACGCGAAAAACTTATTCAGCTTCGATGCTGATAAATTTACGGTTTTTCGGGCCTTTAACtt             >  2:2194377/1‑100 (MQ=255)
                                                                                        tccgatctcGTGTTGCGGGGCTTTCATCCGTTACCGGGACGCGAAAAACTTATTCAGCTTCGATGCTGATAAATTTACGGTTTTTCGGGCCTTTAACTTc            <  1:2933284/92‑1 (MQ=255)
                                                                                        ccccGCAACGTGTTGCGGGGCTTTCATCCGTTACCGGGACGCGAAAAACTTATTCAGCTTCGATGCTGATAAATTTACGGTTTTTCGGGCCTTTAACttn            <  2:3072771/100‑2 (MQ=255)
                                                                                        ccccGCAACGTGTTGCGGGGCTTTCATCCGTTACCGGGACGCGAAAAACTTATTCAGCTTCGATGCTGATAAATTTACGGTTTTTCGGGCCTTTAACTTc            >  1:1005697/1‑100 (MQ=255)
                                                                                        ccccGCAACGTGTTGCGGGGCTTTCATCCGTTACCGGGACGCGAAAAACTTATTCAGCTTCGATGCTGATAAATTTACGGTTTTTCGGGCCTTTAACTTc            >  2:788575/1‑100 (MQ=255)
                                                                                        ccccGCAACGTGTTGCGGGGCTTTCATCCGTTACCGGGACGCGAAAAACTTATTCAGCTTCGATGCTGATAAATTTACGGTTTTTCGGGCCTTTAACTTc            >  1:974061/1‑100 (MQ=255)
                                                                                        ccccGCAACGTGTTGCGGGGCTTTCATCCGTTACCGGGACGCGAAAAACTTATTCAGCTTCGATGCTGATAAATTTACGGTTTTTCGGGCCTTTAACTTc            >  1:3460906/1‑100 (MQ=255)
                                                                                        ccccGCAACACGTTGCGGGGCTTTCATCCGTTACCGGGACGCGAAAAACTTATTCAGCTTCGATGCTGATAAATTTACGGTTTTTCGGGCCTTTAACTTc            >  2:2375222/1‑100 (MQ=255)
                                                                                         cccGCAACGTGTTGCGGGGCTTTCATCCGTTACCGGGACGCGAAAAACTTATTCAGCTTCGATGCTGATAAATTTACGGTTTTTCGGGCCTTTAACTTCg           >  2:2032063/1‑100 (MQ=255)
                                                                                          ccGCAACGTGTTGCGGGGCTTTCATCCGTTACCGGGACGCGAAAAACTTATTCAGCTTCGATGCTGATAAATTTACGGTTTTTCGGGCCTTTAACTTCg           <  1:882960/99‑1 (MQ=255)
                                                                                             cAACGTGTTGCGGGGCTTTCATCCGTTACCGGGACGCGAAAAACTTATTCAGCTTCGATGCTGATAAATTTACGGTTTTTCGGGCCTTTAACTTCGAAtt       <  2:3983530/100‑1 (MQ=255)
                                                                                              aaCGTGTTGCGGGGCTTTCATCCGTTACCGGGACGCGAAAAACTTATTCAGCTTCGATGCTGATAAATTTACGGTTTTTCGGGCCTTTAACTTCGAAttt      <  1:1790463/100‑1 (MQ=255)
                                                                                               aCGTGTTGCGGGGCTTTCATCCGTTACCGGGACGCGAAAAACTTATTCAGCTTCGATGCTGATAAATTTACGGTTTTTCGGGCCTTTAACTTCGAATTTc     <  1:3396227/100‑1 (MQ=255)
                                                                                                cGTGTTGCGGGGCTTTCATCCGTTACCGGGACGCGAAAAACTTATTCAGCTTCGATGCTGATAAATTTACGGTTTTTCGGGCCTTTAACTTcagatcgga    >  2:2933274/1‑92 (MQ=255)
                                                                                                  tgtTGCGGGGCTTTCATCCGTTACCGGGACGCGAAAAACTTATTCAGCTTCGATGCTTATAAATTTACTGTTTTTCGGGCCTTTAATTTCGAATTTCatt  <  1:2832245/100‑3 (MQ=255)
                                                                                                  tgtTGCGGGGCTTTCATCCGTTACCGGGACGCGAAAAACTTATTCAGCTTCGATGCTGATAAATTTACGGTTTTTCGGGCCTTTAACTTCGAATCTCACt  >  1:1117559/1‑100 (MQ=255)
                                                                                                  |                                                                                                   
TGCTGCTTCATACCGTATTCCCTGCAAAGCGCATTTGCGGATTAACGTCGCACCATGCGCATTTTTACCGGTTATCGAATGTAAAAAGCCCCGCAACGTGTTGCGGGGCTTTCATCCGTTACCGGGACGCGAAAAACTTATTCAGCTTCGATGCTGATAAATTTACGGTTTTTCGGGCCTTTAACTTCGAATTTCACT  >  NC_000913/3332725‑3332922

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: