Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 2,173,363 Δ2 bp intergenic (‑1/+1) gatC ← / ← gatC pseudogene, galactitol‑specific enzyme IIC component of PTS;transport; Transport of small molecules: Carbohydrates, organic acids, alcohols; PTS system galactitol‑specific enzyme IIC/pseudogene, galactitol‑specific enzyme IIC component of PTS;transport; Transport of small molecules: Carbohydrates, organic acids, alcohols; PTS system galactitol‑specific enzyme IIC

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,173,3610C.92.5% 156.6 / 3.4 40pseudogene (2/442 nt)gatCpseudogene, galactitol‑specific enzyme IIC component of PTS;transport; Transport of small molecules: Carbohydrates, organic acids, alcohols; PTS system galactitol‑specific enzyme IIC
Reads supporting (aligned to +/- strand):  ref base C (0/3);  new base . (23/14);  total (23/17)
Fisher's exact test for biased strand distribution p-value = 6.88e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
*NC_0009132,173,3620C.92.3% 149.5 / 3.5 40pseudogene (1/442 nt)gatCpseudogene, galactitol‑specific enzyme IIC component of PTS;transport; Transport of small molecules: Carbohydrates, organic acids, alcohols; PTS system galactitol‑specific enzyme IIC
Reads supporting (aligned to +/- strand):  ref base C (0/3);  new base . (23/13);  total (23/17)
Fisher's exact test for biased strand distribution p-value = 6.13e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CGCCAGTTGGGTGTGCAAACCAATAGTTTGCGTCGCGATCCACAGGGTGATGCTCATAATGATGACACCCGAGATTAAGGTGCGGAACAGATTTCCACGATGCACGGCGACCGCCATCGCCACGAAGAAGCCGATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCACCCCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTAATCACCCGTGGCATTAGCAGCATCACTGCCGCCGTTTTTACCGCCAGCTGCAATACACCTTTGACATCGTAACCCGCGAGGATGCCGATAATCAGCCCCATCACAA  >  NC_000913/2173131‑2173593
                                                                                                                                                                                                                                        ||                                                                                                                                                                                                                                       
