Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,285,877 C→G 15.6% P134A (CCG→GCG)  agaD → N‑acetylgalactosamine‑specific enzyme IID component of PTS

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,285,8770CG15.6% 143.5 / 4.8 62P134A (CCG→GCG) agaDN‑acetylgalactosamine‑specific enzyme IID component of PTS
Reads supporting (aligned to +/- strand):  ref base C (14/37);  new base G (8/2);  total (22/40)
Fisher's exact test for biased strand distribution p-value = 2.85e-03
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GCGGCGATGAAAGATAACCTCGAATTTATTAATACCCACCCGAATCTGGTCGGATTCCTGATGGGGTTATTAATTTCGATGGAAGAAAAAGGAGAAAACCGCGACACCATTAAAGGCCTCAAAGTGGCACTGTTTGGCCCAATCGCCGGGATTGGCGATGCGATTTTCTGGTTTACTTTGTTGCCGATTATGGCGGGAATTTGCTCATCATTTGCCAGCCAGGGAAACCTGCTGGGGCCGATTCTATTTTTCGCCGTTTACCTGCTTATCTTTTTCCTGCGCGTCGGCTGGACCCACGTCGGTTATTCAGTCGGCGTGAAGGCGATCGATAAAGTGCGAGAGAACTCGCAGATGATTGCCCGTTCGGCAACCATCCTCGGGATCACGGTAATCGGCGGGCTGATCGCTTCGTATGTGCATATTAACGTGGTGACATCGTTTGCCATCGACAATACCCACAGCGTTGCG  >  NC_000913/3285640‑3286107
                                                                                                                                                                                                                                              |                                                                                                                                                                                                                                      
