Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,332,822 G→A 15.7% 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,8220GA15.7% 548.9 / 71.8 220intergenic (‑87/+40)yhbE/rpmAEamA family inner membrane putative transporter/50S ribosomal subunit protein L27
Reads supporting (aligned to +/- strand):  ref base G (98/86);  new base A (11/24);  total (109/111)
Fisher's exact test for biased strand distribution p-value = 2.59e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 224 total aligned reads displayed.

TGCTGCTTCATACCGTATTCCCTGCAAAGCGCATTTGCGGATTAACGTCGCACCATGCGCATTTTTACCGGTTATCGAATGTAAAAAGCCCCGCAACGTGTTGCGGGGCTTT‑‑‑CATCCGTTACCGGGACGCGAAAAACTTATTCAGCTTCGATGCTGATAAATTTACGGTTTTTCGGGCCTTTAACTTCGAATTTCAC  >  NC_000913/3332725‑3332921
                                                                                                 |                                                                                                       
tgctgcTTCATACCGTATTCCCTGCCATGCGCATTTGCGGATTAACGTCGCACCACGTGCATTGTTACCTGTTATGGAATTTAAACAGCCCCGCAACgtg                                                                                                       >  2:1121036/1‑100 (MQ=255)
tgctgcTTCATACCGTATTCCCTGCAAAGCGCATTTGCGGATTAACGTCGCACCATGCGCATTTTTACCGGTTATCGAATGTAAAAAGCCCCGCAACgtg                                                                                                       >  1:4043088/1‑100 (MQ=255)
 gctgctTCATACCGTATTCCCTGCAAAGCGCATTTGCGGATTAACGTCGCACCATGCGCATTTTTACCGGTTATCGAATGTAAAAAGCCCCGCAACgtgt                                                                                                      <  1:3293837/100‑1 (MQ=255)
    gctTCATACCGTATTCCCTGCAAAGCGCATTTGCGGATTAACGTCGCACCATGCGCATTTTTACCGGTTATCGAATGTAAAAAGCCCCGCAACGTGTTGc                                                                                                   <  1:2655456/100‑1 (MQ=255)
    gctTCATACCGTATTCCCTGCAAAGCGCATTTGCGGATTAACGTCGCACCATGCGCATTTTTACCGGTTATCGAATGTAAAAAGCCCCGCAACGTGTTGc                                                                                                   <  1:1277072/100‑1 (MQ=255)
    gctTCATACCGTATTCCCTGCAAAGCGCATTTGCGGATTAACGTCGCACCATGCGCATTTTTACCGGTTATCGAATGTAAAAAGCCCCGCAACGTGTTGc                                                                                                   >  1:1948209/1‑100 (MQ=255)
    gctTCATACCGTATTCCATGCAAAGCGCATTTGCGGATTAACGTCGCACCATGCGCATTTTTACCGGTAACGGATG‑‑‑AAAGCCCCGCCACACGTTGCgggg                                                                                               <  2:2939204/100‑1 (MQ=21)
    gattCATACCGTATTCCCTGCAAAGCGCATTTGCGGATTAACGTCGCACCATGCGCATTTTTACCGGTTATCGAATGTAAAAAGCCCCGCAACGTGTTg                                                                                                    <  1:2468508/97‑1 (MQ=255)
     ctTCATACCGTATTCCCTGCAAAGCGCATTTGCGGATTAACGTCGCACCATGCGCATTTTTACCGGTTATCGAATGTAAAAAGCCCCGCAACGTGTTGCg                                                                                                  <  1:119979/100‑1 (MQ=255)
       tCATACCGTATTCCCTGCAAAGCGCATTTGCGGATTAACGTCGCACCATGCGCATTTTTACCGGTTATCGAATGTAAAAAGCCCCGCAACGTGTTGCggg                                                                                                >  1:1601183/1‑100 (MQ=255)
        cATCCCGTATCCCCTGCAAAGCTCATTTGCAGATTAACTTCGCACCATGCGCATTTTTACCGGTTATCGAATGTAAAAAACCCCGCAACACGTTGCgggg                                                                                               <  2:2556424/100‑1 (MQ=255)
        cATACCGTATTCCCTGCAAAGCGCATTTGCGGATTAACGTCGCACCATGCGCATTTTTACCGGTTATCTAATGTAAAAAGCCCCGCTACACGTTGCggtg                                                                                               <  1:2824605/100‑3 (MQ=255)
        cATACCGTATTCCCTGCAAAGCGCATTTGCGGATTAACGTCGCACCATGCGCATTTTTACCGGTTATCGAATGTAAAAAGCCCCGCAACGTGTTGCgggg                                                                                               <  1:906797/100‑1 (MQ=255)
        cATACCGTATTCCCTGCAAAGCGCATTTGCGGATTAACGTCGCACCATGCGCATTTTTACCGGTTATCGAATGTAAAAAGCCCCGCAACACGTTGCgggg                                                                                               <  1:621475/100‑1 (MQ=255)
        cATACCGTATTCCCTGCAAAGCGCATTTGCGGATTAACGTCGCACCATGCGCATTTTTACCGGTTATCGAATGTAAAAAGCCCCGCAACACGTTGCgggg                                                                                               <  1:1637839/100‑1 (MQ=255)
        cATACCGTATTCCCTGCAAAGCGCATTTGCGGATTAACGTCGCACCATGCGCATTTTTACCGGTTATCGAATGTAAAAAGCCCCGCAACACGTTGCgggg                                                                                               <  1:228128/100‑1 (MQ=255)
        cATACCGTATTCCCTGCAAAGCGCATTTGCGGATTAACGTCGCACCATGCGCATTTTTACCGGTTATCGAATGTAAAAAGCCCCGCAACACGTTGCgggg                                                                                               <  1:2594898/100‑1 (MQ=255)
        cATACCGTATTCCCTGCAAAGCGCATTTGCGGATTAACGTCCCACCATGCGCATTTTTACCGGTTATCGAATGTAAAAAGCCCCGCAACACGTTGCgggg                                                                                               <  2:4123213/100‑1 (MQ=255)
         aTACCGTATTCCCTGCAAAGCGCATTTGCGGATTAACGTCGCACCATGCGCATTTTTACCGGTTATCGAATGTAAAAAGCCCCGCAACACGTTGCGGGGc                                                                                              <  1:1860096/100‑1 (MQ=255)
         aTACCGTATTCCCTGCAAAGCGCATTTGCGGATTAACGTCGCACCATGCGCATTTTTACCGGTTATCGAATGTAAAAAGCCCCGCAACACGTTGCGGGGc                                                                                              <  1:2124287/100‑1 (MQ=255)
         aTACCGTAGTCCCTGCAAAGCGCATTTGCGGAGTAACGTCCCACCATGCGCATTTTTACCGGTTATCGAATGTAAAAAGCCCCGCAACGTGTTGCGGGGc                                                                                              <  2:3832714/100‑1 (MQ=255)
          tACCGTATTCCCTGCAAAGCGCATTTGCGGATTAACGTCGCACCATGCGCATTTTTACCGGTTATCGAATGTAAAAAGCCCCGCAACGTGTTGCGGGGCt                                                                                             >  1:2076085/1‑100 (MQ=255)
          tACCGTATTCCCTGCAAAGCGCATTTGCGGATTAACGTCGCACCATGCGCATTTTTACCGGTTATCGAATGTAAAAAGCCCCGCAACACGTTGCGGGGCt                                                                                             <  2:1703231/100‑1 (MQ=255)
          tACCGTACTCCCTGCAAAACACATTTTCGCATTAACGTTGCAACATTCGCATTTTTACCGGTTATCGAATGTAAAAAGCCCCGCAACACTTTGCgggg                                                                                               <  2:3890186/98‑1 (MQ=255)
              gTATTCCCTGCAAAGCGCATTTGCGGATTAACGTCGCACCATGCGCATTTTTACCGGTTATCGAATGTAAAAAGCCCCGCAACGTGTTGCGGGGCTTT‑‑‑Ca                                                                                      <  1:2895714/100‑1 (MQ=255)
                aTTCCCTGCAAAGCGCATTTGCGGATTAACGTCGCACCATGCGCATTTTTACCGGTTATCGAATGTAAAAAGCCCCGCAACGTGTTGCGGGGCTTT‑‑‑CATc                                                                                    <  1:451562/100‑1 (MQ=255)
                aTTCCCTGCAAAGCGCATTTGCGGATTAACGTCGCACCATGCGCATTTTTACCGGTTATCGAATGTAAAAAGCCCCGCAACGTGTTGCGGGGCTTT‑‑‑CATc                                                                                    >  2:3703021/1‑100 (MQ=255)
                 ttCCCTGCAAAGCGCATTTGCGTATTAACTTCGCACCATGCGCATTTTTACCGGTTATCGAATGTAAAAAGCCCCGCAACGTGTTGCGGGGCTTT‑‑‑CATcc                                                                                   <  1:3128170/100‑1 (MQ=255)
                   cgcagggagtcgcACTTGCTGGTCAACGTCGCATTATAAAAATTTTTACCGGTTATCGAACGTAAAAACACCCGCAACGTGTTGCGGGGCTTT‑‑‑CATCCGt                                                                                >  2:767600/11‑100 (MQ=255)
                      tGCAAAGCGCTTTTGCGGATTAACGTCGCACCATGCGCATTTTTACCGGTTATCGAATGTAAAAAGTCTCGCAACGTGTTGCGGGGCTTT‑‑‑CATCCg                                                                                  <  2:3485344/96‑1 (MQ=255)
                       gCAAAGCGCATTTGCGGATTAACGTCGCACCATGCGCATTTTTACCGGTTATCGAATGTAAAAAGCCCCGCAACGTGTTGCGGGGCTTT‑‑‑CATCCGTTAcc                                                                            >  1:3032180/1‑100 (MQ=255)
                         aaaGCGCATTTGCGGATTAACGTCGCACCATGCGCATTTTTACCGGTTATCGAATGTAAAAAGCCCCGCAACGTGTTGCGGGGCTTT‑‑‑CATCCGTTACCgg                                                                          <  2:3107870/100‑1 (MQ=255)
                          aaGCGCATTTGCGGATTAACGTCGCACCATGCGCATTTTTACCGGTTATCGAATGTAAAAAGCCCCGCAACGTGTTGCGGGGCTTT‑‑‑CATCCGTTACCggg                                                                         >  1:2872910/1‑100 (MQ=255)
                          aaGCGCATTTGCGGATTAACGTCGCACCATGCGCATTTTTACCGGTTATCGAATGTAAAAAGCCCCGCAACGTGTTGCGGGGCTTT‑‑‑CATCCGTTACCggg                                                                         >  2:385830/1‑100 (MQ=255)
                            tcgcATGTGCGCAGGAACGGCGCACCATGCCGATATTGACCGGGTATCGAATGTAAAAAGTCCCGCAACGTGTTGAAGGGTTTT‑‑‑ATACCGTTTCCgg                                                                          <  2:1812669/96‑1 (MQ=255)
                            gcgcATTTGCGGATTAACGTCGCACCATGCGCATTTTTACCGGTTATCGAATGTAAAAAGCCCCGCAACGTGTTGCGGGGCTTT‑‑‑CATCCGTTACCGGGAc                                                                       <  2:2117144/100‑1 (MQ=255)
                             cgcATTTGCGGATTAACGTCGCACCATGCGCATTTTTACCGGTTATCGAATGTAAAAAGCCCCGCAACGTGTTGCGGGGCTTT‑‑‑CATCCGTTACCGGGAcg                                                                      >  1:3858347/1‑100 (MQ=255)
                                  ttGCGGATTAACGTCGCACCATGCGCATTTTTACCGGTTATCGAATGTAAAAAGCCCCGCAACGTGTTGCGGGGCTTT‑‑‑CATCCGTTACCGGGACGCGaaa                                                                 >  2:3157328/1‑100 (MQ=255)
                                    gCGGATTAACGTCGCACCATGCGCATTTTTACCGGTTATCGAATGTAAAAAGCCCCGCAACGTGTTGCGGGGCTTT‑‑‑CATCCGTTACCGGGACGCGaaaaa                                                               >  2:4053804/1‑100 (MQ=255)
                                     cGGATTAACGTCGCACCATGCGCATTTTTACCGGTTATCGAATGTAAAAAGCCCCGCAACGTGTTGCGGGGCTTT‑‑‑CATCCGTTACCGGGACGCGAAAAAc                                                              >  2:183811/1‑100 (MQ=255)
                                     cGGATTAACGTCGCACCATGCGCATTTTTACCGGTTATCGAATGTAAAAAGCCCCGCAACGTGTTGCGGGGCTTT‑‑‑CATCCGTTACCGGGACGCGAAAAAc                                                              <  2:231867/100‑1 (MQ=255)
                                      ggATTAACGTCGCACCATGCGCATTTTTACCGGTTATCGAATGTAAAAAGCCCCGCAACGTGTTGCGGGGCTTT‑‑‑CATCCGTTACCGGGACGCGAAAAACt                                                             >  2:1008613/1‑100 (MQ=255)
                                        aTTAACGTCGCACCATGCGCATTTTTACCGGTTATCGAATGTAAAAAGCCCCGCAACGTGTTGCGGGGCTTT‑‑‑CATCCGTTACCGGGACGCGAAAAACTTa                                                           <  1:1565259/100‑1 (MQ=255)
                                        aTTAACGTCGCACCATGCGCATTTTTACCGGTTATCGAATGTAAAAAGCCCCGCAACGTGTTGCGGGGCTTT‑‑‑CATCCGTTACCGGGACGCGAAAAACTTa                                                           >  1:4059912/1‑100 (MQ=255)
                                         ttAACGTCGCACCATGCGCATTTTTACCGGTTATCGAATGTAAAAAGCCCCGCAACGTGTTGCGGGGCTTT‑‑‑CATCCGTTACCGGGACGCGAAAAACTTAt                                                          <  1:2438248/100‑1 (MQ=255)
                                           aaCGTCGCACCATGCGGATTTTTACCGGTGATCGAATGTAAAAAGCCCCGCAACGTGTTGCGGGGCTTT‑‑‑CATCCGTTAAAAGGAAAAGAAAAACTTAtta                                                        <  2:2813225/100‑2 (MQ=255)
                                           aaCGTCGCACCATGCGCATTTTTACCGGTTATCGAATGTAAAAAGCCCCGCAACGTGTTGCGGGGCTTT‑‑‑CATCCGTTACCGGGACGCGAAAAACTTATTc                                                        <  1:1793886/100‑1 (MQ=255)
                                           aaCGTCGCACCATGCGCATTTTTACCGGTTATCGAATGTAAAAAGCCCCGCAACGTGTTGCGGGGCTTT‑‑‑CATCCGTTACCGGGACGCGAAAAACTTATTc                                                        >  1:2988174/1‑100 (MQ=255)
                                            aCGTCGCACCATGCGCATTTTTACCGGTTATCGAATGTAAAAAGCCCCGCAACGTGTTGCGGGGCTTT‑‑‑CATCCGTTACCGGGACGCGAAAAACTTATTCa                                                       >  1:640206/1‑100 (MQ=255)
                                            aCGTCGCACCATGCGCATTTTTACCGGTTATCGAATGTAAAAAGCCCCGCAACGTGTTGCGGGGCTTT‑‑‑CATCCGTTACCGGGACGCGAAAAACTTATTCa                                                       >  1:642273/1‑100 (MQ=255)
                                             cGTCGCACCATGCGCATTTTTACCGGTTATCGAATGTAAAAAGCCCCGCAACGTGTTGCGGGGCTTT‑‑‑CATCCGTTACCGGGACGCGAAAAACTTATTCAg                                                      <  2:551323/100‑1 (MQ=255)
                                                 gCACCATGCGCATTTTTACCGGTTATCGAATGTAAAAAGCCCCGCAACGTGTTGCGGGGCTTT‑‑‑CATCCGTTACCGGGACGCGAAAAACTTATTCAGCTTc                                                  >  1:3699519/1‑100 (MQ=255)
                                                 gCACCATGCGCATTTTTACCGGTTATCGAATGTAAAAAGCCCCGCAACGTGTTGCGGGGCTTT‑‑‑CATCCGTTACCGGGACGCGAAAAACTTATTCAGCTTc                                                  >  1:3411483/1‑100 (MQ=255)
                                                  cACCATGCGCATTTTTACCGGTTATCGAATGTAAAAAGCCCCGCAACGTGTTGCGGGGCTTT‑‑‑CATCCGTTACCGGGACGCGAAAAACTTATTCAGCTTCg                                                 <  2:3917375/100‑1 (MQ=255)
                                                   aCCATGCGCATTTTTACCGGTTATCGAATGTAAAAAGCCCCGCAACGTGTTGCGGGGCTTT‑‑‑CATCCGTTACCGGGACGCGAAAAACTTATTCAGCTTCGa                                                <  1:183811/100‑1 (MQ=255)
                                                       tGCGCATTTTTACCGGTTATCGAATGTAAAAAGCCCCGCAACGTGTTGCGGGGCTTT‑‑‑CATCCGTTACCGGGACGCGAAAAACTTATTCAGCTTCGATGcc                                            <  1:2178137/100‑2 (MQ=255)
                                                       tGCGCATTTTTACCGGTTATCGAATGTAAAAAGCCCCGCAACGTGTTGCGGGGCTTT‑‑‑CATCCGTTACCGGGACGCGAAAAACTTATTCAGCTTCGATGCt                                            >  1:736459/1‑100 (MQ=255)
                                                        gcgcATTTTTACCGGTTATCGAATGTAAAAAGCCCCGCAACGTGTTGCGGTGCTTT‑‑‑CATCCGTTACCGGGACACGAAAAACCTATTCAGCTTCGATGCt                                            >  2:3554904/1‑99 (MQ=255)
                                                         cgcATTTTTACCGGTTATCGAATGTAAAAAGCCCCGCAACGTGTTGCGGGGCTTT‑‑‑CATCCGTTACCGGGACGCGAAAAACTTATTCAGCTTCGATGCTGa                                          <  2:347583/100‑1 (MQ=255)
                                                         cgcATTTTTACCGGTTATCGAATGTAAAAAGCCCCGCAACGTGTTGCGGGGCTTT‑‑‑CATCCATTAACGGGACGCGAAAAACTTATTCAGCTTCGCAGCTGa                                          >  2:2953666/1‑100 (MQ=255)
                                                          gcATTTTTACCGGTTATCGAATGTAAAAAGCCCCGCAACGTGTTGCGGGGCTTT‑‑‑CATCCGTTACCGGGACGCGAAAAACTTATTCAGCTTCGATGCTGAt                                         <  1:138262/100‑1 (MQ=255)
                                                           cATTTTTACCGGTTATCGAATGTAAAAAGCCCCGCAACGTGTTGCGGGGCTTT‑‑‑CATCCGTTACCGGGACGCGAAAAACTTATTCAGCTTCGATGCTGATa                                        >  1:2713959/1‑100 (MQ=255)
                                                             tttttACCGGTTATCGAATGTAAAAAGCCCCGCAACGTGTTGCGGGGCTTT‑‑‑CATCCGTTACCGGGACGCGAAAAACTTATTCAGCTTCGATGCTGATaaa                                      >  2:83762/1‑100 (MQ=255)
                                                               tttACCGGTTATCGAATGTAAAAAGCCCCGCAACGTGTTGCGGGGCTTT‑‑‑CATCCGTTACCGGGACGCGAAAAACTTATTCAGCTTCGATGCTGATAAAtt                                    <  1:2851734/100‑1 (MQ=255)
                                                               tttACCGGTTATCGAATGTAAAAAGCCCCGCAACGTGTTGCGGGGCTTT‑‑‑CATCCGTTACCGGAACGCGAAAAACTTATTCAGCTTCGATGCTGATAAAtt                                    <  2:2634829/100‑1 (MQ=255)
                                                                  aCCGGTTATCGAATGTAAAAAGCCCCGCAACGTGTTGCGGGGCTTT‑‑‑CATCCGTTACCGGGACGCGAAAAACTTATTCAGCTTCGATGCTGATAAATTTAc                                 <  1:396256/100‑1 (MQ=255)
                                                                   ccGGTTATCGAATGTAAAAAGCCCCGCAACGTGTTGCGGGGCTTT‑‑‑CATCCGTTACCGGGACGCGAAAAACTTATTCAGCTTCGATGCTGATAAATTTACg                                <  2:4014234/100‑1 (MQ=255)
                                                                   ccGGTTATCGAATGTAAAAAGCCCCGCAACGTGTTGCGGGGCTTT‑‑‑CATCCGTTACCGGGACGCGAAAAACTTATTCAGCTTCGATGCTGATAAATTTACg                                >  1:27500/1‑100 (MQ=255)
                                                                      gTTATCGAATGTAAAAAGCCCCGCAACGTGTTGCGGGGCTTT‑‑‑CATCCGTTACCGGGACGCGAAAAACTTATTCAGCTTCGATGCTGATAAATTTACGGtt                             <  2:3353153/100‑1 (MQ=255)
                                                                        tagcGAATGTAAAAAGCCCCGCAACGTGTCGCGGGGCTTT‑‑‑CATCGGTTGCCGGGACGTGACAACTTTATCCAGCTGCTCTTCTGATTAATGTACggat                             >  2:3586797/4‑96 (MQ=255)
                                                                         aTCGAATGTAAAAAGCCCCGCAACGTGTTGCGGGGCTTT‑‑‑CATCCGTTACCGGGACGCGAAAAACTTATTCAGCTTCGATGCTGATAAATTTACGGttttt                          >  1:811758/1‑100 (MQ=255)
                                                                          tCGAATGTAAAAAGCCCCGCAACGTGTTGCGGGGCTTT‑‑‑CATCCGTTACCGGGACGCGAAAAACTTATTCAGCTTCGATGCTGATAAATTTACGGTTTTTc                         >  1:1065048/1‑100 (MQ=255)
                                                                           cGAATGTAAAAAGCCCCGCAACGTGTTGCGGGGCTTT‑‑‑CATCCGTTACCGGGACGCGAAAAACTTATTCAGCTTCGATGCTGATAAATTTACGGTTTTTCg                        >  2:1584175/1‑100 (MQ=255)
                                                                           cGAATGTAAAAAGCCCCGCAACGTGTTGCGGGGCTTT‑‑‑CATCCGTTACCGGGACGCGAAAAACTTATTCAGCTTCGATGCCGATAAATTTACGGTTTTTCg                        >  1:3505780/1‑100 (MQ=255)
                                                                             aaTGTAAAAAGCCCCGCAACGTGTTGCGGGGCTTT‑‑‑CATCCGTTACCGGGACGCGAAAAACTTATTCAGCTTCGATGCTGATAAATTTACGGTTTTTCggg                      <  1:967484/100‑1 (MQ=255)
                                                                               tGTAAAAAGCCCCGCAACGTGTTGCGGGGCTTT‑‑‑CATCCGTTACCGGGACGCGAAAAACTTATTCAGCTTCGATGCTGATAAATTTACGGTTTTTCGGGcc                    <  1:1134544/100‑1 (MQ=255)
                                                                               tGTAAAAAGCCCCGCAACGTGTTGCGGGGCTTT‑‑‑CATCCGTTACCGGGACGCGAAAAACTTATTCAGCTTCGATGCTGATAAATTTACGGTTTTTCGGGcc                    <  1:2029901/100‑1 (MQ=255)
                                                                                 tAAAAAGCCCCGCAACGTGTTGCGGGGCTTT‑‑‑CATCCGTTACCGGGACGCGAAAAACTTATTCAGCTTCGATGCTGATAAATTTACGGTTTTTCGGGCCtt                  >  2:3378282/1‑100 (MQ=255)
                                                                                 tAAAAAGCCCCGCAACGTGTTGCGGGGCTTT‑‑‑CATCCGTTACCGGGACGCGAAAAACTTATTCAGCTTCGATGCTGATAAATTTACGGTTTTTCGGGCCtt                  >  2:1627961/1‑100 (MQ=255)
                                                                                 tAAAAAGCCCCGCAACGTGTTGCGGGGCTTT‑‑‑CATCCGTTACCGGGACGCGAAAAACTTATTCAGCTTCGATGCTGATAAATTTACGGTTTTTCGGGCCtt                  <  2:1969936/100‑1 (MQ=255)
                                                                                 tAAAAAGCCCCGCAACGTGTTGCGGGGCTTT‑‑‑CATCCGTTACCGGGACGCGAAAAACTTATTCAGCTTCGATGCTGATAAATTTACGGTTTTTCGGGCCtt                  >  1:514714/1‑100 (MQ=255)
                                                                                  aaaaaGCCCCGCAACGTGTTGCGGGGCTTT‑‑‑CATCCGTTACCGGGACGCGAAAAACTTATTCAGCTTCGATGCTGATAAATTTACGGTTTTTCGGGCCttt                 >  2:3687885/1‑100 (MQ=255)
                                                                                  aaaaaGCCCCGCAACGTGTTGCGGGGCTTT‑‑‑CATCCGTTACCGGGACGCGAAAAACTTATTCAGCTTCGATGCTGATAAATTTACGGTTTTTCGGGCCttt                 <  1:1767314/100‑1 (MQ=255)
                                                                                  aaaaaGCCCCGCAACGTGTTGCGGGGCTTT‑‑‑CATCCGTTACCGGGACGCGAAAAACTTATTCAGCTTCGATGCTGATAAATTTACGGTTTTTCGGGCCttt                 >  2:1353301/1‑100 (MQ=255)
                                                                                  aaaaaGCCCCGCAACGTGTTGCGGGGCTTT‑‑‑CATCCGTTACCGGGACGCGAAAAACTTATTCAGCTTCGATGCTGATAAATTTACGGTTTTTCGGGCCtt                  >  1:2610050/1‑99 (MQ=255)
                                                                                  aaaaaGCCCCGCAACGTGTTGCGGGGCTTT‑‑‑CATCCGTCACCGGGACGCGAAAAACTTATTCAGCTTCGATGCTGATAAATTTACGGTTTTTCGTGCCttt                 >  2:2249937/1‑100 (MQ=255)
                                                                                   gaaaGCCCCGCAACACGTTGCGGGGCTTT‑‑‑CATCCGTTACCGGGACGCGAAAAACTTATTCAGCTTCGATGCTGATAAATTTACGGTTTTTCGGGCCTTTa                <  2:2665128/99‑1 (MQ=255)
                                                                                    aaaGCCCCGCAACGTGTTGCGGGGCTTT‑‑‑CATCCGTTACCGGGACGCGAAAAACTTATTCAGCTTCGATGCTGATAAATTTACGGTTTTTCGGGCCTTTaa               >  1:2863526/1‑100 (MQ=255)
                                                                                     aaGCCCCGCAACGTGTTGCGGGGCTTT‑‑‑CATCCGTTACCGGGACGCGAAAAACTTATTCAGCTTCGATGCTGATAAATTTACGGTTTTTCGGGCCTTTAAc              >  1:1715340/1‑100 (MQ=255)
                                                                                      aGCCCCGCAACGTGTTGCGGGGCTTT‑‑‑CATCCGTTACCGGGACGCGAAAAACTTATTCAGCTTCGATGCTGATAAATTTACGGTTTTTCGGGCCTTTAAc              <  1:3550711/99‑1 (MQ=255)
                                                                                       gCCCTGCAACACGTTGCGGGGCTTT‑‑‑CATCCGTTACCGGGACGCGAAAAACTTATTCAGCTTCGATGCTGATAAATTTACGGTTTTTCGGGCCTTTAACtt            >  2:4086379/1‑100 (MQ=255)
                                                                                       gCCCCGCAACACGTTGCGGGGCTTTTTACATTCGATAACCGGGACGCGAAAAACTTATTCAGCTTCGATGCTGATAAATTTACGGTTTTTCGGGCCTTTa                >  1:2111110/1‑100 (MQ=18)
                                                                                        ccccGCAACGTGTTGCGGGGCTTT‑‑‑CATCCGTTACCGGGACGCGAAAAACTTATTCAGCTTCGATGCTGATAAATTTACGGTTTTTCGGGCCTTTAACttn           <  2:3072771/100‑2 (MQ=255)
                                                                                        ccccGCAACGTGTTGCGGGGCTTT‑‑‑CATCCGTTACCGGGACGCGAAAAACTTATTCAGCTTCGATGCTGATAAATTTACGGTTTTTCGGGCCTTTAACTTc           >  2:788575/1‑100 (MQ=255)
                                                                                        ccccGCAACGTGTTGCGGGGCTTT‑‑‑CATCCGTTACCGGGACGCGAAAAACTTATTCAGCTTCGATGCTGATAAATTTACGGTTTTTCGGGCCTTTAACTTc           >  1:1005697/1‑100 (MQ=255)
                                                                                        ccccGCAACACGTTGCGGGGCTTT‑‑‑CATCCGTTACCGGGACGCGAAAAACTTATTCAGCTTCGATGCTGATAAATTTACGGTTTTTCGGGCCTTTAACTTc           >  2:2375222/1‑100 (MQ=255)
                                                                                        ccccGCAACACGTTGCGGGGCTTT‑‑‑CATCCGTTACCGGGACGCGAAAAACTTATTCAGCTTCGATGCTGATAAATTTACGGTTTTTCGGGCCTTTAACTTc           >  1:3634328/1‑100 (MQ=255)
                                                                                         cccGCAACGTGTTGCGGGGCTTT‑‑‑CATCCGTTACCGGGACGCGAAAAACTTATTCAGCTTCGATGCTGATAAATTTACGGTTTTTCGGGCCTTTAACTTCg          >  1:168147/1‑100 (MQ=255)
                                                                                         cccGCAACGTGTTGCGGGGCTTT‑‑‑CATCCGTTACCGGGACGCGAAAAACTTATTCAGCTTCGATGCTGATAAATTTACGGTTTTTCGGGCCTTTAACTTCg          >  2:2032063/1‑100 (MQ=255)
                                                                                                 gtgtTGCGGGGCTTT‑‑‑CATCCGTTACCGGGACGCGAAAAACTTATTCAGCTTCGATGCTGATAAATTTACGGTTTTTCGGGCCTTTAACTTCGAATTTCAc  <  1:1448972/100‑1 (MQ=255)
                                                                                                 |                                                                                                       
TGCTGCTTCATACCGTATTCCCTGCAAAGCGCATTTGCGGATTAACGTCGCACCATGCGCATTTTTACCGGTTATCGAATGTAAAAAGCCCCGCAACGTGTTGCGGGGCTTT‑‑‑CATCCGTTACCGGGACGCGAAAAACTTATTCAGCTTCGATGCTGATAAATTTACGGTTTTTCGGGCCTTTAACTTCGAATTTCAC  >  NC_000913/3332725‑3332921

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