Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,026,9410AC36.0% 98.1 / 11.5 50E392D (GAA→GAC) guaDguanine deaminase
Reads supporting (aligned to +/- strand):  ref base A (27/5);  new base C (0/18);  total (27/23)
Fisher's exact test for biased strand distribution p-value = 1.86e-09
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
Rejected as consensus: Frequency below/above cutoff threshold.

TTAAAGTGGGCATGGGAACGGATATCGGTGCCGGAACCACTTTCAACATGCTGCAAACGCTGAACGAAGCCTACAAAGTATTGCAATTACAAGGCTATCGCCTCTCGGCATATGAAGCGTTTTACCTGGCCACGCTCGGCGGAGCGAAATCTCTGGGCCTTGACGATTTGATTGGCAACTTTTT‑ACCTGGCAAA‑GAGGC‑TGATTT‑CGTGGTGATGGAACCCACCGCCACTCCGCTACAGCAGCTGCGCTATGACAACTCTGTTTCTTTAGTCGACAAATTGTTCGTGATGATGACGTTGGGCGATGACCGTTCGATCTACCGCACCTACGTTGATGGTCGTCTGGTGTACGAACGCAACTAATAATAAAACTTTAACATCCTCGTGAGGACATCATTATGTCTGGAGACATCCTACAAACACCGGACGCACCAAAGC  >  NC_000913/3026724‑3027160
                                                                                                                                                                                                                             |                                                                                                                                                                                                                           
