Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,560,071 Δ2 bp pseudogene (382‑383/606 nt) glpR ← pseudogene, DNA‑binding transcriptional repressor;regulator; Energy metabolism, carbon: Anaerobic respiration; repressor of the glp operon

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,560,0710C.96.4% 234.6 / ‑1.7 55pseudogene (383/606 nt)glpRpseudogene, DNA‑binding transcriptional repressor;regulator; Energy metabolism, carbon: Anaerobic respiration; repressor of the glp operon
Reads supporting (aligned to +/- strand):  ref base C (2/0);  new base . (28/25);  total (30/25)
Fisher's exact test for biased strand distribution p-value = 4.95e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
*NC_0009133,560,0720C.96.4% 234.7 / 0.3 55pseudogene (382/606 nt)glpRpseudogene, DNA‑binding transcriptional repressor;regulator; Energy metabolism, carbon: Anaerobic respiration; repressor of the glp operon
Reads supporting (aligned to +/- strand):  ref base C (0/0);  major base . (28/25);  minor base G (2/0);  total (30/25)
Fisher's exact test for biased strand distribution p-value = 4.95e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CGTCAATGACGCGGATCGGCTACGCAAATACCTGGCGCGTTTTGGTCTGACGTGGGAAGCCGTGCAGGATCAGCACAGCTCCAGTTGAATATGGTGGTCCGTCAGCACCTGCATCACGCTTACTGGCGGCGGGGCGTCGGTGTAGACGGCATCTACCATGCTGATGCTGCCCATATTGACCATCGCGTTACGGCCAAATTTCGAGTGATCGACAACCAGCATAACGTGGCGCGAGTTCTCAATAATGGCGCGTTTGGTGCGAACTTCGTGGTAATCGAACTCCAGCAGCGAGCCGTCGCTATCGATGCCGCTTATCCCCAGAATGCCGAAATCAAGGCGGAACTGGGAGATAAAATCGAGCGTCGCTTCGCCAATGATCCCGCCATCGCGGCTGCGTAATTCGCCACCGGCGAGAATGATGCGAAAATCTTCTTTTACCATCAACGTGTTAGCAACGTTGAGATTGTTGGTGACAATGCGCAAATTGCTGTGATTGAGCAGTGCGTGCGCTACCGCTTCCGGCGTGGTGCCGATATCGATAAACAGCGTCGAGCCATTGGGGATTTGCTCCGCCA  >  NC_000913/3559779‑3560353
                                                                                                                                                                                                                                                                                                    ||                                                                                                                                                                                                                                                                                         
