Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 207,473 T→G 31.5% M783R (ATG→AGG)  dnaE → DNA polymerase III alpha subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913207,4730TG31.5% 53.2 / 3.8 40M783R (ATG→AGG) dnaEDNA polymerase III alpha subunit
Reads supporting (aligned to +/- strand):  ref base T (18/8);  new base G (2/12);  total (20/20)
Fisher's exact test for biased strand distribution p-value = 2.20e-03
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GCGGATATGCTGCGTCGTGCGATGGGTAAGAAAAAGCCGGAAGAGATGGCTAAGCAACGTTCTGTATTTGCTGAAGGTGCAGAAAAGAACGGAATCAACGCTGAACTGGCGATGAAAATC‑TTCGACCTGGTGG‑AGAAATTCGC‑TGGTTACGGATTTAACAAATCGCACTCTGCGGCCTATGCTTTGGTGTCATATCAAACGTTATGGCTGAAAGCGCACTATCCTGCGGAGTTTATGGCGGCGGT‑AATGACCGCCG‑ATATGGACAACACCGAGAAGGTGGTGGGTCTGGTGGATGAGTGCTGGCGGATGGGGCTGAAAATCCTGCCACCAGATATAAACTCCGGTCTTTACCATTTCCACGTCAACGACGACGGCGAAATCGTGTATGGTATTGGCGCGATCAAAGGGGTCGGTGAAGGTCCGATTGAGGCCATCATCGAAGCCCGTAATAAAGGCGGCTACTTCCGCGA  >  NC_000913/207238‑207707
                                                                                                                                                                                                                                              |                                                                                                                                                                                                                                            
