Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 2,173,363 Δ2 bp intergenic (‑490/+918) gatD ← / ← gatB galactitol‑1‑phosphate 5‑dehydrogenase/galactitol‑specific PTS enzyme IIB component

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,173,3630C.91.4% 134.6 / 7.3 35intergenic (‑490/+919)gatD/gatBgalactitol‑1‑phosphate 5‑dehydrogenase/galactitol‑specific PTS enzyme IIB component
Reads supporting (aligned to +/- strand):  ref base C (1/2);  new base . (23/9);  total (24/11)
Fisher's exact test for biased strand distribution p-value = 2.27e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.44e-01
*NC_0009132,173,3640C.94.3% 142.6 / 3.9 35intergenic (‑491/+918)gatD/gatBgalactitol‑1‑phosphate 5‑dehydrogenase/galactitol‑specific PTS enzyme IIB component
Reads supporting (aligned to +/- strand):  ref base C (1/1);  new base . (23/10);  total (24/11)
Fisher's exact test for biased strand distribution p-value = 5.36e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.17e-01

CCAGCATTAGCCGCCAGTTGGGTGTGCAAACCAATAGTTTGCGTCGCGATCCACAGGGTGATGCTCATAATGATGACACCCGAGATTAAGGTGCGGAACAGATTTCCACGATGCACGGCGACCGCCATCGCCACGAAGAAGCCGATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAG‑CAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCACCCCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTAATCACCCGTGGCATTAGCAGCATCACTGCCGCCGTTTTTACCGCCAGCTGCAATACACCTTTGACATCGTAACCCGCGAGGATGCCGATAATCAGCCCCATCACAAAACCCACG  >  NC_000913/2173120‑2173601
                                                                                                                                                                                                                                                    ||                                                                                                                                                                                                                                             