cGTCAATGACGCGGATCGGCTACGCAAATACCTGGCGCGTTTTGGTCTGACGTGGGAAGCCGTGCAGGATCAGCACAGCTCCAGTTGAATATGGTGGTCCGTCAGCACCTGCATCACGCTTACTGGCGGCGGGGCGTCGGTGTAGACGGCATCTACCATGCTGATGCTGCCCATATTGACCATCGCGTTACGGCCAAATTTCGAGTGATCGACAACCAGCATAACGTGGCGCGAGTTCTCAATAATGGCGCGTTTGGTGCGAACTTCGTGGTAATCGAACTCCAGCAGCGAG‑‑GTCGCTAt                                                                                                                                                                                                                                                                                   <  2:209132/300‑1 (MQ=255)
   cAATGACGCGGATCGGCTACGCAAATACCTGGCGCGTTTTGGTCTGACGTGGGAAGCCGTGCAGGATCAGCACAGCTCCAGTTGAATATGGTGGTCCGTCAGCACCTGCATCACGCTTACTGGCGGCGGGGCGTCGGTGTAGACGGCATCTACCATGCTGATGCTGCCCATATTGACCATCGCGTTACGGCCAAATTTCGAGTGATCGACAACCAGCATAACGTGGCGCGAGTTCTCAATAATGGCGCGTTTGGTGCGAACTTCGTGGTAATCGAACTCCAGCAGCGAG‑‑GTCGCTATCGAt                                                                                                                                                                                                                                                                               >  2:211787/1‑301 (MQ=255)
     aTGACGCGGATCGGCTACGCAAATACCTGGCGCGTTTTGGTCTGACGTGGGAAGCCGTGCAGGATCAGCACAGCTCCAGTTGAATATGGTGGTCCGTCAGCACCTGCATCACGCTTACTGGCGGCGGGGCGTCGGTGTAGACGGCATCTACCATGCTGATGCTGCCCATATTGACCATCGCGTTACGGCCAAATTTCGAGTGATCGACAACCAGCATAACGTGGCGCGAGTTCTCAATAATGGCGCGTTTGGTGCGAACTTCGTGGTAATCGAACTCCAGCAGCGAG‑‑GTCGCTATCGATg                                                                                                                                                                                                                                                                              >  1:428405/1‑300 (MQ=255)
       gACGCGGATCGGCTACGCAAATACCTGGCGCGTTTTGGTCTGACGTGGGAAGCCGTGCAGGATCAGCACAGCTCCAGTTGAATATGGTGGTCCGTCAGCACCTGCATCACGCTTACTGGCGGCGGGGCGTCGGTGTAGACGGCATCTACCATGCTGATGCTGCCCATATTGACCCTCGCGTTACGGCCAAATTTCGAGTGATCGACAACCAGCATAACGTGGCGCGAGTTCTCAATAATGGCGCGTTTGGTTCGAACTTCGTGGTAATCGACCTCCAGC‑‑CGCGCGGTCGCTATCgcttcc                                                                                                                                                                                                                                                                            >  1:95107/1‑295 (MQ=255)
        cgggggATCTGCGAAGCAAATACCGAGCGGGGGTCGGGCTGACGGTGGACGGCGGACAGGATAAGCACAGCGCCAGTTGGATCTGGTGGTCCGTCAAGACCTGCATCACGCTTACTGGCGCCGGGGCGTCGGTGTAGACGGCATCTACCATGCTGATGCTGCCCATATTGACCATCGCGTTACGGCCAGATTTCGAGTGATCGACAACGAGCATAACGTGGCGCGAGTTCTCAATAATGGCGCGTTTGGTGCGAACGTCGTGGTAATCGAACTCCAGCAGCGAG‑‑GTCGCTATCGATGCCGc                                                                                                                                                                                                                                                                          <  1:429495/297‑1 (MQ=255)
        aCGCGGATCGGCTACGCAAATACCTGGCGCGTTTTGGTCTGACGTGGGAAGCCGTGCAGGATCAGCACAGCTCCAGTTGAATATGGTGGTCCGTCAGCACCTGCATCACGCTTACTGGCGGCGGGGCGTCGGTGTAGACGGCATCTACCATGCTGATGCTGCCCATATTGACCATCGCGTTACGGCCAAATTTCGAGTGATCGATAACCAGCATAACGTGGCGCGAGTTCTCAATAATGGCGCGTTTGGTGCGAACTTCGTGGTAATCGAACTCCACCAGCGAG‑‑GTCGCTATCGATGCCg                                                                                                                                                                                                                                                                           >  1:145620/1‑300 (MQ=255)
                  atctcGCAAATACCTGGCGCGTTTTGGTCTGACGTGGGAAGCCGTGCAGGATCAGCACAGCTCCAGTTGAATATGGTGGTCCGTCAGCACCTGCATCACGCTTACTGGCGGCGGGGCGTCGGTGTAGACGGCATCTACCATGCTGATGCTGCCCATATTGACCATCGCGTTACGGCCAAATTTCGAGTGATCGACAACCAGCATAACGTGGCGCGAGTTCTCAATAATGGCGCGTTTGGTGCGAACTTCGTGGTAATCGAACTCCAGCAGCGAG‑‑GTCGCTATCGATGCCGCTTATCCCCa                                                                                                                                                                                                                                                                 <  2:371376/296‑1 (MQ=255)
                      cGCAAATACCTGGCGCGTTTTGGTCTGACGTGGGAAGCCGTGCAGGATCAGCACAGCTCCAGTTGAATATGGTGGTCCGTCAGCACCTGCATCACGCTTACTGGCGGCGGGGCGTCGGTGTAGACGGCATCTACCATGCTGATGCTGCCCATATTGACCATCGCGTTACGGCCAAATTTCGAGTGATCGACAACCAGCATAACGTGGCGCGAGTTCTCAATAATGGCGCGTTTGGTGCGAACTTCGTGGTAATCGAACTCCAGCAGCGAG‑‑GCCGCTATCGAGTCCGCTTCTcccccagtc                                                                                                                                                                                                                                                             >  1:370911/1‑295 (MQ=255)
                             aCCTGGCGCGTTTTGGTCTGACGTGGGAAGCCGTGCAGGATCAGCACAGCTCCAGTTGAATATGGTGGTCCGTCAGCACCTGCATCACGCTTACTGGCGGCGGGGCGTCGGTGTAGACGGCATCTACCATGCTGATGCTGCCCATATTGACCATCGCGTTACGGCCAAATTTCGAGTGATCGACAACCAGCATAACGTGGCGCGAGTTCTCAATAATGGCGCGTTTGGTGCGAACTTCGTGGTACTCGAACTCCAGCAGCGAG‑‑GTCGCTATCGATGCCGCTTATCCCCAGAATGCCGaaaa                                                                                                                                                                                                                                                     >  1:209307/1‑300 (MQ=255)
                                           ggTCTGACGTGGGAAGCCGTGCAGGATCAGCACAGCTCCAGTTGAATATGGTGGTCCGTCAGCACCTGCATCACGCTTACTGGCGGCGGGGCGTCGGTGTAGACGGCATCTATCAAGCTGATGCTGCCCATATTGACACTCGCGTTACGGCCAAATCTCGAGTGATCGACAACCAGCATAACGTGTCGAGAGTTCTCACTAATGGAGCGTTTGGTGCGAGTGTGGTGGGACTCGACCTCCAGCAGC‑‑GCGGTCGCTTGAGCATCCGTATACCCCCAGAAAGGCGAAATACAGGCGGaaatcg                                                                                                                                                                                                                                       >  1:175565/1‑297 (MQ=255)
                                             tCTGACGTGGGAAGCCGTGCAGGATCAGCACAGCTCCAGTTGAATATGGTGGTCCGTCAGCACCTGCATCACGCTTACTGGCGGCGGGGCGTCGGTGTAGACGGCATCTACCATGCTGATGCTGCCCATATTGACCATCGCGTTACGGCCAAATTTCGAGTGATCGACAACCAGCATAACGTGGCGCGAGTTCTCAATAATGGCGCGTTTGGTGCGAACTTCGTGGTAATCGAACTCCAGCAGCGAG‑‑GTCGCTATCGATGCCGCTTATCCCCAGAATGCCGAAATCAAGGCGGAACTggg                                                                                                                                                                                                                                      <  2:428952/300‑1 (MQ=255)
                                                        aaGCCGTGCAGGATCAGCACAGCTCCAGTTGAATATGGTGGTCCGTCAGCACCTGCATCACGCTTACTGGCGGCGGGGCGTCGGTGTAGACGGCATCTACCATGCTGATGCTGCCCATATTGACCATCGCGTTACGGCCAAATTTCGAGTGATCGACAACCAGCATAACGTGGCGCGAGTTCTCAATAATGGCGCGTTTGGTCCGAACTTCGTGGTAATCGAAGTCCAGCAGCGAG‑‑GTCGCTATCGATGCCGCTTATCCCCAGAATGCCGAAATAAAGGCGGAACTGGGAGACAAAATCGa                                                                                                                                                                                                                          >  1:466714/1‑301 (MQ=255)
                                                                cAGGATCAGCACAGCTCCAGTTGAATATGGTGGTCCGTCAGCACCTGCATCACGCTTACTGGCGGCGGGGCGTCGGTGTAGACGGCATCTACCATGCTGATGCTGCCCATATTGACCATCGCGTTACGGCCAAATTTCGAGTGATCGACAACCAGCATAACGTGGCGCGAGTTCTCAATAATGGCGCGTTTGGTGCGAACTTCGTGGTAATCGAACTCCAGCAGCGAG‑‑GTCGCTATCGATGCCGCTTATCCCCAGAATGCCGAAATCAAGGCGGAACTGGGAGATAAAATCGAGCGTCGc                                                                                                                                                                                                                   <  1:211532/300‑1 (MQ=255)
                                                                  ggATCAGCACAGCTCCAGTTGAATATGGTGGTCCGTCAGCACCTGCATCACGCTTACTGGCGGCGGGGCGTCGGTGTAGACGGCATCTACCATGCTGATGCTGCCCATATTGACCATCGCGTTACGGCCAAATTTCGAGTGATCGACAACCAGCATAACGTGGCGCGCGTTCTCAATAATGGCGCGTTTGGTGCGAACTTCGTGGTAATCGAACTCCAGCAGCGAG‑‑GTCGCTATCGATGCCGCTTATCCCCAGAATGCCGAAATAAAGGCGGAACTGGGAGATAAAATCGAGTGTCGcat                                                                                                                                                                                                                 >  1:144737/1‑298 (MQ=255)
                                                                             tcttccgatctatatGGTGGTCCGTCAGCACCTGCATCACGCTTACTGGCGGCGGGGCGTCGGTGTAGACGGCATCTACCATGCTGATGCTGCCCATATTGACCATCGCGTTACGGCCAAATTTCGAGTGATCGACAACCAGCATAACGTGGCGCGAGTTCTCAATAATGGCGCGTTTGGTGCGAACTTCGTGGTAATCGAACTCCAGCAGCGAG‑‑GTCGCTATCGATGCCGCTTATCCCCAGAATGCCGAAATCAAGGCGGAACTGGGAGATAAAATCGAGCGTCGTTTCGCCAATGATcc                                                                                                                                                                                                     <  2:118773/290‑1 (MQ=255)
                                                                              cTCCAGTTGAATATGGTGGTCCGTCAGCACCTGCATCACGCTTACTGGCGGCGGGGCGTCGGTGTAGACGGCATCTACCATGCTGATGCTGCCCATATTGACCATCGCGTTACGGCCAAATTTCGAGTGATCGACAACCAGCATAACGTGGCGCGAGTTCTCAATAATGGCGCGTTTGGTGCGAACTTCGTGGTAATCGAACTCCAGCAGCGAG‑‑GTCGCTATCGATGCCGCTTATCCCCAGAATGCCGAAATCAAGGCGGAACTGGGAGATAAAATCGAGCGTCGCTTCGCCAATGATcc                                                                                                                                                                                                     >  2:115566/1‑300 (MQ=255)
                                                                                        atatGGTGGTCCGTCAGCACCTGCATCACGCTTACTGGCGGCGGGGCGTCGGTGTAGACGGCATCTACCATGCTGATGCTGCCCATATTGACCATCGCGTTACGGCCAAATTTCGAGTGATCGACAACCAGCATAACGTGGCGCGAGTTCTCAATAATGGCGCGTTTGGTTCGAACTTCGTGGTAATCGAACTCCAGCAGCGAG‑‑GTCGCTATCGATGCCGCTTATCCCCCGAATGCCGAATTCAAGGCGGACCTGGGAGATAAAATCGAGCGTCGTTTCGCCAATGATccagatcggaaga                                                                                                                                                                                          >  1:118622/1‑290 (MQ=255)
                                                                                         tatGGTGGTCCGTCAGCACCTGCATCACGCTTACTGGCGGCGGGGCGTCGGTGTAGACGGCATCTACCATGCTGATGCTGCCCATATTGACCATCGCGTTACGGCCAAATTTCGAGTGATCGCCAACCAGCATAACGTGGCGCGAGTTCTCAATAATGGCGCGTTTGGTGCGAACTTCGTGGTAATCGAACTCCAGCAGCGAG‑‑GTAGCTATCGATGCCGCTTATCCCCAGAATGCCGAAATCAAGGCGGAACTGGGAGATAAAATCGAGCGTCGCTTCGCCAATGATCCCGCCATCGCGGc                                                                                                                                                                                         <  2:467303/301‑1 (MQ=255)
                                                                                                                        tACTGGCGGCGGGGCGTCGGTGTAGACGGCATCTACCATGCTGATGCTGCCCATATTGACCATCGCGTTACGGCCAAATGTCGAGTGATCGACAACCAGCATAACGTGGCGCGAGTTCTCAATAATGGCGCGTTTGGTGCGAACTTCGTGGTAATCGAACTCCAGCAGCGAG‑‑GTCGCTATCGATGCCGCTTATCCCCAGAATGCCGAAATCAAGGCGGAACTGGGAGATAAAATCGAGCGTCGCTTCGCCAATGATCCCGCCATCGCGGCTGCGTAATTCGCCACCGGCGAGAATGATGCg                                                                                                                                                          <  2:95217/301‑1 (MQ=255)
                                                                                                                          cTGGCGGCGGGGCGTCGGTGTAGACGGCATCTACCATGCTGATGCTGCCCATATTGACCATCGCGTTACGGCCAAATTTCGAGTGATCGACAACCAGCATAACGTGGCGCGAGTTCTCAATAATGGCGCGTTTGGTGCGAACTTCGTGGTAATCGAACTCCAGCAGCGAG‑‑GTCGCTATCGATGCCGCTTATCCCCAGAATGCCGAAATCAAGGCGGAACTGGGAGATAAAATCGAGCGTCGCTTCGCCAATTATCCCGCCATCGCGGCTGCGTAATTCGCCACCGGCGAGAATGATGCGaa                                                                                                                                                        >  1:203492/1‑301 (MQ=255)
                                                                                                                            ggcggcGGGGCGTCGGTGTAGACGGCATCTACCATGCTGATGCTGCCCATATTGACCATCGCGTTACGGCCAAATTTCGAGTGATCGACAACCAGCATAACGTGGCGCGAGTTCTCAATAATGGCGCGTTTGGTGCGAACTTCGTGGTAATCGAACTCCAGCAGCGAG‑‑GTCGCTATCGATGCCGCTTATCCCCAGAATGCCGAAATCAAGGCGGAACTGGGAGATAAAATCGAGCGTCGCTTCGCCAATGATCCCGCCATCGCGGCTGCGTAATTCGCCACCGGCCAGAATGATGCGaaaa                                                                                                                                                      >  1:285603/1‑301 (MQ=255)
                                                                                                                            ggcggcGGGGCGTCGGTGTAGACGGCATCTACCATGCTGATGCTGCCCATATTGACCATCGCGTTACGGCCAAATTTCGAGTGATCGACAACCAGCATAACGTGGCGCGAGTTCTCAATAATGGCGCGTTTGGTGCGAACTTCGTGGTAATCGAACTCCAGCAGCGAG‑‑GTCGCTATCGATGCCGCTTATCCCCAGAATGCCGAAATCAAGGCGGAACTGGGAGATAAAATCGAGCGTCGCTTCGCCAATGAGCCCGCCATCGCTGCTGCGTAATTCGCCAACGGCGAAAATGATGCGaaaa                                                                                                                                                      >  1:62875/1‑301 (MQ=255)
                                                                                                                                 cggGGCGTCGGTGTAGACGGCATCTACCATGCTGATGCTGCCCATATTGACCATCGCGTTACGGCCAAATTTCGAGTGATCGACAACCAGCATAACGTGGCGCGAGTTCTCAATAATGGCGCGTTTGGTGCGAACTTCGTGGTAATCGAACTCCAGCAGCGAG‑‑GTCGCTATCGATGCCGCTTATCCCCAGAATGCCGAAATCAAGGCGGAACTGGGAGATAAAATCGAGCGTCGCTTCGCCAATGATCCCGCCATCGCGGCTGCGTAATTCGCCACCGGCGAGAATGATGCGAAAAtcttc                                                                                                                                                 <  2:62955/301‑1 (MQ=255)
                                                                                                                                    ggCGTCGGTGTAGACGGCATCTACCATGCTGATGCTGCCCATATTGACCATCGCGTTACGGCCAAATTTCGAGTGATCGACAACCAGCATAACGTGGCGCGAGTTCTCAATAATGGCGCGTTTGGTGCGAACTTCGTGGTAATCGAACTCCAGCAGCGAG‑‑GTCGCTATCGATGCCGCTTATCCCCAGAATGCCGAAATCAAGGCGGAACTGGGAGATAAAATCGAGCGTCGCTTCGCCAATGATCCCGCCATCGCGGCTGCGTAATTCGCCACCGGCGAGAATGGTGCGAAAATCTTCttt                                                                                                                                              >  2:15277/1‑301 (MQ=255)
                                                                                                                                     gCGTCGGTGTAGACGGCATCTACCATGCTGATGCTGCCCATATTGACCATCGCGTTACGGCCAAATTTCGAGTGATCGACAACCAGCATAACGTGGCGCGAGTTCTCAATAATGGCGCGTTTGGTGCGAACTTCGTGGTAATCGAACTCCAGCAGCGAG‑‑GTCGCTATCGATGCCGCTTATCCCCAGAATGCCGAAATCAAGGCGGCACTGGGAGATAAAATCGAGCGTCGCTTCGCCCATGATCCCGCCATCGCGGCTGCGTAATTCCCCACCGGCGAGAATGATGCGAAACTcttctt                                                                                                                                               >  1:416722/1‑299 (MQ=255)
                                                                                                                                            tgtAGACGGCATCTACCATGCTGATGCTGCCCATATTGACCATCGCGTTACGGCCAAATTTCGAGTGATCGACAACCAGCATAACGTGGCGCGAGTTCTCAATAATGGCGCGTTTGGTGCGAACTTCGTGGTAATCGAACTCCAGCAGCGAG‑‑GTCGCTATCGATGCCGCTTATCCCCAGAATGCCGAAATCAAGGCGGAACTGGGAGATAAAATCGAGCGTCGCTTCGCCAATGATCCCGCCATCGCGGCTGCGTAATTCGCCACCGGCGAGAATGATGCGAAAATCTTCTTTTACCATc                                                                                                                                       <  2:285949/300‑1 (MQ=255)
                                                                                                                                            tgtAGACGGCATCTACCATGCTGATGCTGCCCATATTGACCATCGCGTTACGGCCAAATTTCGAGTGATCGACAACCAGCATAACGTGGCGCGAGTTCTCAATAATGGCGCGTTTGGTGCGAACTTCGTGGTAATCGAACTCCAGCAGCGAG‑‑GTCGCTATCGATGCCGCTTATCCCCAGAATGCCGAAATCAAGGCGGAACTGGGAGATAAAATCGAGCGTCGCTTCGCCAATGATCCCGCCATCGCGGCTGCGTAATTCGCCACCGGCGAGAATGATGCGAAAATCTTCTTTTACCATc                                                                                                                                       >  1:172791/1‑300 (MQ=255)
                                                                                                                                                gACGGCATCTACCATGCTGATGCTGCCCATATTGACCATCGCGTTACGGCCAAATTTCGAGTGATCGACAACCAGCATAACGTGGCGCGAGTTCTCAATAATGGCGCGTTTGGTGCGAACTTCGTGGTAATCGAACTCCAGCAGCGAG‑‑GTCGCTATCGATGCCGCTTATCCCCAGAATGCCGAAATCAAGGCGGAACTGGGAGATAAAATCGAGCGTCGCTTCGCCAATGATCCCGCCATCGCGGCTGCGTAATTCGCCACCGGCGAGAATGATGCGAAAATCTTCTTTTACCATCAACgt                                                                                                                                  >  2:2677/1‑301 (MQ=255)
                                                                                                                                                 tcttccgatctcATGCTGATGCTGCCCATATTGACCATCGCGTTACGGCCAAATTTCGAGTGATCGACAACCAGCATAACGTGGCGCGAGTTCTCAATAATGGCGCGTTTGGTGCGAACTTCGTGGTAATCGAACTCCAGCAGCGAG‑‑GTCGCTATCGATGCCGCTTATCCCCAGAATGCCGAAATCAAGGCGGAACTGGGAGATAAAATCGAGCGTCGCTTCGCCAATGATCCCGCCATCGCGGCTGCGTAATTCGCCACCGGCGAGAATGATGCGAAAATCTTCTTTTACCATCAACgtg                                                                                                                                 <  2:108133/290‑1 (MQ=255)
                                                                                                                                                 aCGGTATGTACCAGGCTGATGCTGCACATATTGACCATCGAGGTACGGCCAAATGTCGAGTGATCGACAACCAGCATAACGTGGCGCGAGTTCTCAATAATGGCGCGTTTGGTGCGAACTTCGTGGTAATCGAACTCCAGCAGCGAG‑‑GTCGCTATCGATGCCGCTTATCCCCAGAATGCCGAAATCAAGGCGGAACTGGGAGATAAAATCGAGCGTCGCTTCGCCAATGATCCCGCCATCGCGGCTGCGTAATTCGCCACCGGCGAGAATGATGCGAAAATCTTCTTTTACCATCAACg                                                                                                                                   <  1:263906/299‑1 (MQ=255)
                                                                                                                                                           ccATGCTGATGCTGCCCATATTGACCATCGCGTTACGGCCAAATTTCGAGTGATCGACAACCAGCATAACGTGGCGCGAGTTCTCAATAATGGCGCGTTTGGTGCGAACTTCGTGGTAATCGAACTCCAGCAGCGAG‑‑GTCGCTATCGATGCCGCTTATCCCCAGAATGCCGAAATCAAGGCGGAACTGGGAGATAAAATCGAGCGTCGCTTCGCCACTGATCCCGCCATCGCGGCTGCGTAATTCGCCACCGGCGAGAATGATGCGAAAATCTTCTTTTACCATCAACGTGTTAGCAACgg                                                                                                                       >  1:384278/1‑300 (MQ=255)
                                                                                                                                                           ccATGCTGATGCTGCCCATATTGACCATCGCGTTACGGCCAAATTTCGAGTGATCGACAACCAGCATAACGTGGCGCGAGTTCTCAATAATGGCGCGTTTGGTGCGAACTTCGTGGTAATCGAACTCCAGCAGCGAG‑‑GTCGCTATCGATGCCGCTTATCCCCAGAATGCCGAAATCAAGGCGGAACTGGGAGATAAAATCGAGCGTCGCTTCGCCAATGATCCCGCCATCGCGGCTGCGTAATTCGCCACCGGCGAGAATGATGCGAAAATCTTCTTTTACCATCAACGTGTTAGCAACGt                                                                                                                       >  2:281588/1‑301 (MQ=255)
                                                                                                                                                            cATGCTGATGCTGCCCATATTGACCATCGCGTTACGGCCAAATTTCGAGTGATCGACAACCAGCATAACGTGGCGCGAGTTCTCAATAATGGCGCGTTTGGTGCGAACTTCGTGGTAATCGAACTCCAGCAGCGAG‑‑GTCGCTATCGATGCCGCTTATCCCCAGAATGCCGAAATCAAGGCGGAACTGGGAGATAAAATCGAGCGTCGCTTCGCCAATGATCCCGCCATCGCGGCTGCGTAATTCGCCACCGGCGAGAATGATGCGAAAATCTTCTTTTACCATCAACgtgagatcggaaga                                                                                                                      >  1:108003/1‑290 (MQ=255)
                                                                                                                                                                  tatgctgCCCATAATGCCCAGCGCGTTACGTCCAAATTTCGAGTGATCGACAACCAGCATAACGTGGCGCGAGTTCTCAATAATGGCGCGTTTGGTGCTAACTTCGTGGTAATCGAACTCCAGCAGCGAG‑‑GTCGCTATCGATGCCGCTTCTCCCCAGAATGCCGAAATCAAGGCGGAACTGGGAGATAAAATCGAGCGTCGCTTCGCCAATGAGCCCGCCATCGCGGCTGCGTAATTCGCCACCGGCGAGAATGATGCGAAAATCTTCTTTTACCATCAACgttgctaacacgttgatggt                                                                                                                <  1:161426/300‑19 (MQ=255)
                                                                                                                                                                   atgctgCCCATATTGACCATCGCGTTACGGCCAAATTTCGAGTGATCGACAACCAGCATAACGTGGCGCGAGTTCTCAATAATGGCGCGTTTGGTGCGAACTTCGTGGTAATCGAACTCCAGCAGCGAG‑‑GTCGCTATCGATGCCGCTTATCCCCAGAATGCCGAAATCAAGGCGGAACTGGGAGATAAAATCGAGCGTCGCTTCGCCAATGATCCCGCCATCGCGGCTGCGTAATTCGCCACCGGCGAGAATGATGCGAAAATCTTCTTTTACCATCAACGTGTTAGCAACGTTGAGAttt                                                                                                               >  1:59558/1‑300 (MQ=255)
                                                                                                                                                                         cccATATTGACCATCGCGTTACGGCCAAATTTCGAGTGATCGACAACCAGCATAACGTGGCGCGAGTTCTCAATAATGGCGCGTTTGGTGCGAACTTCGTGGTAATCGAACTCCAGCAGCGAG‑‑GTCGCTATCGATGCCGCTTATCCCCAGAATGCCGAAATCAAGGCGGAACTGGGAGATAAAATCGAGCGTCGCTTCGCCAATGATCCCGCCATCGCGGCTGCGTAATTCGCCACCGGCGAGAATGATGCGAAAATCTTCTTTTACCATCAACGTGTTAGCAACGTTGAGATTGTTGGt                                                                                                          <  2:173001/300‑1 (MQ=255)
                                                                                                                                                                               ttGACCATCGCGTTACGGCCAAATTTCGAGTGATCGACAACCAGCATAACGTGGCGCGAGTTCTCAATAATGGCGCGTTTGGTGCGAACTTCGTGGTAATCGAACTCCAGCAGCGAG‑‑GTCGCTATCGATGCCGCTTATCCCCAGAATGCCGAAATCAAGGCGGAACTGGGAGATAAAATCGAGCGTCGCTTCGCCAATGATCCCGCCATCGCGGCTGCGTAATTCGCCACCGGCGAGAATGATGCGAAAATCTTCTTTTACCATCAACGTGTTAGCAACGTTGAGATTGTTGGTGACAATg                                                                                                   <  2:384758/301‑1 (MQ=255)
                                                                                                                                                                                tGACCATCGCGTTACGGCCAAATTTCGAGTGATCGACAACCAGCATAACGTGGCGCGAGTTCTCAATAATGGCGCGTTTGGTGCGAACTTCGTGGTAATCGAACTCCAGCAGCGAG‑‑GTCGCTATCGATGCCGCTTATCCCCAGAATGCCGAAATCAAGGCGGAACTGGGAGATAAAATCGAGCGTCGCTTCGCCAATGATACCGCCATCGCGGCTGCGTAATTCGCCACCGGCGAGAATGATGCGAAAATCTTCTTTTACCATCAACCTGTTAGCAACGTTGAGATTGTTGGTGACAATg                                                                                                   >  2:433605/1‑300 (MQ=255)
                                                                                                                                                                                  aCCATCGCGTTACGGCCAAATTTCGAGTGATCGACAACCAGCATAACGTGGCGCGAGTTCTCAATAATGGCGCGTTTGGTGCGAACTTCGTGGTAATCGAACTCCAGCAGCGAG‑‑GTCGCTATCGATGCCGCTTATCCCCAGAATGCCGAAATCAAGGCGGAACTGGGAGATAAAATCGAGCGTCGCTTCGCCAATGATCCCGCCATCGCGGCTGCGTAATTCGCCACCGGCGAGAATGATGCGAAAATCTTCTTTTACCATCAACGTGTTAGCAACGTTGAGATTGTTGGTGACAATgcg                                                                                                 <  2:203747/300‑1 (MQ=255)
                                                                                                                                                                                           ttACGGCCAAATTTCGAGTGATCGACAACCAGCATAACGTGGCGCGAGTTCTCAATAATGGCGCGTTTGGTGCGAACTTCGTGGTAATCGAACTCCAGCAGCGAG‑‑GTCGCTATCGATGCCGCTTATCCCCAGAATGCCGAAATCAAGGCGGAACTGGGAGATAAAATCGAGCGTCGCTTCGCCAATGATCCCGCCATCGCGGCTGCGTAATTCGCCACCGGCGAGAATGATGCGAAAATCTTCTTTTACCATCAACGTGTTAGCAACGTTGAGATTGTTGGTGACAATGCGCAAATTGCtg                                                                                       <  2:175789/301‑1 (MQ=255)
                                                                                                                                                                                              cGGCCAAATTTCGAGTGATCGACAACCAGCATAACGTGGCGCGAGTTCTCAATAATGGCGCGTTTGGTGCGAACTTCGTGGTAATCGAACTCCAGCAGCGAG‑‑GTCGCTATCGATGCCGCTTATCCCCAGAATGCCGAAATCAAGGCGGAACTGGGAGATAAAATCGAGCGTCGCTTCGCCAATGATCCCGCCATCGCGGCTGCGTAATTCGCCACCGGCGAGAATGATGCGAAAATCTTCTTTTACCATCAACGTGTTAGCAACGTTGAGATTGTTGGTGACAATGCGCAAATTGCTGTGa                                                                                    >  2:253954/1‑301 (MQ=255)
                                                                                                                                                                                                 ccAAATTTCGAGTGATCGATAACCAGCATAACGTGGCGCGAGTTCTCAATAATGGCGCGTTTGGTGCGAACTTCGTGGTAATCGAACTCCAGCAGCGAG‑‑GTCGCTATCGATGCCGCTTATCCCCAGAATGCCGAAATCAAGGCGGAACTGGGAGATAAAATCGAGCGTCGCTTCGCCAATGATCCCGCCATCGCGGCTGCGTAATTCGCCACCGGCGAGAATGATGCGAAAATCTTCTTTTACCATCAACGTGTTAGCAACGTTGAGATTGTTGGTGACAATGCGCAAATTGCTGTGAtt                                                                                  <  2:145800/300‑1 (MQ=255)
                                                                                                                                                                                                             tGATCGACAACCAGAATAACGGGGCGCGAGTTCTCAATAATGGCGCGTTTGGTGCGAACTTCGTGGTAATCGAACTCCAGCAGCGAG‑‑GTCGCTATCGATGCCGCTTATCCCCAGAATGCCGAAATCAAGGCGGAACTGGGAGATAAAATCGAGCGTCGCTTCGCCAATGATCCCGCCATCGCGGCTGCGTAATTCGCCACCGGCGAGAATGATGCGAAAATCTTCTTTTACCATCAACGTGTTAGCAACGTTGAGATTGTTGGTGACAATGCGCAAATTGCTGTGATTGAGCAgtgcgtg                                                                      <  1:15262/300‑1 (MQ=255)
                                                                                                                                                                                                                     aaCCAGCATAACGTGGCGCGAGTTCTCAATAATGGCGCGTTTGGTGCGAACTTCGTGGTAATCGAACTCCAGCAGCGAG‑‑GTCGCTATCGATGCCGCTTATCCCCAGAATGCCGAAATCAAGGCGGAACTGGGAGATAAAATCGAGCGTCGCTTCGCCAATGATCCCGCCATCGCGGCTGCGTAATTCGCCACCGGCGAGAATGATGCGAAAATCGTCTTTTACCATCAACGTGTTAGCAACGTTGAGATTGTTGGTGACAATGCGCAAATTGCTGTGATTGAGCAGTGCGGTCGCTACCg                                                              >  1:177347/1‑300 (MQ=255)
                                                                                                                                                                                                                               aCGTGGCGCGAGTTCTCAATAATGGCGCGTTTGGTGCGAACTTCGTGGTAATCGAACTCCAGCAGCGAG‑‑GTCGCTATCGATGCCGCTTATCCCCAGAATGCCGAAATCAAGGCGGAACTGGGAGATAAAATCGAGCGTCGCTTCGCCAATGATCCCGCCATCGCGGCTGCGTAATTCGCCACCGGCGAGAATGATGCGAAAATCTTCTTTTACCATCAACGTGTTAGCAACGTTGAGATTGTTGGTGACAATGCGCAAATTGCTGTGATTGAGCAGTGCGTGCGCTACCGCTTCCGGCgtg                                                   <  2:59625/301‑1 (MQ=255)
                                                                                                                                                                                                                                           ttCTCAATAATGGCGCGTTTGGTGCGAACTTCGTGGTAATCGAACTCCAGCAGCGAG‑‑GTCGCTATCGATGCCGCTTATCCCCAGAATGCCGAAATCAAGGCGGAACTGGGAGATAAAATCGAGCGTCGCTTCGCCAATGATCCCGCCATCGCGGCTGCGTAATTCGCCACCGGCGAGAATGATGCGAAAATCTTCTTTTACCATCAACGTGTTAGCAACGTTGAGATTGTTGGTGACAATGCGCAAATTGCTGTGATTGAGCAGTGCGTGCGCTACCGCTTCCGGCGTGGTGCCGATATc                                        <  2:209562/300‑1 (MQ=255)
                                                                                                                                                                                                                                           agctcAATAATGGATCGTTTGGTGCGACCTTCTGGGTAATCGAACCCCAGCAGCGAG‑‑GTCGCTATCGATGCCGCATATACCAGGAATGCCGAAATCAAGGCGGAACTGGGAGATAAAATCGCGCGTCGCATGGCCACAGATCCCGGCATCGCGGCTGCGTAATTCGCCACCGGCGAGAATGATGCGAAAATCTTCTTTTACCATCAACGTGTTAGCAACGTTGAGATTGTTGGTGACAATGCGCAAATTGCTGTGATTGAGCAGTGCGTGCGCTACCGCTTCCGGCGTGGTGCCGATATc                                        <  1:433058/298‑1 (MQ=255)
                                                                                                                                                                                                                                              tcAATAATGGCGCGTTTGGTGCGAACTTCGTGGTAATCGAACTCCAGCAGCGAG‑‑GTCGCTATCGATGCCGCTTATCCCCAGAATGCCGAAATCAAGGCGGAACTGGGAGATAAAATCGAGCGTCGCTTCGCCAATGATCCCGCCATCGCGGCGGCGTAATTCGCCACCGGCGAGAATGATGCGGAAATCTTCTTTTACCATCAACGTGTTAGCAACGTTGAGATTGTTGGTGACAATGCGCAAATTGCTGTGATTGAGCAGTGCGTGCGCTACCGCTTCCGGCGTGGTGCCGATATCGATa                                    <  2:166645/301‑1 (MQ=255)
                                                                                                                                                                                                                                               cAATAATGGCGCGTTTGGTGCGAACTTCGTGGTAATCGAACTCCAGCAGCGAG‑‑GTCGCTATCGATGCCGCTTATCCCCAGAATGCCGAAATCAAGGCGGAACTGGGAGATAAAATCGAGCGTCGCTTCGCCAATGATCCCGCCATCGCGGCTGCGTAATTCGCCACCGGCGAGAATGATGCGGAAATCTTCTTTTACCATCAACGTGTTAGCAACGTTGAGATTGTTGGTGACAATGCGCAAATTGCTGTGATTGAGCAGTGCGTGCGCTACCGCTTCCGGCGTGGTGCCGATTTCGATaa                                   >  1:166435/1‑301 (MQ=255)
                                                                                                                                                                                                                                                        cgcgGTTGGGGCGAACTTCGGGGTAATCGAAAGCAAGAAGCGAG‑‑GTCGCTATCGATGCCGCTTATCACCAGAATGCCGAAATCAAGGCGGAACTGGGAGAAAAAATCGAGCGTCGCTTCGCCAATGATCCCGCCATCGCGGCTGCGTAATTCGCCACCGGAGAGAATGATGCGAAAATCTTCTTTTACCATCAACGTGTTAGCAACGTTGAGATTGTTGGTGACAATGCGCAAATTGCTGTGATTGAGCAGTGCGTGCGCTACCGCTTCCGGCGTGGTGCCGATATCGATAAACAGCGTc                           <  1:281242/300‑1 (MQ=255)
                                                                                                                                                                                                                                                                 tGCGAACTTCGTGGTAATCGAACTCCAGCAGCGAG‑‑GTCGCTATCGATGCCGCTTATCCCCAGAATGCCGAAATCAAGGCGGAACTGGGAGATAAAATCGAGCGTCGCTTCGCCAATGATCCCGCCATCGCGGCTGCGTAATTCGCCACCGGCGAGAATGATGCGAAAATCTTCTTTTACCATCAACGTGTTAGCAACGTTGAGATTGTTGGTGACAATGCGCAAATTGCTGTGATTGAGCAGTGCGTGCGCTACCGCTTCCGTCGTGGTGCCGATATCGATAAACAGCGTCGAGCCATAgg                 >  2:163113/1‑301 (MQ=255)
                                                                                                                                                                                                                                                                 tGCGAACTTCGTGGTAATCGAACTCCAGCAGCGAG‑‑GTCGCTATCGATGCCGCTTATCCCCAGAATGCCGAAATCAAGGCGGAACTGGGAGATAAAATCGAGCGTCGCTTCGCCAATGATCCCGCCATCGCGGCTGCGTAATTCGCCACCGGCGAGAATGATGCGAAAATCTTCTTTTACCATCAACGTGTTAGCAACGTTGAGATTGTTGGTGACAATGCGCAAATTGCTGTGATTGAGCAGTGCGTGCGCTACCGCTTCCGGCGTGGTGCCGATATCGATAAACAGCGTCGAGCCATTgg                 >  2:447785/1‑301 (MQ=255)
                                                                                                                                                                                                                                                                       cttccaggTAATCTAAGTCCAGAAGCGAG‑‑GTCGCTATCGATGCCGATTATACCCAGAATGCCGAAATCAAGGCGGAACTGGGAGATAAAATCGAGCGTCGCTTCGCCAATGATACCACCATCGCGGCTGCGTAATTCGCCACCGGCGAGAATGATGCGAAAATCTTCTTTTACCATCAACGTGTTAGCAACGTTGAGATTGTTGGTGACAATGCGCAAATTGCTGTGATTGAGCAGTGCGTGCGCTACCGCTTCCGGCGTGGTGCCGATATCGATAAACAGCGTCGAGCCATTGGGGAttt           <  1:115416/295‑1 (MQ=255)
                                                                                                                                                                                                                                                                                aaTCGAACTCCAGCAGCGAG‑‑GTCGCTATCGATGCCGCTTATCCCCAGAATGCCGAAATCAAGGCGGAACTGGGAGATAAAATCGAGCGTCGCTTCGCCAATGATCCCGCCATCGCGGCTGCGTAATTCGCCACCGGCGAGAATGATGCGAAAATCTTCTTTTACCATCAACGTGTTAGCAACGTAGAGATTGGTGGTGAAAATGCGCAAATTGCTGTGATTTAGCAGTGCGTGCGCTACCGCTTCAGGCGTGGTGCCGGTAGCGCTAAACTGCGTCGAGCACTTGGAGTTTTGATCAGcc   >  1:162306/1‑300 (MQ=255)
                                                                                                                                                                                                                                                                                 aTCGAACTCCAGCAGCGAG‑‑GTCGCTATCGATGCCGCTTATCCCCAGAATGCCGAAATCAAGGCGGAACTGGGAGATAAAATCGAGCGTCGCTTCGCCAATGATCCCGCCATCGCGGCTGCGTAATTCGCCACCGGCGAGAATGATGCGAAAATCTTCTTTTACCATCAACGTGTTAGCAACGTTGAGATTGTTGGTGACAATGCGCAAATTGCTGTGATTGAGCAGTGCGTGCGCTACCGCTTCCGGCGTGGTGCCGATATCGATAAACAGCGTCGAGCCAGTGGGGATTTGCTCCGCCa  >  1:295687/1‑300 (MQ=255)
                                                                                                                                                                                                                                                                                                    ||                                                                                                                                                                                                                                                                                         
