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.85.9% 244.3 / 28.8 71pseudogene (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/10);  new base . (44/17);  total (44/27)
Fisher's exact test for biased strand distribution p-value = 1.83e-05
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.85.9% 238.9 / 24.5 71pseudogene (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/10);  new base . (44/17);  total (44/27)
Fisher's exact test for biased strand distribution p-value = 1.83e-05
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.

CAGTTGGGTGTGCAAACCAATAGTTTGCGTCGCGATCCACAGGGTGATGCTCATAATGATGACACCCGAGATTAAGGTGCGGAACAGATTTCCACGATGCACGGCGACCGCCATCGCCACGAAGAAGCCGATGGTGGCAAGATCGCCAAACGG‑CAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCACCCCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCC‑AATCAGG‑AACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTAATCACCCGTGGCATTAGCAGCATCACTGCCGCCGTTTTTACCGCCAGCTGCAATACACCTTTGACATCGTAACCCGCGAGGATGCCGATAATCAGCCC  >  NC_000913/2173134‑2173585
                                                                                                                                                                                                                                    ||                                                                                                                                                                                                                                 
cAGTTGGGTGTGCAAACCAATAGTTTGCGTCGCGATCCACAGGGTGATGCTCATAATGATGACACCCGAGATTAAGGTGCGGAACAGATTTCCACGATGCACGGCGACCGCCATCGCCACGAAGAAGCCGATGGTGGCAAGATCGCCAAACGG‑CAGCACCTGATTACCCGGCACAAAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCA‑‑CCGCCGTATGt                                                                                                                                                                                                                        >  1:963581/1‑238 (MQ=255)
          tgCAAACCAATAGTTTGCGTCGCGATCCACAGGGTGATGCTCATAATGATGACACCCGAGATTAAGGTGCGGAACAGATTTCCACGATGCACGGCGACCGCCATCGCCACGAAGAAGCCGATGGTGGCAAGATCGCCAAACGG‑CAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAgcagca                                                                                                                                                                                                              >  2:1028834/1‑238 (MQ=255)
             aaaCCAATAGTTTGCGTCGCGACCCCCGGGGTGATGCTAATAATGATGACCCCCGAGATTAAGGTGCGGAACAGATTTCCACGATGCACGGCGACCGCCATCGCCACGAAGAAGCCGATGGTGGCAAGATCGCCAAACGG‑CAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGGGTGGGATAAAAATCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGc                                                                                                                                                                                                          <  2:450376/239‑1 (MQ=255)
             aaaCCAAGAGTTGGCGTCGCGATCAACAGGGTGATGCTAATAATGATGCCACCCGGGATTAGGGTGCGGAACAGATTCCCCGGAGCCACGGCGACCGCCATCGCCACGAAGAAGCCGATGGTGGCAAGATCGCCAAACGG‑CAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGc                                                                                                                                                                                                          <  2:532374/239‑1 (MQ=255)
              acccAATAGTTTGCGTCGCGATCCACAGGGTGATGCTCATAGTGATGACACCCGAGATTAAGGTGCGGAACAGATTCCCACGATGCACGGCGACCGCCATCGCCACGAAGAAGCCGATGGTGGCAAGATCGCCAAACGG‑CAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGcc                                                                                                                                                                                                         <  1:379959/237‑1 (MQ=255)
              aaCCAATAGTTTGCGTCGCGATCCACAGGGTGATGCTCATAATGATGACACCCGAGATTAAGGTGCGGAACAGATTTCCACGATGCACGGCGACCGCCATCGCCACGAAGAAGCCGATGGTGGCAAGATCGCCAAACGG‑CAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGcc                                                                                                                                                                                                         <  1:951661/239‑1 (MQ=255)
                         tGCGTCGCGATCCACAGGGTGATGCTCATAATGATGACACCCGAGATTAAGGTGCGGAACAGATTTCCACGATGCACGGCGACCGCCATCGCCACGAAGAAGCCGATGGTGGCAAGATCGCCAAACGG‑CAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGAGAGAGGGATAAAAATCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGcc                                                                                                                                                                                              >  1:1106778/1‑239 (MQ=255)
                                cgATCCACAGGGTGATGCTCATAATGATGACACCCGAGATTAAGGTGCGGAACAGATTTCCACGATGCACGGCGACCGCCATCGCCACGAAGAAGCCGATGGTGGCAAGATCGCCAAACGG‑CAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCC‑AATCAgg                                                                                                                                                                                      >  2:1020035/1‑239 (MQ=255)
                                       caGGGTGATGCTCATAATGATGACACCCGAGATTAAGGTGCGGAACAGATTTCCACGATGCACGGCGACCGCCATCGCCACGAAGAAGCCGATGGTGGCAAGATCGCCAAACGG‑CAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCC‑AATCAGG‑AACt                                                                                                                                                                                 >  1:549773/1‑236 (MQ=255)
                                                    aTAATGATGACACCCGAGATTAAGGTGCGGAACAGATTTCCACGATGCACGGCGACCGCCATCGCCACGAAGAAGCCGATGGTGGCAAGATCGCCAAACGG‑CAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCA‑‑CCGCCGTCTGTCCCAGCAGCAACGCCGGATCAAGGCC‑AATCAGG‑AACACCTGACCGCCAAACt                                                                                                                                                                  >  1:306245/1‑238 (MQ=255)
                                                    aTAATGATGACACCCGAGATTAAGGTGCGGAACAGATTTCCACGATGCACGGCGACCGCCATCGCCACGAAGAAGCCGATGGTGGCAAGATCGCCAAACGG‑CAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCC‑AATCAGG‑AACTCCTGACCGCCGAACt                                                                                                                                                                  >  1:934908/1‑238 (MQ=255)
                                                          atgaCACCCGAGATTAAGGTGCGGAACAGATTTCCACGATGCACGGCGACCGCCATCGCCACGAAGAAGCCGATGGTGGCAAGATCGCCAAACGG‑CAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCC‑AATCAGG‑AACTCCTGACCGCCGAACTTCg                                                                                                                                                               >  1:211875/1‑235 (MQ=255)
                                                           cgccccgcGAGTTGAAGTTGCGGAATAAGTTCCCGCGATGCACGGCACCGGCAATCGCCACGAAGAAGCCTATGGTGGCAAAATCGCCAAGGGG‑CGGCACCGGATTCCCGGGCACAAAAACAGCAATTAAAATGGTGAGGGGGATAAAAATCAGGCTGGCAGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCC‑AATCAGG‑AACTCCTGACCGCCGAACTTCGCCTGt                                                                                                                                                          <  2:585287/232‑1 (MQ=255)
                                                            gaCACCCGAGATTAAGGTGCGGAACAGATTTCCACGATGCACGGCGACCGCCATCGCCACGAAGAAGCCGATGGTGGCAAGATCGCCAAACGG‑CAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCC‑AATCAGG‑AACTCCTGACCGCCGAACTTCGCCTGTa                                                                                                                                                         >  2:379762/1‑239 (MQ=255)
                                                              cacCCGAGATTAAGGTGCGGAACAGATTTCCACGATGCACGGCGACCGCCATCGCCACGAAGAAGCCGATGGTGGCAAGATCGCCAAACGG‑CAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCC‑AATCAGG‑AACTCCTGACCGCCGAACTTCGCCTGTaa                                                                                                                                                        >  1:437893/1‑238 (MQ=255)
                                                              cacCCGAGATTAAGGTGCGGAACAGATTTCCACGATGCACGGCGACCGCCATCGCCACGAAGAAGCCGATGGTCGCAAGATCGCCAAACGG‑CAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGAGTGATAAAAATCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCAAGCATCAACGCCGGATCAAGGCA‑AATCAGG‑ACATCCTGACCGCCGAACTTCGCCCGaaaa                                                                                                                                                       >  2:801139/1‑239 (MQ=255)
                                                                     gataagggggaggAAAAAATTCCCACGATCCACGGGGACCGCCATCGCAACGAAGAAGCCGATGGGGGCAAGACCGCAAAACGG‑CAGCACCTGATTCCCCGGCACACAAACAGCAATAAAAAGGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCC‑AATCAGG‑AACTCCTGACCGCCGAACTTCGCCTGTAAACGACTa                                                                                                                                                 <  2:437893/227‑1 (MQ=255)
                                                                         agggTGCGGAACAGTTTTCCAGGAGGCACGGCGACCGCCATCCGCACGAAGAGGCGGAGGGGGGCAAGATCGCAAAACGC‑CAGCCCCGGATGCCCCCGCACAAAAACAGCAATTAAAATGGTGGGTGGGATAAAAACCAGGCTTGCCGATAC‑‑CCCCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCC‑AATCAGG‑AACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGc                                                                                                                                            <  2:175186/237‑1 (MQ=255)
                                                                         aaGGTGCGGAACAGATTTCCACGATGCACGGCGACCGCCATCGCCACGAAGAAGCCGATGGTGGCAAGATCGCCAAACGG‑CAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCAGGATCAAGGCC‑ACTCAGG‑AACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGaa                                                                                                                                             >  2:701346/1‑237 (MQ=255)
                                                                               cGGAACAGATTTCCACGATGCACGGCGACCGCCATCGCCACGAAGAAGCCGATGGTGGCAAGATCGCCAAACGG‑CAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCC‑AATCAGG‑AACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCtt                                                                                                                                      >  1:1178340/1‑239 (MQ=255)
                                                                               cGGAACAGATTTCCACGATGCACGGCGACCGCCATCGCCACGAAGAAGCCGATGGTGGCAAGATCGCCAAACGG‑CAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCC‑AATCAGG‑AACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCtt                                                                                                                                      >  1:421052/1‑239 (MQ=255)
                                                                                ggAACAGATTTCCACGATGCACGGCGACCGCCATCGCCACGAAGAAGCCGATGGTGGCAAGATCGCCAAACGG‑CAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCC‑AATCAGG‑AACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTa                                                                                                                                     >  1:10624/1‑239 (MQ=255)
                                                                                 gAACAGATTTCCACGATGCACGGCGACCGCCATCGCCACGAAGAAGCCGATGGTGGCAAGATCGCCAAACGG‑CAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCC‑AATCAGG‑AACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTa                                                                                                                                     >  1:291728/1‑238 (MQ=255)
                                                                                       aTTTCCACGATGCACGGCGACCGCCATCGCCACGAAGAAGCCGATGGTGGCAAGATCGCCAAACGG‑CAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCC‑AATCAGG‑AACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGACgg                                                                                                                              >  2:1141592/1‑239 (MQ=255)
                                                                                          cccACGATGCACGGCGACCGCCATCGCCACGAAGAAGCCGATGGTGGCAAGATCGCCAAACGG‑CAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCC‑AATCAGG‑AACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCg                                                                                                                           <  1:379762/238‑1 (MQ=255)
                                                                                            cACGATGCACGGCGACCGCCATCGCCACGAAGAAGCCGATGGTGGCAAGATCGCCAAACGG‑CAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCC‑AATCAGG‑AACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGtt                                                                                                                         >  1:122542/1‑239 (MQ=255)
                                                                                               gtagCACGGCGACGCCCATCGCCACAAAGAAGCCGATGTTGGCAAGTGCGCCACACGG‑AGGCACCTGATTACCCGGCCCACAAACAGCAATAAAATGGGTGAGGGGGATAAAAACCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCC‑AATCAGG‑AACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTaaa                                                                                                                      <  2:835028/236‑1 (MQ=255)
                                                                                                     cGGCGACCGCCATCGCCACGAAGAAGCCGATGGTGGCAAGATCGCCAAACGG‑CAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCC‑AATCAGG‑AACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAAccatcc                                                                                                                >  2:26347/1‑239 (MQ=255)
                                                                                                         gcccGCCATCGCCACGAAGAAGCCGATGGTGGCAAGATCGCCAAACGG‑CAGCACCTGATTACCCGGCACACAAACAGCAATTAAAAGGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCC‑AATCAGG‑AACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAACCATCCatga                                                                                                            <  1:1066292/237‑1 (MQ=255)
                                                                                                                  cGCCACGAAGAAGCGGATGTTGGCAAGATCGCCAAACGG‑CAGCACCTGATAACCCGGCACACAAACAGCAATTAAAAGGGTGAGTGGGAAAAAAATCAGGCTTGCCGAGACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCC‑AATCAGG‑AACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTaa                                                                                                   <  1:615949/239‑1 (MQ=255)
                                                                                                                          gaagcccGATGGGGGCAAGAGCGCAAACCGGACCGC‑CCGTATTACCCGCCCAAAAAACAGCAATTAAAAGGGGGGGTGGGATAAAAACGAGCGCCGCCCTGACCA‑‑CCGCCGTAGGTCCCAGCAACAACGCCGGATCAAGTCC‑AAGCAGG‑AACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTAATCAccc                                                                                             <  2:10624/232‑1 (MQ=255)
                                                                                                                               ccGATGGTGGCAAGATCGCCAAACGG‑CAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCC‑AATCAGG‑AACTCCTGa                                                                                                                                                                            >  2:942226/1‑153 (MQ=255)
                                                                                                                               ccGATGGTGGCAAGATCGCCAAACGG‑CAGCACCTGATTACCCGGCACACAAACAGAAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCC‑AATCAGG‑AACTCCTGa                                                                                                                                                                            <  1:942226/153‑1 (MQ=255)
                                                                                                                                   tggtggCAAGATCGCCAAACGG‑CAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCC‑AATCAGG‑AACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTAATCACCCGTGGCATTagc                                                                                  >  1:830025/1‑239 (MQ=255)
                                                                                                                                        gCAAGATCGCCAAACGG‑CAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCC‑AATCAGG‑AACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTaaa                                                                                                                      <  1:575805/198‑1 (MQ=255)
                                                                                                                                        gCAAGATCGCCAAACGG‑CAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCC‑AATCAGG‑AACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTaaa                                                                                                                      >  2:575805/1‑198 (MQ=255)
                                                                                                                                         cAAGATCGCCAAACGG‑CAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCC‑AATCAGG‑AACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTAATCACCCGTGGCATTAGCAGCATc                                                                            >  1:270504/1‑239 (MQ=255)
                                                                                                                                         cAAGATCGCCAAACGG‑CAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCC‑AATCAGG‑AACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTAATCACCCGTGGCATTAGCAGCATc                                                                            >  1:599007/1‑239 (MQ=255)
                                                                                                                                         cAAGATCGCCAAACGG‑CAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCC‑AATCAGG‑AACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGATTAATCACCCGTGGCATTAGCAGCATc                                                                            >  1:616595/1‑239 (MQ=255)
                                                                                                                                             aTCGCCAAACGG‑CAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCC‑AATCAGG‑AGCTCCTGACCGCCGAACTTAGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAACCATCCCTGATGGGCTTAATCACCCGTGGCATCAGCAGCATctct                                                                         >  1:438332/1‑238 (MQ=255)
                                                                                                                                                gCCAAACGG‑CAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCC‑AATCAGG‑AACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTaaa                                                                                                                      >  1:533902/1‑190 (MQ=255)
                                                                                                                                                gCCAAACGG‑CAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCC‑AATCAGG‑AACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTAATCACCCGTGGCATTAGCAGCATCACTgcc                                                                      >  2:788109/1‑238 (MQ=255)
                                                                                                                                                               ccTGATCACCCGACCCAAAAACAGAAATAAAATGGGTGAGTGGGAAAAAACCCAGCCTTGCCGAT‑CC‑CCCCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCC‑AATCAGG‑AACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTAATCACCCGTGGCATTAGCAGCATCACTGCCGCCGTTTTTACCGcc                                                       <  2:599007/239‑1 (MQ=255)
                                                                                                                                                                  gtttccccGACACACAAACAGCATTAAAAAGGGGGAGGGGGAAAAAATACAGTCTGGCGG‑‑ACCCCCCCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCC‑AACCAGG‑AACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTAATCACCCGTGGCATTAGCAGCATCACTGCCGCCGTTTTTACCGCCAGc                                                    <  2:1130819/234‑1 (MQ=255)
                                                                                                                                                                  gATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGTTTCCCGATAC‑‑CCCCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCC‑AATCAGG‑AACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTAATCACCCGTGGCATTAGCAGCATCACTGCCGCCGTTTTTACCGCCAGc                                                    <  1:349288/239‑1 (MQ=255)
                                                                                                                                                                  gATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCC‑AATCAGG‑AACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTAATCACCCGTGGCATTAGCAGCATCACTGCCGCCGTTTTTACCGCCa                                                      >  1:1130819/1‑237 (MQ=255)
                                                                                                                                                                     tACCCCGCACCCAAAATGCCATTAAAATTGCGAATTGGAAAAAAAACCAGCCCGGCGAGTCC‑CCCCTCCGTGCGTCCCACCACAAACGCGAGA‑CAACGCCAAAGAAGG‑ACCGCCTGACCGCCGAACTTCGCCTGTAAACGACAACGAGCCTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTAATCCCCCGTGGCATTAGCAGCATCACTGCCGCCGTTTTTACCGCCAGc                                                    <  2:211875/237‑1 (MQ=255)
                                                                                                                                                                        ccGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCC‑AATCAGG‑AACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTAATCACCCGTGGCATTAGCAGCATCACTGCCGCCGTTTTTACCGCCAGCTGCAAt                                              >  2:670566/1‑239 (MQ=255)
                                                                                                                                                                        ccGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCC‑AATCAGG‑AACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTAATCACCCGTGGCATTAGCAGCATCACTGCCGCCGTTTTTACCGCCAGCTGCAAt                                              >  1:634360/1‑239 (MQ=255)
                                                                                                                                                                          ggCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCC‑AATCAGG‑AACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTaaa                                                                                                                      >  1:692067/1‑165 (MQ=255)
                                                                                                                                                                          ggCACAAAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATCCCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCC‑AATCAGG‑AACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTAATCACCCGTGGCATTAGCAGCATCACTGCCGCCGTTTTTACCGCCAGCTGCAATac                                            <  1:57844/239‑1 (MQ=255)
                                                                                                                                                                            cacacaAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCC‑AATCAGG‑AACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTAATCACCCGTGGCATTAGCAGCATCACTGCCGCCGTTTTTACCGCCCGCTGCAATaca                                           >  1:464050/1‑238 (MQ=255)
                                                                                                                                                                             acacaAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCC‑AATCAGG‑AACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGAt                                                                                                                                >  2:263446/1‑152 (MQ=255)
                                                                                                                                                                             acacaAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCC‑AATCAGG‑AACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGAt                                                                                                                                <  1:263446/152‑1 (MQ=255)
                                                                                                                                                                                caAACAGCATTTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATAC‑‑CCCCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCC‑AATCAGG‑AACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTAATCACCCGTGGCATTAGCAGCATCACTGCCGCCGTTTTTACCGCCAGCTGCAATACACCttt                                      <  1:159509/239‑1 (MQ=255)
                                                                                                                                                                                 aaaCAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCC‑AATCAGG‑AACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTAATCACCCGTGGCATTAGCAGCATCACTGCCGCCGTTTTTACCGCCAGCTGCAATACACCTTTg                                     >  2:1019130/1‑239 (MQ=255)
                                                                                                                                                                                      gCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCAAGCAGCAACGCCGGATCAAGGCC‑AATCAGG‑AACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGAGCGTTAAACCATCCATGATGGGTTTAATCACCCGTGGCATTAGCAGCATCAATGACGCCGTTTTTACCGCCAGCTGCAATACACCTTTGACa                                  >  1:229311/1‑237 (MQ=255)
                                                                                                                                                                                                    gAGTGGGATAAAAATCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCC‑AATCAGG‑AACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTAATCACCCGTGGCATTAGCAGCATCACTGCCGCCGTTTTTACCGCCAGCTGCAATACACCTTTGACATCGTAACCCGCGAGGa                  >  1:931629/1‑239 (MQ=255)
                                                                                                                                                                                                     aGTGGGATAAAAATCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCC‑AATCAGG‑AACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTAATCACCCGTGGCATTAGCAGCATCACTGCCGCCGTTTTTACCGCCAGCTGCAATACACCTTTGACATCGTAACCCGCGAGGAt                 >  2:1030583/1‑239 (MQ=255)
                                                                                                                                                                                                     aGTGGGATAAAAATCAGCCTGGCCGATCCCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCC‑AATCAGG‑AACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTAATCACCCGTGGCATTAGCAGCATCACTGCCGCCGTTTTTACCGCCAGCTGCAATACACCTTTGACATCGTAACCCGCGAGGa                  <  1:881192/238‑1 (MQ=255)
                                                                                                                                                                                                     aGTGGAATAAAAACCGGCCTGCCCAATAC‑‑CCCCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCC‑AATCAGG‑AACTCCTGCCCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTAATCACCCGTGGCATTAGCAGCATCACTGCCGCCGTTTTTACCGCCAGCTGCAATACACCTTTGACATCGTAACCCGCGAGGa                  <  2:1178340/238‑1 (MQ=255)
                                                                                                                                                                                                      tggggAAAAAAAACAGGCGTGCCGAATCCACCCC‑CCGTATGGCCCAGCAGCAACGCCGGATCAAGGCC‑AAACAGA‑AACCCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTAATCACCCGTGGCATTAGCAGCATCACTGCCGCCGTTTTTACCGCCAGCTGCAATACACCTTTGACATCGTAACCCGCGAGGa                  <  2:549773/236‑1 (MQ=255)
                                                                                                                                                                                                       tGGGATAAAAATCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCC‑AATCAGG‑AACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTAATCAccc                                                                                             >  2:476801/1‑161 (MQ=255)
                                                                                                                                                                                                       tGGGATAAAAATCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCC‑AATCAGG‑AACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTAATCAccc                                                                                             <  1:476801/161‑1 (MQ=255)
                                                                                                                                                                                                       tGGGATAAAAATCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCC‑AATCAGG‑AACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTAATCACCCGTGGCATTAGCAGCATCACTGCCGCCGTTTTTACCGCCAGCTGCAATACACCTTTGACATCGTAACCCGCGAGGATg                >  1:109127/1‑238 (MQ=255)
                                                                                                                                                                                                         ggAAAAAAATCAGGCTTGCCGAT‑CC‑CCCCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCC‑AATCAGG‑AACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTAATCACCCGTGGCATTAGCAGCATCACTGCCGCCGTTTTTACCGCCAGCTGCAATACACCTTTGACATCGTAACCCGCGAGGATGcc              <  1:125475/238‑1 (MQ=255)
                                                                                                                                                                                                          gATAAAAATCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCC‑AATCAGG‑AACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTAATCACCCGTGGCATTAGCAGCATCACTGCCGCCGTTTTTACCGCCAGCTGCAATACACCTTTGACATCGTAACCAGCGAGGATGCCGa            >  2:1046543/1‑239 (MQ=255)
                                                                                                                                                                                                           aTAAAAATCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCC‑AATCAGG‑AACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTAATCACCCGTGGCATTAGCAGCATCACTGCCGCCGTTTTTACCGCCAGCTGCAATACACCTTTGACATCGTAACCCGCGAGGATGCCGAt           >  1:909739/1‑239 (MQ=255)
                                                                                                                                                                                                             aaaaaCCGGGCTTGCCGAT‑CC‑CCCCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCC‑AATCAGG‑AACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTAATCACCCGTGGCATTAGCAGCATCACTGCCGCCGTTTTTACCGCCAGCTGCAATACACCTTTGACATCGTAACCCGCGAGGATGCCGATaa         <  1:145213/239‑1 (MQ=255)
                                                                                                                                                                                                                    aGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCC‑AATCAGG‑AACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTAATCACCCGTGGCATTAGCAGCATCACTGCCGCCGTTTTTACCGCCAGCTGCAATACACCTTTGACATCGTAACCCGCGAGGATGCCGATAATCAGcc   >  2:754397/1‑238 (MQ=255)
                                                                                                                                                                                                                     ggCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCC‑AATCAGG‑AACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTAATCACCCGTGGCATTAGCAGCATCACTGCCGCCGTTTTTACCGCCAGCTGCAATACACCTTTGACATCGTAACCCGCGAGGATGCCGATAATCAGccc  >  1:209974/1‑238 (MQ=255)
                                                                                                                                                                                                                         tcgcccgtaccccctGCCGTATCGCCCCGCACCAACGCCGG‑TCAAGGGC‑CAACAGGTACCTCCTACCCGCCGAACTTCGCCTATAACCGACTACGAGCCTGCGTAGCGATGGGCGTAAAACCATCCATGATGGGTTTAATCTCCCGTGGCATTAGCAGCATCACTGCCGCCGTTTTTACCGCCAGCTGCAATACACCTTTGACATCGTACCCCGCGAGGATGCCGATAATCAGccc  <  2:830025/225‑1 (MQ=255)
                                                                                                                                                                                                                                    ||                                                                                                                                                                                                                                 
CAGTTGGGTGTGCAAACCAATAGTTTGCGTCGCGATCCACAGGGTGATGCTCATAATGATGACACCCGAGATTAAGGTGCGGAACAGATTTCCACGATGCACGGCGACCGCCATCGCCACGAAGAAGCCGATGGTGGCAAGATCGCCAAACGG‑CAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCACCCCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCC‑AATCAGG‑AACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTAATCACCCGTGGCATTAGCAGCATCACTGCCGCCGTTTTTACCGCCAGCTGCAATACACCTTTGACATCGTAACCCGCGAGGATGCCGATAATCAGCCC  >  NC_000913/2173134‑2173585

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: ‑