Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,173,361 Δ2 bp 100% pseudogene (1‑2/442 nt) gatC ← 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.1% 291.2 / 11.5 76pseudogene (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/6);  new base . (51/19);  total (51/25)
Fisher's exact test for biased strand distribution p-value = 8.10e-04
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
Rejected as polymorphism: Variant not supported by required number of reads on each strand.
Rejected as polymorphism: Polymorphic indel expands or contracts a homopolymer stretch.
*NC_0009132,173,3620C.92.3% 287.4 / 5.0 76pseudogene (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 . (51/19);  total (51/25)
Fisher's exact test for biased strand distribution p-value = 2.48e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
Rejected as polymorphism: Variant not supported by required number of reads on each strand.
Rejected as polymorphism: Polymorphic indel expands or contracts a homopolymer stretch.

ATTAGCCGCCAGTTGGGTGTGCAAACCAATAGTTTGCGTCGCGATCCACAGGGTGATGCTCATAATGA‑TGACACCCGAGATTAAGGTGCGGAACAGATTTCCACGATGCACGGCGACCGCCATCGCCACGAAGAAGCC‑GATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCACCCCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTAATCACCCGTGGCATTAGCAGCATCACTGCCGCCGTTTTTACCGCCAGCTGCAATACACCTTTGACATCGTAACCCGCGAGGATGCCGATAATCAGCCCCATCAC  >  NC_000913/2173125‑2173591
                                                                                                                                                                                                                                              ||                                                                                                                                                                                                                                     
aTTAGCCGCCAGTTGGGTGGGCAAACCAATAGTTTGCGTCGCGATCCACAGGGTGATGCTCATAATGA‑GGACCCCCGAGATTAAGGTGCGGAACAGATTTCCACGATGCACGGCGACCGCCATCGCCACGAAGAAGCC‑GATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATaccaccg                                                                                                                                                                                                                                      <  1:318121/239‑2 (MQ=255)
         cAGTTGGGTGTGCAAACCAATAGTTTGCGTCGCGATCCACAGGGTGATGCTCATAATGA‑TGACACCCGAGATTAAGGTGCGGAACAGATTTCCACGATGCACGGCGACCGCCATCGCCACGAAGAAGCC‑GATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCA‑‑CCGCCGTATGt                                                                                                                                                                                                                            >  1:1043797/1‑238 (MQ=255)
               ggTGTGCAAACCAATAGTTTGCGTCGCGATCCACAGGGTGATGCTCATAATGA‑TGACACCCGAGATTAAGGTGCGGAACAGATTTCCACGATGCACGGCGACCGCCATCGCCACGAAGAAGCC‑GATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCcagca                                                                                                                                                                                                                     >  1:910943/1‑239 (MQ=255)
                   tgCAAACCAATAGTTTGCGTCGCGATCCACAGGGTGATGCTCATAATGA‑TGACACCCGAGATTAAGGTGCGGAACAGATTTCCACGATGCACGGCGACCGCCATCGCCACGAAGAAGCC‑GATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAGTGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAgcagca                                                                                                                                                                                                                  >  1:299097/1‑238 (MQ=255)
                             tAGTTTGCGTCGCGATCCACAGGGTGATGCTCATAATGA‑TGACACCCGAGATTAAGGTGCGGAACAGATTTCCACGATGCACGGCGACCGCCATCGCCACGAAGAAGCC‑GATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCa                                                                                                                                                                                                       >  1:413786/1‑239 (MQ=255)
                                 ttGCGTCGCGATCCACAGGGTGATGCTCATAATGA‑TGACACCCGAGATTAAGGTGCGGAACAGATTTCCACGATGCACGGCGACCGCCATCGCCACGAAGAAGCC‑GATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGc                                                                                                                                                                                                   <  1:147214/239‑1 (MQ=255)
                                  tGCGTCGCGATCCACAGGGTGATGCTCATAATGC‑TGACACCCGAGATTAAGGTGCGGAACAGATTTCCACGCTGCACGGCGACCGCCATCGCCACGAAGAAGCC‑GATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAACCAGCAATTAAAATGGTGAGTGGGCTAAAAATCAAGCTTGCCGATACCA‑‑CCGCCGTATGTCACAGCAGCAACGCCGGATCAAGGcc                                                                                                                                                                                                  >  2:230544/1‑239 (MQ=255)
                                  tGCGTCGCGATCCACAGGGTGATGCTCATAATGA‑TGACACCCGAGATTAAGGTGCGGAACAGATTTCCACGATGCACGGCGACCGCCATCGCCACGAAGAAGCC‑GATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGcc                                                                                                                                                                                                  >  1:317096/1‑239 (MQ=255)
                                  tGCGTCGCGATCCACAGGGTGATGCTCATAATGA‑TGACACCCGAGATTAAGGTGCGGAACAGATTTCCACGATGCACGGCGACCGCCATCGCCACGAAGAAGCC‑GATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGc                                                                                                                                                                                                   >  2:785623/1‑238 (MQ=255)
                                    cGTCGCGATCCACAGGGTGATGCTCATAATGA‑TGACACCCGAGATTAAGGTGCGGAACAGATTTCCACGATGCACGGCGACCGCCATCGCCACGAAGAAGCC‑GATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCA‑‑CCGCCGTATGACCCAGCAGCAACGCCGGATCAAGGCCaa                                                                                                                                                                                                >  1:633417/1‑239 (MQ=255)
                                        gcgATCCACAGGGGGATGCTCATAATGA‑TGACCCCCGAGATTAAGGTGCGGAACAGATTTCCACGATGCACGGCGACCGCCATCGCCACGAAGAAGCC‑GATGGTGGAAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACAAAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAg                                                                                                                                                                                            <  1:401306/239‑1 (MQ=255)
                                         acagcaacaGGGTGATGCTCAAAATGA‑TGACACCCGAGATTAAGGTGCGGAACAGTTTTCCACGATGCACGGCGCCCCCCATCGCCACGAAGAAGCC‑GATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATCCCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAgg                                                                                                                                                                                           <  1:806369/233‑1 (MQ=255)
                                                   ggTGATGCTC‑TAAAGAGTGACCCCCGACATGAAGGGGCGAAAAAAATTCCCACGATGCACGTCGACCGCCATCGCAACGAAGAAGCC‑GATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGCCACACAAACAGCAATTAAAATGGGGAGTGGGATAAAAATAAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAgg                                                                                                                                                                                           <  2:738097/229‑1 (MQ=255)
                                                       aTGCTCATAATGA‑TGACACCCGAGATTAAGGTGCGGAACAGATTTCCACGATGCACGGCGACCGCCATCGCCACGAAGAAGCC‑GATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGAccgcc                                                                                                                                                                             >  1:124802/1‑239 (MQ=255)
                                                        tGCTCATAATGA‑TGACACCCGAGATTAAGGTGCGGAACAGATTTCCACGATGCACGGCGACCGCCATCGCCACGAAGAAGCC‑GATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGAccgccg                                                                                                                                                                            >  1:946745/1‑239 (MQ=255)
                                                        tGCTCATAATGA‑TGACACCCGAGATTAAGGTGCGGAACAGATTTCCACGATGCACGGCGACCGCCATCGCCACGAAGAAGCC‑GATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGAccgccg                                                                                                                                                                            >  1:877541/1‑239 (MQ=255)
                                                        tGCTCATAATGA‑TGACACCCGAGATTAAGGTGCGGAACAGATTTCCACGATGCACGGCGACCGCCATCGCCACGAAGAAGCC‑GATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGAccgccg                                                                                                                                                                            >  1:77787/1‑239 (MQ=255)
                                                        tGCTCATAATGA‑TGACACCCGAGATTAAGGTGCGGAACAGATTTCAACGATGCACGGCGACCGCCATCGCCACGAAGAAGCC‑GATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGAACACAAAGAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGACGATACCA‑‑CCGCCGTATGTCCCAGCAGCACAGCCGGATCCAGGACACTAAGGAACTCCTCAcccccc                                                                                                                                                                            >  2:806230/1‑238 (MQ=255)
                                                             aTAATGA‑TGACACCCGAGATTAAGGTGCGGAACAGATTTCCACGATGCACGGCGACCGCCATCGCCACGAAGAAGCC‑GATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGAc                                                                                                                                                                                 >  1:738097/1‑229 (MQ=255)
                                                               aaTGA‑TGACACCCGAGATTAAGGTGCGGAACAGATTTCCACGATGCACGGCGACCGCCATCGCCACGAAGAAGCC‑GATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAACTGATGAGTGGGATAAAAATCAGGCTTGCCGATACCA‑‑CCGCAGTATGTCCCAGCAGCAACGTCGGATCAAGGACAATCAGGAACTCCTGACCGCCGTACTTc                                                                                                                                                                      >  2:892948/1‑238 (MQ=255)
                                                               aaTGA‑TGACACCCGAGATTAAGGTGCGGAACAGATTTCCACGATGCACGGCGACCGCCATCGCCACGAAGAAGCC‑GATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTc                                                                                                                                                                      >  1:94884/1‑238 (MQ=255)
                                                               aaTGA‑TGACACCCGAGATTAAGGTGCGGAACAGATTTCCACGATGCACGGCGACCGCCATCGCCACGAAGAAGCC‑GATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAAAAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTc                                                                                                                                                                      >  1:635686/1‑238 (MQ=255)
                                                               aaTGA‑TGACACCCGAGATTAAGGAGCGGAACAGATTTCCACGCTGCACGGCGACCGCCATCGCCACGAAGAAGCC‑GATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTTCCGATACCA‑‑CCGCCGAATGTCCCAGCAGCAAGGCCGGATCAAGGCCCATCAGGAACTCCTGACCGCATAACTTc                                                                                                                                                                      >  2:851477/1‑238 (MQ=255)
                                                                         acCCGAGATTAAGGTGCGAAACAGATTCCCACGATGCGCGCCGACCGCGATCCCAACGAAGAAGCC‑GAGGGTGGAAAGATCGCCAAACGGCATCACCTGATTACCCGGCACACAAACAACAATAAAAATGGTGAGTGGGATAAAAATCGGGCTTGCCGCTCCCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAAc                                                                                                                                                            <  2:1084829/239‑1 (MQ=255)
                                                                             gagaTTAAGGTGCGGAACAGATTTCCACGATGCACGGCGACCGCCATCGCCACGAAGAAGCC‑GATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGCACAAAGCCAATCAGGAAATCCTGCCCGCCGACCTTCGCCTGTAAACGACt                                                                                                                                                        >  2:981360/1‑239 (MQ=255)
                                                                                         cGGAACAGATTTCCACGATGCACGGCGACCGCCATCGCCACGAAGAAGCC‑GATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATc                                                                                                                                                                                                        >  1:1018032/1‑179 (MQ=255)
                                                                                         cGGAACAGATTTCCACGATGCACGGCGACCGCCATCGCCACGAAGAAGCC‑GATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATc                                                                                                                                                                                                        <  2:1018032/179‑1 (MQ=255)
                                                                                         cGGAACAGATTTCCACGATGCACGGCGACCGCCATCGCCACGAAGAAGCC‑GATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCtt                                                                                                                                            >  1:170874/1‑239 (MQ=255)
                                                                                         cGGAACAGATTTCCACGATGCACGGCGACCGCCATCGCCACGAAGAAGCC‑GATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCtt                                                                                                                                            >  2:455909/1‑239 (MQ=255)
                                                                                          ggAACAGATTTCCACGATGCCCGGCGACCGCCATCGCCACGAAGAAGCC‑GATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCAGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTa                                                                                                                                           >  2:731316/1‑239 (MQ=255)
                                                                                          ggAACAGATTTCCACGATGCACGGCGACCGCCATCGCCACGAAGAAGCC‑GATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTa                                                                                                                                           >  1:894789/1‑239 (MQ=255)
                                                                                            aaCAGATTTCCACGATGCACGGCGACCGCCATCGCCACGAAGAAGCC‑GATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGAGAGTGGGATAAAAATCAGGCTTGCCGATACAA‑‑ACGCCGTATGTCCCAGCAGCAACGCCGGACCAAGGCCCATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGAATCCGACCCTTCTTAGc                                                                                                                                         >  2:746471/1‑239 (MQ=255)
                                                                                               caagtccagcGAGGCACGGCGCCCACCACCGCCCAGAAGAAGCC‑GATGGTGGAAAGACCGCCAAACGGCGCCACCTAATTCCCCGCCCCACAACCAGAAATAAAAATGTTGGGGGGAAAAAAAATAAGGCTGGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGAt                                                                                                                                      <  2:254655/230‑1 (MQ=255)
                                                                                               aGATTTCCACGATGCACGGCGACCGCCATCGCCACGAAGAAGCC‑GATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGAt                                                                                                                                      <  1:455909/239‑1 (MQ=255)
                                                                                               aGATTTCCACGATGCACGGCGACCGCCATCGCCACGAAGAAGCC‑GATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGAt                                                                                                                                      <  1:785623/239‑1 (MQ=255)
                                                                                                   ttCCACGATGCACGGCGACCGCCATCGCCACGAAGAAGCC‑GATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGc                                                                                                                                  >  1:495794/1‑239 (MQ=255)
                                                                                                      cACGATGCACGGCGACCGCCATCGCCACGAAGAAGCC‑GATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAAt                                                                                                                                                                                               >  1:851865/1‑175 (MQ=255)
                                                                                                           gggcgtcgagccgccaccgccaGGAAGAGACC‑GATGGGGGAAAGGTCGGAAAAGGGCAGCCCCTGATTCCCCGGCACACAAACATCAATTAAATGTGGGAGGGGGAAAAAAACCAGGCTTGAC‑AAACC‑CACCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTaaa                                                                                                                            <  2:1043797/227‑1 (MQ=255)
                                                                                                                ggCGACCGCCATCGCCACGAAGAAGCC‑GATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAAccatcca                                                                                                                     >  1:36084/1‑239 (MQ=255)
                                                                                                                 gCGACCGCCATCGCCACGAAGAAGCC‑GATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAAccatccat                                                                                                                    >  2:50034/1‑239 (MQ=255)
                                                                                                                     ccGCCATCGCCACGAAGAAGCC‑GATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAgg                                                                                                                                                                                           >  1:1047628/1‑164 (MQ=255)
                                                                                                                     ccGCCATCGCCACGAAGAAGCC‑GATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGGACTCCTGACCGCCGAACTTCCCCTGTAAACGACTACGAGCCTGCTTAGCGATGGCCGTTAAAACATCCatgatg                                                                                                                >  2:844334/1‑239 (MQ=255)
                                                                                                                     ccGCCATCGCCACGAAGAAGCC‑GATGGAGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAgg                                                                                                                                                                                           <  2:1047628/164‑1 (MQ=255)
                                                                                                                       gCCATCGCCACGAAGAAGCC‑GATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAACCATCCATGATggg                                                                                                              >  1:425522/1‑239 (MQ=255)
                                                                                                                        cctccGCCACGAA‑AAGCCTCATGTGGGTAAGTTCCCAAAACGGCAGCACCGAGTTACTCCACACACAACCGAAAATTAAAATGGGGAGGGGTAAAAAACACAGGTTGGCGGATAC‑‑CTCCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGAAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGt                                                                                                             <  2:2641/235‑1 (MQ=255)
                                                                                                                             gCCACGAAGAGGCC‑GATGGTGGCAAGACCGCCAAACGGCGGCACCGAATACCCCGGCACAAAAACACCAATAAAAATGGGGAGTGGGATAAAAACCAGGCGGGCCGAT‑CC‑CCCCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTAAt                                                                                                        <  2:36084/239‑1 (MQ=255)
                                                                                                                             gCCACGAAGAAGCC‑GATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATCCCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTAAt                                                                                                        <  1:167948/239‑1 (MQ=255)
                                                                                                                             gCCACGAAGAAGCC‑GATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGAT‑CC‑CCCCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTAAt                                                                                                        <  1:160579/239‑1 (MQ=255)
                                                                                                                                 cGAAGAAGCC‑AATGGTGCCAAAATCGCCAAACGGCAGACCCTG‑TTCCCCGCCACAAAAAAAGCAATAAAAATGGGTAGGGGGAAAAAAACCAGGCTTGCCAATAC‑‑CCCCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTAATCAccc                                                                                                   <  2:282876/239‑1 (MQ=255)
                                                                                                                                   aagaagCC‑GATGGTGGCAAGATCGCCAAACGGCAGCACCTGATGACCCGGCACAAAAACAGCAATTAAAAGGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTAATCAccc                                                                                                   <  1:196821/238‑1 (MQ=255)
                                                                                                                                    agaggCC‑GAGGTTGGCAAGACCGCAAAGCGTCAGCACCGGATTACCCGGCAAACAAACAGCAATAAAAATGGTGTGTGGGAAAAAAATCAGGCTTGCCAATCCCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTAATCACCCGt                                                                                                 <  1:782169/239‑1 (MQ=255)
                                                                                                                                    agaagCC‑GATGGTGGCAAGATCGCCAAACGGCAGCCCCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTAATCACCCGt                                                                                                 <  1:261118/239‑1 (MQ=255)
                                                                                                                                             aTGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTAATCAc                                                                                                     >  1:99229/1‑227 (MQ=255)
                                                                                                                                                           gCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCCACGCCGGATCAAGGCCAATCAGGAAATCCTGACCGACGAACTCCGCCTGTAAACGACTCCGAGCATACTTAGCGATGGGCGTTAAACGATCCAGGATGGGTTTAATCACCCGTAGCATTAGCAGATTCACTgccg                                                                           >  2:671629/1‑239 (MQ=255)
                                                                                                                                                                 cGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTAATCACCCGTGGCATTAGCAGCATCACTGCCGCCGttt                                                                     >  2:437444/1‑239 (MQ=255)
                                                                                                                                                                       caCCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTAATCACCCGTGGCATTAGCAGCATCACTGCCGCCGTTTTTAcc                                                                >  2:409662/1‑238 (MQ=255)
                                                                                                                                                                            gATTCCCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTAATCACCCGTGGCATTAGCAGCATCACTGCCGCCGTTTTTACCGCCAGc                                                          <  1:836615/239‑1 (MQ=255)
                                                                                                                                                                              ttccccGGCACACAAACAGCAATTAAAATGGTGAGGGGAATAAAAACCAGGCTTCCCAATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTAATCACCCGTGGCATTAGCAGCATCACTGCCGCCGTTTTTACCGCCAGCTg                                                        <  1:919440/236‑1 (MQ=255)
                                                                                                                                                                                   cGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTAATCACCCGTGGCATTAGCAGCATCACTGCCGCCGTTTTTACCGCCAGCTGCAATa                                                   <  1:409662/239‑1 (MQ=255)
                                                                                                                                                                                    ggCCCACAAACAGCACTTAAAATGGTGATTGGGATAAAAACCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCACTCAGGAACTCCTTACCGCCCACCTTCGCCTGTAAACGACTACGAGCCTCCTTAGCGATGGGCGTTAAACAATCCATGATGGGTTTAATCACCCGTGGCATTAACAGCATACCTGCCGCCGTTTTTAGCGCCAGCTGCattca                                                  >  2:318417/1‑235 (MQ=255)
                                                                                                                                                                                    ggCACACAACCAGAAATTAAAAGGTTGAGGGGAATAAAAAGCAGGTTTGCTGATAC‑‑CTCCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTAATCACCCGTGGCATTAGCAGCATCACTGCCGCCGTTTTTACCGCCAGCTGCAATac                                                  <  1:230544/239‑1 (MQ=255)
                                                                                                                                                                                    ggCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTAATCAccc                                                                                                   >  2:771856/1‑190 (MQ=255)
                                                                                                                                                                                    ggCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTAATCAccc                                                                                                   <  1:771856/190‑1 (MQ=255)
                                                                                                                                                                                    ggCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTAATCACCCGTGGCATTAGCAGCATCACTGCCGCCGTTTTTACCGCCAGCTGCAATac                                                  >  1:747035/1‑239 (MQ=255)
                                                                                                                                                                                    ggCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTAATCACCCGTGGCATTAGCAGCATCACTGCCGCCGTTTTTACCGCCAGCTGCAATa                                                   >  2:729731/1‑238 (MQ=255)
                                                                                                                                                                                                  aaTTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTAATCACCCGTGGCATTAGCAGCATCACTGCCGCCGTTTTTACCGCCAGCTGCAATACACCTTTGACATCGt                                    <  1:437444/239‑1 (MQ=255)
                                                                                                                                                                                                   aTTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTAATCACCCGTGGCATTAGCAGCATCACTGCCGCCGTTTTTACCGCCAGCTGCAATaca                                                 >  1:996684/1‑225 (MQ=255)
                                                                                                                                                                                                   aTTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTAATCACCCGTGGCATTAGCAGCATCACTGCCGCCGTTTTTACCGCCAGCTGCAATACACCTTTGACATCGTa                                   >  1:1132080/1‑239 (MQ=255)
                                                                                                                                                                                                         aTGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTAATCACCCGTGGCATTAGCAGCATCACTGCCGCCGTTTTTACCGCCAGCTGCAATACACCTTTGACATCGTAACCcgc                             >  1:582332/1‑239 (MQ=255)
                                                                                                                                                                                                            gTGAGTGGGATAAAAATCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTAATCACCCGTGGCATTAGCAGCATCACTGCCGCCGTTTTTACCGCCAGCTGCAATACACCTTTGACATCGTAACCCGCGAg                          >  1:696556/1‑239 (MQ=255)
                                                                                                                                                                                                              gAGTGGGATAAAAATCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTAATCACCCGTGGCATTAGCAGCATCACTGCCGCCGTTTTTACCGCCAGCTGCAATACACCTTTGACATCGTAACCCGCGAGGa                        >  1:938112/1‑239 (MQ=255)
                                                                                                                                                                                                                    gATAAAAATCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGCACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTAATCACCCGTGGCATTAGCAGCATCACTGCCGCCGTTTTTACCGCCAGCTGCAATACACCTTTGACATCGTAACCCGCGAGGATGCCGa                  >  1:358282/1‑239 (MQ=255)
                                                                                                                                                                                                                               ggCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTAATCACCCGTGGCATTAGCAGCATCACTGCCGCCGTTTTTACCGCCAGCTGCAATACACCTTTGACATCGTAACCCGCGAGGATGCCGATAATCAGcccc       >  2:564164/1‑239 (MQ=255)
                                                                                                                                                                                                                                   tGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTAATCACCCGTGGCATTAGCAGCATCACTGCCGCCGTTTTTACCGCCAGCTGCAATACACCTTTGACATCGTAACCCGCGAGGATGCCGATAATCAGCCCCATc    >  1:416101/1‑238 (MQ=255)
                                                                                                                                                                                                                                   tGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTAATCACCCGTGGCATTAGCAGCATCACTGCCGCCGTTTTTACCGCCAGCTGCAATACACCTTTGACATCGTAACCCGCGAGGATGCCGATAATCAGCCCCATc    >  1:48562/1‑238 (MQ=255)
                                                                                                                                                                                                                                    gCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGCACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTAATCACCCGTGGCATTAGCAGCATCACTGCCGCCGTTTTTACCGCCAGCTGCAATACCCCTTTGACATCTTAACCCGCGAGGATGCCGAGAATCAGCCCcaccac  >  2:645118/1‑239 (MQ=255)
                                                                                                                                                                                                                                    gCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTAATCACCCGTGGCATTAGCAGCATCACTGCCGCCGTTTTTACCGCCAGCTGCAATACACCTTTGACATCGTAACCCGCGAGGATGCCGATAATCAGCCCCATcac  >  1:966573/1‑239 (MQ=255)
                                                                                                                                                                                                                                              ||                                                                                                                                                                                                                                     
ATTAGCCGCCAGTTGGGTGTGCAAACCAATAGTTTGCGTCGCGATCCACAGGGTGATGCTCATAATGA‑TGACACCCGAGATTAAGGTGCGGAACAGATTTCCACGATGCACGGCGACCGCCATCGCCACGAAGAAGCC‑GATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCACCCCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTAATCACCCGTGGCATTAGCAGCATCACTGCCGCCGTTTTTACCGCCAGCTGCAATACACCTTTGACATCGTAACCCGCGAGGATGCCGATAATCAGCCCCATCAC  >  NC_000913/2173125‑2173591

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