gcgGATATGCTGCGTCGTGCGATGGGTAAGAAAAAGCCGGAAGAGATGGCTAAGCAACGTTCTGTATTTGCTGAAGGTGCAGAAAAGAACGGAATCAACGCTGAACTGGCGATGAAAATC‑TTCGACCTGGTGG‑AGAAATTCGC‑TGGTTACGGATTTAACAAATCGCACTCTGCGGCCTATGCTTTGGTGTCATATCAAACGTTATGGCTGAAAGGGCACTATCTTTCGGAGCTTATggc                                                                                                                                                                                                                                           >  1:923063/1‑239 (MQ=255)
 cgGATATGCTGCGTCGTGCGATGGGTAAGAAAAAGCCGGAAGAGATGGCTAAGCAACGTTCTGTATTTGCTGAAGGTGCAGAAAAGAACGGAATCAACGCTGAACTGGCGATGAAAATC‑TTCGACCTGGTGG‑AGAAATTCGC‑TGGTTACGGATTTAACAAATCGCACTCTGCGGCCTATGCTTTGGTTTCATATCAAACGTTATGGCTGAAAGCGCCCCATCCTTCGGGGGTTAGggcg                                                                                                                                                                                                                                          >  1:843911/1‑239 (MQ=255)
  gaatatGCTGCGGCGGGCGGTGGGTAAGGAAAAGCCGGAAGAGATGGCTAAGCAGCGTTCTGTATTTGCTGAAGGTGCAGAAAAGAACGGAATCAACGCTGAACTGGCGATGAAAATC‑TTCGACCTGGTGG‑AGAAATTCGC‑TGGTTACGGATTTAACAAATCGCACGCTGCGGCCTATGCTTTGGTGTCATATCAAACGTTATGGCTGAAAGCGCACTATCCTGCGGAGTTTATggcg                                                                                                                                                                                                                                          <  2:116761/236‑1 (MQ=255)
   gATAGGGTGCGTCGTGCGATGGGTAAGAAAAAGCCGGAATGAATGGCTAAGCAACGTTCTGTATTTGCTGAATGTGCAGAAAAGAACGGAATCAACGCTGAACTGGCGATGAAAATC‑TTCGACCTGGTGG‑AGAAATTCGC‑TGGTTACGGATTTAACAAATCGCACTCTGCGGCCTATGCTTTGGTGTCATATCAAACGTTATGGCTGAAAGCGCACTATCCTGCGGAGTTTATggcgg                                                                                                                                                                                                                                         <  1:1022336/238‑1 (MQ=255)
            cgtcgtGCGATGGGTAAGAAAAAGCCGGAAGAGATGGCTAAGCAACGTTCTGTATTTGCTGAAGGTGCAGAAAAGAACGGAATCAACGCTGAACTGGCGATGAAAATC‑TTCGACCTGGTGG‑AGAAATTCGC‑TGGTTACGGATTTAACAAATCGCACTCTGCGGCCTATGCTTTGGTGTCATATCAAACGTTATGGCTGAAAGCGCCCTATCCTGCGGAGTTTATGGCGGCGGT‑AAtta                                                                                                                                                                                                                               >  1:774207/1‑236 (MQ=255)
             gtcgtGCGATGGGTAAGAAAAAGCCGGAAGAGATGGCTAAGCAACGTTCTGTATTTGCTGAAGGTGCAGAAAAGAACGGAATCAACGCTGAACTGGCGATGAAAATC‑TTCGACCTGGTGG‑AGAAATTCGC‑TGGTTACGGATTTAACAAATCGCACTCTGCGGCCTATGCTTTGGTGTCATATCAAACGTTATGGCTGAAAGCGCACTATCCTGCGGAGTTTATGGCGGCGGT‑AATGAc                                                                                                                                                                                                                              >  1:706665/1‑238 (MQ=255)
              tcgtGCGATGGGTAAGAAAAAGCCGGAAGAGATGGCTAAGCAACGTTCTGTATTTGCTGAAGGTGCAGAAAAGAACGGAATCAACGCTGAACTGGCGATGAAAATC‑TTCGACCTGGTGG‑AGAAATTCGC‑TGGTTACGGATTTAACAAATCGCACTCTGCGGCCTATGCTTTGGTGTCATATCAACCGTTATGGCTGAAAGGGCACTATCCTGCGGAGTTTAGGGCGG‑GGT‑aaataaca                                                                                                                                                                                                                            >  1:651150/1‑232 (MQ=255)
              tcgtGCGATGGGTAAGAAAAAGCCGGAAGAGATGGCTAAGCAACGTTCTGTATTTGCTGAAGGTGCAGAAAAGAACGGAATCAACGCTGAACTGGCGATGAAAATC‑TTCGACCTGGTGG‑AGAAATTCGC‑TGGTTACGGATTTAACAAATCGCACTCTGCGGCCTATGCTTTGGTGTCATATCAAACGTTATGGCTGAAAGCGCACTATCCTGCGGAGTTTATGGCGGCGGC‑AAAGAcc                                                                                                                                                                                                                             >  1:290680/1‑238 (MQ=255)
              tcgtGCGATGGGTAAGAAAAAGCCGGAAGAGATGGCTAAGCAACGTTCTGTATTTGCTGAAGGTGCAGAAAAGAACGGAATCAACGCTGAACTGGCGATGAAAATC‑TTCGACCTGGTGG‑AGAAATTCGC‑TGGTTACGGATTTAACAAATCGCACTCTGCGGCCTATGCTTTGGTGTCATATCAAACGTTATGGCTGAAAGCGCACTATCCTCCGGAGTTTATGGCGGCGGt‑gatgaac                                                                                                                                                                                                                             >  1:586103/1‑236 (MQ=255)
                 tGCGATGGGTAAGAAAAAGCCGGAAGAGATGGCTAAGCAACGTTCTGTATTTGCTGAAGGTGCAGAAAAGAACGGAATCAACGCTGAACTGGCGATGAAAATC‑TTCGACCTGGTGG‑AGAAATTCGC‑TGGTTACGGATTTAACAAATCGCACTCTGCGGCCTATGCTTTGGTGTCATATCAAACGTTATGGCTGAAAGCGCACTATCCTGCGGGGTTTATGCCGGCGGT‑AATggaccccc                                                                                                                                                                                                                         >  1:635114/1‑232 (MQ=255)
                                  aGCCGGAAGAGATGGCTAAGCAACGTTCTGTATTTGCTGAAGGTGCAGAAAAGAACGGAATCAACGCTGAACTGGCGATGAAAATC‑TTCGACCTGGTGG‑AGAAATTCGC‑TGGTTACGGATTTAACAAATCGCACTCTGCGGCCTATGCTTTGGTGTCATATCAAACGTTATGGCTGAAAGCGCACTATCCTGCGGGGTTTATTGCGGCGGTAAATGA‑CGCCG‑GCATGGACCAAAcccgg                                                                                                                                                                                                       >  2:1061549/1‑236 (MQ=255)
                                  aGCCGGAAGAGATGGCTAAGCAACGTTCTGTATTTGCTGAAGGTGCAGAAAAGAACGGAATCAACGCTGAACTGGCGATGAAAATC‑TTCGACCTGGTGG‑AGAAATTCGC‑TGGTTACGGATTTAACAAATCGCACTCTGCGGCCTATGCTTTGGTGTCATATCAAACGTTATGGCTGAAAGCGCACTATCCTGCGGAGTTTATGGCGGCGGG‑AAGGACCGCCG‑ATAAGGACAAAACCggg                                                                                                                                                                                                       >  1:529788/1‑237 (MQ=255)
                                  aGCCGGAAGAGATGGCTAAGCAACGTTCTGTATTTGCTGAAGGTGCAGAAAAGAACGGAATCAACGCTGAACTGGCGATGAAAATC‑TTCGACCTGGTGG‑AGAAATTCGC‑TGGTTACGGATTTAACAAATCGCACTCTGCGGCCTATGCTTTGGGGTCATATAAAACGTTATGGCTGAAAGCGCACTATCCTGCGGGGGTTATGGCGGCGGG‑AAACACCCCCGAATATGGAAAAACCCga                                                                                                                                                                                                        >  2:881315/1‑239 (MQ=255)
                                   cgcGGAAAAGAGGGCAAAGCAACGTTTTGTATTGGCGGAAGGGCAAAAAAAGAACGAAATCAACGCTGAACTGGCGATGAAAATC‑TTCGACCTGGTGG‑AGAAATTCGC‑TGGTTACGGATTTAACAAATCGCACTCGGCGGCCTATGCTTGGTTGTCATATCAACCGTTATGGCTGAAAGCGCACTATCCTGGGGAGTTTATGGCGGCGGT‑AATGACCGCCG‑ATATGGACAACACCgaga                                                                                                                                                                                                      <  1:604393/237‑1 (MQ=255)
                                    cgggAAGATGTGGAAAAGCAAGTTTTTGTATTTGCGGAAGGGGAAAAAAAGAACGGAATCAACGCTGAACGGGAGATGAAAATC‑TTCGCCGTGGTGG‑AGAAATTCGC‑TGGTTACGAATTAAACAAATCGCACTCGGCGGCTTATGCTTTGGTGTAATATCAAACTTTATGGCTGAAAGCGCACTATCCTGCGGAGTTTATGGCGGCGGT‑AATGACCGCCG‑ATATGGACAACACCgag                                                                                                                                                                                                       <  1:408258/235‑1 (MQ=255)
                                         agagATGGCTAAGGAATGTTCTTTATTTGCTGAAGGTGAAGAAAAGAACGGAACCACCGCTGAACTGGCAATGAAAACC‑TTCTCCCTGGTGG‑AAAAATTCGC‑TGGTTACGGATTTAACAAATCGCACTCGGCGGCCTATGCTTTGGTGTCATATCAAACGTTATGGCTGAAAGCGCACTATCCGGCGGAGTTTATGGCGGCGGT‑AATGACCGCCG‑ATATGGACAACACCGAGAAggtgg                                                                                                                                                                                                <  2:269548/239‑1 (MQ=255)
                                               gggtAAGCAATGTTCTTTATTTGCGGAAGGTGAAAAAAAGCACGGAAATAACGTTGACTGGGCGATAAAAATC‑TTCAACCTGGTGG‑AAAAATCCGC‑TGGTTACGGATTTAACAAATCGCACTCGGCGCCCTATGCTTTGGTGTCATATCAAACTTTATGGCTGAAAGCGCACTATCCTGCGGAGTTTATGGGGGCGGT‑AATGACCGCCG‑ATATGGACAACACCGAGAAGGTGGTGGGTc                                                                                                                                                                                          <  1:237247/236‑1 (MQ=255)
                                                    aGCAACGTTCTGTATTTGCTGAAGGTGCAGAAAAGAACGGAATCAACGCTGAACTGGCGATGAAAATC‑TTCGACCTGGTGG‑AGAAATTCGC‑TGGTTACGGATTTAACAAATCGCACTCTGCGGCCTATGCTTTGGTGTCATATCAAACGTTATGGCTGAACGCGCACTATCCTGCGGAGTTTATGGCGGCGGG‑AATGACCGCCG‑TTATGGACACCACCGGGAAGGTGGGGGGTCtgggt                                                                                                                                                                                     >  1:639567/1‑237 (MQ=255)
                                                    aGCAACGTTCTGTATTTGCTGAAGGTGCAGAAAAGAACGGAATCAACGCTGAACTGGCGATGAAAATC‑TTCGACCTGGTGG‑AGAAATTCGC‑TGGTTACGGATTTAACAAATCGCACTCTGCGGCCTATGCTTTGGTGTCATATCAAACGTTATGGCTGAAAGCGCACTATCCTGCGGAGTTTATGGCGGCGGG‑AATGAACGCCG‑ATATGGGCAAAACCGCGAAGGTGGGGGGTCtcgtg                                                                                                                                                                                     >  1:311921/1‑239 (MQ=255)
                                                          gTTCTGTATTTGCTGAAGGTGCAGAAAAGAACGGAATCAACGCTGAACTGGCGATGAAAATC‑TTCGACCTGGTGG‑AGAAATTCGC‑TGGTTACGGATTTAACAAATCGCACTCTGCGGCCTATGCTTTGGTGTCATATCAAACGTTATGGCTGAAAGCGCACTATCCTGCGGAGTTTATGGCGGCGGT‑AATGACCGCCG‑ATATTGACAAAACCGAGAAGGAGGTGGGTGtggtggggagg                                                                                                                                                                               >  1:447021/1‑234 (MQ=255)
                                                                        gAAGGTGCAGAAAAGAACGGAATCAACGCTGAACTGGCGATGAAAATC‑TTCGACCTGGTGG‑AGAAATTCGC‑TGGTTACGGATTTAACAAATCGCACTCTGCGGCCTATGCTTTGGTGTCATATCAAACGTTATGGCTGAAAGCGCACTATCCTGCGGAGTTTATGGCGGCGGG‑AATGACCGCCC‑GTATGGACAAACCCGCAAGAGTGGGGGGTCTGGGGGGGGGGGGGTGGCGGAgggg                                                                                                                                                                 >  1:462788/1‑239 (MQ=255)
                                                                                       aaCGGAATCAACGCTGAACTGGCGATGAAAATC‑TTCGACCTGGTGG‑AGAAATTCGC‑TGGTTACGGATTTAACAAATCGCACTCTGCGGCCTATGCTTTGGTGCCATATCAAACGTTATGGCTGAAAGCGCACTATCCTGCGGAGTTTATGGCGGCGGT‑AATGACCGCCG‑ATATTGACAACAACGAGAAGGTGGGGGGTCTGGTGGGTGAGGGCTGGCGGATGGGGCTGGAAAACCTGc                                                                                                                                                   >  1:26969/1‑238 (MQ=255)
                                                                                         cGGAATCAACGCTGAACTGGCGATGAAAATC‑TTCGACCTGGTGG‑AGAAATTCGC‑TGGTTACGGATTTAACAAATCGCACTCTGCGGCCTATGCTTTGGTGTCATATCAAACGTTATGGCTGAAAGCGCACTATCCTGCGGAGTTTATGGCGGCGGT‑AATGAACGCCG‑CTATGGGCAAAAACGAGAAGGGGGCGGGTCTGGTGGATGAGGGCTGGCGGGTGGGGCTGAAAATCCTGtcac                                                                                                                                                >  1:331340/1‑239 (MQ=255)
                                                                                         cGGAATAACCGCGGAGCGGGCGATGAAATTC‑TTCGCCCGGGTGG‑AAAAATTCGC‑TGGTTACGTATTTAAAAAATCGCCCTCTGCGGCCTATGCTTTGGTGTCATACCAACCTTTATGGCTGAAAGCGCACTACCCTGGGGAGTTTAGGGCGGCGGT‑AAGGACCGCCG‑ATATGGACAACACCGAGAAGGTGGTGGGTCTGGTGGATGAGTGCTGGCGGATGGGGCTGAAAATCCTGccac                                                                                                                                                <  1:926488/239‑1 (MQ=255)
                                                                                          ggatacaagtgtGAACTGGCGA‑GAAAATT‑GTCGACCTGGGGGAAAAAATTAGT‑TGGTTGGGGTTTTAACAAACCGCAGTCGCCGCTCTCTGGTTTGGGGTGAGTTCAGTCGTTAGGGCTGAAGGCGCTCGCTCCGGGGGAGTTTAGGGGGGCGTT‑AATGACCGTCG‑ATATGGACAACACCGAGAAGGTGGTGGGTCTGGTGGATGAGTGCTGGCGGATGGGGCTGAAAATCCTGccac                                                                                                                                                <  2:737571/227‑1 (MQ=255)
                                                                                           gAATCAACGGTAAAGTGGCGA‑GTAAATCGTTCGACTTGG‑GG‑GGAAATTCCCTTGTTTGCGGTTTTAACAATTCGCTCCTGGCGCCTTATGTTTTGGTGTCATATCAAACGTTAGGGGTAAAAGCGCACTACCCTGGGGATTTTAGGGCGGGGGT‑AATGACCGCCG‑ATATGGACAACACCGAGAAGGTGGTGGGTCTGGTGGATGAGTGCTGGCGGATGGGGCTGAAAATCCTGccac                                                                                                                                                <  1:881315/237‑1 (MQ=255)
                                                                                            aaTCACCGCTGACTGGGAAATAAAACTC‑TTCGGCGTGGTGG‑AAAAGTTCGT‑GGGTAGCGGATTACACAACCCGCCTTGTGCGGCTTTGGTTTGGGTGTCAATACAAGCTTTAGGGCGAAAGGCGCACTACCCTGCGGAGTTTGGGGCGGCGTT‑AATGACCGCCG‑ATATGGACAACACCGAGAAGGTGGTGGGTCTGGTGGATGAGTGCTGGCGGATGGGGCTGAAAATCCTGccacca                                                                                                                                              <  2:706665/238‑1 (MQ=255)
                                                                                                      gAACTGGCGATGAAAATC‑TTCGACCTGGTGG‑AGAAATTCGC‑TGGTTACGGATTTAACAAATCGCACTCTGCGGCCTATGCTTTGGTGTCATATCAAACGTTATGGCTGAAAGCGCACTATCCTGCGGAGTTTATGGCGGCGGT‑AATGACCGCCG‑ATATGGACAAAACCGGGAAGGTGGGGGGTCTGGTGGGTGAGTGCTGGCGGATGGGGCCGAAAAACCTTccccccagttaaaaat                                                                                                                                    >  1:716342/1‑227 (MQ=255)
                                                                                                                                                      tcgggATTTAACAAATCGCCCCCTGGGGCCTTGGCTTGGGTTTAATATAAACTGTTATGGGTGAAGCCGCCTTACCCTGGGGAGTTTTGGGGGGCGGT‑AATGACCGCCG‑ATATGGACAACACCGAGAAGGTGGTGGGTCTGGTGGATGAGTGCTGGCGGATGGGGCTGAAAATCCTGCCACCAGATATAAACTCCGGTCTTTACCATTTCCACGTCAACGACGACGGCGAAATCGTGTa                                                                                      <  1:586815/236‑1 (MQ=255)
                                                                                                                                                      tagggATTAAACAAAGCACACTTGGGGGCTTTTGTTTGGGTAACATAAAAAACGTTATGGCAGAAAGCCCACTTTCCTGGGGGTTTTAGGGCGGCGGT‑AATGACCGCCG‑ATATGGACAACACCGAGAAGGTGGTGGGTCTGGTGGATGAGTGCTGGCGGATGGGGCTGAAAATCCTGCCACCAGATATAAACTCCGGTCTTTACCATTTCCACGTCAACGACGACGGCGAAATCGTGTa                                                                                      <  1:601483/236‑1 (MQ=255)
                                                                                                                                                      taggaaTTTAAAAAATCCCCCGCTGCGGCCATTGTTTTGGTTGTATACAAAACGTTAGGGATGAAAGCGCACTACCCTGCGGAGTTTAGGGCGGCGGT‑AATGACCGCCG‑ATATGGACAACACCGAGAAGGTGGTGGGTCTGGTGGATGAGTGCTGGCGGATGGGGCTGAAAATCCTGCCACCAGATATAAACTCCGGTCTTTACCATTTCCACGTCAACGACGACGGCGAAATCGTGTa                                                                                      <  1:762252/234‑1 (MQ=255)
                                                                                                                                                              tAACAAATCGCACTCTGCGGCCTATGCTTTGGTGTCATATCAAACGTTATGGCTGAAAGCGCACTATCCTGCGGAGTTTATGGCGGCGGT‑AATGACCGCCG‑ATATGGACAACACCGAGAAGGGGGTGGGTCTGGTGGAGGAGTGCTGGCGGGTGGGGGTGGAAATCCTGCCACCAGATATAAAATCCGGTCTTTTACATTTCCACGGCCACGACGACCGGGAAAACTGGTATGGGAttt                                                                              >  1:916260/1‑238 (MQ=255)
                                                                                                                                                                aCAAACCGCATTTGGCGGCTTGTGTTTTGGGGTCATATCAAACGTTAGGGGTAAAAGCGCACTGTCCTGGGGATTTTAGGGCGGGGTA‑AATGACCGCCG‑ATATGGACAACACCGAGAAGGTGGTGGGTCTGGTGGATGAGTGCTGGCGGATGGGGCTGAAAATCCTGCCACCAGATATAAACTCCGGTCTTTACCATTTCCACGTCAACGACGACGGCGAAATCGTGTATGGTATTGgc                                                                            <  1:92067/239‑1 (MQ=255)
                                                                                                                                                                               cGGCTTATGTTTGGGGTTAATATCAAAGGTTTGGGGTGAAAGCGCACTCTCCTGCGGATTTTATGGGGGGGGT‑AATGACCGCCG‑ATATGGACAACACCGAGAAGGTGGTGGGTCTGGTGGATGAGTGCTGGCGGATGGGGCTGAAAATCCTGCCACCAGATATAAACTCCGGTCTTTACCATTTCCACGTCAACGACGACGGCGAAATCGTGTATGGTATTGGCGCGATCAAAGGGGt                                                              <  1:590082/238‑1 (MQ=255)
                                                                                                                                                                                         ttgggtttaatatCAAGCGTTATGGCTGAAGGCGCACTATCCGGCGGGGTTTAGGGGGGCGGT‑AATGACCGCCG‑ATATGGACAACACCGAGAAGGTGGTGGGTCTGGTGGATGAGTGCTGGCGGATGGGGCTGAAAATCCTGCCACCAGATATAAACTCCGGTCTTTACCATTTCCACGTCAACGACGACGGCGAAATCGTGTATGGTATTGGCGCGATCAAAGGGGTCGGTGAAGGTc                                                   <  1:655832/230‑1 (MQ=255)
                                                                                                                                                                                                   tataAAACGTTATGGCTAAAAGCGCACTATCCGGCGGAGTTTATGGCGGCGGT‑AATGACCGCCG‑ATATGGACAACACCGAGAAGGTGGTGGGTCTGGTGGATGAGTGCTGGCGGATGGGGCTGAAAATCCTGccac                                                                                                                                                <  1:619568/136‑1 (MQ=255)
                                                                                                                                                                                                                               tcgcgggggtgTTTAGGGGGGGGGT‑AATGGCCGCCG‑ATATGGACAACACCGAGAAGGTGGTGGGTCTGGTGGATGAGTGCTGGCGGATGGGGCTGAAAATCCTGCCACCAGATATAAACTCCGGTCTTTACCATTTCCACGTCAACGACGACGGCGAAATCGTGTATGGTATTGGCGCGATCAAAGGGGTCGGTGAAGGTCCGATTGAGGCCATCATCGAAGCCCGTAATAAAggcggc             <  1:765437/229‑1 (MQ=255)
                                                                                                                                                                                                                                 ccTGCGGTTTTTAGGGCGGCTGA‑TAAGCCCGCCG‑ATATGGACAACACCGAGAAGGTGGTGGGTCTGGTGGATGAGTTCTGGCGGATGGGGCTGAAAATCCTGCCACCAGATATAAACTCCGGTCTTTACCATTTCCACGTCAACGACGACGGCGAAATCGTGTATGGTATTGGCGCGATCAAAGGGGTCGGTGAAGGTCCGATTGAGGCCATCATCGAAGCCCGTAATAAAGGCGGCTa           <  1:683781/239‑1 (MQ=255)
                                                                                                                                                                                                                                          tttATGGCGGCGGT‑AATGACCGCCG‑ATATGGACAACACCGAGAAGGGGGTGGGTCTGGTGGATGTGGGCTGGCGGATGGGGCTGAAAATCCTGCCACCAGGTATAAACTCCGGTCTTTACCATTTCCACGTCAACGACGACGGCGAAATCGTGTATTGGATTGgc                                                                            >  1:133405/1‑165 (MQ=255)
                                                                                                                                                                                                                                           ttAGGGGGGCGGT‑AATGACCGCCG‑ATATGGACAACACCGAGAAGGTGGTGGGTCTGGTGGATGAGTGCTGGCGGATGGGGCTGAAAATCCTGCCACCAGATATAAACTCCGGTCTTTACCATTTCCACGTCAACGACGACGGCGAAATCGTGTATGGTATTGGCGCGATCAAAGGGGTCGGTGAAGGTCCGATTGAGGCCATCATCGAAGCCCGTAATAAAGGCGGCTACTTCCGCGa  <  1:4049/238‑1 (MQ=255)
                                                                                                                                                                                                                                              |                                                                                                                                                                                                                                            
GCGGATATGCTGCGTCGTGCGATGGGTAAGAAAAAGCCGGAAGAGATGGCTAAGCAACGTTCTGTATTTGCTGAAGGTGCAGAAAAGAACGGAATCAACGCTGAACTGGCGATGAAAATC‑TTCGACCTGGTGG‑AGAAATTCGC‑TGGTTACGGATTTAACAAATCGCACTCTGCGGCCTATGCTTTGGTGTCATATCAAACGTTATGGCTGAAAGCGCACTATCCTGCGGAGTTTATGGCGGCGGT‑AATGACCGCCG‑ATATGGACAACACCGAGAAGGTGGTGGGTCTGGTGGATGAGTGCTGGCGGATGGGGCTGAAAATCCTGCCACCAGATATAAACTCCGGTCTTTACCATTTCCACGTCAACGACGACGGCGAAATCGTGTATGGTATTGGCGCGATCAAAGGGGTCGGTGAAGGTCCGATTGAGGCCATCATCGAAGCCCGTAATAAAGGCGGCTACTTCCGCGA  >  NC_000913/207238‑207707

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