CGTCAATGACGCGGATCGGCTACGCAAATACCTGGCGCGTTTTGGTCTGACGTGGGAAGCCGTGCAGGATCAGCACAGCTCCAGTTGAATATGGTGGTCCGTCAGCACCTGCATCACGCTTACTGGCGGCGGGGCGTCGGTGTAGACGGCATCTACCATGCTGATGCTGCCCATATTGACCATCGCGTTACGGCCAAATTTCGAGTGATCGACAACCAGCATAACGTGGCGCGAGTTCTCAATAATGGCGCGTTTGGTGCGAACTTCGTGGTAATCGAACTCCAGCAGCGAGCCGTCGCTATCGATGCCGCTTATCCCCAGAATGCCGAAATCAAGGCGGAACTGGGAGATAAAATCGAGCGTCGCTTCGCCAATGATCCCGCCATCGCGGCTGCGTAATTCGCCACCGGCGAGAATGATGCGAAAATCTTCTTTTACCATCAACGTGTTAGCAACGTTGAGATTGTTGGTGACAATGCGCAAATTGCTGTGATTGAGCAGTGCGTGCGCTACCGCTTCCGGCGTGGTGCCGATATCGATAAACAGCGTCGAGCCATTGGGGATTTGCTCCGCCA  >  NC_000913/3559779‑3560353

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