ccAGCATTAGCCGCCGGTTGGTTGTCCAAACAAATATTTGGCGCCGCGATCCACAGGGTGATGCTCATAATGATGACCCCCGAGATTAAGGTGCGGAACAGATTTCCACGATGCACGGCGACCGCCATCGCCACGAAGAAGCCGATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAG‑CAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCACC‑‑GCCGTATg                                                                                                                                                                                                                                       <  2:251048/251‑1 (MQ=255)
     aTTAGCCGCCAGTTGGGTGTGCAAACCAATAGTTTGCGTCGCGATCCACAGGGTGATGCTCATAATGATGACACCCGAGATTAAGGTGCGGAACAGATTTCCACGATGCACGGCGACCGCCATCGCCACGAAGAAGCCGATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAG‑CAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCACC‑‑GCCGTATGTccc                                                                                                                                                                                                                                   >  2:52366/1‑250 (MQ=255)
     aTTAGCAGCCAGTTGGGTGTGCAAACCAATAGTTTGCGTCGCGATCCACAGGGTGATGCTCATAATGATGACACCCGAGCTTCAGGCGCGGAACAGATTTCCACGATGCACGGCGACCGCCATCGCCACGAAGAAGCCGATGGGGGCAAGATCGCCAAACGGCCGCACCTGATTACCCGGCACACCACCAG‑CACTTAAAATGGTGTGAGGGATAACAATCAGGCTTGCCGATACCAcccacgtatgtccc                                                                                                                                                                                                                                     >  2:397732/1‑239 (MQ=255)
                       gtgCAAACCAATAGTTTGCGTCGCGATCCACAGGGTGATGCTCATAATGATGACACCCGAGATTAAGGTGCGGAACAGATTTCCACGATGCACGGCGACCGCCATCGCCACGAAGAAGCCGATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAG‑CAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCACC‑‑GCCGTATGTCCCAGCAGCAACGCCGGATCa                                                                                                                                                                                                                 <  1:408168/250‑1 (MQ=255)
                              ccactAGTTTGCGTCGCGATCCACAGGGTGATGCTCATAATGATGACACCCGAGATTAAGGTGCGGAACAGATTTCCACGATGCACGGCGACCGCCATCGCCACGAAGAAGCCGCTGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAG‑CAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCACC‑‑GCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAAt                                                                                                                                                                                                         >  2:40786/5‑251 (MQ=255)
                                     tttGCGTCGCGATCCACAGGGTGATGCTCATAATGATGACACCCGAGATTAAGGTGCGGAACAGATTTCCACGATGCACGGCGACCGCCATCGCCACGAAGAAGCCGATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAG‑CAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCACC‑‑GCCGTATGTCCCAGCAGCAACGCCGGATCACGGCCAATCAGGAAc                                                                                                                                                                                                  >  1:117397/1‑251 (MQ=255)
                                                       gggTGATGCTCATAATGATGACCCCCGAGATTAAGGTGCGGAACAGATTTCCACGATGCACGGCGACCGCCATCGCCACGAAGAAGCCGATGGTGGAAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAG‑CAATTAAAATGGTGAGTGGGATAAAAATCAGTCTTGCCGATACCACC‑‑GCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTc                                                                                                                                                                                <  1:119510/251‑1 (MQ=255)
                                                        gtggtgccccaaaaTGATGACCCCGGAGATTAAGGTCCGAAAAAGATTGCCCCGAGGCCGGGCTACCGCCACCGCCACGAAGAAGCCGATGGTGGCAAGATCGCAAAACGGCAGCACCTGATTCCCCGGCACACAACCAG‑CAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCACC‑‑GCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCg                                                                                                                                                                               <  2:227865/239‑1 (MQ=255)
                                                        gggaaTGCTAATAAGGATGACACCCGAGTTTAAGGTGCGGAACAGATTCCCACGAGGCACGCCGACCCCAAACGCCACGAAGAAGCCGAGGGTGGCAAGATCGCAAAACGGCGGCACCTGATTACCCGGCACACAAACAG‑CAATAAAAAGGGTGGGTGGGATAAAAATCAGGCGTGCCGATACCACC‑‑GCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCg                                                                                                                                                                               <  2:74613/247‑1 (MQ=255)
                                                        ggTGATGCTCATAATGATGACACCCGAGATTAAGGTGCGGAACAGATTTCCACGATGCACGGCGACCGCCATCGCCACGAAGAAGCCGATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAG‑CAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCACC‑‑GCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCg                                                                                                                                                                               >  1:397492/1‑251 (MQ=255)
                                                              gCTCATAATGCTGACACCCGAGATTAAGGTGCGGAACCGATTTCCACGATGCACGGCGACCGCCATCGCCACGAAGAAGCCGATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAG‑CAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCACC‑‑GCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGt                                                                                                                                                                          >  2:128824/1‑250 (MQ=255)
                                                                     atgatgACACCCGAGATTAAGGTGCGGAACAGATTTCCACGATGCACGGCGACCGCCATCGCCACGAAGAAGCCGATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAG‑CAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCACC‑‑GCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACt                                                                                                                                                                  >  1:299190/1‑251 (MQ=255)
                                                                                 gagaTTAAGGTGCGGAACAGATTTCCACGATGCACGGCGACCGCCATCGCCACGAAGAAGCCGATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAG‑CAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCACC‑‑GCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCtt                                                                                                                                                      >  1:434796/1‑251 (MQ=255)
                                                                                           tGCGGAACAGATTTCCACGATGCACGGCGACCGCCATCGCCACGAAGAAGCCGATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAG‑CAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCACC‑‑GCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGc                                                                                                                                            >  2:96952/1‑251 (MQ=255)
                                                                                                aaCAGATTTCCACGATGCACGGCGACCGCCATCGCCACGAAGAAGCCGATGGCGGCAAGAGCGCCAAACGGCAGCACCTGATGACCCGGCACACAAACAG‑CACTTAACATGGGGAGTGGGATAAAAATCAGGCTTGCCGATACCACC‑‑GCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCACTCAGGAACTCCTAACAGCCGCACTTCACCTGTGAACGACTACAAGCCGGCTTAGCGATGGGCGtga                                                                                                                                        >  2:139847/1‑248 (MQ=255)
                                                                                                              aTGCACGGCGACCGCCATCGCCACGAAGAAGCCGATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAG‑CAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCACC‑‑GCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAACCATCCATGATgg                                                                                                                         >  2:361781/1‑251 (MQ=255)
                                                                                                                          cGCCATCGCCACGAAGAAGCCGATGGTGGCAAGATCGCCAACCGGCAGCACCTGATTAACCGGCACACAAACAG‑CAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCACC‑‑GCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAAATTCGCCTGTAAACGCCTACGAGACTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTAATCAccc                                                                                                             >  2:239637/1‑251 (MQ=255)
                                                                                                                          cGCCATCGCCACGAAGAAGCCGATGGTGGCAAGATCGCCAAACGGCACCACCTGATTACCCGGCACACAAACAG‑CAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATCCCACC‑‑GCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTAATCAccc                                                                                                             <  1:169264/251‑1 (MQ=255)
                                                                                                                                 cccACGAAAAACCCGATGGTGGAAAGAGCCCCAACCGCCAGCCCATTATAACCCTGCAC‑CAAAAAGGCAATAAAAATGGTGGTTGGGATAAAAATCAGGCTTGCCGATCCCACC‑‑GCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAACCCTCCATGATGGGTTTAATCACCCGTGGCAt                                                                                                      <  2:1675/250‑1 (MQ=255)
                                                                                                                                       aagaagCCGATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAG‑CAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCACC‑‑GCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTAATCACCCGTGGCATTagc                                                                                                  >  1:16074/1‑249 (MQ=255)
                                                                                                                                         gaagCCGATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAG‑CAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCACC‑‑GCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTAATCACCCGTGGCATTagcagc                                                                                               >  2:75006/1‑250 (MQ=255)
                                                                                                                                                             cGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAG‑CAATTAAAATGGTGAGTGGGATAAAAATCAGGCTGGCCGATACCACC‑‑GCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTAATCACCCGTGGCATTAGCAGCATCACTGCCGCCGTTTTTAcc                                                                          <  1:402338/251‑1 (MQ=255)
                                                                                                                                                                       cagcaCCTGATTACCCGGCACCCAAACAG‑CAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCACC‑‑GCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAACCGACTACGAGCCTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTAATCACCCGTGGCATTAGCAGCATCACTGCCGCCGTTTTTACCGCCAGCTGc                                                                 >  1:183204/1‑250 (MQ=255)
                                                                                                                                                                                aTTACCCGGCACACAAACAG‑CAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCACC‑‑GCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGCACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTAATCACCCGTGGCATTAGCAGCATCCCTGCCGCCGTTTTTACCGCCAGCTGCAATAAACCtt                                                       >  2:405344/1‑251 (MQ=255)
                                                                                                                                                                                    cgcggtacacaAACAG‑CAATTAACATGGAGAGTGGGATAAAAATCAGGCTTGCCGATA‑‑ACACCGCCGTATGTAACAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAACCTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTAATCACACGTGGCATTAGCAGCATCACTGCCGCCGTTTTTACCGCCAGCTGCAATACACCTTTGAc                                                   >  2:303830/7‑251 (MQ=255)
                                                                                                                                                                                       ggCACACAAACAG‑CAATTAAAAGGGTGAGTGGGAAAAAAACCAGGCTTGCCGTTACCACC‑‑GCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTAATCACCCGTGGCATTAGCAGCATCACTGCCGCCGTTTTTACCGCCAGCTGCAATACACCTTTGACATc                                                <  1:262552/251‑1 (MQ=255)
                                                                                                                                                                                        gCACACAAACAG‑CAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCACC‑‑GCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTAATCACCCGTGGCATTAGCAGCATCACTGCCGCCGTTTTTACCGCCAGCTGCAATACACCTTTGACATCg                                               >  1:112805/1‑251 (MQ=255)
                                                                                                                                                                                         cacacaAACAG‑CAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCACC‑‑GCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTAATCACCCGTGGCATTAGCAGCATCACTGCCGCCGTTTTTACCGCCAGCTGCAATACACCTTTGACATCGt                                              >  2:82206/1‑251 (MQ=255)
                                                                                                                                                                                              aaaCGG‑CAATTAAAATGGTGAGTGGGATAAAAATCAGGCTGCCCGAT‑CC‑CCCCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACCCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTAATCACCCGTGGCATTAGCAGCATCACTGCCGCCGTTTTTACCGCCAGCTGCAATACACCTTTGACATCGTAAccc                                         <  1:286840/251‑1 (MQ=255)
                                                                                                                                                                                                     cAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCACC‑‑GCCGTATGTCCCAGCAGAAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTAATCACCAGTGGCATTAGCAGCATCACTGCCGCCGTTTTTACCGCCAGCTGCAATACACCTTTGACATCGTAACCCGCGAgg                                   >  2:413502/1‑251 (MQ=255)
                                                                                                                                                                                                                   agtgggctccaaaacAGGCTTGCAGCTACCACC‑‑GGCGTATGTCCCAGCAGCAAAGCCGGATCACGGCAAATCAGGAACTCCTGACCGCCGAACTTAGCCTGTAAACGACTACGAGCCTGCTTAGCGAAGGGCGTTAAACCATCGATGATGGGTTTAATCACCCGTGGCATTAGCAGCATCACTGCCGCCGTTTTTACCGCCAGCGGCAATACACCTTTGACATCGTAACCCGCGCGGATGCCGAGAATc                       >  2:323382/15‑249 (MQ=255)
                                                                                                                                                                                                                           aaaaaTCAGGCTTGCCGATACCACC‑‑GCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTAATCACCCGTGGCATTAGCAGCATCACTGCCGCCGTTTTTACCGCCAGCTGCAATACACCTTTGACATCGTAACCCGCGAGGATGCCGATAATCAGCCCCATc              >  2:3646/1‑250 (MQ=255)
                                                                                                                                                                                                                              aaTCAGGCTTCCCGATACCACC‑‑GCCGTATGTCCCAGCAGCAACGCCGGATCCAGGCCAATCAGGAACTCCTGACCGCCGACCTTCGCCTGTAAACGACTACGAACCTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTAATCACCCGTGGCATTAGCAGCATCACTGCCGCCGTTTTTAACGCCAGCTGCAATACACCTTTGACATCGTAACCCGCGAGGATGCCGATAATCAGCCCCATcaca           >  2:312518/1‑250 (MQ=255)
                                                                                                                                                                                                                                  aGGCTTGCCGATACCACC‑‑GCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTAATCACCCGTGGCATTAGCAGCATCACTGCCGCCGTTTTTACCGCCAGCTGCAATACACCTTTGACATCGTAACCCGCGAGGATGCCGATAATCAGCCCCATCACAAAAcc      >  2:254611/1‑251 (MQ=255)
                                                                                                                                                                                                                                    gCTTGCCGATACCACC‑‑GCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTAATCACCCGTGGCATTAGCAGCATCACTGCCGCCGTTTTTACCGCCAGCTGCAATACACCTTTGACATCGTAACCCGCGAGGATGCCGATAATCAGCCCCATCACAAAAccc     >  1:337307/1‑250 (MQ=255)
                                                                                                                                                                                                                                     ggtGCCGCTACCACCC‑GCCGTATGTCCCAGCAGCAACGCCGGATAAAGGCCAATCAGGAACTCCTGCCCGCGGAACGTCGCCTGTAAACGACAACGAGCCTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTAATCACCCGTGGCATTAGCAGCATCACTGCCGCCGTTTGTACCGCCAGCTGCAATACACCTTTGACATCGTAACCCGCGAGGATGCCGATAATCAGCCCCATCACAAAACCCa    <  2:201553/249‑1 (MQ=255)
                                                                                                                                                                                                                                         tgccgttccccccGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTAATCACCCGTGGCATTAGCAGCATCACTGCCGCCGTTTTTACCGCCAGCTGCAATACACCTTTGACATCGTAACCCGCGAGGATGCCGATAATCAGCCCCATCACAAAACCCACg  <  1:96952/241‑1 (MQ=255)
                                                                                                                                                                                                                                                    ||                                                                                                                                                                                                                                             
CCAGCATTAGCCGCCAGTTGGGTGTGCAAACCAATAGTTTGCGTCGCGATCCACAGGGTGATGCTCATAATGATGACACCCGAGATTAAGGTGCGGAACAGATTTCCACGATGCACGGCGACCGCCATCGCCACGAAGAAGCCGATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAG‑CAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCACCCCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTAATCACCCGTGGCATTAGCAGCATCACTGCCGCCGTTTTTACCGCCAGCTGCAATACACCTTTGACATCGTAACCCGCGAGGATGCCGATAATCAGCCCCATCACAAAACCCACG  >  NC_000913/2173120‑2173601

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