taaaaGTGGGAATGGGAACGGATATGGGTGCCGAAACCACTTACAACATGCTGCAAACGCTGAACGAAGCATAAAAATTATTGCAATTACAAGGCTACCCCCTCTCGGAATATGAAGCTTTTTACCTGCCCACGCTCGGCGGAGCGAAATCTCTGGGCCTTGACGTTTTGATGGGCAATTTTTT‑CCCTGGAAAA‑GAGGC‑TGTTTT‑CGTGGTGAGGAAACCCACCGCCACTCCGCTAc                                                                                                                                                                                                          <  1:598185/235‑1 (MQ=255)
      gggtcaggagAAGGGATATCGGTGCCGAAACCCATTTCAACATGCTCCAAACGCTAAACGAGGCCTAAAAATTATTGCAATTACAGGGCTTCCGCCTCTCGGCATATGAAGCGTTTTACCTGGCCACGCTGGGCGGAGCGAAATCTCTGGGCTTTGACGATTTGATTGGCAATTTTTT‑CCCTGGAAAA‑GAGGC‑TGATTT‑CGTGGTGAGGGACCCCACCGCCACTCCGCTACagcagc                                                                                                                                                                                                    <  1:775929/228‑1 (MQ=255)
             gggAACGGATATCGGTGCCGGAACCACTTTCAACATGCTGCAAACGCTGAACGAAGCCTACAAAGTATTGCAATTACAAGGCTATCGCCTCTCGGCATATGAAGCGTTTTACCTGGCCACGCTCGGCGGAGCGAAATCTCTGGGCCTTGACGATTTGATTGGCAACTTTTT‑ACCTGGCAAA‑GAGGC‑TGATTT‑CGTGGTGATGGAACCCACCGCCACTCCGCTACAGCAGCTGCGCt                                                                                                                                                                                              >  1:1131767/1‑236 (MQ=255)
                   ggATATCGGTGCCGGAACCACTTTCAACATGCTGCAAACGCTGAACGCAGCCTACAAAGTATTGCAATTACAAGGCTATCGCCTCTCGGCATATGAAGCGTTTTACCTGGCCACGCTCGGCGGAGCGAAATCTCTGGGCCTAGACGATTTGATCGGCAACTTTTT‑ACCTGCCACA‑GAGGC‑TGATTT‑CGTGGTGATGGAACCCACCGCCACTCCTCTACAGCAGCTGCGCTATGACAAc                                                                                                                                                                                      >  2:534909/1‑238 (MQ=255)
                   gaattccGGTCCCGAAACCATTTCAAACATGCTGAAAACGTTAAACGAAGCCAAAAAAGTATTGCAATTACAGGGCTATCGCCTCTGGGAATAGGAGGTGTTTTACCGGGCCACGCTCGGGGGAGCGAAATCTCGGGGCCTTGACGTTTTGATTGGCATCTTTTT‑CCCTGGAAAA‑GGGGC‑TGTTTT‑GGGGGTGATGGACCCCACCGCCACTCCGCTACAGCAGCTGCGCTAt                                                                                                                                                                                            <  2:961409/226‑1 (MQ=255)
                    gATATCGGGGCCGGACCCCTTTTCAAACTGCTCCAATCGCTGAACAAAGCCAAAAAATTATTGCAATTAAAAGCATCTCGCCTCGGGGCTTTTGAAGCTTTTTACCGGGCCACGCTCGGCGGAGCGAATCTTCTGGGCCTTGGCGATTTGATTGCCACCTTTTT‑CCCTGGCAAA‑GGGGG‑TGATTG‑CGTGGTGATGCACCCCACCGCCACTCCGCTACAGCAGCTGCGCTATGACAAc                                                                                                                                                                                      <  2:501492/237‑1 (MQ=255)
                           gTGCCGGAACCACTTTCAACATGCTGCAAACGCTGAACGAAGCCTACAAAGTATTGCAATTACAAGGCTATCGCCTCTCGGCATATGAAGCGTTTTACCTGGCCACGCTCGGCGGAGCGAAATCTCTGGGCCTTGACGATTTGATTGGCAACTTTTT‑ACCTGGCAAA‑GAGGC‑TGATTT‑CGTGGTGATGGAACCCACCGCCACTCCGCTACAGCAGCTGCGCTATGACAAc                                                                                                                                                                                      >  1:863772/1‑230 (MQ=255)
                             gCCGGAACCACTTTCAACATGCTGCAAACGCTGAACGAAGCCTACAAAGTATTGCAATTACAAGGCTATCGCCTCTCGGCATATGAAGCGTTTTACCTGGCCACGCTCGGCGGAGCGAAATCTCTGGGCCTTGACGATTTGATTGGCAACTTTTT‑ACCTGGCAAA‑GAGGC‑TGATTT‑CGTGGTGATGGAACCCACCGCCACTCCGCTACAGCAGCTGCGCTATGACAACTCTg                                                                                                                                                                                  >  1:961409/1‑232 (MQ=255)
                                                tgctgcAAACGCTGAACGAAGCCTACAAAGTATTGCAATTACAAGGCTATCGCCTCTCGGCATATGAAGCGTTTTACCTGGCCACGCTCGGCGGAGCGAAATCTCTGGGCCTTGACGATTTGATTGGCAACTTTTT‑ACCTGGCAAC‑GAGGC‑TGATTT‑CGTGGTGATGGAACCCACCGCCACTCCGCTACAGCAGCTGCGCTATGACAACTCTGTTTCTTTAGTCGACAAATTGTTCGtg                                                                                                                                                        >  2:969275/1‑239 (MQ=255)
                                                          gataaaCGAAGCCAAAAAATTATTCCATTTACAAGGCTATCGCCTCTCGGCATAGGAAGGTTTTTCCCTGGCCCCCCTCGGCGGAGCAAAAGCTCGGGGCCTGGACGATTTGATTGGCAACTTTTT‑CCCGGGAAAA‑GAGGC‑TGATTT‑CGTGGTGATGGACCCCACCGCCACTCCGCTACAGCAGCTGCGCTATGACAACTCTGTTTCTTTAGTCGACAAATTGTTCGTGATGATGACGt                                                                                                                                              <  1:430020/235‑1 (MQ=255)
                                                             gAACGAAGCCTACAAAGTATTGCAATTACAAGGCTATCGCCTCTCGGCATATGAAGCGTTTTACCTGGCCACGCTCGGCGGAGCGAAATCTCTGGGCCTTGACGATTTGATTGGCAACTTTTT‑ACCTGGCAAA‑GAGGC‑TGATTT‑CGTGGTGATGGAACCCACCGCCACTCCGCTACAGCAGCTGCGCTATGACAACTCTGTTTCTTTAGTCGACAAATTGTTCGTGATGATGACGTTgg                                                                                                                                           >  1:823406/1‑239 (MQ=255)
                                                                        aCAAAGTATTGCAATTACAAGGCTATCGCCTCTCGGCATATGAAGCGTTTTACCTGGCCACGCTCGGCGGAGCGAAATCTCTGGGCCTTGACGATTTGATTGGCAACTTTTT‑ACCTGGCAAA‑GAGGC‑TGATTT‑CGTGGTGATGGAACCCACCGCCACTCCGCTACAGCAGCTGCGCTATGACAACTCTGTTTCTTTAGTCGACAAATTGTTCGTGATGATGACGTTGGGCGATGACCGt                                                                                                                                >  1:55115/1‑239 (MQ=255)
                                                                        aCAAAGTATTGCAATTACAAGGCTATCGCCTCTCGGCATATGAAGCGTTTTACCTGGCCACGCTCGGCGGAGCGAAATCTCTGGGCCTTGACGATTTGATTGGCAACTTTTT‑ACCTGGCAAA‑GAGGC‑TGATTT‑CGTGGTGATGGAACCCACCGCCACTCCGCTACAGCAGCTGCGCTATGACAACTCTGTTTCTTTAGTCGACAAATTGTTCGTGATGATGACGTTGGGCGATGACCGt                                                                                                                                >  1:903291/1‑239 (MQ=255)
                                                                         cAAATTTTTGCAATTACAAGGCTGTCGCCTCTCGGCATATGAATCTTTCTCCCTGCCCACCCTCGGCGGAGCGAAATCTCTGGCCTTTGCCGTTTTGTTTGGCACCTTTTT‑CCCTGCAAAA‑GAGCT‑TTTTTT‑CGTGGTGATGGACCCCACCGCCACTCCGCTACagcagc                                                                                                                                                                                                    <  2:682309/170‑1 (MQ=255)
                                                                             gTATTGCAATTACAAGGCTATCGCCTCTCGGCATATGAAGCTTTTTACCGGGCCACGCTCGGCGGAGCGAAATCTCGGGGCTTTGACGATTTGATTGGCATTTTTTT‑ACCTGGCAAA‑GAGGC‑TTATTT‑GGTGGGGATGAACCCCACCGCCACTCCGCTACAGCAGCTGCGCTATGACAACTCTGTTTCTTTAGTCGACAAATTGTTCGTGATGATGACGTTGGGCGATGACCGtt                                                                                                                               <  1:963182/235‑1 (MQ=255)
                                                                              ttttGCAATAAAAAGGCTATCGCCTCTCGGCATATGAAGTTTTTTCCGTGGCCCCCCTCGGCGGAGCGAAATCTCGGGCCCTTGACGTTTTGTTGGGCACTTTTTT‑CCCTGGAAAA‑GGGGC‑TGTTTT‑GGTGGTGAGGAACCCCCCCGCCACTCCGCTACAGCAGCTGCGCTATGAc                                                                                                                                                                                         <  2:50834/174‑1 (MQ=255)
                                                                                 tGCAATTACAAGGCTATCGCCTCTCGGCATATGAAGCGTTTTACCTGGCCACGCTCGGCGGAGCGAAATCTCTGGGCCTTGACGATTTGATTGGCAACTTTTT‑ACCTGGCAAA‑GAGGC‑TGATTT‑CGTGGTGATGGAACCCACCGCCACTCCGCTa                                                                                                                                                                                                           >  1:1201497/1‑155 (MQ=255)
                                                                                   cAATTACAAGGCTATCGCCTCTCGGCATATGAAGCGTTCTACCTGGCCACGCTCGGCGGAGCGAAATCTCTGGGCCTTGACGATTTGATTGGCAACTTTTT‑ACCTGGCAAA‑GAGGC‑TGATTT‑CGTGGTGATGGAACCCACCGCCACTCCGCTACAGCAGCTGCGCTATGa                                                                                                                                                                                          >  1:682309/1‑170 (MQ=255)
                                                                                      ttACAAGGCTATCGCCTCTCGGCATATGAAGCGTTTTACCTGGCCACGCTCGGCGGAGCGAAATCTCTGGGCCTTGACGATTTGATTGGCAACTTTTT‑ACCTGGCAAA‑GAGGC‑TGATTT‑CGTGGTGATGGAACCCACCGCCACTCCGCTACAGCAGCTGCGCTATGACAACTCTGTTTCTTTAGTCGACAAATTGTTCGTGATGATGACGTTGGGCGATGACCGTTCGATCTACCGCa                                                                                                                   >  2:569895/1‑238 (MQ=255)
                                                                                        aCAAGGCTATCGCCTCTCGGCATATGAAGCGTTTTACCTGGCCACGCTCGGCGGAGCGAAATCTCTGGGCCTTGACGATTTGATTGGCAACTTTTT‑ACCTGGCAAA‑GAGGC‑TGATTT‑CGTGGTGATGGAACCCACCGCCACTCCGCTACAGCAGCTGCGCTATGACAACTCTGttt                                                                                                                                                                               >  1:50834/1‑176 (MQ=255)
                                                                                        aCAAGGCTATCGCCTCTCGGCATATGAAGCGTTTTACCTGGCCACGCTCGGCGGAGCGAAATCTCTGGGCCTTGACGATTTGATTGGCAACTTTTT‑ACCTGGCAAA‑GAGGC‑TGATTT‑CGTGGTGATGGAACCCACCGCCACTCCGCTACAGCAGCTGCGCTATGACAACTCTGTTTCTTTAGTCGACAAATTGTTCGTGATGATGACGTTGGGCGATGACCGTTCGATCTACCGCACCt                                                                                                                >  1:168297/1‑239 (MQ=255)
                                                                                               tATCGCCTCTCGGCATATGAAGCTTTTTACCTGGCCACGCTCGGCGGAGCGAATTCTCTGGGCCTTGACGATTTGATTGGCAACTTTTT‑ACCTGGCAAA‑GAGGC‑TGATTT‑CGTGGTGATGGAACCCACCGCCACTCCGCTACAGCAGCTGCGCTATGACAACTCTGTTTCTTTAGTCGACAAATTGTTCGTGATGATGACGTTGGGCGATGACCGTTCGATCTACCGCACCTACGTTGa                                                                                                         <  1:569895/239‑1 (MQ=255)
                                                                                                 tCGCCTCTCGGCATATGAAGCGCTTTACCTGGCCACGCTCGGCGGAGCGAAATCTCTTGGCCTTGACGATTTGATTGGCAACTTTTT‑ACCTGGCAAA‑GAGGC‑TGATTT‑CGTGGTGATGGAACCCACCGCCACTCCGTTa                                                                                                                                                                                                           >  1:1073498/1‑139 (MQ=255)
                                                                                                     ctctCGGCATATGAAGCGTTTTCCCTGGCCCCGCTGGGGGGAGCGAAATCTCTGGGCTTTGCCGATTTGATTGGCAACTTTTT‑CCCTGGCAAA‑GAGGC‑TGTTTT‑GGTGGTGAGGGAACCCACCGCCACTCCGCTACAGCAGCTgcgc                                                                                                                                                                                               <  1:905584/147‑1 (MQ=255)
                                                                                                        tcGGCATATGAAGCGTTTTACCTGGCCACGCTCGGCGGAGCGAAATCTCTGGGCCTTGACGATTTGATTGGCAACTTTTT‑ACCTGGCAAA‑GAGGC‑TGATTT‑CGTGGTGATGGAACCCACCGCCACTCCGCTACAGCAGCTGCGCTATGACAACTCTGTTTCTTTAGTCGACAAATTGTTCGTGATGATGACGTTGGGCGATGACCGTTCGATCTACCGCACCTACGTTGATGgtcg                                                                                                   >  1:810809/1‑236 (MQ=255)
                                                                                                        tcGGCATATGAAGCGTTTTACCTGGCCACGCTCGGCGGAGCGAAATCTCTGGGCCTTGACGATCTGATTGGCAACTTTTT‑ACCTGGCAAA‑GAGGC‑TGATTT‑CGTGGTGATGGAACCCACCGCCACTCCGCTACAGCAGCTGCGCTATGACAAc                                                                                                                                                                                      >  1:679064/1‑153 (MQ=255)
                                                                                                               atGAAGCGTTTTACCTGGCCACGCTCGGCGGAGCGAAATCTCTGGGCCTTGACGATTTGATTGGCAACTTTTT‑ACCTGGCAAA‑GAGGC‑TGATTT‑CGTGGTGATGGAACCCACCGCCACTCCGCTACAGCAGCTGCGCTATGACAAct                                                                                                                                                                                     >  2:905584/1‑147 (MQ=255)
                                                                                                               atGAAGCGTTTTACCTGGCCACGCTCGGCGGAGCGAAATCTCTGGGCCTTGACGATTTGATTGGCAACTTTGT‑ACCTGGCAAA‑GAGGC‑TGATTT‑CGTGGTGATGGAACCCACCGCCACTCCGCTACAGCAGCt                                                                                                                                                                                                   >  1:1312059/1‑133 (MQ=255)
                                                                                                                       ttttACCTGGCCCCGCGCGGCGGAGCGAAATCTCTGGGCCTTGGTTATTTGATTGGCAACTTTTA‑ACCTGAAAAA‑GAGCC‑TTTTTT‑GGTGATGATGGACCCAGCCACCACTCCGCTACAGCAGCTGCGCTATGACAACTCTGTTTCTTTAGTCGACAAATTGTTCGTGATGATGACGTTGGGCGATGACCGTTCGATCTAc                                                                                                                       <  2:825450/201‑1 (MQ=255)
                                                                                                                                gCCACGCTCGGCGGAGAAAAACCTCTGGGCCTGGACGATTTGATGGGCAACTTTTT‑ACCGGGAAAA‑GAGGC‑GGTTTT‑CGGGGTGAGGAACCCCACCGCCACTCCGCTACAGCAGCTGCGCTATGACAACTCTGTTTCTTTAg                                                                                                                                                                         <  2:661704/142‑1 (MQ=255)
                                                                                                                                 ccACGCTCGGCGGAGCGAAATCTCTGGGCCTTGACGATTTGATTGGCAACTTTTT‑ACCTGGCAAA‑GAGGC‑TGATTT‑CGTGGTGATGGAACCCACCGCCACTCCGCTACAGCAGCTGCGCTATGACAACTCTGTTTCTTTAGTCGACAAATTGTTCGTGATGATGACGTTGGGCGATGACCGTTCGATCTACCGCACCTACg                                                                                                             >  1:825450/1‑201 (MQ=255)
                                                                                                                                          gcggAGCGAAATCTCTGGGCCTTGACGATTTGATTGGCAACTTTTT‑ACCTGGCAAA‑GAGGC‑TGATTT‑CGTGGTGATGGAACCCACCGCCACTCCGCTACAGCAGCTGCGCTATGACAACTCTGTTTCTTTAGTCGACAAAtt                                                                                                                                                               >  1:661704/1‑142 (MQ=255)
                                                                                                                                            ggggcaaattctcGGGGCCTGGGCGTTTTGA‑TGGCAACTTTTT‑CGCTGGGAAAGGGGGC‑TGTTTT‑GGTGGTGATGGACCCCCCCGCCACTCCGCTACAGCAGCTGCGCTATGACAACTCTGTTTCTTTAGTCGACAAATTGTTCGTGATGATGACGTTGGGCGATGACCGTTCGATCTACCGCACCTACGTTGATGGTCGTCTGGTGTACGAACGCAACTAATAATAAAACTTTAACa                                                             <  2:55115/229‑1 (MQ=255)
                                                                                                                                              aGCGAAATCGTGGCGCCTTGACGTTTTGGTT‑GCCATTTTTT‑TCCTGGGCAA‑GGAGC‑TGTTTTGCTGGGAGATGAACCCCACCGCCACTCCGCTACAGCAGCTGCGCTATGACAACTCTGTTTCTTTAGTCGACAAATTGTTCGTGATGATGACGTTGGGCGATGACCGTTCGATCTACCGCACCTACGTTGATGGTCGTCTGGTGTACGAACGCAACTAATAATAAAACTTTAACATc                                                           <  2:772443/238‑1 (MQ=255)
                                                                                                                                               gCGAAATCTCTGGGCCTTGACGATTTGATTGGCAACTTTTT‑ACCTGGCAAA‑GAGGC‑TGATTT‑CGTGGTGATTGAACCCACCGCCACTCCGCTACAGCAGCTGCGCTATGACAACTCTGTTTCTTTAGTCGACAAATTGTTCGTGATGATGACGTTGGGCGATGACCGTTCGATCTACCGCACCTACGTTGATGGTCGTCTGGTGTACGAACGCAACTAATAATAAAACTTTAACATcc                                                          >  1:1354519/1‑238 (MQ=255)
                                                                                                                                                      cttggggCGTTGACGTTTTGATTGGCATCTTTTT‑ACCGGGAAAA‑GGGGT‑TGATTT‑CGTGGTGAGGAACCCCACCGCCACTCCGCTACAGCAGCTGCGCTATGACAACTCTGTTTCTTTAGTCGACAAATTGTTCGTGATGATGACGTTGGGCGATGACCGTTCGATCTACCGCACCTACGTTGATGGTCGTCTGGTGTACGAACGCAACTAATAATAAAACTTTAACATCCTCGTGAgg                                                  <  1:429186/235‑1 (MQ=255)
                                                                                                                                                        ctGGGCTTTGACGTTTGGATTGG‑AACTTTTTCACCTGGAAAG‑GAGGC‑TGTTTT‑CGTGGTGATGGACCCCACCGCCACTCCGCTACAGCAGCTGCGCTATGACAACTCTGTTTCTTTAGTCGACAAATTGTTCGTGATGATGACGTTGGGCGATGACCGTTCGATCTACCGCACCTACGTTGATGGTCGTCTGGTGTACGAACGCAACTAATAATAAAACTTTAACATCCTCGTGAGGAc                                                <  2:1096352/239‑1 (MQ=255)
                                                                                                                                                             ccTGGACTATTTGATTGGCATCTTTTT‑CCCT‑GAAAA‑GGGGCTTGTTTT‑CGGGGTGATGGACCCCCCCGCCACTCCGCTACAGCAGCTGCGCTATGACAACTCTGTTTCTTTAGTCGACAAATTGTTCGTGATGATGACGTTGGGCGATGACCGTTCGATCTACCGCACCTACGTTGATGGTCGTCTGGTGTACGAACGCAACTAATAATAAAACTTTAACATCCTCGTGAGGAcatcat                                           <  2:168297/239‑1 (MQ=255)
                                                                                                                                                              cTTGACGTTTTTATGGGCAATTTTTT‑CCCGGGAAAA‑GAGGC‑TGATTT‑GGTGTTGATGAAACCCACCGCCACTCCGCTACAGCAGCTGCGCTATGACAACTCTGTTTCTTTAGTCGACAAATTGTTCGTGATGATGACGTTGGGCGATGACCGTTCGATCTACCGCACCTACGTTGATGGTCGTCTGGTGTACGAACGCAACTAATAATAAAACTTTAACATCCTCGTGAGGAcatcat                                           <  1:901925/238‑1 (MQ=255)
                                                                                                                                                               tttgaggattGTTTGGCCATTTTTTCCCCTGGAAAA‑GGGTC‑TGTTTG‑CGTGGTAGAGCACCCCCCCGCCACTCCGCTACAGCAGCTGCGCTATGACAACTCTGTTTCTTTAGTCGACAAATTGTTCGTGATGATGACGTTGGGCGATGACCGTTCGATCTACCGCACCTACGTTGATGGTCGTCTGGTGTACGAACGCAACTAATAATAAAACTTTAACATCCTCGTGAGGACATCAtt                                          <  2:903291/231‑1 (MQ=255)
                                                                                                                                                                                   ttttt‑CCCTGGAAAG‑GGGGT‑TGTTTT‑GGGGGTGATGAACCCCACCGCCACTCCGCTACAGCAGCTGCGCTATGa                                                                                                                                                                                          <  1:604442/74‑1 (MQ=255)
                                                                                                                                                                                             ggCAAA‑GAGGC‑TGATTT‑CGTGGTGATGGACCCCACCGCCACTCCGCTACAGCAGCTGCGCTATGACAACTCTGTTTCTTTAGTCGACAAATTGTTCGTGATGATGACGTTGGGCGATGACCGTTCGATCTACCGCACCTACGTTGATGGTCGTCTGGTGTACGAACGCAACTAATAATAAAACTTTAACATCCTCGTGAGGACATCATTATGTCTGGAGACATCCTACAAACACCGGa             <  1:68111/238‑1 (MQ=255)
                                                                                                                                                                                             ggCAAA‑GAGGC‑TGATTT‑CGTGGTGATGGAACCCACCGCCACTCCGCTACAGCAGCTGCGCTATGACAACTCTGTTTCTTTAGTCGACAAATTGTTCGTGATGATGACGTTGGGCGATGACCGTTCGATCTACCGCACCTACGTTGATGGTCGTCTGGTGTACGAACGCAACTAATAATAAAACTTTAACATCCTCGTGAGGACATCATTATGTCTGGAGACATCCTACAAACACCGGa             >  2:68111/1‑238 (MQ=255)
                                                                                                                                                                                              gCAAA‑GAGGC‑TGATTT‑CGTGGTGATGGAACCCACCGCCACTCCGCTACAGCAGCTGCGCTATGACAACTCTGtt                                                                                                                                                                                >  2:604442/1‑74 (MQ=255)
                                                                                                                                                                                              gCAAA‑GAGGC‑TGATTT‑CGTGGTGATGGAACCCACCGCCACTCCGCTACAGCAGCTGCGCTATGACAACTCTGTTTCTTTAGTCGACAAATTGTTCGTGATGATGACGTTGGGCGATGACCGTTCGATCTACCGCACCTACGTTGATGGTCGTCTGGTGTACGAACGCAACTAATAATAAAACTTTAACATCCTCGTGAGGACATCATTATGTCTGGAGACATCCTACAAACACCGGAc            >  2:1371942/1‑238 (MQ=255)
                                                                                                                                                                                              aaaaa‑GGGGC‑TTTTTT‑GGGGGTGTGGGACCCCCCCGCCACTCCGCTACAGCAGCTGCGCTATGACAACTCTGTTTCTTTAGTCGACAAATTGTTCGTGATGATGACGTTGGGCGATGACCGTTCGATCTACCGCACCTACGTTGATGGTCGTCTGGTGTACGAACGCAACt                                                                               <  2:51903/169‑1 (MQ=255)
                                                                                                                                                                                                     agGT‑TGATTT‑CGTGGTGAGGGAACCCACCGCCACTCCGCTACAGCAGCTGCGCTATGACAACTCTGTTTCTTTAGTCGACAAATTGTTCGTGATGATGACGTTGGGCGATGACCGTTCGATCTACCGCACCTACGTTGATGgtcgtc                                                                                                 <  1:377957/147‑1 (MQ=255)
                                                                                                                                                                                                     agGC‑TGATTT‑CGTGGTGATGGAACCCACCGCCACTCCGCTACAGCAGCTGCGCTATGACAACTCTGTTTCTTTAGTCGACAAATTGTTCGTGATGATGACGTTGGGCGATGACCGTTCGATCTACCGCACCTACGTTGATGgtcgtc                                                                                                 >  2:377957/1‑147 (MQ=255)
                                                                                                                                                                                                          tGATTT‑CGTGGTGATGGAACCCACCGCCACTCCGCTACAGCAGCTGCGCTATGACAACTCTGTTTCTTTAGTCGACAAATTGTTCGTGATGATGACGTTGGGCGATGACCGTTCGATCTACCGCACCTACGTTGATGGTCGTCTGGTGTACGAACGCAACTAATAATaaaa                                                                     >  1:51903/1‑171 (MQ=255)
                                                                                                                                                                                                                       gATGGAACCCACCGCCACTCCGCTACAGCAGCTGCGCTATGACAACTCTGTTTCTTTAGTCGACAAATTGTTCGTGATGATGACGTTGGGCGATGACCGTTCGATCTACCGCACCTACGTTGATGGTCGTCTGGTGTACGAACGCAACTAATAATAAAACTTTAACATCCTCGTGAGGACATCATTATGTCTGGAGACATCCTACAAACACCGGACGCACCAAAGc  >  1:439945/1‑226 (MQ=255)
                                                                                                                                                                                                                             |                                                                                                                                                                                                                           
TTAAAGTGGGCATGGGAACGGATATCGGTGCCGGAACCACTTTCAACATGCTGCAAACGCTGAACGAAGCCTACAAAGTATTGCAATTACAAGGCTATCGCCTCTCGGCATATGAAGCGTTTTACCTGGCCACGCTCGGCGGAGCGAAATCTCTGGGCCTTGACGATTTGATTGGCAACTTTTT‑ACCTGGCAAA‑GAGGC‑TGATTT‑CGTGGTGATGGAACCCACCGCCACTCCGCTACAGCAGCTGCGCTATGACAACTCTGTTTCTTTAGTCGACAAATTGTTCGTGATGATGACGTTGGGCGATGACCGTTCGATCTACCGCACCTACGTTGATGGTCGTCTGGTGTACGAACGCAACTAATAATAAAACTTTAACATCCTCGTGAGGACATCATTATGTCTGGAGACATCCTACAAACACCGGACGCACCAAAGC  >  NC_000913/3026724‑3027160

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 13 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 33 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