Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,638,735 T→C P186P (CCT→CCC 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,7350TC93.5% 129.7 / ‑1.7 46P186P (CCT→CCCcreDputative inner membrane protein CreD
Reads supporting (aligned to +/- strand):  ref base T (2/1);  new base C (14/29);  total (16/30)
Fisher's exact test for biased strand distribution p-value = 2.74e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.97e-01

ATTTATACCGGTCAGGTCTGGCACAGTGATTTAACGTTAAAAGCCGATTTCGATGTTTCGCGTCTTAGCGAACTCAACGCGCCAAATATCACCTTAGGCAAGCCATTTATTGTGATTAGCGTCGGGGATGCGCGTGGTATTGGTGTGGTGAAAGCGCCTGAAGTTAACGGAACGGCGCTGACCATTGAACCCGGCACCGGGTTAGAGCAAGGCGGGCAGGGCGTGCATATCCCTTTACCTGAAGGGGACTGGCGGAAGCAGAACCTGAAGCTGAATATGGCCCTGAATTTAAGCGGTACCGGCGATCTTTCTGTGGTGCCTGGCGGGCGTAATAGCGAAATGACCTTAACCAGCAACTGGCCGCATCCCAGTTTTTTAGGTGATTTTCTACCAGCCAAACGGGAAGTTAGCGAGTCAGGTTTTCAGGCGCACTGGCAAAGCAGCTGGTTTGCTAATAATCTCGGTGAGCGTTTTGCTTCAGGCAATG  >  NC_000913/4638496‑4638982
                                                                                                                                                                                                                                               |                                                                                                                                                                                                                                                       
aTTTATACCGGTCAGGTCTGGCACAGTGATTTAACGCTAAAAGCCGATTTCGATGTTTCGCGTCTTAGCGAGCTCAACGCGCCAAATATCATCTTAGGCAAGCCCTTTATTGTGATTAGCGTCGGGGATGCGCGTGGTCTTGGTGTGGTGAAAGCGCCTGAAGTTACAGGAACGGCGCTGACCATTGAACCCGGCCACGGGGTAGAGCAAGGCGGGCAAGGCGTGCATAACCCTTTACCCGCAGGGGGCTg                                                                                                                                                                                                                                              >  2:379258/1‑251 (MQ=255)
   tataCCGGTCAGGTCTGGCACAGTGATTTAACGCTAAAAGCCGATTTCGATGTTTCGCGTCTTAGCGAGCTCAACGCGCCAAATATCATCTTAGGCAAGCCCTTTATTGTGATTAGCGTCGGGGATGCGCGTGGTATTGGTGTGGTGAAAGCGCCTGAAGTTAACGGAACGGCGCTGACCATTGAACCCGGCACCGGGTTAGAGCAAGGCGGGCAGGGCGTGCATATCCCTTTACCCGAAGGGGACTGGc                                                                                                                                                                                                                                            <  1:180003/250‑1 (MQ=255)
      aCCGGTCAGGTCTGGCACAGTTATTTAACGCTAAAAGCCGATTTCGATGTTTCGCGTCTTAGAGAGCTCAACGCGCAAAATATCATCTTAGGCAAGCCCTTTATTGAGATTAGCGTCGGGGAGGCGCGGGGTATTGGTGGGGTGAAAGCGCCTGAAGTCAACGTAAAGGCGCTGACCATTGAACACGGCCCCGGGTTAGAGCAAGGCGGGGCGGGGGGGCATATCCCTTCACCCGAAaggggcatgagggc                                                                                                                                                                                                                                        >  2:4996/1‑241 (MQ=255)
         ggTCAGGTCTGGCACAGTGATTTAACGCTAAAAGCCGATTTCGATGTTTCGCGTCTTAGCGAGCTCAACGCGCCAAATATCATCTTAGGCAAGCCCTTTATTGTGATTAGCGTCGGGGATGCGCGTGGTATTGGTGTGGTGAAAGCGCCTGAAGTTAACGGAACGGCGCTGACCATTGAACCCGGCACCGGGTTAGAGCAAGGCGGGCAGGGCGTGCATATCCCTTTACCCGAAGGGGACTGGCGGGAGc                                                                                                                                                                                                                                      >  1:14395/1‑250 (MQ=255)
             aGGTCTGGCACAGTGATTTAACGCTAAAAGCCGATTTCGATGTTTCGCGTCTTAGCGAGCTCAACGCGCCAAATATCATCTTAGGCAAGCCCTTTATTGTGATTAGCGTCGGGGATGCGCGTGGTATTGGTGTGGTGAAAGCGCCTGAAGTTAACGGAACGGCGCTGACCATTGAACCCGGCACCGGGTTAGAGCAAGGCGGGCAGGGCGTGCATATCCCTTTACCCGAAGGGGACTGGCGGAAGCAGAAc                                                                                                                                                                                                                                 <  1:380662/251‑1 (MQ=255)
                   ggCACAGTGATTTAACGCTAAAAGCCGATTTCGATGTTTCGCGTCTTAGCGAGCTCAACGCGCCAAATATCATCTTAGGCAAGCCCTTTATTGTGATTAGCGTCGGGGATGCGCGTGGTATTGGTGTGGTGAAAGCGCCTGAAGATAACGGAACGGCGCTGACCATTGAACCCGGCCCCGGGGTAGAGCAAGGCGGGCAGGGCGTGCATATCCCTTTACCCGAAGGGGGCTGGCGGAAGCAAAACCTg                                                                                                                                                                                                                              >  2:219591/1‑248 (MQ=255)
                            aTTTAACGCTAAAAGCCGATTTCGATGTTTCGCGTCTTAGCGAGCTCAACGCGCCAAATATCATCTTAGGCAAGCCCTTTATTGTGATTAGCGTCGGGGATGCGCGTGGTATTGGTGTGGTGAAAGCGCCTGAAGTTAACGGAACGGCGCTGACCATTGAACCCGGCACCGGGTTAGAGCAAGGCGGGCGGGGCGTGCATATCCCTTTACCCGCAGGGGACTGGCgggagaggaaacttaaccgtaatatg                                                                                                                                                                                                                  >  2:477872/1‑227 (MQ=255)
                                 acgctAAAAGCCGATTTCGATGTTTCGCGTCTTAGCGAGCTCAACGCGCCAAATATCATCTTAGGCAAGCCCTTTATTGTGATTAGCGTCGGGGATGCGCGTGGTATTGGTGTGGTGAAAGCGCCTGAAGTTAACGGAACGGCGCTGACCATTGAACCCGGCACCGGGTTAGAGCAAGGCGGGCAGGGCGTGCATATCCCTTTACCCGAAGGGGACTGGCGGAAGCAGAACCTGAAGCTGAATATGGCCCt                                                                                                                                                                                                             <  2:406772/247‑1 (MQ=255)
                                        aaGCCTATTGCGATGGTACGCGTCTGAGCGAGCTCAACGCTCCAAATTTCATTTGAGGCAAGCCCTTGATTGTGATTAGCGTCGGGGATGCGCGTGGTATTGGTGTGGTGAAAGCGCCTGAAGTTAACGGAACGGCGCTGACCATTGAACCCGGCACCGGGTTAGAGCAAGGCGGGCAGGGCGTGCATAGCCCGTGAACCGAAGGTGACTGGCGGAAGCAGAACCTGAAGAGGAATATGGCCCTGAATTTa                                                                                                                                                                                                      <  2:673996/251‑1 (MQ=255)
                                        aaGCCGATTTCGATGTTTCGCGTCTTAGCGAGCTCAACGCGCCAAATATCATCTTAGGCAAGCCCTTTATTGTGATTAACGTCGGGGAAGCGCGTGGTATTGGTGTGGTGAAAGCGCCTGAAGTTAACGGAAAGGCGCTGACCATTGAACCCGGCACCGGGTTAGAGCAAGGCGGGCAGGGCGTGCATAACCCTTTACCCGAAGGGGACTGGCGGCAGAAGAACAGGAAGCTGAATATGGCCATGGAttt                                                                                                                                                                                                       >  1:415980/1‑250 (MQ=255)
                                           ccGATTTCGATGTTTCGCGTCTTAGCGAGCTCAACGCGCCAAATATCATCTTAGGCAAGCCCTTTATTGTGATTAGCGTCGGGGATGCGCGTGGTATTGGTGTGGTGAAAGCGCCTGAAGTTAACGGAACGGCGCTGACCATTGAACCCGGCACCGGGTTAGAGCAAGGCGGGCAGGGCGTGCATATCCCTTTACCCGAAGGGGACTGGCGGAAGCAGAACCTGAAGCTGAATATGGCCCTGAATTTAAGc                                                                                                                                                                                                   <  1:453262/251‑1 (MQ=255)
                                                ttccaagTTTCGCCTCTTTGCGAGCTCAACTCGCCAAATATCTTCTTAGGCAAGCCCTTTATTTTGATTAGCGTCGGTGATGCGCGTGGTATTGGTGTGGTGAAAGCTCCTGAAGTTAACGGCACGGCGCGGACCATTGAACCCGGCACCGGGTTAGAGCAAGGCGGGCAGGGCTTGCATATCCCTTTACCCGAAGGTGACTGGCGGAAGCATATCCGGAAGCTGAATAGGGCCCTGAATTTAAGCGGTa                                                                                                                                                                                               <  2:335489/244‑1 (MQ=255)
                                                   gATGTTTCGCGTCTTAGCGAGCTCAACGCGCCAAATATCATCTTAGGCAAGCCCTTTATTGTGATTAGCGTCGGGGATGCGCGTGGTATTGGTGTGGTGAAAGCGCCTGAAGTTAACGGAACGGCGCTGACCATTGAACCCGGCACCGGGTTAGAGCAAGGCGTGCAGGGCGTGCATATCCCTTTACCCGAAGGGGACTGGCGGAAGCAGAACCTGAAGCTGAATATGGCCCTGAATTTAAGCGGTACCgg                                                                                                                                                                                           <  1:4996/251‑1 (MQ=255)
                                                            cgTCTTAGCGAGCTCAACGCGCCAAATATCATCTTAGGCAAGCCCTTTATTGTGATTAGCGTCGGGGATGCGCGTGGTATTGGTGTGGTGAAAGCGCCTGAAGTTAACGGAACGGCGCTGACCATTGAACCCGGCACCGGGTTAGAGCAAGGCGGGCAGGGCGTGCATATCCCTTTACCCGAAGGGGACTGGCGGACGAAGAACCTGCAGCTGAATATGGCCCTGAATTTTAGCGGGACCGGCGAttcttt                                                                                                                                                                                  >  1:384571/1‑246 (MQ=255)
                                                            cgTCTTAGCGAGCTCAACGCGCCAAATATCATCTTAGGCAAGCCCGTTATTGTGATTAGCGTCGGGGATGCGCGTGGTATTGGTGTGGTGAAAGCGCCTGAAGTTAACGGAACGGCGCTGACCATTGAACCCGGCACCGGGTTAGAGCAAGGCGGGCAGGGCGTGCATATCCCTTTACCCGAAGGGGACTGGCGGAAGCAGAACCTGAAGCTGAATATGGCCCTGAATTTAAGCGGTACCGGCGATCTTTc                                                                                                                                                                                  <  1:647741/251‑1 (MQ=255)
                                                                 tAGCGAGCTCAACGCGCCAAATATCATCTTAGGCAAGCCCTTTATTGTGATTAGCGTCGGGGATGCGCGTGGTATTTGTGTGGTGAAAGCGCCTGAAGTTAACGGAACGGCGCTGACCATTGAACCCGGCACCGGGTTAGAGCAAGGCGGGCAGGGCGTGCATATCCCTTTACCCGAAGGGGACTGGCGGAAGCAGAACCTGAAGCTGAATATGGCCCTGAATTTAAGCGGTACCGGCGATCTTTCtg                                                                                                                                                                                <  1:369662/248‑1 (MQ=255)
                                                                       gcttaaCGCGCAAAATATCATCTTATGCAAACCCTTTATTGTGATTAGCGTAGGGGATGCGCGGTGTATTGGTGGGGTGAAAGCGCCTGAAGTTAACGGAACGGCGCTGACCATTGAACCCGGACCCGGGTTAGAGCAAGGCGGGCAGGGCGTGCATATCCCTTTACCCGAAGGGGACTGGCGGAAGCAGAACCTGAAGCTGAATATTGCCCTGAATTTAAGCGGTACCGGCGATCTTTCTGTGTTGCCTg                                                                                                                                                                       <  2:536789/247‑1 (MQ=255)
                                                                           aaCGCGCCAAATATCATCTTAGGCAAGCCCTTTATTGTGATTAGCGTCGGGGATGCGCGTGGTATTGGTGTGGTGAAAGCGCCTGAAGTTAACGGAACGGCGCTGACCATTGAACCCGGCACCGGGTTAGAGCAAGGCGGGCAGGGCGTGCATATCCCTTTACCCGAAGGGGACTGGCGGGATCAGAAACTTGAGCTGAATAGGACCCCGAATTTAAGCGGGACCGGCgctcgtttttttggtgcagccgg                                                                                                                                                                   >  2:222857/1‑229 (MQ=255)
                                                                             cgcgCCAAATATCATCTTAGGCAAGCCCTTTATTGTGATTAGCGTCGGGGATGCGCGTGGTATTGGTGTGGTGAAAGCGCCTGAAGTTAACGGAACGGCGCTGACCATTGAACCCGGCACCGGGTTAGAGCAAGGCGGGCAGGGCGTGCATATCCCTTTACCCGAAGGGGACTGGCGGAAGCAGAACCTGAAGCTGAATATGGCCCTGGATTTAAGCGGTACCGGCGATCTTTCTGTGGGGACTgccgggc                                                                                                                                                                 >  1:510358/1‑251 (MQ=255)
                                                                                 ccAAATATCATCTTAGGCAAGCCCTTTATTGTGACTAGCGTCGGGGATGCGCGTGGTATTGGTGTGGTGAAAGCGCCTGAAGTTAACGGAACGGCGCTGACCATTGAACCCGGCACCGGGTTAGAGCAAGGCGGGCAGGGCGTGCATATCCCTTTACCCGaa                                                                                                                                                                                                                                                      <  1:151548/162‑1 (MQ=255)
                                                                                                   aaGCCATTTATTGTGATTAGCGCCGGGGATGCGCGTGGTATTGGTGTGGTGAAAGCGCCTGAAGTTAACGGAACGGCGCTGACCATTGAACCCGGCACCGGGTTAGAGCAAGGCGGGCAGGGCGTGCATATCCCTTTACCTGAAGGGGACTGGCGGAAGCAGAACCTGAAGCTGAATATGGCCCTGAATTTAAACGGTACCGGCGATCTTTCTGTGGGGCCTGGCGGGCGTAATAGCGAAATTACCTTAAc                                                                                                                                           >  1:283632/1‑251 (MQ=255)
                                                                                                                tgATTAGCGTCGGGGATGCGCGTGGTATTGGTGTGGTGAAAGCGCCTGAAGTTAACGGAACGGCGCTGACCATTGAACCCGGCACCGGGTTAGAGCAAGGCGGGCAGGGCGTGCATATCCCTTTACCCGAAGGGGACTGGCGGAAGCAGGAACTGAAGCTGAATATGGCCCTGAATTTAAGCGGGACCGGCGGTCTTTTTGTGGTGCCTGGCGGGCGTAATAGCGAAATGACCTTAACCAGCAACTGGCCg                                                                                                                              >  1:524681/1‑251 (MQ=255)
                                                                                                                      gcggcGGGGATGCGCGTGGTATTGGGGTGGTGAAAGCGCCGGAAGTTAACGGAACGGCGCTGACCATTGAACCCGGCACCGGGTTAGAGCAAGGCGGGCAGGGCGTGCATATCCCTTTACCCGAAGGGGACTGGCGGAAGCAGAACCTGAAGCTGAATATGGCCCTGAATTTAAGCGGTACCGGCGATCTTTCTGTGGTGCCTGCCGGGCGTAATAGCGAAATGACCTTAACCAGCAACTGGCCGCATccc                                                                                                                        <  2:153704/247‑1 (MQ=255)
                                                                                                                      gCGTCGGGGATGCGCGTGGTATTGGTGTGGTGAAAGCGCCTGAAGTTAACGGAACGGCGCTGACCATTGAACCCGGCACCGGGTTAGAGCAAGGCGGGCAGGGCGTGCATATCCCTTTACCCGAAGGGGACTGGCGGAAGCAGAACCTGAAGCTGAATATGGCCCTGAATTTAAGCGGTACCGGCGATCTTTCTGTGGTGCCTGCCGGGCGTAATAGCGAAATGACCTTAACCAGCAACTGGCCGCATccc                                                                                                                        <  1:222857/251‑1 (MQ=255)
                                                                                                                        gTCGGTGATGCGCGTGGGATGTGTGTGTTGAAAGCGCCTGAAGTTAACGGAACGGCGCTGACCATTTAACCCGGCACCGGGTTAGAGCATGGCGGGCATGGCGTGTATATCCCTTTACCCGAAGGGTACTGGCGGAAGCAGAACCTGAAGCTGAATATGGCCCTGAATTTAAGCGTTACCGGCGATCTTTCTGTGGTGCCTGCCGGGCGTAATAGCGAAATGACCTTAACCAGCAACTTGCCGCATCCCAg                                                                                                                      <  2:471144/251‑1 (MQ=255)
                                                                                                                                      tGGTATTGGTGTGGTGAAAGCGCCTGAAGTTAACGGAACGGCGCTGAACATTGAACCCGGCACCGGGTTAGAGCAAGGAGGGCAGGGCGTGCATATCCCTTTACCCGAAGGGGACTGGCGGAAGCAGAAACTGGAGCTTAATATTGACCAGAATTTAAGCGGGAACGGGGATCTTTTCGTGGTGCCTTCCGGGAGTTAAAGCGAAATGAACTTAACAATAAACTGGCCGCATCCCAGTTTTTTAgggggtt                                                                                                        >  2:430129/1‑246 (MQ=255)
                                                                                                                                      gggTATTGGTGTGGTGAAAGCGCCTGAAGTTAACGGAACGGCGCTGACCATTGAACCCGGCACCGGGTTAGAGCAAGGCGGGCAGGGCGTGCATATCCCTTTACCCGAAGGGGACTGGCGGAAGCAGAACCTGAAGCTGAATATGGCCCTGAATTTAAGCGGTACCGGCGATCTTTCTGTGGTGCCTGCCGGGCGTAATAGCGAAATGACCTTAACCAGCAACTGGCCGCATCCCAGTTTTTTAGGTGAtt                                                                                                        <  1:357736/250‑1 (MQ=255)
                                                                                                                                          aTTGGTGTGGTGAAAGCGCCTGAAGTTTACGGAACGGCGCTGACCATTGAACCCGGTACCGGGTTAGAGCAAGGCGGGCAGGGCGTTCATATCCCTTTTCCCGAAGGGGACTGGCGGAAGCAGAACCTGAAGCTGAATATGGCCCTGAATTTAAGCGGTACCGGCTATCTTTCTGTGGTGCCTGCCGGGCGTAATAGCGAAATGACCTTAACCAGCAACTGGCCGCATCCCAGTTTTTTAGGTGAttttt                                                                                                     <  1:260541/250‑2 (MQ=255)
                                                                                                                                          aTTGGTGTGGTGAAAGCGCCTGAAGTTAACGGAACGGCGCTGACCATTGAACCCGGCACCGGGTTAGAGCAAGGCGGGCAGGGCGTGCATATCCCTTTACCCGAAGGGGACTGGCGGAAGCAGAACCTGAAGCTGAATATGGCCCTGAATTTAAGCGGTACCGGCGATCTTTCTGTGGTGCCTGCCGGGCGTAATAGCGAAATGACCTTAACCAGCAACTGGCCGCATACCAGTTTTTTAGGTGAttttt                                                                                                     <  1:379258/250‑2 (MQ=255)
                                                                                                                                             ggtgtggtgAAAGCGCCTGAAGTTAACGGAACGGCGCTGACCAGTGAACCCGGCATCGGGTTAGAGAAAGGCGGGCAGGGCGTGCATATCCCTTTACCCGAAGGGGACTGGCGGAAGCAGAACCTGAAGCTGAATATGGCCCTGAATTTAAGCGGTACCGGCGATCTTTCTGTGGTGCCTGCCGGGCGTAATAGCGAAATGACCTTAACCATCAACTGTCCGCAGCCCAGTTTTTTAGGTGAttttttgc                                                                                                  <  1:543198/250‑5 (MQ=255)
                                                                                                                                              gtgtggtgAAAGCGCCTGAAGTTAACGGAACGGCGCTGACCATTGAACCCGGCACCGGGTTAGAGCAAGGCGGGCAGGGCGTGCATATCCCTTTACCCGAAGGGGACTGGCGGAAGCAGAACCTGAAGCTGAATATGGCCCTGAATTTAAGCGGTACCGGCGATCTTTCTGTGGTGCCTGCCGGGCGTAATAGCGAAATGACCTTAACCAGCAACTGGCCGCATCCCAGTTTTTTAGGTGAttttttgcca                                                                                                <  1:38523/251‑7 (MQ=255)
                                                                                                                                                gtggtgAAAGCGCCTGAAGTTAACGGAACGGCGCTGACCATTGAACCCGGCACCGGGTTAGAGCAAGGCGGGCAGGGCGTGCATATCCCTTTACCCGAAGGGGACTGGCGGAAGCAGAACCTg                                                                                                                                                                                                                              <  2:620649/123‑1 (MQ=255)
                                                                                                                                                gtggtgAAAGCGCCTGAAGTTAACGGAACGGCGCTGACCATTGAACCCGGCACCGGGTTAGAGCAAGGCGGGCAGGGCGTGCATATCCCTTTACCCGAAGGGGACTGGCGGAAGCAGAAACTg                                                                                                                                                                                                                              >  1:620649/1‑123 (MQ=255)
                                                                                                                                                              ggaatttaccggaacggCGCTGACCATTGAACCCGGCACCGGGTTAGAGCAAGGCGGGCAGGGCGTGCATATCCCTTTACCCGAAGGGGACTGGCGGAAGCAGAACCTGAAGCTGAATATGGCCCTGAATTTAAGCGGTACCGGCGATCTTTCTGTGGTGCCTGCCGGGCGTAATAGCGAAATGACCTTAACCAGCAACTGGCCGCATCCCAGTTTTTTAGGTGATTTTTTGCCAGCCAACCGGGAAGTTa                                                                                <  2:14395/242‑1 (MQ=255)
                                                                                                                                                                     aacggaacggCGCTGACCATTGAACCCGGCACCGGGTTAGAGCAAGGCGGGCAGGGCGTGCATATCCCTTTACCCGAAGGGGACTGGCGGAAGCAGAACCTGAAGCTGAATATGGCCCTGAATTTAAGCGGTACCGGCGATCTTTCTGTGGTGCCTGCCGGGCGGAATAGCGAAATGACCTTAACCAGCCACTGGCCGCATCCCAGTTTTTTAGGTGATTTTTTGCCAGCCAAACGGGAAGTTAGCGAGTc                                                                         >  1:234905/1‑251 (MQ=255)
                                                                                                                                                                       cggaacggCGCTGACCATTGAACCCGGCACCGGGTTAGAGCAAGGCGGGCAGGGCGTGCATATCCCTTTACCCGAAGGGGACTGGCGGAAGCAGAACCTGAAGCTGAATATGGCCCTGAATTTAAGCGGTACCGGCGATCTTTCTGTGGTGCCTGCCGGGCGTAATAGCGAAATGACCTTAACCAGCAACTGGCCGCATCCCAGTTTTTTAGGTGATTTTTTGCCAGCCAAACGGGAAGTTAGCGaatc                                                                         >  1:614657/1‑246 (MQ=255)
                                                                                                                                                                        ggaacggCGCTGACCATTGAACCCGGCACCGGGTTAGAGCAAGGCGGGCAGGGCGTGCATATCCCTTTACCCGAAGGGGACTGGCGGAAGCAGAACCTGAAGCTGAATATGGCCCTGAATTTAAGCGGTACCGGCGATCTTTCTGTGGTGCCTGCCGGGCGTAATAGCGAAATGACCTTAACCAGCAACTGGCCGCATCCCAGTTTTTTAGGTGATTGTTTGCCAGCCAAACGGGAAGTTAGCGAGTCAgg                                                                      <  1:191940/251‑1 (MQ=255)
                                                                                                                                                                                      cATTGAACCCGGCACCGGGTTAGAGCAAGGCGGGCAGGGCGTGCATATCCCTTTACCCGAAGGGGACTGGCGGAAGCAGAACCTGAAGCTGAATATGGCCCTGAATTTAAGCGGTACCGGCGATCTTTCTGTGGTGCCTGCCGGGCGTAATAGCGAAATGACCTTAACCAGCAACTGGCCGCATCCCAGTTTTTTAGGTGATTTTTTGCCAGCCAAACGGGAAGTTAGCGAGTCAGGTTTTCAGGCGCaat                                                        <  1:577543/251‑3 (MQ=255)
                                                                                                                                                                                      cATTGAACCCGGCACCGGGTTAGAGCAAGGCGGGCAGGGCGTGCATATCCCTTTACCCGAAGGGGACTGGCGGAAGCAGAACCTGAAGCTGAATATGGCCCTGAATTTAAGCGGTACCGGCGATCTTTCTGTGGTGCCTGCCGGGCGTAATAGCGAAATGACCTTAACCAGCAACTGGCCGCATCCCAGTTTTTTAGGTGATTTTTTGCCAGCCAAACGGGAAGTTAGCGAGTCAGGTTTTCAGGCGCaat                                                        <  1:440251/251‑3 (MQ=255)
                                                                                                                                                                                                            gagCAAGGCGGGCAGGGCGTGCATATCCCTTTACCCGAAGGGGACTGGCGGAAGCAGAACCTGAAGCTGAATATGGCCCTGAATTTAAGCGGTACCGGCGATCTTTCTGTGGTGCCTGCCGGGCGTAATAGCGAAATGACCTTAACCAGCAACTGGCCGCATCCCAGTTTTTTAGGTGATTTTTTGCCAGCCAAACGGGAAGTTAGCGAGTCAGGTTTTCAGGCGCAATGGCAAAGCAGCTGGTTTGCtaa                                  <  1:132353/251‑1 (MQ=255)
                                                                                                                                                                                                                 aGGCGGGCAGGGAGGGCATATCCCTTTACCTGAAGGGGAATGGCGGAAGCAGAACCTGAAGCTGAATATGGCCCTGAATTTAAACGGTACCGGCGATCTTTCTGGGGTGCCTGGCGGGCGTAATAGCGAAATGACCTTAACCAGCAACTGGCCGCAGCCCAGTTTTTTAGGTGATTTTCTACCAGCCAAACGGTAAGTTAGCGAGTCAGGTTTTCAGGCGCAATGGCAAAGCAGCTGGTTTGCTAATAAtc                             <  2:283632/251‑1 (MQ=255)
                                                                                                                                                                                                                   gCGGGCAGGGCGTGCATATCCCTTTACCCGAAGGGGACTGGCGGAAGCAGAACCTGAAGCTGAATATGGCCCTGAATTTAAGCGGTACCGGCGATCTTTCTGTGGTGCCTGCCGGGCGTAATAGCGAAATGACCTTAACCAGCAACTGGCCGCATCCCAGTTTTTTAGGTGATTTTTTGCCAGCCAAACGGGAAGTTAGCGAGTCAGGTTTTCAGGCGCAATGGCAAAGCAGCTGGTTTGCTAATAAtctc                           <  1:386817/251‑1 (MQ=255)
                                                                                                                                                                                                                           ggCGTGCATATCCCTTTACCCGAAGGGGACTGGCGGAAGCAGAACCTGAAGCTGAATATGGCCCTGAATTTAAGCGGTACCGGCGATCTTTCTGTGGTGCCTGCCGGGCGTAATAGCGAAATGACCTTAACCAGCAACTGGCCGCATCCCAGTTTTTTAGGTGATTTTTTGCCAGCCAAACGGGAAGTTAGCGAGTCAGGTTTTCAGGCGCAATGGCAAAGCAGCTGGTTTGCTAATAATCTCGGTGAGCg                   <  1:112640/251‑1 (MQ=255)
                                                                                                                                                                                                                                gCATATCCCTTTACCCGAAGGGGACTGGCGTAAGCAGAACCTGAAGCTGAATATGGCCCTGAATTTAAGCGGTACCGGCGATCTTTCTGTGGTGCCTGCCGGGCGTAATAGCGAAATGACCTTAACCAGCAACTGGCCGCATCCCAGTTTTTTAGGTGATTTTTTGCCAGCCAAACGGGAAGTTAGCGAGTCAGGTTTTCAGGCGCAATGGCAAAGCAGCTGGTTTGCTAATAATCTCGGTGAGCGTTTTg              <  1:658643/251‑1 (MQ=255)
                                                                                                                                                                                                                                          ttACCCGAAGGGGACTGGCTGAAGCATAACCTGAAGCTGAATATGGCCCTGAATTTAAGCGGTACCGGCTATCTTTCTGTGGTGCCTGCCGGGCGTAATAGCGAAATGACCTTAACCAGCAACTGGCCGCATCCCAGTTTTTTAGGTGATTTTTTGCCAGCCAAACGGGAAGTTAGCGAGTCAGGTTTTCAGGCGCAATGGCAAAGCAGCTGGTTTGCTAATAATCTCGGTGAGCGTTTTGCTTCAGGCa     <  1:555033/250‑1 (MQ=255)
                                                                                                                                                                                                                                            aCCTGAAGGGGACTGGCGGAAGCAGAACCTGAAGCTGAATATGGCCCTGAATTTAAGCGGTACCGGCGCTCTTTCTGTGGTGCCTGGCGGGCGTAATAGCGAAATGACCTTAACCAGCAACTGGCCGCATCCCAGTTTTTTAGGTGATTTTCTACCAGCCAAACGGGAAGTTAGCGAGTCAGGTTTTCAGGCGCACTGGCAAAGCAGCTGGTTTGCTAATAATCTCGGTGAGCGTTTTGCTTCAGGCAATg  >  1:649176/1‑251 (MQ=255)
                                                                                                                                                                                                                                               |                                                                                                                                                                                                                                                       
ATTTATACCGGTCAGGTCTGGCACAGTGATTTAACGTTAAAAGCCGATTTCGATGTTTCGCGTCTTAGCGAACTCAACGCGCCAAATATCACCTTAGGCAAGCCATTTATTGTGATTAGCGTCGGGGATGCGCGTGGTATTGGTGTGGTGAAAGCGCCTGAAGTTAACGGAACGGCGCTGACCATTGAACCCGGCACCGGGTTAGAGCAAGGCGGGCAGGGCGTGCATATCCCTTTACCTGAAGGGGACTGGCGGAAGCAGAACCTGAAGCTGAATATGGCCCTGAATTTAAGCGGTACCGGCGATCTTTCTGTGGTGCCTGGCGGGCGTAATAGCGAAATGACCTTAACCAGCAACTGGCCGCATCCCAGTTTTTTAGGTGATTTTCTACCAGCCAAACGGGAAGTTAGCGAGTCAGGTTTTCAGGCGCACTGGCAAAGCAGCTGGTTTGCTAATAATCTCGGTGAGCGTTTTGCTTCAGGCAATG  >  NC_000913/4638496‑4638982

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