Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,223,934 A→C 18.9% intergenic (+36/+2) dld → / ← pbpG D‑lactate dehydrogenase, FAD‑binding, NADH independent/D‑alanyl‑D‑alanine endopeptidase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,223,9340AC18.9% 75.5 / 2.4 49intergenic (+36/+2)dld/pbpGD‑lactate dehydrogenase, FAD‑binding, NADH independent/D‑alanyl‑D‑alanine endopeptidase
Reads supporting (aligned to +/- strand):  ref base A (18/19);  new base C (2/8);  total (20/29)
Fisher's exact test for biased strand distribution p-value = 1.54e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GCGTTAAAAGAACAGATGCTGGAACTGCTACAGCAGCGCGGCGCGCAGTACCCTGCCGAGCATAACGTCGGTCATTTGTATAAAGCACCGGAGACGTTGCAGAAGTTCTATCGCGAGAACGATCCGACCAACAGCATGAATCCGGGGATCGGTAAAACCAGTAAACGGAAAAACTGGCAGGAAGTGGAGTAAAAATTACGGATGGCAGAGTATCGCCATCCGAATTCACTTAATCGTTCTGTGCCGTCTGCCCCGCCGCCGCCATTTGGGCGGCTTTTTGTTTTTTATAGCTCAACGCTGCTGCCGGCACAGGCATCACTTTACCGGTTTCAATCCAGGTACGCAGGCGGCTGGCATCGGCAAAATGGGTATATTTGCCAAACGCGTCCATCACTACCAGCGCCACCGGTTTATTATTGATAACCGTACGCATCACCAGACAATGGCCCGC  >  NC_000913/2223707‑2224157
                                                                                                                                                                                                                                      |                                                                                                                                                                                                                                
