Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,638,818 G→C G214A (GGC→GCC)  creD → putative inner membrane protein CreD

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,638,8180GC94.0% 131.2 / 3.5 50G214A (GGC→GCC) creDputative inner membrane protein CreD
Reads supporting (aligned to +/- strand):  ref base G (2/1);  new base C (20/27);  total (22/28)
Fisher's exact test for biased strand distribution p-value = 5.76e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.03e-01

TAGGCAAGCCATTTATTGTGATTAGCGTCGGGGATGCGCGTGGTATTGGTGTGGTGAAAGCGCCTGAAGTTAACGGAACGGCGCTGACCATTGAACCCGGCACCGGGTTAGAGCAAGGCGGGCAGGGCGTGCATATCCCTTTACCTGAAGGGGACTGGCGGAAGCAGAACCTGAAGCTGAATATGGCCCTGAATTTAAGCGGTACCGGCGATCTTTCTGTGGTGCCTGGCGGGCGTAATAGCGAAATGACCTTAACCAGCAACTGGCCGCATCCCAGTTTTTTAGGTGATTTTCTACCAGCCAAACGGGAAGTTAGCGAG‑TCAGGTTTTCAGGCG‑CACTGGCAAAGCAGCTGGTTTGCTAATAATCTCGGTGAGCGTTTTGCTTCAGGCAATGATACCGGCTGGGAAAACTTCCCGGCGTTTAGCGTCGCAGTAACGACGCCAGCCGATCAATACCAATTAACTGACCGG  >  NC_000913/4638590‑4639059
                                                                                                                                                                                                                                    |                                                                                                                                                                                                                                                   
tAGGCACGCCCTTTATTGTGATTTGTGACGGGGATGCGCGTTGGATTGGTGTGGTGAAAGCGCCTGAAGTTAACGGAACGGTGCTGGCCATTGAACCCGGCACCGGGTTAGAGCTAGGCGGGCAGGGCGGGCATATCCCTTTACCCGAAGGTGACTGGCGGAAGTATAACCTGAAGTTGAATATGGCCCTGAATTTAATCGGTACCGTCGATCTTTCTATGGTGCCTGCCGGGCGTAATAGCGAAAt                                                                                                                                                                                                                                   <  2:434649/247‑1 (MQ=255)
 aGGCAAGCCCTTTATTGGGATTAGCGTCGGGGAGGCGCGTGGTATTGGTGTGGTGAAAGCGCCTGAAGTTAACGGAACGGCGCTGACCTTTGAACCCGGCACCGGGTTAGAGAAAGGCGGGCAGGGCGTGCATATCCCTTTACCCGAAGGGGACTGGCGGAAGCAGAACCTGAAGCTGAATATGGCCCTGAATTTAAGCGGTACCGGCGATCTTTCTGTGGTGCATGCCGGGCGTAATAGCGAAATGCCCt                                                                                                                                                                                                                              <  2:357375/251‑1 (MQ=255)
  ggCAAGCCCTTTATTGTGATTAGCGTCGGGGATGCGCGTGGTATTGGTGTGGTGAAAGCGCCTGAAGTTAACGGAACGGCGCTGACCATTGAACCCGGCACCGGGTTAGAGCAAGGCGGGCAGGGCGTGCATATCCCTTTACCCGAAGGGGACTGGCGGAAGCAGAACCTGAAGCTGAATATGGCCCTGAATTTAAGCGGTACCGGCGATCTTTCTGTGGTGCCTGCCGGGCGTAATAGCGAAATGACCtt                                                                                                                                                                                                                             <  2:177072/251‑1 (MQ=255)
  ggCAAGCCCTTTATTGTGATTAGCGTCGGGGATGCGCGAGGTATTCGTGTGGTGAAAGCGCCTGAAGTTAACGGAACGGCGCTGACCATTGAACCCGGCACCGGGGTAGCGCAAGGCGGGCCGGGCGTGCATATCACTTTACCCGAAGGGGACTGGCGGA‑GCAGAACCTGAAGCTGAATATTGCCCTGTCTTTTAGCGGTACCTGGGATCATTATGTGGTGACTGCCGGGCGTAATAGCGaaaagccctt                                                                                                                                                                                                                             >  1:138318/1‑243 (MQ=255)
      aGCCCTTTATTGGGATTAGCGTCGGGGATGCGCGTGGTATTGGTGTGGTGAAAGCGCCTGAAGTTAACGGAACGGCGCTGACCATTGAACCCGGCACCGGGTTAGAGCAAGGCGGGCAGGGCGTGCATATCCCTTTACCCGAAGGGGACTGGCGGAAGCAGAACCTGAAGCTGAATATGGCCCTGAATTTAAGCGGTACCGGCGATCTTTCTGTGGTGCCTGCCGGGCGTAATAGCGAAATGACCTTAAcc                                                                                                                                                                                                                         <  2:82346/251‑1 (MQ=255)
        ccctttAATGTGATTAGCGTCGTGGATGCGCGTGGTATTGGTTTGGTGAAAGCGCCTGAAGGTAACGGAACGGCTCTGACCATTGAACCCGGCACCGGGTTAGAGCAAGGCGGGCAGGGCGTGCATATCCCTTTACCCGAAGGGGACTGGCGGAAGCAGAACCTGAAGCTGAATATGGCCCTGAATTTAAGCGTGACCGGCTATCTTTCGGTGGTGCCTGCCGGGCGTAAGAGCGAAATTACCTTAACCAg                                                                                                                                                                                                                       <  2:284805/248‑1 (MQ=255)
                 gtgATTAGCGTCGGGGATGCGCGTGGTATTGGTGTGGTGAAAGCGCCTGAAGTTAACGGAACGGCGCTGACCATTGAACCCGGCACCGGGTTAGAGCAAGGCGGGCAGGGCGTGCATATCCCTTTACCCGAAGGGGACTGGCGGAAGCAGAACCTGAAGCTGAATATGGCCCTGAATTTAAGCGGTACCGGCGATCTTTCTGTGGTGCCTGCCGGGCGTAATAGCGAAATGACCTTAACCAGCAACTGGcc                                                                                                                                                                                                              <  2:134438/251‑1 (MQ=255)
                         cGTCGGGGATGCGCGTGGTATTGGTGTGGTGAAAGCGCCTGAAGTTAACGGAACGGCGCTGACCATTGAACCCGGCACCGGGTTAGAGCAAGGCGGGCAGGGCGTGCATATCCCTTTACCCGAAGGGGACTGGCGGAAGCAGAACCTGAAGCTGAATATGGCCCTGAATTTAAGCGGTACCGGCGATCTTTCTGTTGTGCCTGCCGGGCGTAATAGCGAAATGACCTTAACCAGCAACTGGCCGCATccc                                                                                                                                                                                                       >  1:274168/1‑250 (MQ=255)
                            cGGGGATGCGCGTGGTATTGGTGTGGTGAAAGCGCCTGAAGTTAACGGAACGGCGCTGACCATTGAACCCGGCACCGGGTTAGAGCAAGGCGGGCAGGGCGTGCATATCCCTTTACCCGAAGGGGACTGGCGGAAGCAGAACCTGGAGCTGAATATGGCCCTGAATTTAAGCGGTACCGGCGAACTTTATGAGGTGCCTTCCGGGCGGAATAGAGAAATGACCTTAACCAGCAACTGGCCGCATCCCAGtt                                                                                                                                                                                                   >  1:51566/1‑251 (MQ=255)
                                         ggTATTGGTGTGGTGAAAGCGCCTGAAGTTAACGGAACGGCGCTGACCATTGAACCCGGCACCGGGTTAGAGCAAGGCGGGCAGGGCGTGCATATCCCTTTACCCGAAGGGGACTGGCGGAAGCAAAACCTGAAGCTGAATATGGCCCTGAATTTAAGCGGTACCGGCGATCTTTCTGTGGTGCCTGCCGGGCGTAATAGCGAAATGACCTTAACCAGCAACTGGCCGCATCCCAGTTTTTTAGGGGAttt                                                                                                                                                                                      >  1:28664/1‑251 (MQ=255)
                                           tATTGGTGTGGTGAAAGCGCCTGAAGTTAACGGAACGGCGCTGACCATTGAACCCGGCACCGGGTTAGAGCAAGGCGGGCAGGGCGTGCATATCCCTTTACCCGAAGGGGACTGGCGGAAGCAGAAACTGAAGCAGAAAATGGCCCGGGATTTAAGGGGGACCGGAGATCTTTCTGGGGTGCCTGCCGGGCGTAATAGCGAAATGACCTTAACCAGCAACTGGGCGCATCCCAGTTTTTTTGGAGAttttt                                                                                                                                                                                    >  1:200790/1‑250 (MQ=255)
                                            aTTGGTGTGGTGAAAGCGCCTGAAGTTAACGGAACGGCGCTGACCATTGAACCCGGCACCGGGTTAGAGCAAGGCGGGCAGGGCGTGCATATCCCTTTACCCGAAGGGGACTGGCGGAAGCAGAACCTGAAGCTGAATATGGCCCTGAATTTAAGCGGTACCGGCGATCTTTCTGTGGTGCCTGCCGGGCGTAATAGCGAAATGACCTTAACCAGCAACTGGCCGCATCCCAGTTTTTTAGGTGAttttt                                                                                                                                                                                    <  2:105091/250‑2 (MQ=255)
                                                gtgtggtgAAAGCGCCTGAAGTTAACGGAACGGCGCTGACCATTGAACCCGGCACCGGGTTAGAGCAAGGCGGGCAGGGCGTGCATATCCCTTTACCCGAAGGGGACTGGCGGAAGCAGAACCTGAAGCTGAATATGGCCCTGAATTTAAGCGGTACCGGCGATCTTTCTGTGGTGCCTGCCGGGCGTAATAGCGAAATGACCTTAACCAGCAACTGGCCGCATCCCAGTTTTTTAGGTGAttttt                                                                                                                                                                                    >  1:254876/1‑245 (MQ=255)
                                                gtgtggtgAAAGCGCCTGAAGTTAACGGAACGGCGCTGACCATTGAACCCGGCACCGGGTTAGAGCAAGGCGGGCAGGGCGTGCATATCCCTTTACCCGAAGGGGACTGGCGGAAGCAGAACCTGAAGCTGAATATGGCCCTGAATTTAAGCGGTACCGGCGATCTTTCTGTGGTGCCTGCCGGGCGTAATAGCGAAATGACCTTAACCAGCAACTGGCCGCATCCCAGTTTTTTAGGTGAttttt                                                                                                                                                                                    <  2:254876/246‑2 (MQ=255)
                                                 ggtggtgAAAGCGCCTGAAGTTAACGGAACGGCGCTGACCATTGAACCCGGCACCGGGTTAGAGCAAGGCGGGCAGGGCGTGCATATCCCTTTACCCGAAGGGGACTGGCGGAAGCAGAACCTGAAGCTGAATATGGCCCTGAATTTAAGCGGTACCGGCGATCTTTCTGTGGTGCCTGCCGGGCGTAATAGCGAAATGACCTTAACCAGCAACTGGCCGCATCCCAGTTTTTTAGGTGAttttttgccag                                                                                                                                                                              <  2:28664/250‑8 (MQ=255)
                                                  gtggtgAAAGCGCCTGAAGTTAACGGAACGGCGCTGACCATTGAACCCGGCACCGGGTTAGAGCAAGGCGGGCAGGGCGTGCATATCCCTTTACCCGAAGGGGACTGGCGGAAGCAGAACCTGAAGCTGAATATGGCCCTGAATTTAAGCGGTACCGGCGATCTTTCTGTGGTGCCTGCCGGGCGTAATAGCGAAATGACCTTAACCAGCAACTGGCCGCATCCCAGTTTTTTAGGTGATTTTTTGCCAGc                                                                                                                                                                             >  1:210857/1‑251 (MQ=255)
                                                  gtggtgAAAGCGCCTGAAGTTAACGGAACGGCGCTGACCATTGAACCCGGCACCGGGTTAGAGCAAGGCGGGCAGGGCGCGCATATCCCTTTACCCGAAGGGGACTGGCGGAAGCAAAACCTGAAGCAGAAAATGGACCTGAATTTAAGCGGTCCCGGCGATCTTTCTGTGGTGCCTGCCGGGCGTCAAAGCGAAATGACCTTAACCAGGAACTGGCCGCATGCCAGGTTTTTAGGGGATTTTTTGCCAGc                                                                                                                                                                             >  1:22394/1‑251 (MQ=255)
                                                                 gAAGTTAACGGAACGGCGCTGACCATTGAACCCGGCACCGGGTTAGAGCAAGGCGGGCAGGGCGTGCATATCCCTTTACCCGAAGGGGACTGGCGGAAGCAGAACCTGAAGCTGAATATGGCCCTGAATTTAAGCGGTACCGGCGATCTTTCTGTGGTGCCTGCCGGGCGTAATAGCGAAATGACCTTAACCAGCAACTGGCCGCATCCCAGTTTTTTAGGTGATTTTTTGCCAGCCCAACGGGAAGtt                                                                                                                                                                >  1:285231/1‑249 (MQ=255)
                                                                      tAACGGAACGGCGCTGACCATTGAACCCGGCACCGGGTTAGAGCAAGGCGGGCAGGGCGTGCATATCCCTTTACCCGAAGGGGACTGGCGGAAGCAGAACCTGAAGCTGAATATGGCCCTGAATTTAAGCGGTACCGGCGATCTTTCTGTGGTGCCTGCCGGGCGTAATAGCGAAATGACCTTAACCAGCAACTGGCCGCATCCCAGTTTTTTAGGTGATTTTTTGCCAGCCAAACGGGAAGTTAGCGAg                                                                                                                                                          <  1:319578/250‑1 (MQ=255)
                                                                        acggaactgCGCCGACCATTGAACCCGGCCTCGGGTTTGAGCAAGGCGGGCAGGGCGTGCATATCCCTTTACCCGAAGGTGACTGGCGTAATCAGAACCTGAAGCTGAATATGGCCCGGAATTTAAGCGGTGCCGTCGATCTTTCTGTGGGGCCTGCCGGGCGTAATGGCGAAAGGACCTTAACCAGCAACTGGGCGCATCCCAGTTTTTTGGGTGATTTTTTGCCAGCCAAACGGGAAGTTAGCGAggca                                                                                                                                                       <  2:284386/251‑4 (MQ=255)
                                                                               ggCGCTGACAATTGAACCCGGCACCGGGATAGAGCAAGGCGGGCAGGGCGTGCATATCCCTTTACCCGAAGGGGACTGGCGGAAGCAGAACCTGAAGCTGAATATGGCCCTGAATTTAAGCGGTACCGGCGATCTTTCTGTGGTGCCTGCCGGGCGTAATAGCTAAATGACCTTAACCAGCAACTGGCCGAATCCCAGTTTTTTAGGTGATTTTTTGCCAGCCAAACGGGAAGTTAGCGAG‑TCAGGTTTTc                                                                                                                                               <  2:34728/251‑1 (MQ=255)
                                                                                 cgcTGACCATTGAACCCGGCACCGGGCTAGAGCAAGGCGGGCAGGGCGTGCATATCCCTTTACCCGAAGGGGACTGGCGGAAGCAGAACCTGAAGCTGAATATGGCCCTGAATTTAAGCGGTACCGGCGATCTTTCTGTGGTGCCTGCCGGGCGTAATAGCGAAATGACCTTAACCAGCAACTGGCCGCATCCCAGTTTTTTAGGTGATTTTTTGCCAGCCAAACGGGAAGTTAGCGAG‑TCAGGTTTTc                                                                                                                                               >  1:33660/1‑249 (MQ=255)
                                                                                  gcTGACCATTGAACCCGGCACCGGGTTAGAGCAAGGCGGGCAGGGCGTGCATATCCCTTTACCCGAAGGGGACTGGCGGAAGCAGAACCTGAAGCTGAATATGGCCCTGAATTTAAGCGGTACCGGCGATCTTTCTGTGGTGCCTGCCGGGCGTAATAGCGAAATGACCTTAACCAGCAACTGGCCGCATCCCAGTTTTTTAGGTGATTTTTTGCCAGCCAAACGGGAAGTTAGCGAG‑TCAGGTTTTCAgg                                                                                                                                            <  1:288207/251‑1 (MQ=255)
                                                                                     gACCATTGAACCCGGCACCGGGTTAGAGCAAGGCGGGCAGGGCGTGCATATCCCTTTACCCGAAGGGGACTGGCGGAAGCAGAACCTGAAGCTGAATATGGCCCTGAATTTAAGCGGTACCGGCGATCTTTCTGTGGTGCCTGCCGGGCGTAATAGCGAAATGACCTTAACCAGCAACTGGCCGCATCCCAGTTTTTTAGGTGATTTTTTGCCAGCCAAACGGGAAGTTAGCGAG‑TCAGGTTTTCAGgc                                                                                                                                           <  1:4354/249‑1 (MQ=255)
                                                                                                    cACCGGGTTAGAGCAAGGCGGGCAGGGCGTGCATATCCCTTTACCCGAAGGGGACTGGCGGAAGCAGAACCTGAAGCTGAATATGGCCCTGAATTTAAGCGGTACCGGCGATCTTTCTGTGGTGCCTGCCGGGCGTAATAGCGAAATGACCTTAACCAGCAACTGGCCGCATCCCAGTTTTTTAGGTGATTTTTTGCCAGCCAAACGGGAAGTTAGCGAG‑TCAGGTTTTCAGGCG‑CAATGGCAAAGCAGCt                                                                                                                         <  1:103881/251‑1 (MQ=255)
                                                                                                      ccGGGTTAGAGCAAGGCGGGCAGGGCGTGCATATCCCTTTACCCGAAGGGGACTGGCGGAAGCAGATCCTGAAGCTGAATATGGCCCTGAATTTAAGCGGTACCGGCGATCTTTCTGTGGTGCCTGCCGGGCGTAATAGCGAAATGACCTTAACCAGCAACTGGCCGCATCCCAGTTTTTTAGGTGATTTTTTGCCAGCCAAACGGGAAGTTAGCGAG‑TCAGGTTTTCAGGCG‑CAATGGCAAAGCAGCTgg                                                                                                                       <  1:338871/251‑1 (MQ=255)
                                                                                                                       gggCAGGGCGTGCATATCCCTTTACCCGAAGGGGACTGGCGGAAGCAGAACCTGAAGCTGAATATGGCCCTGAATTTAAGCGGTACCGGCGATCTTTCTGTGGTGCCTGCCGGGCGTAATAGCGAAATGACCTTAACCAGCAACTGGCCGCATCCCAGTTTTTTAGGTGATTTTTTGCCAGCCAAACGGGAAGTTAGCGAG‑TCAGGTTTTCAGGCG‑CAATGGCAAAGCAGCTGGTTTGCTAATAATCTCgg                                                                                                      >  1:98773/1‑251 (MQ=255)
                                                                                                                           agggcggtcATATCCCTTTACCCGAAGGGGACTGGCGGAAGCAGAACCTGAAGCTGAATATGGCCCTGAATTTAAGCGGTACCGGCGATCTTTCTGTGGTGCCTGCCGGGCGTAATAGCGAAATGACCTTAACCAGCAACTGGCCGCATCCCAGTTTTTTAGGTGATTTTTTGCCAGCCAAACGGGAAGTTAGCGAG‑TCAGGTTTTCAGGCG‑CAATGGCAAAGCAGCTGGTTTGCTAATAATCTCGGTGAg                                                                                                  <  2:285231/243‑1 (MQ=255)
                                                                                                                           aGGGCGTGCATATCACTTTACCCGAAGGGGACTGGCGGAAGCAGAAACTGAAGCTGAATATGGCCCTGAATTTAAGCGGTACCGGCGATCTTTCTGTGGGGCCTGCCGGGCGTAATAGCGAAATGACCTTAACCAGCAACTGGCAGCATCCAAGTTTTTTAGGTGATTTTTTGCCAGCCAAACGGGAAGTTAGCGAG‑TCAGGTTTTCAGGAGTCA‑TGGCAAAGCAGCTGGTTTGCTAATAATCTCGGTg                                                                                                    >  1:118254/1‑249 (MQ=255)
                                                                                                                                gTGCATATCCCTTTACCTGAAGGGGACTGGCGGAAGCAGAACCTGAAGCTGAATATGGCCCTGAATTTAAGCGGTACCGGCGATCTTTCTGTGGTGCCTGGCGGGCGTAATAGCGAAATGACCTTAACCAGCAACTGGcc                                                                                                                                                                                                              <  1:101647/140‑1 (MQ=255)
                                                                                                                                gTGCATATCCCTTTACCCGAAGGGGACTGGCGGAAGCAGAACCTGAAGCTGAATATGGCCCTGAATTTAAGCGGTACCGGCGATCTTTCTGTGGTGCCTGCCGGGCGTAATAGCGAAATGACCTTAACCAGCAACTGGCCGCATCCCAGTTTTTTAGGTGATTTTTTGCCAGCCAAACGGGAAGTTAGCGAG‑TCAGGTTTTCAGGCG‑CAATGGCAAAGCAGCTGGTTTGCTAATAATCTCGGTGAGCGttt                                                                                             <  1:24101/251‑1 (MQ=255)
                                                                                                                                 tGCATATCCCTTTACCCGAAGGGGACTGGCGGCAGCAGAACATGAAGCTGAATATGGCCCTGAATTTAAGCGGTACCGGCGATCTTTCTGGGGTGCCTGCCGGGCGGACTAGCGAAATGGCCTTAACCAGCAACTGGGCGAATCCCAGTTTTTTAGGTGATTTTTTGCCAGCCACACGGGAAGTTAGGGAGATC‑GGGTTTCAGGCG‑GAATGGGAAAGAAGCTGGTTTGCtattaaatccggtgggggtttt                                                                                            >  2:58897/1‑235 (MQ=255)
                                                                                                                                  gCATATCCCTTTACCTGAAGGGGACTGGCGGAAGCAGAACCTGAAGCCGCATATGGACCTGAATTTAAGCGGTACCGGCGATCTTTCTGTGGTGCCTGGCGGGCGTAATAGCGAAATGACCTTAACCAGCAACTGGCCGCATCCCAGTTTTTTAGGTGATTTTCTACCAGCCAAACGGGAAGTTAGCGAG‑TCAGGTTTTCAGGCG‑CACTGGCAAAGCAGCTGGTTTGCTAATAATCTCGGTGAGCGTTTTg                                                                                           >  1:174892/1‑251 (MQ=255)
                                                                                                                                      atCCCTTTACCCGAAGGGGACTGGCGGAAGCAGAACCTGAAGCTGAATATGGCCCTGAATTTAAGCGGTACCGGCGATCTTTCTGTGGTGCCTGCCGGGCGTAATAGCGAAATGACCTTAACCAGCAACTGGCCGCATCCCAGTTTTTTAGGTGATTTTTTGCCAGCCAAACGGGAAGTTAGCGAG‑TCAGGTTTTCAGGCG‑CAATGGCAAAGCAGCTGGTTTGCTAATAATCTCGGTGAGCGTTTTGCTTc                                                                                       >  1:390238/1‑251 (MQ=255)
                                                                                                                                                ccgAAGGGGACTGGCGGAAGCAGAACCTGAAGCTGAATAAGGGCCTGAATTTAAGCGGTACCGGCGATCTTTCTGTGGGGCCTGCCGGGCGTAATAGCGAAATGACCTTAACCAGCAACTGGCCGCATCCCCGTTTTTTAGGTGATTTTTTGCCAGCCAAACGGGAAGTTAGCGAG‑TCAGGTTTTCAGGCG‑CAATGGCAAAGCAGCTGGTTTGCTAATAATCTCGGGGAGCGTTTTGCTTCAGGCAATGa                                                                              >  1:377795/3‑250 (MQ=255)
                                                                                                                                                 cgAAGGGGACTGGCGGAAGCAGAACCTGAAGCTGAATATGGCCCTGAATTTAAGCGGTACCGGCGATCTTTCTGTGGTGCCTGCCGGGCGTAATAGCGAAATGACCTTAACCAGCAACTGGCCGCATCCCAGTTTTTTAGGTGATTTTTTGCCAGCCAAACGGGAAGTTAGCGAG‑TCAGGTTTTCAGGCG‑CAATGGCAAAGCAGCTGGTTTGCTAATAATCTCGGTGAGCGTTTTGCTTCAGGCAATGAt                                                                             >  1:283157/2‑250 (MQ=255)
                                                                                                                                                  gAAGGGGACTGGCGGAAGCAGAACCTGAAGCTGAATATGGCCCTGGATTTAAGCGGTACCGGCGATCTTTCTTTGGTGCCTGCCGGGCGGAATATCGAAATGACCTTAACCCGCAACTGGCCGCATCCCCGTTTTTTAGGAGATTTTTTGCCCGCCAAACGGGAAGTTAGCGAG‑TCAGGTTTTCAGACG‑CAATAGCAACACCGCTGGTTTGCTAATAATCTCGGTGAGCGGTTTGCTTCAGGCAATgtga                                                                            >  2:128913/1‑247 (MQ=255)
                                                                                                                                                       ggACTGGCGGAAGCGGAACTTGAAGCTGAATATGGCCCTGAAATTAAGCGGTACCGGCGATCTTTCTGTGGTGCCTGCCGGGCGTAATAGCGAAATGACCTTAACCAGCAACTGGCCGCATCCCAGTTTTTTAGGTGATTTTTTGCCAGCCAAACGGGAAGTTAGCTAG‑TCAGGTTTTCAGGCG‑CAATGGCAAAGCAGCTGGTTTGCTAATAATCTCGGTGAGCGTTTTGCTTCAGGCAATGATACCGGCt                                                                      <  2:51116/251‑1 (MQ=255)
                                                                                                                                                                     acaaCCTGAAGCTGAATATGGCCCTGAATTTAAGCGGTACCGTCGATCTTTCTGTGGTGCCTGCCGGGCGTAATAGCGAAATGACCTTAACCAGCAACTGGCCGCATCCCAGTTTTTTAGGTGATTTTTTGCCAGCCAAACGGGAAGTTAGCGAG‑TCAGGTTTTCAGGCG‑CAATGGCAAAGCAGCTGGTTTGCTAATAATCTCGGTGAGCGTTTTGCTTCAGGCAATGATACCGGCTGGGAAAACTTcctc                                                        <  2:200790/249‑3 (MQ=255)
                                                                                                                                                                     aGAACCTGAAGCTGAATATGGCCCTGAATTTACGCGGTACCGGCGATCTTTCTGTGGTGCCTGGCGGGCGTAATAGCGAAATGAACTTAACCAGCAACTGGCCGCATCCCAGTTTTTTAGGTGATTTTCTACCAGCCAAACGGGAAGTTAGCGAG‑TCAGGTTTTCAGGCG‑CAATGGCATAGCAGCTGGTTTGCTAATAATCTCGATGAGCGATTTGCTTCAGGCAATGATAACGGCTGGGAAAACtttccc                                                        >  1:146923/1‑247 (MQ=255)
                                                                                                                                                                      gAACCTGAAGCTGAATATGGCCCTGAATTTAAGCGGTACCGGCGATCTTTCTGTGGTGCCTGCCGGGCGTAATAGCGAAATGACCTTAACCAGCAACTGGCCGCATCCCAGTTTTTTAGGTGATTTTTTGCCAGCCAAACGGGAAGTTAGCGAG‑TCAGGGTTTCAGGCG‑CAATGGCAAAGCAGCTGGTTTGCTAATAATCTCGGTGAGCGTTTTGCTTCAGGCAATGATACCGGCTGGGAAAAATTccccg                                                       >  1:406393/1‑249 (MQ=255)
                                                                                                                                                                        aCCTGAAACTGAATATGGCCCTGAATTTAAGCGGTACCGGCGATCTTTCTGTGGTGCCTGCCGGGCGTAAATGCGAAATGACCTTAACCAGCAACTGGCCGCATCCCAATTTTGTAGGTGATTTTTTGCCAGCCAAACGGGAAGTTAGCGAG‑TCAGGTTTTCAGGCG‑CAATGGCAAAGCAGCTGGTTTGCTAATAATCTCGGTGAGCGTTTTGCTTCAGGCAATGATACCGGCTGGGAAAACTTCCCcgcg                                                     <  2:377795/251‑1 (MQ=255)
                                                                                                                                                                                              gAATTTAAGCGGTACCGGCGATCTTTCTGTGGTGCCTGCCGGGCGTAATAGCGAAATGACCTTAACCAGCAACTGGCCGCATCCCCGTTTTTTAGGTGATTTTTTGCCAGCCAAACGGGAAGTTAGCGAG‑TCAGGTTTTCAGGCG‑CAATGGCAACGCAGCTGGTTTGCTAATAATCTCTGTGGGCGGTTTGCTTCAGGCAATGATCCCGGCTGGGAAAAATTACCCGCGTTTAAGGTCGCAgagacaaccg                               >  1:269500/1‑242 (MQ=255)
                                                                                                                                                                                                aTTTAAGCGGTACCGGCGATCTTTCTGTGGTGCCTGCCGGGCGTAATAGCGAAATGACCTTAACCAGCAACTGGCCGCATCCCAGTTTTTTAGGTGATTTTTTGCCAGCCAAACGGGAAGTTAGCGAG‑TCAGGTTTTCAGGCG‑CAATGGCAAAGCAGCTGGTTTGCTAATAATCTCGGTGAGCGTTTTGCTGCAGGCAATGATACCGGCTGGGAAAACTTCCCCGCGGTTAGCGTCGCAGTGACGACGCCa                             <  1:378572/251‑1 (MQ=255)
                                                                                                                                                                                                    aaGCGGTCCCGGCGATCTTTCTGTGGTGCCTGCCGGGCGTAATAGCGAAATGACCTTAACCAGCAACTGGCCGCATCCCAGTTTTTTAGGTGATTTTTTGCCAGCCAAACGGGAAGTTAGCGAG‑TCAGGTTTTCAGGCG‑CAATGGCAAAGCAGCTGGTTTGCTAATAATCTCGGTGAGCGTTTTGCTTCAGGCAATGATACCGGCTGGGAAAACTTCCCCGCGTTTAGCGTCGCAGTGACGACGCCAGcc                          <  2:283157/250‑1 (MQ=255)
                                                                                                                                                                                                     aGCGGTCCCGGCGATCTTGCTGTGGTGCCGGCCGGGCGTAATAGCGAAATGACCTTAACCAGCAACTGGCCGCATCCCAGTTTGTTAGGTGATTTTTTGCCAGCCAAACGGGAAGTTAGCGAG‑TCAGGTTTTCAGGCG‑CAATGGCAAAGCAGCTGGTTTGCTAATAATCTCGGTGAGCGTTTTGCTTCAGGCAATGATACCGGCTGGGAAAACTTCCCCGCGTTTAGCGTCGCAGTGACGACGCCAGCCGa                        <  2:445126/251‑1 (MQ=255)
                                                                                                                                                                                                         gttccGGCGATCTTTCTGGGGTGCCTGCCGGGCGTAATAGCGAAATGACCTTAACCAGCAACTGGCCGCATCCCAGTTTTTTAGGTGATTTTTTGCCAGCCAAACGGGAAGTTAGCGAG‑GCAGGTTTTCAGGCG‑CAATGGCAAAGCAGCTGGTTTGCTAATAATCTCGGTGAGCGTTTTGCTTCAGGCAATGATTCCGGCTGGGAAAACTTCCCCGCGTTTAGCGTCGCAGTGACGACGCCAGCCGATCaa                    <  1:418483/248‑1 (MQ=255)
                                                                                                                                                                                                                      ttCTGTGGTGCCTGCCGGGCGTAATAGCGAAATGACCTTAACCAGCAACTGGCCGCATCCCAGTTTTTTAGGTGATTTTTTGCCAGCCAAACGGGAAGTTAGCGAG‑TCAGGTTTTCAGGCG‑CAATGGCAAAGCAGCTGGTTTGCTAATAATCTCGGTGAGCGTTTTGCTTCAGGCAATGATACCGGCTGGGAAAACTTCCCCGCGTTTTGCGTCGCCGTGACGCCGGCAGCCGGTCAAttcacattaa          >  2:165486/1‑239 (MQ=255)
                                                                                                                                                                                                                      ttCTGTGGTGCCTGCCGGGCGTAATAGCGAAATGACCTTAACCAGCAACTGGCCGCATCCCAGTTTTTTAGGTGATTTTTTGCCAGCCAAACGGGAAGTTAGCGAG‑TCAGGTTTTCAGGCG‑CAATGGCAAAGCAGCTGGTTTGCTAATAATCTCGGTGAGCGTTTTGCTTCAGGCAATGATACCGGCTGGGAAAACTTCCCCGCGTTTAGCGTCGCAGTGACGACGCCAGCCGATCAATACCAATTaa          <  1:165486/248‑1 (MQ=255)
                                                                                                                                                                                                                           gggggCCTGCCGGGCGTAAAAGCGAAATGACCTTAACCAGCAACTGGCCGCATCCCAGTTTTTTAGGTGATTTTTTGCCAGCCAAACGGGAAGTTAGCGAT‑TCAGGTTTTCAGGCG‑CAATGGAAAAGCAGCTGGTTTGCTAATAATCTCGGTGAGCGTTTTGCTTCAGGCAATGATACCGGCTGGGAAAACTTCCCCGCGTTTAGCGTCGCAGGGACGACGCCAGCCGAGCAATACCAATTAAcggtccgg  <  2:367211/247‑8 (MQ=255)
                                                                                                                                                                                                                                    |                                                                                                                                                                                                                                                   
TAGGCAAGCCATTTATTGTGATTAGCGTCGGGGATGCGCGTGGTATTGGTGTGGTGAAAGCGCCTGAAGTTAACGGAACGGCGCTGACCATTGAACCCGGCACCGGGTTAGAGCAAGGCGGGCAGGGCGTGCATATCCCTTTACCTGAAGGGGACTGGCGGAAGCAGAACCTGAAGCTGAATATGGCCCTGAATTTAAGCGGTACCGGCGATCTTTCTGTGGTGCCTGGCGGGCGTAATAGCGAAATGACCTTAACCAGCAACTGGCCGCATCCCAGTTTTTTAGGTGATTTTCTACCAGCCAAACGGGAAGTTAGCGAG‑TCAGGTTTTCAGGCG‑CACTGGCAAAGCAGCTGGTTTGCTAATAATCTCGGTGAGCGTTTTGCTTCAGGCAATGATACCGGCTGGGAAAACTTCCCGGCGTTTAGCGTCGCAGTAACGACGCCAGCCGATCAATACCAATTAACTGACCGG  >  NC_000913/4638590‑4639059

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