gcggcgATGAAAGATAACCTCGAATTTATTAATACCCACCCGAATCTGGTCGGATTCCTGATGGGGTTATTAATTTCGATGGAAGAAAAAGGAGAAAACCGCGACACCATTAAAGGCCTCAAAGTGGCACTGTTTGGCCCAATCGCCGGGATTGGCGATGCGATTTTCTGGTTTACTTTGTTGCCGATTATGGCGGGAATTTGCTCATCATTTGCCAGCCAGGGAAACCTGCTGGGGc                                                                                                                                                                                                                                        <  1:317352/238‑1 (MQ=255)
              tAACCTCGAATTTATTAATACCCACCCGAATCTGGTCGGATTCCTGATGGGGTTATTAATTTCGATGGAAGAAAAAGGAGAAAACCGCGACACCATTAAAGGCCTCAAAGTGGCACTGTTTGGCCCAATCGCCGGGATTGGCGATGCGATTTTCTGGTTTACTTTGTTGCCGATTATGGCGGGAATTTGCTCATCATTTGCCAGCCAGGGAAACCTGCTGGGGCCGATTCTAttttt                                                                                                                                                                                                                           <  1:416279/237‑1 (MQ=255)
              tAACCCCGAATTTATTAAAACCCACCCGAATCTGGTCGGATTCCTGATGGGGTTAGTAAGTTCGATGGAAGAAAAAGGAGAAAACCGCGACACCATTAAAGGCCTCAAAGTGGCACTGTTTGGCCCAATCGCCGGGATTGGCGATGCGATTTTCTGGTTTACTTTGTTGCCGATTATGGCGGGAATTTGCTCATCATTTGCCAGCCAGGGAAACCTGCTGGGGCCGATTCTAttttt                                                                                                                                                                                                                           <  1:984156/237‑1 (MQ=255)
                aCCTCGAATTTATTAATACCCACCCGAATCTGGTCGGATTCCTGATGGGGTTATTAATTTCGATGGAAGAAAAAGGAGAAAACCGCGACACCATTAAAGGCCTCAAAGTGGCACTGTTTGGCCCAATCGCCGGGATTGGCGATGCGATTTTCTGGTTTACTTTGTTGCCGATTATGGCGGGGATTTGGTCATTATTTTTCAGCCCGGGGAAACTTCTGGGGGCGGTTTTTTTTTTTGcc                                                                                                                                                                                                                       >  1:172346/1‑239 (MQ=255)
                    cGAATTTATTAATACCCACCCGAATCTGGTCGGATTCCTGATGGGGTTATTAATTTCGATGGAAGAAAAAGGAGAAAACCGCGACACCATTAAAGGCCTCAAAGTGGCACTGTTTGGCCCAATCGCCGGGATTGGCGATGCGATTTTCTGGTTTACTTTGTTGCCGATTATGGCGGGAATTTGCTCATCATTTGCCAGCCAGGGAAACCTGCTGGGGCCGATTCTATTTTTCGCCGttt                                                                                                                                                                                                                   <  1:157058/239‑1 (MQ=255)
                      aaTTTATTAATACCCACCCGAATCTGGTCGGATTCCTGATGGGGTTATTAATTTCGATGGAAGAAAAAGGAGAAAACCGCGACACCATTAAAGGCCTCAAAGTGGCACTGTTTGGCCCAATCGCCGGGATTGGCGATGCGATTTTCTGGTTTACTTTGTTGCCGATTATGGCGGGAAATTTCTCATCGTTTTGCATCCAAGGGAACACGCTGGGGGCGGTTCTtttttttttcccttt                                                                                                                                                                                                                  >  1:199853/1‑229 (MQ=255)
                                     acccGAATCTGGTCGGATTCCTGATGGGGTTATTAATTTCGATGGAAGAAAAAGGAGAAAACCGCGACACCATTAAAGGCCTCAAAGTGGCACTGTTTGGCCCAATCGCCGGGATTGGCGATGCGATTTTCTGGTTTACTTTGTTGCCGATTATGGCGGGGATTTGTTTTTCTTTTTCCCTCCCCGGGAAACTGCTTGGGGCGGGTTTTTTTTTTTCTGTTTTCCTTCTTTTCtttttt                                                                                                                                                                                                  >  2:572863/1‑238 (MQ=255)
                                     acccGAATCTGGTCGGATTCCTGATGGGGTTATTAATTTCGATGGAAGAAAAAGGAGAAAACCGCGACACCATTAAAGGCCTCAAAGTGGCACTGTTTGGCCCAATCGCCGGGATTGGCGATGCGATTTTCTGGTTTACTTTGTTGCCGATTATGGCGGGAGTTTGTTCATCCTTTGCCACCCAGGGAAACATGCTGGGGCCGCTTATCTTTTTTTCCGGTTATCTGCTTATCtttttt                                                                                                                                                                                                  >  2:1104478/1‑238 (MQ=255)
                                     acccGAATCTGGTCGGATTCCTGATGGGGTTATTAATTTCGATGGAAGAAAAAGGAGAAAACCGCGACACCATTAAAGGCCTCAAAGTGGCACTGTTTGGCCCAATCGCCGGGATTGGCGATGCGATTTTCTGGTTTACTTTGTTGCCGATTATGGCGGGAATATGCTCCTCCTTTTCCTACCACCGGAACATGCTGGGGGCGCTTCTATTTTTTCCCGTTTTGCCTCTTCTCtttttt                                                                                                                                                                                                  >  2:335367/1‑238 (MQ=255)
                                      cccGAATCTGGTCGGATTCCTGATGGGGTTATTAATTTCGATGGAAGAAAAAGGAGAAAACCGCGACACCATTAAAGGCCTCAAAGTGGCACTGTTTGGCCCAATCGCCGGGATTGGCGATGCGATTTTCTGGTTTACTTTGTTGCCGATTATGGCGGGGATTTGCTTATCTTTTTTCCCCCAGGGAAACCTTCTGGGGGCGTTTCTATTTTTTGCCGCTTTTCCGCtttttttttt                                                                                                                                                                                                   >  1:453465/1‑237 (MQ=255)
                                      cccGAATCTGGTCGGATTCCTGATGGGGTTATTAATTTCGATGGAAGAAAAAGGAGAAAACCGCGACACCATTAAAGGCCTCAAAGTGGCACTGTTTGGCCCAATCGCCGGGATTGGCGATGCGATTTTCTGGTTTACTTTGTTGCCGATTATGGCGGCGACTTTTTCAACTCTTTCCATTCCGGGGAAACTGCTAGGGGCGGTTCTTTTTTTTGTCCTTTATCCGCCTGTCttttttt                                                                                                                                                                                                 >  1:238489/1‑237 (MQ=255)
                                          aaTCTGGTCGGATTCCTGATGGGGTTATTAATTTCGATGGAAGAAAAAGGAGAAAACCGCGACACCATTAAAGGCCTCAAAGTGGCACTGTTTGGCCCAATCGCCGGGATTGGCGATGCGATTTTCTGGTTTACTTTGTTGCCGATTATGGCGGGGATTTGCTCATCATTTTCTAGCCCTGGGAGCCTGCTTGGGGCCCGTGTATTTTTTTCCCTTTTCCTTCTTATCTTTTTTCtccc                                                                                                                                                                                             >  2:597869/1‑236 (MQ=255)
                                           aTCTGGTCGGATTCCTGATGGGGTTATTAATTTCGATGGAAGAAAAAGGAGAAAACCGCGACACCATTAAAGGCCTCAAAGTGGCACTGTTTGGCCCAATCGCCGGGATTGGCGATGCGATTTTCTGGTTTACTTTGTTGCCGATTATGGCGGGAATTTGCTCATCATTTGCCAGCCAGGGAAACCTGCTGGGGCCGATTCTATTTTTCGCCGTTTACCTGCTTATCTTTTTCCTgcgc                                                                                                                                                                                            <  1:901008/239‑1 (MQ=255)
                                             cTGGGCGGATTCCGGATGGGGTTATTAATTGCGATGGAAGAAAAAGGAGAAATCCGCGACACCATTAAAGGCCGAAAAGTGGCACTGTTTGGCCCAATCGCCGGGATTGGCGATGCGATTTTCTGGTTTACTTTGTTGCCGATTATGGCGGGAATTTGCTCATCATTTGCCAGCCAGGGAAACCTGCTGGGGCCGATTCTATTTTTCGCCGTTTACCTGCTTATCTTTTTCCTGCGCGt                                                                                                                                                                                          <  2:450217/239‑1 (MQ=255)
                                                          tGATGGGGTTATTAATTTCGATGGAAGAAAAAGGAGAAAACCGCGACACCATTAAAGGCCTCAAAGTGGCACTGTTTGGCCCAATCGCCGGGATTGGCGATGCGATTTTCTGGTTTACTTTGTTGCCGATTATGGCGGGAATTTGCTCATCATTTGCCAGCCAGGGAAACCTGCTGGGGCCGATTCTATTTTTCGCCGTTTACCTGCTTATCTTTTTCCTGCGCGTCGGCTGGACCCAc                                                                                                                                                                             <  1:572863/239‑1 (MQ=255)
                                                                    attaATTTCGATGGAAGAAAAAGGAGAAAACCGCGACACCATTAAAGGCCTCAAAGTGGCACTGTTTGGCCCAATCGCCGGGATTGGCGATGCGATTTTCTGGTTTACTTTGTTGCCGATTATGGCGGTGATTTTCTCATCACTTGCCCGCCCGGGGAACCTGCTTGGGGCCGTTCTATTTTTTCCCGTTTAACTCTCTATCTTTTTTCTGCGCGGGGGTGGGAACCACGTCCGTtttt                                                                                                                                                                   >  1:1084748/1‑239 (MQ=255)
                                                                          ttCGATGGAAGAAAAAGGAGAAAACCGCGACACCATTAAAGGCCTCAAAGTGGCACTGTTTGGCCCAATCGCCGGGATTGGCGATGCGATTTTCTGGTTTACTTTGTTGCCGATTATGGCGGGAATTTTCTCATCCTTTTCCAGCCACGGAAACCTGCTGGGGCCGGTTCTTTTTTTTGTCGCTTTCCTGCTTATTTTTTTTTTGCGCGCTGGCTGGACCCACGTCGGGTTTTCAGTCg                                                                                                                                                             >  1:865003/1‑239 (MQ=255)
                                                                            cGATGGAAGAAAAAGGAGAAAACCGCGACACCATTAAAGGCCTCAAAGTGGCACTGTTTGGCCCAATCGCCGGGATTGGCGATGCGATTTTCTGGTTTACTTTGTTGCCGATTATGGCGGGAAATTGCTCATTATTTTCCCCCCCGGGAAACCTGCTTGGGGCGGTTTTATTTTTTGCCGTTTTCCCTCTTATTTTTTTTCTGCGCGTCGGCTGGGCCCCCGCCCGTttttttctcgt                                                                                                                                                            >  1:289789/1‑231 (MQ=255)
                                                                               tGGAAGAAAAAGGAGAAAACCGCGACACCATTAAAGGCCTCAAAGTGGCACTGTTTGGCCCAATCGCCGGGATTGGCGATGCGATTTTCTGGTTTACTTTGTTGCCGATTATGGCGGGGAATTGTTCATCCTTTGCCCGCCATGGGAACCCGCTGGGGCCGAGTCTATTTTTTGCCGCTTTCCTGCCTTTTTTTTTTCTTCTCGTCGGCTGGGACCACGCCGGTTTTTCTCTTGGCGTg                                                                                                                                                        >  1:767665/1‑239 (MQ=255)
                                                                                 gaagaaAAAGGCGAAAACCGCGACACCATTAAAGGCCTCAAAGTGGCACTGTTTGGCCCAATCGCCGGGATTGGCGATGCGATTTTCTGGTTTACTTTGTTGCCGATTATGGCGGGAATTTGCTCATCATTTGCCAGCCAGGGAAACCTGCTGGGGCCGATTCTATTTTTCGCCGTTTACCTGCTTATCTTTTTCCTGCGCGTCGGCTGGACCCACGTCGGTTATTCAGTCGGCGTGa                                                                                                                                                       <  2:126848/238‑1 (MQ=255)
                                                                                      aaaaGGAGAAAACCGCGACACCATTAAAGGCCTCAAAGTGACACTGTTTGGCCCAATCGCCGGGATTGGCGATGCGATTTGCTGGTTTACTTTGTTGCCGATTATGGCGGGAATTTGCTCATCATTTGCCAGCCAGGGAAACCTGCTGGGGCCGATTCTATTTTTCGCCGTTTCCCTGCTTATCTTTTTCCTGCGCGTCGGCTGGACCCACGTCGGTTATTCAGTCGGCGTGAAGGc                                                                                                                                                   <  1:1054598/237‑1 (MQ=255)
                                                                                             gaAAACCGCGACACCATTAAAGGCCTCAAAGTGGCACTGTTTGGCCCAATCGCCGGGATTGGCGATGCGATTTTCTGGTTTACTTTGTTGCCGATTATGGCGGGAATTTGCTCATCATTTGCCAGCCAGGGAAACCTGCTGGGGCCGATTCTATTTTTCGCCGTTTACCTGCTTATCTTTTTCCTGCGCGTCGGCTGGACCCACGTCGGTTATTCAGTCGGCGTGAAGGCGATCGATa                                                                                                                                           <  1:668328/238‑1 (MQ=255)
                                                                                                              tAAAGGCCTCAAAGTGGCTCTGTTTCGCTCAATCGCCGGGATTGGCGATGCGATTTTCTGGTTTACTTTGTTTCCGATTATGGCGGGAATTTGCTCATCATTTGCCAGCCAGGGAAACCTGCTGGTGCCGATTCTATTTTTCGCCGTTTACCTGCTTATCTTTTTCCTGCGCGTCGGCTGGACCCACGTCGGTTATTCAGTCGGCGTGAAGGCGATCGATAAAGTGCGAGAGAACTCGc                                                                                                                         <  2:664892/239‑1 (MQ=255)
                                                                                                               aaaGGCCTCAAAGTGGCGCTGTTTGGCCCAATCGCCTGGATTGGCGATGCGATTTTCTGGTTTGCGTTGGTGCCGATTATGGCGGGAATTTGCTCATCATTTGCCAGCCAGGGAAACCTGCTGGGGCCGATTCTATTTTTCGCCGTTTACCTGCTTATCTTTTTCCTGCGCGTCGGCTGGACCCACGTCGGTTATTCAGTCGGCGTGAAGGCGATCGATAAAGTGCGAGAGAACTCGCa                                                                                                                        <  2:721168/239‑1 (MQ=255)
                                                                                                               aaaGGCCTCAAAGTGGCACTGTTTGGCCCAATCGCCGGGATTGGCGATGCGATTTTCTGGTTTACTTTGTTGCCGATTATGGCGGGAATTTGCTCATCATTTGCCAGCCAGGGAAACCTGCTGGGGCCGATTCTATTTTTCGCCGTTTACCTGCTTATCTTTTTCCTGCGCGTCGGCTGGACCCACGTCGGTTATTCAGTCGGCGTGAAGGCGATCGATAAAGTGCGAGAGAACTc                                                                                                                           <  1:609755/236‑1 (MQ=255)
                                                                                                                     cTCAAAGTGGCACTGTTTGGCCCAATCGCCGGGATTGGCGATGCGATTTTCTGGTTTACTTTGTTGCCGATTATGGCGGGAATTTTCTCATCATTTTTCCGCCCGGGGAACCTGCTTGGGCCGATTCTATTTTTTTCCGTTTTTCTGCTTTTCTTTTTTCTTCGCGTCGGCTGGACCCCCCTCGGTTATTCTGTCGGCGTGAAGGCGATCGATAAAATGCGGGAGAACTCGCAgatgat                                                                                                                  >  1:1023607/1‑239 (MQ=255)
                                                                                                                        aaaGTGGCACTGTTTGGCCCAATCGCCGGGATTGGCGATGCGATTTTCTGGTTTACTTTGTTGCCGATTATGGCGGGAATTTGCTCATCATTTGCCAGCCAGGGAAACCTGCTGGGGCCGATTCTATTTTTCGCCGTTTACCTGCTTATCTTTTTCCTGCGCGTCGGCTGGACCCACGTCGGTTATTCAGTCGGCGTGAAGGCGATCGATAAAGTGCGAGAGAACTCGCAGATGATTGc                                                                                                               <  1:756322/239‑1 (MQ=255)
                                                                                                                          aGTGGCACTGTTTGGCCCAATCGCCGGGATTGGCGATGCGATTTTCTGGTTTACTTTGTTGCCGATTATGGCGGGAATTTGCTCATCATTTGCCAGCCAGGGAAACCTGCTGGGGCCGATTCTATTTTTCGCCGTTTACCTGCTTATCTTTTTCCTGCGCGTCGGCTGGACCCAc                                                                                                                                                                             <  1:119103/175‑1 (MQ=255)
                                                                                                                                   gTTTGGCCCAATCGCCGGGATTGGCGATGCGATTTTCTGGTTTACTTTGTTGCCGATTATGGCGGGAATTTGCTCATCATTTGCCAGCCAGGGAAACCTGCTGGGGCCGATTCTATTTTTCGCCGTTTACCTGCTTATCTTTTTCCTGCGCGTCGGCTGGACCCACGTCGGTTATTCAGTCGGCGTGAAGGCGATCGATAAAGTGCGAGAGAACTCGCAGATGATTGCCCGTTCGGCaa                                                                                                    <  1:259864/239‑1 (MQ=255)
                                                                                                                                     ttGGCCCAATCGCCGGGATTGGCGATGCGATTTTCTGGTTTACTTTGTTGCCGATTATGGCGGGAATTTGCTCATCATTTGCCAGCCAGGGAAACCTGCTGGGGCCGATTCTATTTTTCGCCGTTTACCTGCTTATCTTTTTCCTGCGCGTCGGCTGGACCCACGTCGGTTATTCAGTCGGCGTGAAGGCGATCGATAAAGTGCGAGAGAACTCGCAGATGATTGCCCGTTCGGCAAcc                                                                                                  <  1:1104478/239‑1 (MQ=255)
                                                                                                                                                              tgcgatTTTCTGGTTTACTTTGTTGCCGATTATGGCGGGAATTTGCTCATCATTTGCCAGCCAGGGAAACCTGCTGGGGCCGATTCTATTTTTCGCCGTTTACCTGCTTATCTTTTTCCTGCGCGTCGGCTGGACCCACGGCGGTTATTCAGTCGGCGTGAAGGCGATCGATAAAGTGCGAGAGAACTCGCAGATGATTGCCCGTTCGGCAACCATCCTCGGGATCACGGTAATcggc                                                                          <  1:446358/238‑1 (MQ=255)
                                                                                                                                                                    tttCTGGTTTACTGTGTTGCCGATTATGGCGGGAATTTGCTCACCATTTGCCGGCCAGTGAAACCTGCTGGGGCAGATTCTATTTTTCGCCGTTTACCTGCTTATCTTTTTGCTGCGCGTCGGCTGGACCCACGTCGGTTATTCAGTCGGCGTGAAGGCGATCGATAAAGTGCGAGAGAACTCGCAGATGATTGCCCGTTCGGCAACCATCCTCGGGATCACGGTAATCGGCGGGCTg                                                                    <  2:453465/238‑1 (MQ=255)
                                                                                                                                                                     ttCTGCTTTACTTGGTTGCCGATTATGGCGGGAATGTGCTTATCATTTCCCAGGCAGGGAACCCTGCTGGTGTAGATTCGATTTTTCGCCGTTTACCTGCTTATCTCTTTCCTTCGCGTCGGCTGGACCCACGCCGGTTATTCAGTCGGCGTGAAGGCGATCGATTAAGTGCGAGAGAACTCGCAGATGATTGCGCGTTCGGCAACCATCCGCGAGAACACAGTAATCCGCGGGCTCAt                                                                  <  2:879544/239‑1 (MQ=255)
                                                                                                                                                                       gtgttttACTTTGTTGCCGATTATGGCGAGAATTTGCTCATCATTTGCCAGCCAGTGAATCCTGCTGGGGGCGATTCTATTTTTCGCCGTTTACCTGCGTATCTTTTTCCTGCGCTTCGGCTGGACCCACGTCGGTTATTCAGTCGGCGTGAAGGCGATCGATAAAGTGCGAGAGAACTCGCAGATGATTGCCCGTTCGGCAACCATCCTCGGGATCACGGTAATCGGCGGGCTGATCg                                                                <  2:430205/235‑1 (MQ=255)
                                                                                                                                                                         ggTTTACTTTGTTGCCGATTATGGCGGGAATTTGCTCATCATTTGCCAGCCAGGGAAACCTGCTGGGGCCGATTCTATTTTTCGCCGTTTACCTGCTTATCTTTTTCCTGCGCGTCGGCTGGACCCACGTCGGTTATTCAGTCGGCGTGAAGGCGATCGATAAAGTGCGAGAGAACTCGCAGATGATTGCCCGTTCGGCAACCATCCTCGGGATCACGGTAATCGGCGGGCTGATCGCt                                                              <  1:973120/239‑1 (MQ=255)
                                                                                                                                                                           tttACTTTGTTGCCGATTATGGCGGGAATTTGCTCATCATTTGCCAGCCAGGGAAACCTGCTGGGGCCGATTCTATTTTTCGCCGTTTACCTGCTTATCTTTTTCCTGCGCGTCGGCTGGACCCACGTCGGTTATTCAGTCGGCGTGAAGGCGATCGATAAAGTGCGAGAGAACTCGCAGATGATTGCCCGTTCGGCAACCATCCTCGGGATCACGGTAATCGGCGGGCTGATCGCtt                                                             <  1:141301/238‑1 (MQ=255)
                                                                                                                                                                             tACTTTGTTGCCGATTATGGCGGGAATTTGCTCATCATTTGCCAGCCAGGGAAACCTGCTGGGGCCGATTCTATTTTTCGCCGTTTACCTGCTTATCTTTTTCCTGCGCGTCGGCTGGACCCACGTCGGTTATTCAGTCGGCGTGAAGGCGATCGATAAAGTGCGAGAGAACTCGCAGATGATTGCCCGTTCGGCAACCATCCTCGGGATCACGGTAATCGGCGGGCTGATCGCTTCGt                                                          <  1:597869/239‑1 (MQ=255)
                                                                                                                                                                             gacttgtgtgCCGTTTATGGCGGAAATGTGCTCATCTTTGGCCAGCCAGGGGAACCTGCTTGGGCCGATTCTATTTTTCGCCGTTTTCCTGCTTATCTTTGTCCGGCCCGTCGGCTGGACCCACGTCGGTTATTCAGTCGGCGTGAAGGCGATCGATAAAGTGCGAGAGAACTCGCAGATGATTGCCCGTTCGGCAACCATCCTCGGGATCACGGTAAGCGGCGGGCTGAGCGCTGCGt                                                          <  2:1083752/231‑1 (MQ=255)
                                                                                                                                                                               cTTTGTTGCCGATTATGGCGGGAATTTGCTCATCATTTGCCAGCCAGGGAAACCTGCTGGGGCCGATTCTATTTTTCGCCGTTTACCTGCTTATCTTTTTCCTGCGCGTCGGCTGGACCCACGTCGGTTATTCAGTCGGCGTGAAGGCGATCGATAAAGTGCGAGAGAACTCGCAGATGATTGCCCGTTCGGCAACCATCCTCGGGATCACGGTAATCGGCGGGCTGATCGCTTCGt                                                          <  1:1024574/237‑1 (MQ=255)
                                                                                                                                                                                tttGTTGCCGATTATGGCGGGAATTTGCTCATCATTTGCCAGCCAGGGAAACCTGCTGGGGCCGATTCTATTTTTCGCCGTTTACCTGCTTATCTTTTTCCTGCGCGTCGGCTGGACCCACGTCGGTTATTCAGTCGGCGTGAAGGCGATCGATAAAGTGCGAGAGAACTCGCAGATGATTGCCCGTTCGGCAACCATCCTCGGGATCACGGTAATCGGCGGGCTGATCGCTTCGTAtg                                                       <  1:471160/239‑1 (MQ=255)
                                                                                                                                                                                 ttgttgCCGATTATGGCGGGAATTTGCTCATCATTTGCCAGCCAGGGAAACCTGCTGGGGCCGATTCTATTTTTCGCCGTTTACCTGCTTATCTTTTTCCTGCGCGTCGGCTGGACCCACGTCGGTTATTCAGTCGGCGTGAAGGCGATCGATAAAGTGCGAGAGAACTCGCAGATGATTGCCCGTTCGGCAACCATCCTCGGGATCACGGTAATCGGCGGGCTGATCGCTTCGTAtgt                                                      <  1:1015497/239‑1 (MQ=255)
                                                                                                                                                                                  tgttgCCGATTATGGCGGGAATTTGCTCATCATTTGCCAGCCAGGGAAACCTGCTGGGGCCGATTCTATTTTTCGCCGTTTACCTGCTTATCTTTTTCCTGCGCGTCGGCTGGACCCACGTCGGTTATTCAGTCGGCGTGAAGGCGATCGATAAAGTGCGAGAGAACTCGCAGATGATTGCCCGTTCGGCAACCATCCTCGGGATCACGGTAATCGGCGGGCTGATCGCTTCGTAtgtg                                                     <  1:1146835/239‑1 (MQ=255)
                                                                                                                                                                                     tgCCGATTATGGCGGGAATTTGCTCATCATTTGCCAGCCAGGGAAACCTGCTGGGGCCGATTCTATTTTTCGCCGTTTACCTGCTTATCTTTTTCCTGCGCGTCGGCTGGACCCACGTCGGTTATTCAGTCGGCGTGAAGGCGATCGATAAAGTGCGAGAGAACTCGCAGATGATTGCCCGTTCGGCAACCATCCTCGGGATCACGGTAATCGGCGGGCTGATCGCTTCGTATGTGCat                                                  <  1:1023347/239‑1 (MQ=255)
                                                                                                                                                                                      gCCGATTATGGCGGGAATTTGCTCATCATTTGCCAGCCAGGGAAACCTGCTGGGGCCGATTCTATTTTTCGCCGTTTGCCTGCTTATCTTTTTCCTGCGCGTCGGCTGGACCCACGTCGGTTATTCAGTCGGCGTGAAGGCGATCGATAAAGTGCGAGAGAACTCGCAGATGATTGCCCGTTCGGCAACCATCCTCGGGATCACGGTAATCGGCGGGCTGATCGCTTCGTATGTGCata                                                 <  1:992362/239‑1 (MQ=255)
                                                                                                                                                                                        cGATTATGGCGGGAATTTGCTCATCATTTGCCAGCCAGGGAAACCTGCTGGGGCCGATTCTATTTTTCGCCGTTTACCTGCTTATCTTTTTCCTGCGCGTCGGCTGGACCCACGTCGGTTATTCAGTCGGCGTGAAGGCGATCGATAAAGTGCGAGAGAACTCGCAGATGATTGCCCGTTCGGCAACCATCCTCGGGATCACGGTAATCGGCGGGCTGATCGCTTCGTATGTGCATAtt                                               <  1:819912/239‑1 (MQ=255)
                                                                                                                                                                                           ttATGGCGGGAATTTGGTCATCATTTGCCAGCCAGGGAACCCTGCTGGGGGCGATTCTATTTTTCGCCGTTTACCTGCTTATCTTTTTCCTGCGCGTCGGCTGGACCCACGTCGGTTATTCAGTCGGCGTGAAGGCGATCGATAAAGTGCGAGAGAACTCGCAGATGATTGCCCGTTCGGCAACCATCCTCGGGATCACGGTAATCGGCGGGCTGATCGCTTCGTATGTGCATATTaa                                             <  2:767665/238‑1 (MQ=255)
                                                                                                                                                                                                  gggAATTTGCTCATCATTTTCCAGCCAGGGAATCCTGCTGGTGCCGATTCTATTTTTCGCAGTTTACCTGCTTATCTTTTTCCTGCGCGTCGGCGGGACCCACGTCGGTTATTCAGTCGGCGTGAAGGCGATCGATAAAGTGCGAGAGAACTCGCAGATGATTGCCCGTTCGGCAACCATCCTCGGGATCACGGTAATCGGCGGGCTGATCGCTTCGTATGTGCATATTAACGTGGTGa                                     <  2:880709/239‑1 (MQ=255)
                                                                                                                                                                                                        ttGCTCATCATTTGCCAGCCAGGGAAACCTGCTGGGGCCGATTCTATTTTTCGCCGTTTACCTGCTTATCTTTTTCCTGCGCGTCGGCTGGACCCACGTCGGTTATTCAGTCGGCGTGAAGGCGATCGATAAAGTGCGAGAGAACTCGCAGATGATTGCCCGTTCGGCAACCATCCTCGGGATCACGGTAATCGGCGGGCTGATCGCTTCGTATGTGCATATTAACGTGGTGACATCGt                               <  1:335367/239‑1 (MQ=255)
                                                                                                                                                                                                             catcctTTGCCAGCCAGGGAAACCTGCTGGGGCCGATTCTATTTTTCGCCGTTTACCTGCTTATCTTTTTCCTGCGCGTCGGCTGGACCCACGTCGGTTATTCAGTCGTCGTGAAGGCGATCGATACAGTGCGAGAGAACTCGCAGATGATTGCCCGTTCGGCAACCATCCTCGGGATCACGTTAATCGGCGGGCTGATCGCTTCGTATGTGCATATTTACTGGGTGACATCTTGTGc                           >  2:720967/1‑238 (MQ=255)
                                                                                                                                                                                                               tcatTTGCCAGCCAGGGAAACCTGCTGGGGCCGATTCTATTTTTCGCCGTTTACCTGCTTATCTTTTTCCTGCGCGTCGGCTGGACCCACGTCGGTTATTCAGTCGGCGTGAAGGCGATCGATAAAGTGCGAGAGAACTCGCAGATGATTGCCCGTTCGGCAACCATCCTCGGGATCACGGTAATCGGCGGGCTGATCGCTTCGTATGTGCATATTAACGTGGTGACATCGTTTGCCAt                        <  1:445312/239‑1 (MQ=255)
                                                                                                                                                                                                               tcatTTGCCAGCCAGGGAAACCTGCTGGGGCCGATTCTATTTTTCGCCGTTTACCTGCTTATCTTTTTCCTGCGCGTCGGCTGGACCCACGTCGGTTATTCAGTCGGCGTGAAGGCGATCGATAAAGTGCGAGAGAACTCGCAGATGATTGCCCGTTCGGCAACCATCCTCGGGATCACGGTAATCGGCGGGCTGATCGCTTCGTATGTGCATATTAACGTGGTGACATCGTTTGCCAt                        <  1:949865/239‑1 (MQ=255)
                                                                                                                                                                                                                     gccagccaGGGAAACCTGCTGGGGCCGATTCTATTTTTCGCCGTTTACCTGCTTATCTTTTTCCTGCGCGTCGGCTGGACCCACGTCGGTTATTCAGTCGGCGTGAAGGCGATCGATAAAGTGCGAGAGAACTCGCAGATGATTGCCCGTTCGGCAACCATCCTCGGGATCACGGTAATCGGCGGGCTGAt                                                                  <  1:566656/191‑1 (MQ=255)
                                                                                                                                                                                                                     gccagccaGGGAAACCTGCTGGGGCCGATTCTATTTTTCGCCGTTTACCTGCTTATCTTTTTCCTGCGCGTCGGCTGGACCCACGTCGGTTATTCAGTCGGCGTGAAGGCGATCGATAAAGTGCGAGAGAACTCGCAGATGATTGCCCGTTCGGCAACCATCCTCGGGATCACGGTAATCGGCGGGCTGAt                                                                  >  2:566656/1‑191 (MQ=255)
                                                                                                                                                                                                                     gccagccaGGGAAACCTGCTGGGGCCGATTCTATTTTTCGCCGTTTACCTGCTTATCTTTTTCCTGCGCGTCGGCTGGACCCACGTCGGTTATTCAGTCGGCGTGAAGGCGATCGATAAAGTGCGAGAGAACTCGCAGATGATTGCCCGTTCGGCAACCATCCTCGGGATCACGGTAATCGGCGGGCTGATCGCTTCGTATGTGCCTATTAACGCGGTGACATCGTTTTCCATcaacaa                  >  1:701277/1‑239 (MQ=255)
                                                                                                                                                                                                                       cagccagGGAAACCTGCTGGGGCCGATTCTATTTTTCGCCGTTTACCTGCTTATCTTTTTCCTGCGCGTCGGCTGGACCCACGTCGGTTATTCAGTCGGCGTGAAGGCGATCGATAAAGTGCGAGAGAACTCGCAGATGATTGCCCGTTCGGCAACCATCCTCGGGATCACGGTAATCGGCGGGCTGATCGCTTCGTATGTGCATATTAACGTGGTGACATCGTTTGCCATCGACAAt                 <  1:319092/238‑1 (MQ=255)
                                                                                                                                                                                                                        agccagGGAAACCTGCTGGGGCCGATTCTATTTTTCGCCGTTTACCTGCTTATCTTTTTCCTGCGCGTCGGCTGGACCCACGTCGGTTATTCAGTCGGCGTGAAGGCGATCGATAAAGTGCGAGAGAACTCGCAGATGATTGCCCGTTCGGCAACCATCCTCGGGATCACGGTAATCGGCGGGCTGATCGCTTCGTATGTGCATATTAACGTGGTGACATCGTTTGCCATCGACAATa                >  1:333276/1‑238 (MQ=255)
                                                                                                                                                                                                                        agccagGGAAACCTGCTGGGGCCGATTCTATTTTTCGCCGTTTACCTGCTTATCTTTTTCCTGCGCGTCGGCTGGACCCACGTCGGTTATTCAGTCGGCGTGAAGGCGATCGATAAAGTGCGAGAGAACTCGCAGATGATTGCCCGTTCGGCAACCATCCTCGGGATCACGGTAATCGGCGGGCTGATCGCTTCGTATGTGCATATTAACGTGGTGACATCGTTTGCCATCGACAATa                >  1:374297/1‑238 (MQ=255)
                                                                                                                                                                                                                          ccagGGAAACCTGCTGGGGCCGATTCTATTTTTCGCCGTTTACCTGCTTATCTTTTTCCTGCGCGTCGGCTGGACCCACGTCGGTTATTCAGTCGGCGTGAAGGCGATCGATAAAGTGCGAGAGAACTCGCAGATGATTGCCCGTTCGGCAACCATCCTCGGGATCACGGTAATCGGCGGGCTGATCGCTTCGTATGTGCATATTAACGTGGTGACATCGTTTGCCATCGAc                    <  1:1147726/232‑1 (MQ=255)
                                                                                                                                                                                                                          ccagGGAAACCTGCTGGGGCCGATTCTATTTTTCGCCGTTTACCTGCTTATCTTTTTCCTGCGCGTCGGCTGGACCCACGTCGGTTATTCAGTCGGCGTGAAGGCGATCGATAAAGTGCGAGAGAACTCGCAGATGATTGCCCGTTCGGCAACCATCCTCGGGATCACGGTAATCGGCGGGCTGATCGCTTCGTATGTGCATATTAACGTGGTGACATCGTTTGCCATCGAc                    >  2:1147726/1‑232 (MQ=255)
                                                                                                                                                                                                                          ccagGGAAACCTGCTGGGGCCGATTCTATTTTTCGCCGTTTACCTGCTTATCTTTTTCCTGCGCGTCGGCTGGACCCACGTCGGTTATTCAGTCGGCGTGAAGGCGATCGATAAAGTGCGAGAGAACTCGCAGATGATTGCCCGTTCGGCAACCATCCTCGGGATCACGGTAATCGGCGGGCTGATCGCTTCGTATGTGCATATTAACGTGGTGACATCGTTTGCCATCGACAATAccc             >  1:694073/1‑239 (MQ=255)
                                                                                                                                                                                                                                   ccTGCTGGGGCCGCTTCTATTTTTCGCCGTTTACCTGCTTATCTTTTTCCTGCGCGTCGGCTGGACCCACGTCGGTTATTCAGTCGGCGTGAAGGCGATCGATAAAGTGCGAGAGAACTTGCAGATGATTGCCCGTTCGGCAACCATCCTCGGGATCAGGGTAATCGGCGGGATGATCGCTTCGCTTGTGCATATTATCGTGTGGACATCGGTTGTCCTCGACAATAccccacgccttg    >  2:902434/1‑230 (MQ=255)
                                                                                                                                                                                                                                     tgctgGGGCCGATTCTATTTTTCGCCGTTTACCTGCTTATCTTTTTCCTGCGCGTCGGCTGGACCCACGTCGGTTATTCAGTCGGCGTGAAGGCGATCGATAAAGTGCGAGAGAACTCGCAGATGATTGCCCGTTCGGCAACCATCCTCGGGATCACGGTAATCGGCGGGCTGATCGCTTCGTATGTGCATATTAACGTGGTGACATCGTTTGCCATCGACAATACCCACAGCGTTgcg  <  1:613108/239‑1 (MQ=255)
                                                                                                                                                                                                                                      gctgGGGCCGATTGTATTTATCGCCGTTTGCCTGCTTATCTTTTTCCTGCGCGTCCGCTGGACCCACGTCGGTTATTCAGTCGGCGTCAAGGCGATCGATAAAGTGCGAGCGAACTCGCAGATGATTGCCCGTTCGGCAACCATCCTAGGGATCTCGGCAATCGGCGGGCTGACTCCTTCGTATGAGAATCTTATCGAGCTGACATCGTGTGCCATCGACAATACCCACAGCGTTgc   >  2:366787/1‑237 (MQ=255)
                                                                                                                                                                                                                                              |                                                                                                                                                                                                                                      
GCGGCGATGAAAGATAACCTCGAATTTATTAATACCCACCCGAATCTGGTCGGATTCCTGATGGGGTTATTAATTTCGATGGAAGAAAAAGGAGAAAACCGCGACACCATTAAAGGCCTCAAAGTGGCACTGTTTGGCCCAATCGCCGGGATTGGCGATGCGATTTTCTGGTTTACTTTGTTGCCGATTATGGCGGGAATTTGCTCATCATTTGCCAGCCAGGGAAACCTGCTGGGGCCGATTCTATTTTTCGCCGTTTACCTGCTTATCTTTTTCCTGCGCGTCGGCTGGACCCACGTCGGTTATTCAGTCGGCGTGAAGGCGATCGATAAAGTGCGAGAGAACTCGCAGATGATTGCCCGTTCGGCAACCATCCTCGGGATCACGGTAATCGGCGGGCTGATCGCTTCGTATGTGCATATTAACGTGGTGACATCGTTTGCCATCGACAATACCCACAGCGTTGCG  >  NC_000913/3285640‑3286107

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: