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.88.9% 129.6 / 6.0 36pseudogene (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 (1/3);  new base . (21/11);  total (22/14)
Fisher's exact test for biased strand distribution p-value = 2.77e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
*NC_0009132,173,3620C.94.1% 130.5 / 2.2 36pseudogene (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 (1/1);  new base . (21/11);  total (22/14)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CCTGATCCATTGAAGCCACCATACCCCCGGCTTTTAACGCCCCAGCATTAGCCGCCAGTTGGGTGTGCAAACCAATAGTTTGCGTCGCGATCCACAGGGTGATGCTCATAATGATGACACCCGAGATTAAGGTGCGGAACAGATTTCCACGATGCACGGCGACCGCCATCGCCACGAAGAAGCCGATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCACCCCGCCGTATGTCCCAGCAGCAACGCCGGATC‑AAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCT‑TAGCGATGGGCGTTAAACCATCCATGATGGGTTTAATCACCCGTGGCATTAGCAGCATCACTGCCGCCGTTTTTACCGCCAGCTGCAATACACCTTTGACATCGTAACCCGCGAGGATGCCGATAATCAGCCCCATCACAAAACCCACGGTGACAGGCTCGCCAAATGGACCAAATTTGCGCTGAATATCGTCGGCG  >  NC_000913/2173079‑2173649
                                                                                                                                                                                                                                                                                          ||                                                                                                                                                                                                                                                                                                 
ccTGATCCATTGAAGCCACCATACCCCCGGCTTTTAACGCCCCAGCATTAGCCGCCAGTTGGGTGTGCAAACCAATAGTTTGCGTCGCGATCCACAGGGTGATGCTCATAATGATGACACCCGAGATTAAGGTGCGGAACAGATTTCCACGATGCACGGCGACCGCCATCGCCACGAAGAAGCCGATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCA‑‑CCGCCGTATCTCCcagcag                                                                                                                                                                                                                                                                                >  2:236174/1‑301 (MQ=255)
                       cccccGGCTTTTAACGCCCCAGCATTAGCCGCCAGTTGGGTGTGCAAACCAATAGTTTGCGTCGCGATCCACAGGGTGATGCTCATAATGATGACACCCGAGATTAAGGTGCGGAACAGATTTCCACGATGCACGGCGACCGCCATCGCCACGAAGAAGCCGATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACACCCAGCAATTAAACTGGTGAGTGGGATAAAACTCAGCCTTGCCGATACCA‑‑CCACCGTCTGGCCAAGCCGCACCGCCGGCAC‑AAGGCCACTc                                                                                                                                                                                                                                                         >  1:3906/1‑300 (MQ=255)
                          cccggttttAAAGACCCGGCATTAGGCGCCAGTTGGGTGTGCAAACCAATAGGTTGGGTCGCGATCCACAGGGTGATGCTCATAATGATGACACCCGAGATTAAGGTGCGGAACAGATTTCCACGATGCACGGCGACCGCCATCGCCACGAAGAAGCCGATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATC‑AAGGCCAATCAGGa                                                                                                                                                                                                                                                     <  1:256177/296‑1 (MQ=255)
                                     cGCCCCAGCATTAGCCGCCAGTTGGGTGTGCAAACCAATAGTTTGCGTCGCGATCCACAGGGTGATGCTCATAATGATGACACCCGAGCTTAAGGTGCGGAACAGATTTCCACGATGCACGGCGACCGCCATCGCCACGAAGCAGCCGACGGTGGCAAGATCGCCAAACGGCAGCACCTGCTTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAACTCATGCTTGCCGATCCCACCC‑‑CCGTTTTTCCCAGCAGCAACCCCGGCAC‑CAGCCCTACCACGCACTCCTGAcc                                                                                                                                                                                                                                           >  1:14001/1‑300 (MQ=255)
                                      gCCCCAGCATTAGCCGCCAGTTGGGTGTGCAAACCAATAGTTTGCGTCGCGATCCACAGGGTGATGCTCATAATGATGACACCCGAGATTAAGGTGCGGAACAGATTTCCACGATGCACGGCGACCGCCATCGCCACGAAGAAGCCGATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGGGGGATAAAAATCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCACAGCCGGATC‑AAGGCCAATCAGGAACCCCTGAccg                                                                                                                                                                                                                                          >  1:65266/1‑300 (MQ=255)
                                         cgagaagtagaagagaGTTGTGTGTGCAAACCAAATGTTTGCGTAGGGATCGCAAGGGAGATGCTGATAATGATGACACCCGGGATTAAGCGGCGGAACAGAATGCCACGATGCACGGCGACCGCCATCGCCACGAAGAAGCCGATGGTGGCAAGAACGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGGTACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATC‑AAGGCCAATCAGGAACTCCTGACCGCCGa                                                                                                                                                                                                                                      <  1:183716/286‑1 (MQ=255)
                                                  gccgccAGTTGGGTGTGCAAACCAATAGTTTGCGTCGCGATCCACAGGGTGATGCTCATAATGATGACACCCGAGATTAAGGTGCGGAACAGATTTCCACGATGCACGGCGACCGCCATCGCCACGAAGAAGCCGATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATC‑AAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGc                                                                                                                                                                                                                               >  2:291976/1‑299 (MQ=255)
                                                             ggTGTTCAAACCAATAGTGTGCGTCGCGATCCACAGGGTGATGCTCAGAATGATGACAACCGAGATTAAGGTGGGGAACAGATTTCCACGATGAACGGCTACCGCCATCGCCACGAAGAAGCCGATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATC‑AAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACt                                                                                                                                                                                                                   <  1:236166/300‑1 (MQ=255)
                                                                                        gATCCACAGGGTGATGCTCATAATGATGACACCCGAGATTAAGGTGCGGAACAGATTTCCACGATGCACGGCGACCGCCATCGCCACGAAGAAGCCGATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATC‑AAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCT‑TAGCGATGGGCGTTaaa                                                                                                                                                                                      >  2:53642/1‑301 (MQ=255)
                                                                                                           aTAATGATGACACCCGAGATTTAGGTGCGGAACAGATTTCCACGATGCACGGCGACCGCCATCGCCACGAAGAAGCCGATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAACATCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATC‑AAGGCCACTCAGGAACCCCTGACCGCCGCACTTCGCCTGTAAACGAACACGAGCC‑GCTCTAGCAATCGACGTTAAAcaacacactctgggttt                                                                                                                                                                     >  1:222668/1‑283 (MQ=255)
                                                                                                             aaTGATGACACCCGAGATTAAGGTGCGGAACAGATTTCCACGATGCACGGCGACCGCCATCGCCACGAAGAAGCCGATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATC‑AAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCT‑TAGCGATGGGCGTTAAACCATCCATGATGGGTTTAAt                                                                                                                                                                  >  2:264147/1‑300 (MQ=255)
                                                                                                                gatggCACCCCAGATTAAGGGTCGAAACAGATTTCCACGATGCACGGCGACCGCCATCGCCACGAAGAAGCCGATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATC‑AAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCT‑TAGCGATGGGCGTTAAAcccatcatggatggtttaacgcc                                                                                                                                                               <  1:183294/300‑22 (MQ=255)
                                                                                                                       cccGAGATTAAGGGGCGGAACAGATTTCCACGATGCACGGCGACCGCCATCGCCACGAAGAAGCCGATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATC‑AAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCT‑TAGCGATGGGCGTTAAACCATCCATGATGGGTTTAATCACCCGTGGCa                                                                                                                                                       <  1:264139/301‑1 (MQ=255)
                                                                                                                                     gCGGAACAGATTTCCACGATGCACGGCGACCGCCATCGCCACGAAGAAGCCGATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATC‑AAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCT‑TAGCGATGGGCGTTAAACCATCCATGATGGGTTTAATCACCCGTGGCATTAGCAGCATCAc                                                                                                                                          <  2:206787/300‑1 (MQ=255)
                                                                                                                                                       aTGCACGGCGACCGCCATCGCCACGAAGAAGCCGATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATC‑AAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCT‑TAGCGATGGGCGTTAAACCATCCATGATGGGTTTAATCACCCGTGGCATTAGCAGCATCACTGCCGCCGTTTTTACCGcc                                                                                                                       <  2:65268/301‑1 (MQ=255)
                                                                                                                                                            cGGCGACCGCCATCGCCACGAAGAAGCCGATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATC‑AAGGCCAATCAGGAACTCCTGCCCTCCGCACTTCGCCTATAAACGACTCCGAGCCTGCT‑TAGCGCTTGCCGTTAAACCACCCATGATGGGTTTAATCACCCGTTGCATTAGCAGCATCTCTTACGCCGTTTTTACCGCCaccc                                                                                                                   >  1:223966/1‑297 (MQ=255)
                                                                                                                                                                          gCCACGAAGACGCCGATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATC‑AAGGCCAATCAGGAACTCCTGACCGCCGTACTTCGCCTGTAAACGACTACGAACCTGCT‑TAGCGCTGGTCGTTAAACCATCCATGATTGGTTTAATCTCTCGTGGCATTCGATGCATCACTGCCCCCGTTTTTACGTGCAGCCGCAATAACCCTTTgt                                                                                                    >  1:285531/1‑300 (MQ=255)
                                                                                                                                                                           ccACGAAGAAGCCGATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATC‑AAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCT‑TAGCGATGGGCGTTAAACCATCCATGATGGGTTTAATCACCCGTGGCATTAGCAGCATCACTGCCGCCGTTTTTACCGCCAGCTGCAATACACCttcgc                                                                                                    >  2:5773/1‑297 (MQ=255)
                                                                                                                                                                           ccACGAAGAAGCCGATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATC‑AAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACAACTACGAGCCTGCT‑TAGCGATGGGCGTTAAACCATCCATGATGGGTTTAATCACCCGTGGCATTAGCAGCATCACTGCCGCCGTTTTTACCGCCAGCTGCAATACACCTTTGa                                                                                                    <  2:316702/300‑1 (MQ=255)
                                                                                                                                                                            cACGAAGAAGCCGATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATC‑AAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCT‑TAGCGATGGGCGTTAAACCATCCATGATGGGTTTAATCACCCGTGGCATTAGCAGCATCACTGCCGCCGTTTTTACCGCCAGCTGCAATACACCTTTGACa                                                                                                  >  2:281026/1‑301 (MQ=255)
                                                                                                                                                                              cGAAGAAGCCGATGGTGGCAAGATCGCCAAACGGCAGCACCTGATCACCCGGCACACAAACAGCAATCAAAATGGTGAGTGGGATAAAAATCAGGCTTGCTGTTACGA‑‑CCGCCGTATGTTCCAGCAGCAACACCGGATC‑CAGGCCAATCAGGAACTCATGACCTCCGAACTTCGCCTGTAAACGACTACGAGCCTGCT‑TACCGTTGTGCGTTAAACCATACATCATGGTTTTAATCCCTCGTGGTTTTGGCTGCAGCACTGCCCCCATTTTTACCCTCAGcacaacaccactttttacac                                                                                                 >  1:204171/1‑281 (MQ=255)
                                                                                                                                                                                       cGATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATC‑AAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCT‑TAGCGATGGGCGTTAAACCATCCATGATGGGTTTAATCACCCGTGGCATTAGCAGCATCACTGCCGCCGTTTTTACCGCCAGCTGCAATACACCTTTGACATCGTAACCcg                                                                                        >  2:307356/1‑300 (MQ=255)
                                                                                                                                                                                              tggcgtaagatcAAACGGCAGGACCTGTTGACCCGAGACGCAAACGGCGCTTAATAAGTGGAGGGGGATAAAGGTAAGGCTGGCC‑ATAAC‑CACCACAGTAGGTCCCAGAAGCAACTCCGCATCAAACGCCAGT‑AGGATCGCCTGACCGCCGAGCTTCGATTGTACACGTCTCCAAGCGTGCT‑TGGCGATGGGCTTAAGACCATCCCTGATGGGTTTAATCACCCGTTGCATTAGCAGCATCGCTGCCACCGTTTTTACCGCCAGCTGCAATACACCTTTGACATCGTAACCCGCAAgggt                                                                                 <  1:275777/289‑3 (MQ=255)
                                                                                                                                                                                                   gATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATC‑AAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCT‑TAGCGATGGGCGTTAAACCATCCATGATGGGTTTAATCACCCGTGGCATTAGCAGCATCACTGCCGCCGTTTTTACCGCCAGCTGCAATACACCTTTGACATCGTAACCCGCGATGATGCCg                                                                             >  2:274582/1‑299 (MQ=255)
                                                                                                                                                                                                           cttacgacttgcgaaTTAACCGGCCAAAAACCAGCAATGTAAAGGGTGAGGTGGATAAAAACCAGGCTTGCCGATAC‑‑CGCCGCCGTAAGTCCCAGCAGCAACGCCGGATC‑AAGGCCAATCAGGAACTCCTGACCCCCGAACTTCGCCTGTAAACGACTACGAGCCTGCT‑TAGCGATGGGCGTTAAACCATCCATGATGGGTTTAATCACCCGTGGCATTAGCAGCATCACTGCCGCCGTTTTTACCGCCAGCTGCAATACACCTTTGACATCGTAACCCGCGAGGATGCCGATAATCAGcc                                                                   <  1:288293/287‑1 (MQ=255)
                                                                                                                                                                                                                   caCCTGATTAACCGGCACACAAACAGCAATTAAAAGGTTGAGTGGGATAAAAATCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATC‑AAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCT‑TAGCGATGGGCGTTAAACCATCCATGATGGGTTTAATCACCCGTGGCATTAGCAGCATCACTGCCGCCGTTTTTACCGCCAGCTGCAATACACCTTTGACATCGTAACCCGCGAGGATGCCGATAATCAGCCCCATca                                                             <  2:272529/299‑1 (MQ=255)
                                                                                                                                                                                                                     ccTGATTACCCGGCACACAAACAGTAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATC‑AAGGCCAATCAGGAACTCCTGACTGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCT‑TAGCGATGGGCGTTAAACCATCCATGATGGGTTTAATCACCCGTGGCATTAGCAGCATCACTGCCGCCGATTGTACCGCCAGCAGCAATACACATTTGACATCGTAACCCGCGAGGATGCCGATAAACAGcaagatcggaag                                                         >  2:288301/1‑290 (MQ=255)
                                                                                                                                                                                                                        gATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATC‑AAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCT‑TAGCGATGGGCGTTAAACCATCCATGATGGGTTTAATCACCCGTGGCATTAGCAGCATCACTGCCGCCGTTTTTACCGCCAGCTGCAATACACCTTTGACATCGTAACCCGCGAGGATGCCGATAATCAGCCCCATCACAAAAc                                                       <  2:3907/300‑1 (MQ=255)
                                                                                                                                                                                                                                   acaaaAACAGCAATTAAAATGGTGAGGGGGATAAAAATCAGGCCTGCCAATAC‑‑CCCCGCCGTATGTCCCAGCAGCAACGCCGGATC‑AAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCT‑TAGCGATGGGCGTTAAACCATCCATGATGGGTTTAATCACCCGTGGCATTAGCAGCATCACTGCCGCCGTTTTTACCGCCAGCTGCAATACACCTTTGACATCGTAACCCGCGAGGATGCCGATAATCAGCCCCATCACAAAACCCACGGTGACAg                                           <  1:281018/301‑1 (MQ=255)
                                                                                                                                                                                                                                    cacaAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATC‑AAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCT‑TAGCGATGGGCGCTAAACCATCCATGATGGGTTTAATCACCCGTGGCACTAGCAGCATCACTGCCGCCGTTTTTACCGCCAGCTGCAATACACCTTTGACATCGTAACCCCCGAGGATGCCGATAATCAGCCCCATCACAAAACCCACGGTGACAgg                                          >  2:21169/1‑301 (MQ=255)
                                                                                                                                                                                                                                        aaCAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATC‑AAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCT‑TAGCGATGGGCGTTAAACCATCCATGATGGGTTTAATCACCCGTGGCATTAGCAGCATCACTGCCGCCGTTTTTACCGCCAGCTGCAATACACCTTTGACATCGTAACCCGCGAGCATGCCGATAATCAGCCCCATCACACAACCCACGGTCACAGCCTc                                       >  2:33760/1‑300 (MQ=255)
                                                                                                                                                                                                                                               aTTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATC‑AAGTCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCT‑TAGCGATGGGCGTTAAACCATCCATGATGGGTTTAATCACCCGTGGCATTAGCAGCATTACTGCCGCCGTTTCTACCGCCAGCTGCAATACACCTTTGACATCGTAACCCGCGAGGATGCCGATAATCAGTCTCATCACCAAACCCACGGTTACAGGCCCGCCAAATg                               >  2:257901/1‑301 (MQ=255)
                                                                                                                                                                                                                                                               ggATAAAAATCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATC‑AAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCT‑TAGCGATGGGCGTTAAACCATCCATGATGGGTTTAATCACCCGTGGCATTAGCAGCATCACTGCCGCCGCTTTTACCGCCAGAGGCAATACAACTTTGCCATCGTAACCCGCAAGGATGCCGACACTCAGACGCATCACAAAAACCAAGGGAACAGACCGCCCAAATGTACCAAATTTGCGCta               >  1:309101/1‑300 (MQ=255)
                                                                                                                                                                                                                                                                 aTAAAAATCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATC‑AAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCT‑TAGCGATGGGCGTTAAACCATCCATGATGGGTTTAATCACCCGTGGCCTTAGCAGCATCACTGCCGCCGTTTTTACCGCCAGCTGCAATACACCTTTGACATCGTAACCCGCGAGGATGCCGATAATCAGCCCCATCAAAAAACCCCCGGTGACAGGCCCGCACAATGGACACAAttgtcgccgat             >  1:259756/1‑292 (MQ=255)
                                                                                                                                                                                                                                                                      aaTCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATC‑AAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCT‑TAGCGATGGGCGTTAAACCATCCATGATGGGTTTAATCACCCGTTGCATTAGCAGCATCACTGCCGCCGTTTTTACCGCCAGCCGCCATACACCTTTGACATCGTAACCCGCGAGGATGCCGATCATCAGCCCCATCACAAACCCCCCAGTGACAGGCTCCCCAAATGAACAAAATTTCCCCCAAATAtca        >  1:156544/1‑300 (MQ=255)
                                                                                                                                                                                                                                                                             cTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATC‑AAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCT‑TAGCGATGGGCGTTAAACCATCAATGATGGGTTTAATAACCAGTGGCATTAGCAGCATCACTGCCGCCGTTTTTACCGCCAGCTGCAATACACCTTTGACATCGTAACCCGCGAGGATGCCGATAATCAGCCCCATCACAAAACCCACGGTGACAGGCTCGCCAAATGGACCAAATTTGCGCTGAATATCGTCGgcg  <  2:14002/300‑1 (MQ=255)
                                                                                                                                                                                                                                                                                          ||                                                                                                                                                                                                                                                                                                 
CCTGATCCATTGAAGCCACCATACCCCCGGCTTTTAACGCCCCAGCATTAGCCGCCAGTTGGGTGTGCAAACCAATAGTTTGCGTCGCGATCCACAGGGTGATGCTCATAATGATGACACCCGAGATTAAGGTGCGGAACAGATTTCCACGATGCACGGCGACCGCCATCGCCACGAAGAAGCCGATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCACCCCGCCGTATGTCCCAGCAGCAACGCCGGATC‑AAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCT‑TAGCGATGGGCGTTAAACCATCCATGATGGGTTTAATCACCCGTGGCATTAGCAGCATCACTGCCGCCGTTTTTACCGCCAGCTGCAATACACCTTTGACATCGTAACCCGCGAGGATGCCGATAATCAGCCCCATCACAAAACCCACGGTGACAGGCTCGCCAAATGGACCAAATTTGCGCTGAATATCGTCGGCG  >  NC_000913/2173079‑2173649

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