gCGTTAAAAGAACAGATGCTGGAACTGCTACAGCAGCTCGGCGCGCAGTACCCTGCCGAGCATAACGTCGGTCATTTGTATAAAGCCACGGAGACGTTGCAGAAGTTCTATCGCGAGAACGATCCGACCAACAGCATGAATCCGGGGATCGGTAAAACCAGTAAACGGAAAAACTGGCAGGAAGTGGAGTAAAAATTACGGATGGCAGAGTATCGCCATCCGAATTCa                                                                                                                                                                                                                                  <  1:532819/228‑1 (MQ=255)
  gTTAAAAGAACAGATGCTGGAACTGCTACAGCAGCGCGGCGCGCAGTACCCTGCCGAGCATAACGTCGGTCATTTGTATAAAGCACCGGAGACGTTGCAGAAGTTCTATCGCGAGAACGATCCGACCAACAGCATGAATCCGGGGATCGGTAAAACCAGTAAACGGAAAAACTGGCAGGAAGTGGAGTAAAAATTACGGATGGCAGGGTATCGCCATCCGAATTCACt                                                                                                                                                                                                                                >  2:415382/1‑228 (MQ=255)
            cAGATGCTGGAACTGCTACAGCAGCGCGGCGCGCAGTACCCTGCCGAGCATAACGTCGGTCATTTGTATAAAGCACCGGAGACGTTGCAGAAGTTCTATCGCGAGAACGATCCGACCAACAGCATGAATCCGGGGATCGGTAAAACCAGTAAACGGAAAAACTGGCAGGAAGTGGAGTAAAAATTACGGACGGCAGAGTATAGCCATCCGAATTAACTTAATCGttt                                                                                                                                                                                                                       >  2:377899/1‑226 (MQ=255)
             aggtgcgtgaccTGCCAAAGCAGCGCGGCGCGCAGTTCCCGGCCGAGCATAACGTCGGTCATTTGTATAAAGCACCGGAGACGTTGCAGAAGTTCTATCGCGAGAACGATCCGACCAACAGCATGAATCCGGGGATCGGTAAAACCAGTAAACGGAAAAACTGGCAGGAAGTGGAGTAAAAATTACGGATGGCAGAGTATCGCCATCCGAATTCACTTAATCGTTCt                                                                                                                                                                                                                      <  2:83864/216‑1 (MQ=255)
               aTGCTGGAACTGATACAGCAGCGCGGCGCGCAGTTCCCTGCCGAGCATAACGTCGGTCATTTGTATAAAGCCCCGGAGACGTTGCAGAAGTTCTATCGCGAGAACGATCCGACCAACAGCATGAATCCGGGGATCGGTAAAACCAGTAAACGGAAAAACTGGCAGGAAGTGGAGTAAAAATTACGGATGGCAGAGTATCGCCATCCGAATTCACTTAATCGTTCtgtg                                                                                                                                                                                                                   <  2:485205/228‑1 (MQ=255)
                    ggAACTGCTACAGCAGCGCGGCGCGCAGTACCCTGCCGAGCATAACGTCGGTCATTTGTATAAAGCACCGGAGACGTTGCAGAAGTTCTATCGCGAGAACGATCCGACCAACAGCATGAATCCGGGGATCGGTAAAACCAGTAAACGGAAAAACTGGCAGGAAGTGGAGTAAAAATTACGGATGGCAGAGTATCGCCATCCGAATTCACTTAATCGTTATGTGCCGt                                                                                                                                                                                                               >  2:533891/1‑227 (MQ=255)
                         aagtgtaaaaaaagacagcccAGTACCCGGCCGATCATAACGTGCGTCATTGTGAAAAAGCCCCGGAGACGTTGGAGAAGTTCTAGCGAGACAACGATCCGACCAACAGCAAGAATCCGGGGATCGGTAAAACCAGTAAACGGAAAAACTGGCAGGAAGTGGAGTAAAAATTACGGATGGCAGAGTCTCGCCATCCGAATTCACTTAATCGTTCTGTGCCGTCGGcc                                                                                                                                                                                                          <  1:27997/207‑1 (MQ=255)
                            taaagcaggcgggcgcgCCGTACCCTGCCGAGCATAACGGCGGTCATTTGTATAAAGCACCGGAGACGTTGCAGAAGTTCTATCGCGAGAACGATCCGACCAAAAGCATGAATCCGGGGATCGGTAAAACCAGTAAACGGAAAAACTGGCAGGAAGTGGAGTAAAAATTACGGATGGCAGAGTATCGCCATCCGAATTCACTTAATCGTTCTGTGCCGTCTGCCccgc                                                                                                                                                                                                      <  2:479143/217‑1 (MQ=255)
                                 cagcGGGGCGCGCAGGACCCTGCCGAGCATAACGTCGGTCATTTGTATAAAGCACCGGAGACGGTGCAGAAGTTCTATCGCAAGAACGATCCGACCAACAGCAGGAATCCGGGGATCGGTAAAACCAGTAAACGGAAAAACTGGCAGGAAGTGGAGTAAAAATTACGGATGGCAGAGTATCGCCATCCGAATTCACTTAATCGTTCTGTGCCGTCTGCCccgccgcc                                                                                                                                                                                                  <  2:289417/227‑1 (MQ=255)
                                     gcggcgcgCAGTACCCTGCCGAGCATAACGTCGGTCATTTGTATAAAGCACCGGAGACGTTGCAGAAGTTCTATCGCGAGAACGATCCGACCAACAGCATGAATCCGGGGATCGGTAAAACCAGTAAACGGAAAAACTGGCAGGAAGTGGAGTAAAAATTACGGATGGCAGAGTATCGCCATCCGAATTCACTTATTCGTTCTGTGCCGTCTGCACCGCCGCGGCCAt                                                                                                                                                                                             >  2:94513/1‑228 (MQ=255)
                                                      gCCGAGCATAACGTCGGTCATTTGTATAAAGCACCGGAGACGTTGCAGAAGTTCTATCGCGAGAACGATCCGACCAACAGCATGAATCCGGGGATCGGTAAAACCAGTAAACGGAAAAACTGGCAGGAAGTGGAGTAAAAATTACGGATGGCAGAGAATCGCCATCCGAATTCACTTAATCGTTCTGTGCCGTTTGCCCCgccgccgccctgcctcttatacccatct                                                                                                                                                                            >  2:27997/1‑209 (MQ=255)
                                                        cGAGCATAACGTCGGTCATTTGTATAAAGCACCGGAGACGTTGCAGAAGTTCTATCGCGAGAACGATCCGACCAACAGCATGAATCCGGGGATCGGTAAACCAGTAAACGGAAAAACTGGCAGGAAGTGGAGTAAAAATTACGGATGGCAGAGTATCGCCATCCGAATTCACTTAATCGTTCTGTGCCGTCTGCCCCCCCGCCGCCATTTGGGCGGCTTTTTGtttt                                                                                                                                                                          >  2:61024/1‑227 (MQ=255)
                                                          aGCATAACGTCGGTCATTTGTATAAAGCACCGGAGACGTTGCAGAAGTTCTATCGCGAGAACGATCCGACCAACAGCATGAATCCGGGGATCGGTAAAACCAGTAAACGGAAAAACTGGCAGGAAGTGGAGTAAAAATTACGGATGGCAGAGTATCGCCATCCGAATTCACTTAATCGTTCTGTGCCGTTTGCCCCACCGCCGCCATTTGGGCGGctttttcttttt                                                                                                                                                                         >  2:599653/1‑227 (MQ=255)
                                                          aGCATAACGTCGGTCATTTGTATAAAGCACCGGAGACGTTGCAGAAGTTCTATCGCGAGAACGATCCGACCAACAGCATGAATCCGGGGATCGGTAAAACCAGTAAACGGAAAAACTGGCAGGAAGTGGAGTAAAAATTACGGATGGCAGAGAATCGCCATCCGAATTCACTTAATCGTTCTGTGCCGTCCGCCCCGCCGCCGCCATTTGGGGCGCTtttttttttt                                                                                                                                                                         >  1:377438/1‑227 (MQ=255)
                                                            cATAACGTCGGTCATTTGTATAAAGCACCGGAGACGTTGCAGAAGTTCTATCGCGAGAACGATCCGACCAACAGCATGAATCCGGGGATCGGTAAAACCAGAAAACGGAAAACCAGGCAGGAAGGGGAATATAAATTACGGAAGGCAGCGGAACGCCAACCGGATTCACTTAATAGTTCTGTGCCGTCTGCCCCGCCCCCGCCATTTGGGCGGCTTTTTGttttt                                                                                                                                                                         >  1:699751/1‑225 (MQ=255)
                                                                                  aaGCACCGGGGCCGTTGCAGAAGTTCTATCGCGAGAACGATCCGACCAACAGCATGAATCCGGGGATCGGTAAAACCAGTAAACGGAAAAACTGGCAGGAAGTGGAGTAAAAATTACGGATGGCAGAGTATCGCCATCCGAATTCACTTAATCGTTCTGTGCCGTCTGCCCCGCCGCCGCCATTTGGGCGGCTTTTTGTTTTTTATAGCTCAACGCTGCTGCCGGca                                                                                                                                                 <  2:530714/227‑1 (MQ=255)
                                                                                       acGGAGACGTTGCAGGAGTTCGATCGCGAGAACGATCCGACCAACAGCATGAATCCGGGGATCGGTAAAACCAGTAAACGGAAAAACTGGCAGGAAGTGGAGTAAAAATTACGGATGGCAGAGTATCGCCATCCGAATTCACTTAATCGTTCTGTGCCGTCTGCCCCGCCGCCGCCATTTGGGCGGCTTTTTGTTTTTTATAGCTCAACGCTGCTGCCGGCACAGGCa                                                                                                                                           <  2:401587/227‑1 (MQ=255)
                                                                                        cGGAGACGTTGCAGAAGTTCTATCGCGAGAACGATCCGACCAACAGCATGAATCCGGGGATCGGTAAAACCAGTAAACGGAAAAACTGGCAGGAAGTGGAGTAAAAATTACGGATGGCAGAGTATCGCCATCCGAATTCCCTTAATCGTTCTGTGCCGCCGTCCCCGCCCCCGCCCTTTGGGCGGCTTTTTTTTTTTTAATGCTCacacccgctaccgggccaggcct                                                                                                                                          >  1:532211/1‑206 (MQ=255)
                                                                                             aCGTTGCAGAAGTTCTATCGCGAGAACGATCCGACCAACAGCATGAATCCGGGGATCGGTAAAACCAGTAAACGGAAAAACTGGCAGGAAGTGGAGTAAAAATTACGGATGGCAGAGTATCGCCATCCGAATTACCTTAATCGTTCTGTGCCGTCTGCCCCGCCGCCGCCATTTGGGCGGCTTTTTGTTTTTTTTAGCTCAAACCTGCTGCCGGCACAGGCAtacctt                                                                                                                                     >  2:354843/1‑223 (MQ=255)
                                                                                                 tGCAGAAGTTCTACAGCGAAAAAGATCCGACCAACAGCATGAATCCGGGGATCGGTAAAACCAGTAAACGGAAAAACTGTCAGGAAGTGGAGTAAAAATTCCGGATGGCAGAGTATCGCCATCCGAATTCACTTAATCTTTCTGTGCCGTCTGCCCCGCCGCCGCCATTTGGGCGGCTTTTTGTTTTTTATAGCTCAACGCTGCTGCCGGCACAGGCATCACTTTAcc                                                                                                                                 <  1:430476/228‑1 (MQ=255)
                                                                                                 tGCAGAAGTGCTATCGCGAGAACGATCCGGCCAACAGCATGAATCCGGTGATCGGTAAAACCAGTAAACGGAAAAACTGGCAGGAAGTGGAGTAAAAATTACGGATGGCAGAGTATCGCCCTCCGAATTCACTTAATCGTTCTGTGCCGTCTGCCCCGCCGCCGCCATTTGGGCGGCTTTTTGTTTTTTATAGCTCAACGCTGCTGCCGGCACAGGCATCACTTTAc                                                                                                                                  <  2:679290/227‑1 (MQ=255)
                                                                                                                                aaaaaGCACGAAACCGGGGACAGGAAAACCCAGTAAACGAAAAACCTGGCAGGAAGTGGTGTAAAAATTGCGGATGGCAGAGTATCGCCATCCTAATCCCCTTAATCGTTCTGTTCCGTCTGCCCCGCCGCCGCCATTTGGGCGGCTTTTTGGTTTTTATAGCTCAACGCTGCTGCCGGCACAGGCATCACTTTACCGGTTTCAATCCAGGTACGCAGGCGGCTGGCa                                                                                                  <  1:415382/224‑1 (MQ=255)
                                                                                                                                       tggagtcgggggCTCGGTAACCAAGGAAACGGAAAAACGGGCAGATAGTGGAGAAAAAATACGGGGTGGCAGACTTTCGCCACCCGATTTCCCTTAATCGTTCTGGGCCGTCTGCCCCGCCGCCGCCATTTGGGCGGCTTTTTGTTTTTTATAGCTCAACGCTGCTGCCGGCACAGGCATCACTTTACCGGTTTCAATCCAGGTACGCAGGCGGCTGGCATCGGCaaa                                                                                          <  1:61024/220‑1 (MQ=255)
                                                                                                                                        tGAATCCGGGGATCGGTAAAACCAGTAAACGGAAAAACTGGCAGGAAGTGGAGTAAAAATTACGGATGGCAGAGTATCGCCATCCGAATTCACTTAATCGTTCTGTGCCGTCTGCCCCGCCGCCGCCATTTGGGCGGCTTTTTGTTTTTTATAGCTCAACGCTGCTGCCGGGCCAGGCATCACTTTACCGGTTTCAATCCAGGTACCCCGGCGCCTGGCacaggaaaa                                                                                          >  1:172319/1‑220 (MQ=255)
                                                                                                                                        tGAATCCGGGGATCGGTAAAACCAGTAAACGGAAAAACTGGCAGGAAGTGGAGAAAAACTTACGGATGACAGAGTATAGCCAACCGAATTCACTTAATCGTTTGCGGCCGGCTGCCCCGCCGCCGCCATTTGGGCTGCTTTTTTTTTTTTCTTGTTCACCGCCGACGCCGGCACAGTGTTCACTTTTCCGGTTTTCATCTAGCTCCGCAAGCGCGTGGcccccgcaac                                                                                          >  1:186556/1‑219 (MQ=255)
                                                                                                                                         aggaccGGGGGTCCGGAAAACCCAGTAAACAGAAAACCTGGAAGGAATGTGATCAAAAATTGCGGATGGCGTAGCATCGCCCTCCGAATGCACTTAATCGTTCTGTGCCGTCGGCCCCGCCGCCGCCATTTGGGCGGCTTTTTGTTTTTTATAGCTCAACGCTGCTGCCGGCACAGGCATCACTTTACCGGTTTCAATCCAGGTACGCAGGCGGCTTGCATCGGCaaa                                                                                          <  1:239060/224‑1 (MQ=255)
                                                                                                                                                gggATCGGTAAAACCAGTAAACGGAAAAACTGGCAGGAAGTGGAGTAAAAATTACGGATGGCAGAGTATCGCCATCCGAATTCACTTAATCGTTCTGTGCCGTCTGCCCCGCCGCCGCCATTTGGGCGGCTTTTTGTTTTTTATAGCTCAACGCTGCTGCCGGCCCAGGCATCCCTTTACCGGTTTCAATCCAGGGACGCAGGCGGCTGGCATCTGTAAAATGGGcat                                                                                  >  2:239060/1‑228 (MQ=255)
                                                                                                                                                    aagggTAAAAACCGGAAACAGGAAAAACTGGAAGGAAGTGGAGTAAAAATTCCGGATGGCGAGGTTACCCCACCCGAATTCACTTATTCGTTCGTGGCCGTTTGCCCCGCCGCCGCCATTTGGGCGGCTTTTTGTTTTTTATAGCTCAACGCTGCTGCCGGCACAGGCATCACTTTACCGGTTTCAATCCAGGTACGCAGGCGGCTGGCATCGGCAAAATGGGtata                                                                                 <  1:199247/224‑1 (MQ=255)
                                                                                                                                                        gTAAAACCAGTAAACGGAAAAACTGGCAGGAAGTGGAGTAAAAATTTCGGATGGCAGAGTATCGCCATCCGAATTCACTTAACCGTTCTGTGCCGTCTGCCCCGCCGCCGCCATTTGGGCGGCTTTTTGTTTTTTATAGCTCAACGCTGCTGCCGGCACAGGCATCACTTTACCGGTTTCAATCCAgg                                                                                                                   <  2:620127/188‑1 (MQ=255)
                                                                                                                                                              ccAGTAAACGGAAAAACTGGCAGGAAGTGGAGTAAAAATTACGGAAGGCAGGGTATCGCCATCCGAATTCACTTAATCGTTCTGGGCCGGCTGCCCCGCCGCCGCCATTTGGGCGGCTTTTTGTTTTTTATAGTTCCACGCCGCGGCCGCCCCAGGCATCCCTTTGCCGTTTTTAATCCAGGTACGCAGGCGGCCTGGACCGGCAAAATGGGGATATTTTCCAAAccc                                                                     >  1:308528/1‑226 (MQ=255)
                                                                                                                                                                gtgaaaGGGAAAAACTGGCAGGAGGTGGAGTAAAAATTACGGATGGCAGGGTCACGCCATCCGATTTCTCTTATTCGTTTTTGGCGGTCTGCCCCGCCGCCGCCATTTGGGCGGCTTTTTGTTTGTTATAGCTCAACGCTGCTGCCGGCACAGGCATCACTTTACCGGTTTCAATCCAGGTACGCAGGCGGCTGGCATCGGCAAAATGGGTATATTTGCCAAAcgcg                                                                    <  1:591919/224‑1 (MQ=255)
                                                                                                                                                                ggaaaaCGGAAAACTGGGAGGGAAGGGGGTTAAAATTTCCGGAGGGCAGCGTAGCGCCAGCCGATTTCCCTGATCTGTTCTGTGCCGGCTGCCCCGCCGCCGCCATTTGGGCGGCTTTTTGTTTTTTATAGCTCAACGCTGCTGCCGGCACAGGCATCACTTTACCGGTTTCAATCCAGGTACGCAGGCGGCTGGCATCGGCAAAATGGGTATATTTGCCAAACGCGt                                                                   <  1:354843/225‑1 (MQ=255)
                                                                                                                                                                 gTAAACGGAAAAACTGGCAGGAAGTGGAGTAAAAATTACGGATGGCAGAGTATCGCCATCCGAATTCACTTAATCGTTCTGTGCCGTCTGCCCCGCCGCCGCCATTTGGGCGGCTTTTTGTTTTTTATAGCTCAACGCTGCTGCCGGCACAGGCATCACTTTACCGGTTTCAATCCAGGTACGCAGGCGGCTGGCATCGGCAAAATGGGTATATTTGCCAAACGCGt                                                                   <  2:452172/227‑1 (MQ=255)
                                                                                                                                                                 gTAAACGGAAAAACTGGCAGGAAGTGGAGTAAAAATTACGGATGGCAGAGTATCGCCATACGAATTCACTTAATCGTTCTGTGCCGTCTGCCCCGCCGCCGCCATTTGGGCGGCTTTTTGTTTTTTATAGCTCAACGCTGCTGCCGGCACAGGCATCACTTTACCGGTTTCATTCCAGGTACGCAGGGGGCTGGCATCGGCAAAATGGGTATATTTGCAAACCGCGt                                                                   >  2:289964/1‑227 (MQ=255)
                                                                                                                                                                  tAAACGGAAAAACTGGCAGGAAGTGGAGTAAAAATTACGGATGGCAGAGTATCGCCATCCGAATTCACTTAATCGTTCTGTGCCGTCTGCCCCGCCGCCGCCATTTGGGCGGCTTTTTGTTTTTTATAGCTCAACGCTGCTGCCGGCCCAGGCATCACTTTACCGGTTTCAATCCAGGTACGCAGGCGGCTGGCATAGGCAAAATggggac                                                                                  >  2:199247/1‑208 (MQ=255)
                                                                                                                                                                     aCGGAAAACCTGGCAGGAAGTGGAGTAAAAATTACGGATGGCAGAGTATCGCCATCCGAATTCCCTTAATCGTTCTGTGCCGTCTGCCCCGCCGCCGCCATTTGGGCGGCTTTTTGTTTTTTATAGCTCAACGCTGCTGCCGGCACAGGCATCACTTTACCGGTTTCAATCCAGGTACGCAGGCGGCTGGCATCGGCAAAATGGGTATATTTGCCAAACGCGTCCAt                                                               <  2:308528/227‑1 (MQ=255)
                                                                                                                                                                           aaaaCTGGCAGGAAGTGGAGTAAAAATTACGGATGGCAGAGTATCGCCATCCGAATTCACTTAATCGTTCTGTGCCGTCTGCCCCGCCGCCGCCATTTGGGCGGCTTTTTGTTTTTTATAGCTCAACGCTGCTGCCGGCACAGGCATCACTTTACCGGTTTCAATCCAGGTACGCAGGCGGCTGGCATCGGCAAAATGGGTATATTTGCCAAACGCGTCCATCACTAc                                                         <  2:377438/228‑1 (MQ=255)
                                                                                                                                                                                     cggaatgtgggTAAAAATTACGGAAGCACGAGTACCCCCTCCCAAATTCCTTTACTCTTTTTGGGCCTTCTGCCCCGCCGCCGCCATTTGGGCGGCTTTTTGTTTTTTATAGCTCAACGCTGCTGCCGGCACAGGCATCACTTTACCGGTTTCAATCCAGGTACGCAGGCGGCTGGCATCGGCAAAATGGGTATATTTGCCAAACGCGTCCATCACTACCAGCGCCAc                                               <  1:674019/218‑1 (MQ=255)
                                                                                                                                                                                       gaatgggAGTAAAAAGTCCGGATGGCAGAGTATCGCCATCCGAATTCACTTAATCGTTCTGTGCCGTCTGCCCCGCCGCCGCCATTTGGGCGGCTTTTTGTTTTTTATAGCTCAACGCTGCTGCCGGCACAGGCATCACTTTACCGGTTTCAATCCAGGTACGCAGGCGGCTGGCATCGGCAAAATGGGTATATTTGCCAAACGCGTCCATCACTACCAGCGCCACCg                                             <  2:35934/223‑1 (MQ=255)
                                                                                                                                                                                         aggggATTAAAAATTACGGATGGCAGGGGTTCGCCATCCGAATTCCCTTAATCGTTCTGTGCCGTCTGCCCCGCCGCCGCCATTTGGGCGGCTTTTTGTTTTTTATAGCTCAACGCTGCTGCCGGCACAGGCATCACTTTACCGGTTTCAATCCAGGTACGCAGGCGGCTGGCATCGGCAAAATGGGTATATTTGCCAAACGCGTCCATCACTACCAGCGCCACCGGt                                           <  2:532211/225‑1 (MQ=255)
                                                                                                                                                                                         aGTGCATTAAAAATTACGGGTGGGAAAGTATCGGCATCGGAATTTCATTAATCGTTCTGTTCCGTCTGCCACGGCGCCGCCATTTGGGCGGCTTTTTGTTTTTTATAGCTCAACGCTGCTGCCGGCACAGGCATCACTTTACCGGTTTCAATCCAGGTACGCAGGCGGCTGGCATCGGCAAAATGGGTATATTTGCCAAACGCGTCCATCACTACCAGCGCCACCGGt                                           <  2:234595/228‑1 (MQ=255)
                                                                                                                                                                                            ggcagaaaaaTTGCGGATGGCAGTGTAGCGCCAACCGAATTCACTTAATCGTTCTGTGCCGTCGGCCCCGCCGCCGCCATCTGGGCGGCTTTTTGTTTTTTATAGCTCAACGCTGCTGCCGGCACAGGCATCACTTTACCGGTTTCAATCCAGGTACGCAGGCGGCTGGCATCGGCAAAATGGGTATATTTGCCAAACGCGTCCATCACTACCAGCGCCACCGGTtta                                        <  2:1404/223‑1 (MQ=255)
                                                                                                                                                                                             gAGTAAAAATTACGAATGGCAGAGTATCGCCATCCGAATTCACTTAATCGTTCTGTGCCGGCTGCCCCGCCGCCGCCATTTGGGCGGCTTTTTGTTTTTTATAGCTCACCGCTGCTGCCGGCCCGGGCATCACTTTTCATGTTTCTATCCACGTACGCAGCCGTCTGGCATCGGCAAAATGGGTATATTTTCCACACCCGGCCATCACTCCCATCCCACACCGGttaa                                        >  1:332343/1‑226 (MQ=255)
                                                                                                                                                                                                  aaaaTTACGGATGGCAGAGTATCGCCATCCGAATTCACTTAATCGTTCTGTGCCGTCTGCCCCGCCGCCGCCATTTGGGCGGCTTTTTGTTTTTTATAGCTCAACGCTGCTGCCGGCACAGGCATCACTTTACCGGTTTCAATCCAGGTACGCCGGCGGCTGGCATCGGCAAACTGGGTATCTTTGCCAAACCCGTCCATACCTACCCGCGCCACCGGTTTattatt                                   >  2:222017/1‑227 (MQ=255)
                                                                                                                                                                                                    aaTTGCGGATGGCAGAGTATCCCCATCCGAATTCACTTAATCGTTCTGTGCCGTCTGCCCCGCCGCCGCCATTTGGGCGGCTTTTTGTTTTTTATAGCTCAACGCTGCTGCCGGCACAGGCATCACTTTACCGGTTTCAATCCAGGTACGCAGGCGGCTGGCATCGGCAAAATGGGTATATTTGCCAAACGCGTCCATCACTACCAGCGCCACCGGTTTATTATTGa                                 <  2:324112/227‑1 (MQ=255)
                                                                                                                                                                                                                cggattgcccccctCCGATTTCCCTTATTCGTTTTTTTCCGTCTGCCCCGCCGCCGCCATTTGGGCGGCTTTTTGTTGTTTATAGCTCAACGCTGGCGCCGGCACAGGCATAACTTTACCGGTTTCAATCCAGGTACGCAGGCGGCTAGCATCGGCAAAATGGGTATATTTGCCAAACGCGTCCATCACTACCAGCGCCACCGGTTTATTATTGATAACCGTACGCAt                    <  1:533891/215‑1 (MQ=255)
                                                                                                                                                                                                                cAGAGTCTCGCCATCCGATTTCACTTAATCGTTCTGGCCCGGCTGCCCCGCCGCCGCCATTTTGGCGGCTTTTTGTTTTTTATAGCTCAACGCTGCTGCCGGCACAGGCATCACTTTACCGGTTTCAATCCAGGTACGCAGGCGGCTGGCATCGGCAAAATGGGTATATTTGCCAAACGCGTCCATCACTACCAGCGCCACCGGTTTATTATTGATAACCGTACGCAt                    <  2:188569/228‑1 (MQ=255)
                                                                                                                                                                                                                              ccgttttCAGTTAATAGTTCTTCGCCGTCTGCCCCCCCCCCGCCCTTTGGGCGGCTTTTTTTTTTTTAAAGCTCACCGCTGCTGGCGGCACAGGCATCACTTTACCGGTTTCaagc                                                                                                                       <  2:615022/111‑3 (MQ=255)
                                                                                                                                                                                                                               cGAATTCACTTAATCGTTCTGTGCCGTCTGCCCCGCCGCCGCCATTTGGGCGGCTTTTTGTTTTTTATAGCTCAACGCTGCTGCCGGCACAGGCATCACTTTACCGGTTTCAATCCAGGTACGCAGGCGGCTGGCATCGGCAAAATGGGGATATTTGCCAAACGCGTCCATCACTCCCAGCGCCCCCGGGTTATTATTGATAACCGTACGCATCACCAGACaaggtcc     >  2:670284/1‑223 (MQ=255)
                                                                                                                                                                                                                                  aTTCTCTTAAGCGTTCTTGCCCGCCTGCCCCCCCGTCGCCATTTGGGCGGTCTTTTGTTTTGTATAGCTCAACGGTGCTGCCGGCACAGGCATCACTTTACCGGTTGCAATCCAGGTACGCAGGCGGCTGGCATCGGCAAAATGGGTATATTTGCCAAACGCGTCCATCACTACCAGCGCCACCGGTTTATTATTGATAACCGTACGCATCACCAGACAATGGCccgc  <  2:618194/228‑1 (MQ=255)
                                                                                                                                                                                                                                      |                                                                                                                                                                                                                                
GCGTTAAAAGAACAGATGCTGGAACTGCTACAGCAGCGCGGCGCGCAGTACCCTGCCGAGCATAACGTCGGTCATTTGTATAAAGCACCGGAGACGTTGCAGAAGTTCTATCGCGAGAACGATCCGACCAACAGCATGAATCCGGGGATCGGTAAAACCAGTAAACGGAAAAACTGGCAGGAAGTGGAGTAAAAATTACGGATGGCAGAGTATCGCCATCCGAATTCACTTAATCGTTCTGTGCCGTCTGCCCCGCCGCCGCCATTTGGGCGGCTTTTTGTTTTTTATAGCTCAACGCTGCTGCCGGCACAGGCATCACTTTACCGGTTTCAATCCAGGTACGCAGGCGGCTGGCATCGGCAAAATGGGTATATTTGCCAAACGCGTCCATCACTACCAGCGCCACCGGTTTATTATTGATAACCGTACGCATCACCAGACAATGGCCCGC  >  NC_000913/2223707‑2224157

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