cgcgAGTTGGGTGTACAAACGAATAGTTGGGGTGGGAGTACACGGTGTGATGCTAAAAATGAGGACACCCGAGATTAAGGGGGGAAACAGTTTCCCAGGATGCCCGGCGCCCGCCATCGCCACGAAGAAGCCGAGGGTGGCAAGATCGCAAAACGGCAGCACCTGATTCCCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATaccaccgccgtat                                                                                                                                                                                                                                  <  2:244565/239‑8 (MQ=255)
   cAGTTGGGTGTGCAAACCAATAGTTTGCGTCGCGATCCACAGGGTGATGCTCATAATGATGACACCCGAGATTAAGGTGCGGAACAGATTTCCACGATGCACGGCGACCGCCATCGCCACGAAGAAGCCGATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGCTACCA‑‑CCGCCGTATGt                                                                                                                                                                                                                              >  2:226221/1‑238 (MQ=255)
             tgCAAACCAATAGTTTGCGTCGCGATCCACAGGGTGATGCTCATAATGATGACACCCGAGATTAAGGTGCGGAACAGATTTCCACGATGCACGGCGACCGCCATCGCCACGAAGAAGCCGATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAgcagca                                                                                                                                                                                                                    >  1:317162/1‑238 (MQ=255)
                     aaTGGTTGGCGGCGCGATCCACAGGGTGATGCTCATAATGATGACCCCGGAGATTAAGGTGCGGAACAGATTTCCACGATGCACGGCGACCGCCATCGCCAAGAAGAAGCCGATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGAt                                                                                                                                                                                                           <  1:71597/239‑1 (MQ=255)
                           ttGCGTCGCGATCCACAGGGTGATGCTCATAATGATGACACCCGAGATTAAGGTGCGGAACAGATTTCCACGATGCACGGCGACCGCCATCGCCACGAAGAAGCCGATGGTGGCAAGATCGCCAAACTGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCA‑‑CCGCCGGATGTCCCAGCAGCAACGCCGGATCAAgg                                                                                                                                                                                                      >  2:698336/1‑238 (MQ=255)
                                       ccACGGGGTGATGCTCAAAAGGATGCCCCCCGAAATTAGGGTGGGGAACAGTTTCCCACGAGGCCCGGCGACCGCCAGCGCCACGAAGAAGCCGATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATaccaccg                                                                                                                                                                                                                                        <  2:117132/194‑2 (MQ=255)
                                                 aTGCTCATAATGATGACACCCGAGATTAAGGTGCGGAACAGATTTCCACGATGCACGGCGACCGCCATCGCCACGAAGAAGCCGATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCA‑‑CCGCCGTATGTcc                                                                                                                                                                                                                            >  1:117132/1‑194 (MQ=255)
                                                  tGCTCATAATGATGACACCCGAGATTAAGGTGCGGAACAGATTTCCACGATGCACGGCGACCGCCATCGCCACGAAGAAGCCGATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGAt                                                                                                                                                                                                           >  1:617003/1‑210 (MQ=255)
                                                              tgaCACCCGAGATTAAGGTGCGGAACAGATTTCCACGATGCACGGCGACCGCCATCGCCACGAAGAAGCCGATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCg                                                                                                                                                                       >  2:72284/1‑234 (MQ=255)
                                                                                tggggAAAAGATTTCCAGGATGCGCGGCGACCCCCATCGCCACGAAGACGCCGATGGTGGAAAGATCGCCAAACGGCAGCACCTGATTCCCCGGCCCACAAACACCAATAAAAATGGTGGGTGGGATAAAAATCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACg                                                                                                                                                       <  2:40264/229‑1 (MQ=255)
                                                                                  cGGAACAGATTTCCACGATGCACGGCGACCGCCATCGCCACGAAGAAGCCGATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATaccac                                                                                                                                                                                                                                          >  2:365806/1‑149 (MQ=255)
                                                                                  cGGAACAGATTTCCACGATGCACGGCGACCGCCATCGCCACGAAGAAGCCGATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATaccac                                                                                                                                                                                                                                          <  1:365806/149‑1 (MQ=255)
                                                                                          aTTTCCACGATGCACGGCGACCGCCATCGCCACGAAGAAGCCGATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTa                                                                                                                                             >  1:40264/1‑232 (MQ=255)
                                                                                           gtcccACGGAGCGCGGCGCCCGCAGTCGCCACAGAGAACCGGGAGGGTGCAGGATCGACAAGCTGCACACCTCGTTAGCCGGGCAGACAACCAGCATATAAAGTGTGGAGTGGTATAAAACCCCGCCTGGCGGAACCAACCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAAGCAGGAACCCCTGGCCGCCGAACTTCGCCTGTAAACGACTACGAGCCGGCTTAGCGAt                                                                                                                                        <  2:7039/233‑1 (MQ=255)
                                                                                              ccACGATGCACGGCGACCGCCATCGCCACGAAGAAGCCGATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAg                                                                                                                                                                                                       >  1:76398/1‑170 (MQ=255)
                                                                                                 cGATGCACGGCGACCGCCATCGCCACGAAGAAGCCGATGGTGGCAAGATCGCCAAACGGCAGCACCTGATGACCCGGCATACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACg                                                                                                                                                       >  1:849709/1‑215 (MQ=255)
                                                                                                             cccccattccccAGGAAGAGCCGGATGGTGGCAAGACCGCAAACCGGAGGCCCTTGTTACCCCGGCAACAAACAAGCATTAAAAAGGTGGGGGGGAAAAAAAATAGGGCTTCCGGTTCCCCCCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAACCATCCatgat                                                                                                                   <  2:156647/229‑1 (MQ=255)
                                                                                                              ccGCCATCGCCACGAAGAAGCCGATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGCACTCCTGACCGCCGAAATTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAACCATCCatgatg                                                                                                                  >  2:583596/1‑239 (MQ=255)
                                                                                                              ccGCCATCGCCACGAAGAAGCCGATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAACCATCCatgatg                                                                                                                  >  2:84335/1‑239 (MQ=255)
                                                                                                                      accggaaaaaagCCGAGGGTGGAAAAACCCCCAAACGGCGCCCCCTGATTCCCCGGCACACAAACAGCAATAAAATAGTGGGGGGGAAAAAAAACCGGCTTGGCCGATCCCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTAAt                                                                                                          <  2:316544/230‑1 (MQ=255)
                                                                                                                                     aTGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTAATCACCCGTGGc                                                                                                >  1:316544/1‑234 (MQ=255)
                                                                                                                                      tggtggCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTAATCACCCGTGGCATTagc                                                                                          >  1:211631/1‑239 (MQ=255)
                                                                                                                                        gtggCAAGATCGCCAAACGGCAGCACCTGATAACCCGGCACACAAACAGCAATTAAAAGGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAgg                                                                                                                                                                                                      <  2:837423/129‑1 (MQ=255)
                                                                                                                                                  cccaaaaCGGCAGCCCCTGATTACCCGGCACACAACCAGCAATTAAAAGGGTGGGTGGAAAAAAAACCGGCCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGc                                                                                                                                                                                                     <  2:432881/116‑1 (MQ=255)
                                                                                                                                                  cGCCAAACGGCTGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGa                                                                                                                                                                                            >  1:837423/1‑129 (MQ=255)
                                                                                                                                                   gCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTAATCACCCGTGGCATTAGCAGCATCACTgccg                                                                             <  1:84335/239‑1 (MQ=255)
                                                                                                                                                            cagcaCCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGaa                                                                                                                                                                                           >  1:432881/1‑120 (MQ=255)
                                                                                                                                                               caCCTGATTACCCGGCACACAAACAGCAATAAAATTGGTGAGTGGGATAAAAAGCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCcagcag                                                                                                                                                                                                                      <  2:344632/90‑1 (MQ=255)
                                                                                                                                                                    gATTCCCCGACACAAAAACAGCAATTAAAATGGGGAGGGGAAAAAAAACCAGGCTTGCGAATAC‑‑CCCCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTAATCACCCGTGGCATTAGCAGCATCACTGCCGCCGTTTTTACCGCCAGc                                                            <  1:544246/239‑1 (MQ=255)
                                                                                                                                                                          cccGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGa                                                                                                                                                                                                            >  1:344632/1‑90 (MQ=255)
                                                                                                                                                                           ccGGCACCAAAACAGAAATTAAAATGGTGGGTGGGATAAAAAGCCGGCTTGCCGGTACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTAATCACCCGTGGCATTAGCAGCATCACTGCCGCCGTTTTTACCGCCAGCTGCAAt                                                      <  1:226221/239‑1 (MQ=255)
                                                                                                                                                                               cacacaAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATCCCA‑‑CCGCCGTATGTCCcagcag                                                                                                                                                                                                                      <  2:109038/75‑1 (MQ=255)
                                                                                                                                                                               cacacaAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCcagcag                                                                                                                                                                                                                      >  1:109038/1‑75 (MQ=255)
                                                                                                                                                                                          cAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTAATCACCCGTGGCATTAGCAGCATCACTGTCGCCGTTTTTACTGCCAGCTGCACTACACCCTTGACATCg                                       >  2:758589/1‑239 (MQ=255)
                                                                                                                                                                                                      tGAGTGGGATAAAAATCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTAATCACCCGTGGCATTAGCAGCATCACTGCCGCCGTTTTTACCGCCAGCTGCAATACACCTTTGACATCGTAACCCGCGAgg                           >  1:182370/1‑239 (MQ=255)
                                                                                                                                                                                                        aGTGGGATAAAAATCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAAccat                                                                                                                          >  2:248538/1‑142 (MQ=255)
                                                                                                                                                                                                        aGTGGGATAAAAATCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAAccat                                                                                                                          <  1:248538/142‑1 (MQ=255)
                                                                                                                                                                                                             gaaaaaaaTCAGGCTTCCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTAATCACCCGTGGCATTAGCAGCATCACTGCCGCCGTTTTTACCGCCAGCTGCAATACACCTTTGACATCGTAACCCGCGAGGATGCCg                     <  1:319878/235‑1 (MQ=255)
                                                                                                                                                                                                                aaaaaTCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTaa                                                                                                                                                                >  2:3435/1‑96 (MQ=255)
                                                                                                                                                                                                                aaaaaTCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTaa                                                                                                                                                                <  1:3435/96‑1 (MQ=255)
                                                                                                                                                                                                                  aaaTCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTAATCACCCGTGGCATTAGCAGCATCACTGCCGCCGTTTTTACCGCCAGCTGCAATACACCTTTGACATCGTAACCCGCGAGGATGCCGATAATc               <  1:698336/239‑1 (MQ=255)
                                                                                                                                                                                                                        ggCTGGGCGATACCCCCCGCCTTATGCCCCAACAGCAACGCCGGATCAAGCCCAATGAGGACCCCCGGCCCGCCGAACTTCGCCTTTAAACGACTGCGAGCCTGCTTAGCGATGGGCGTTAACCCATCCATGATGGGTTTAATCGCCCGTGGCATTGGCAGCATCACTGCCGCCGTTTTTACCGCCAGCTGCAATACACCTTTGACGTCGTAACCCGCAAGGATGCCGATAATCAGccc          <  2:211631/239‑1 (MQ=255)
                                                                                                                                                                                                                              gCCGATACCA‑‑CTGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTAATCACCCGTGGCATTAGCAGCATCACTgccg                                                                             >  1:824006/1‑166 (MQ=255)
                                                                                                                                                                                                                              gCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTaaa                                                                                                                              >  1:704004/1‑117 (MQ=255)
                                                                                                                                                                                                                                  gaatccccccGCCGTCTGACCAGGCAACAGCGCGGGAGCAAGGCCAACAAGAACCTCCTGCCCGCCGAACTTCGCCTGTAAACGACTACTAGCCTGCTTAACGGTGGACGTTAAACCATCCATGATGGGTTTAATCCCCCATAACAATAGCAGCATCACTGCCGCCGTTTTTACCGCCAGCTGCAATACACCTTTGACATCGTAACCCGCGAGGATGCCGATAATCAGCCCCATCACaa  <  2:182370/233‑1 (MQ=255)
                                                                                                                                                                                                                                        ||                                                                                                                                                                                                                                       
CGCCAGTTGGGTGTGCAAACCAATAGTTTGCGTCGCGATCCACAGGGTGATGCTCATAATGATGACACCCGAGATTAAGGTGCGGAACAGATTTCCACGATGCACGGCGACCGCCATCGCCACGAAGAAGCCGATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCACCCCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTAATCACCCGTGGCATTAGCAGCATCACTGCCGCCGTTTTTACCGCCAGCTGCAATACACCTTTGACATCGTAACCCGCGAGGATGCCGATAATCAGCCCCATCACAA  >  NC_000913/2173131‑2173593

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