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

GCGGCGATGAAAGATAACCTCGAATTTATTAATACCCACCCGAATCTGGTCGGATTCCTGATGGGGTTATTAATTTCGATGGAAGAAAAAGGAGAAAACCGCGACACCATTAAAGGCCTCAAAGTGGCACTGTTTGGCCCAATCGCCGGGATTGGCGATGCGATTTTCTGGTTTACTTTGTTGCCGATTATGGCGGGAATTTGCTCATCATTTGCCAGCCAGGGAAACCTGC‑TGGGGCCGATTCTATTTTTCGCCGTTTACCTGCTTATCTTTTTCCTGCGCGTCGGCTGGACCCACGTCGGTTATTCAGTCGGCGTGAAGGCGATCGATAAAGTGCGAGAGAACTCGCAGATGATTGCCCGTTCGGCAACCATCCTCGGGATCACGGTAATCGGCGGGCTGATCGCTTCGTATGTGCATATTAACGTGGTGACATCGTTTGCCATCGACAATACCCACAGCGTTGCG  >  NC_000913/3285640‑3286107
                                                                                                                                                                                                                                              |                                                                                                                                                                                                                                      
gcggcgATGAAAGATAACCTCGAATTTATTAATACCCACCCGAATCTGGTCGGATTCCTGATGGGGTTATTAATTTCGATGGAAGAAAAAGGAGAAAACCGCGACACCATTAAAGGCCTCAAAGTGGCACTGTTTGGCCCAATCGCCGGGATTGGCGATGCGATTTTCTGGTTTACTTTGTTGCCGATTATGGCGGGAATTTGCTCATCATTTGCCAGCCAGGGAAACCTGC‑TGGGGc                                                                                                                                                                                                                                        <  1:317352/238‑1 (MQ=255)
              tAACCTCGAATTTATTAATACCCACCCGAATCTGGTCGGATTCCTGATGGGGTTATTAATTTCGATGGAAGAAAAAGGAGAAAACCGCGACACCATTAAAGGCCTCAAAGTGGCACTGTTTGGCCCAATCGCCGGGATTGGCGATGCGATTTTCTGGTTTACTTTGTTGCCGATTATGGCGGGAATTTGCTCATCATTTGCCAGCCAGGGAAACCTGC‑TGGGGCCGATTCTAttttt                                                                                                                                                                                                                           <  1:416279/237‑1 (MQ=255)
              tAACCCCGAATTTATTAAAACCCACCCGAATCTGGTCGGATTCCTGATGGGGTTAGTAAGTTCGATGGAAGAAAAAGGAGAAAACCGCGACACCATTAAAGGCCTCAAAGTGGCACTGTTTGGCCCAATCGCCGGGATTGGCGATGCGATTTTCTGGTTTACTTTGTTGCCGATTATGGCGGGAATTTGCTCATCATTTGCCAGCCAGGGAAACCTGC‑TGGGGCCGATTCTAttttt                                                                                                                                                                                                                           <  1:984156/237‑1 (MQ=255)
                aCCTCGAATTTATTAATACCCACCCGAATCTGGTCGGATTCCTGATGGGGTTATTAATTTCGATGGAAGAAAAAGGAGAAAACCGCGACACCATTAAAGGCCTCAAAGTGGCACTGTTTGGCCCAATCGCCGGGATTGGCGATGCGATTTTCTGGTTTACTTTGTTGCCGATTATGGCGGGGATTTGGTCATTATTTTTCAGCCCGGGGAAACTTC‑TGGGGGCGGTTTTTTTTTTTGcc                                                                                                                                                                                                                       >  1:172346/1‑239 (MQ=255)
                    cGAATTTATTAATACCCACCCGAATCTGGTCGGATTCCTGATGGGGTTATTAATTTCGATGGAAGAAAAAGGAGAAAACCGCGACACCATTAAAGGCCTCAAAGTGGCACTGTTTGGCCCAATCGCCGGGATTGGCGATGCGATTTTCTGGTTTACTTTGTTGCCGATTATGGCGGGAATTTGCTCATCATTTGCCAGCCAGGGAAACCTGC‑TGGGGCCGATTCTATTTTTCGCCGttt                                                                                                                                                                                                                   <  1:157058/239‑1 (MQ=255)
                      aaTTTATTAATACCCACCCGAATCTGGTCGGATTCCTGATGGGGTTATTAATTTCGATGGAAGAAAAAGGAGAAAACCGCGACACCATTAAAGGCCTCAAAGTGGCACTGTTTGGCCCAATCGCCGGGATTGGCGATGCGATTTTCTGGTTTACTTTGTTGCCGATTATGGCGGGAAATTTCTCATCGTTTTGCATCCAAGGGAACACGC‑TGGGGGCGGTTCTtttttttttcccttt                                                                                                                                                                                                                  >  1:199853/1‑229 (MQ=255)
                                     acccGAATCTGGTCGGATTCCTGATGGGGTTATTAATTTCGATGGAAGAAAAAGGAGAAAACCGCGACACCATTAAAGGCCTCAAAGTGGCACTGTTTGGCCCAATCGCCGGGATTGGCGATGCGATTTTCTGGTTTACTTTGTTGCCGATTATGGCGGGGATTTGTTTTTCTTTTTCCCTCCCCGGGAAACTGC‑TTGGGGCGGGTTTTTTTTTTTCTGTTTTCCTTCTTTTCtttttt                                                                                                                                                                                                  >  2:572863/1‑238 (MQ=255)
                                     acccGAATCTGGTCGGATTCCTGATGGGGTTATTAATTTCGATGGAAGAAAAAGGAGAAAACCGCGACACCATTAAAGGCCTCAAAGTGGCACTGTTTGGCCCAATCGCCGGGATTGGCGATGCGATTTTCTGGTTTACTTTGTTGCCGATTATGGCGGGAGTTTGTTCATCCTTTGCCACCCAGGGAAACATGC‑TGGGGCCGCTTATCTTTTTTTCCGGTTATCTGCTTATCtttttt                                                                                                                                                                                                  >  2:1104478/1‑238 (MQ=255)
                                     acccGAATCTGGTCGGATTCCTGATGGGGTTATTAATTTCGATGGAAGAAAAAGGAGAAAACCGCGACACCATTAAAGGCCTCAAAGTGGCACTGTTTGGCCCAATCGCCGGGATTGGCGATGCGATTTTCTGGTTTACTTTGTTGCCGATTATGGCGGGAATATGCTCCTCCTTTTCCTACCACCGGAACATGC‑TGGGGGCGCTTCTATTTTTTCCCGTTTTGCCTCTTCTCtttttt                                                                                                                                                                                                  >  2:335367/1‑238 (MQ=255)
                                      cccGAATCTGGTCGGATTCCTGATGGGGTTATTAATTTCGATGGAAGAAAAAGGAGAAAACCGCGACACCATTAAAGGCCTCAAAGTGGCACTGTTTGGCCCAATCGCCGGGATTGGCGATGCGATTTTCTGGTTTACTTTGTTGCCGATTATGGCGGGGATTTGCTTATCTTTTTTCCCCCAGGGAAACCTTC‑TGGGGGCGTTTCTATTTTTTGCCGCTTTTCCGCtttttttttt                                                                                                                                                                                                   >  1:453465/1‑237 (MQ=255)
                                      cccGAATCTGGTCGGATTCCTGATGGGGTTATTAATTTCGATGGAAGAAAAAGGAGAAAACCGCGACACCATTAAAGGCCTCAAAGTGGCACTGTTTGGCCCAATCGCCGGGATTGGCGATGCGATTTTCTGGTTTACTTTGTTGCCGATTATGGCGGCGACTTTTTCAACTCTTTCCATTCCGGGGAAACTGC‑TAGGGGCGGTTCTTTTTTTTGTCCTTTATCCGCCTGTCttttttt                                                                                                                                                                                                 >  1:238489/1‑237 (MQ=255)
                                          aaTCTGGTCGGATTCCTGATGGGGTTATTAATTTCGATGGAAGAAAAAGGAGAAAACCGCGACACCATTAAAGGCCTCAAAGTGGCACTGTTTGGCCCAATCGCCGGGATTGGCGATGCGATTTTCTGGTTTACTTTGTTGCCGATTATGGCGGGGATTTGCTCATCATTTTCTAGCCCTGGGAGCCTGC‑TTGGGGCCCGTGTATTTTTTTCCCTTTTCCTTCTTATCTTTTTTCtccc                                                                                                                                                                                             >  2:597869/1‑236 (MQ=255)
                                           aTCTGGTCGGATTCCTGATGGGGTTATTAATTTCGATGGAAGAAAAAGGAGAAAACCGCGACACCATTAAAGGCCTCAAAGTGGCACTGTTTGGCCCAATCGCCGGGATTGGCGATGCGATTTTCTGGTTTACTTTGTTGCCGATTATGGCGGGAATTTGCTCATCATTTGCCAGCCAGGGAAACCTGC‑TGGGGCCGATTCTATTTTTCGCCGTTTACCTGCTTATCTTTTTCCTgcgc                                                                                                                                                                                            <  1:901008/239‑1 (MQ=255)
                                             cTGGGCGGATTCCGGATGGGGTTATTAATTGCGATGGAAGAAAAAGGAGAAATCCGCGACACCATTAAAGGCCGAAAAGTGGCACTGTTTGGCCCAATCGCCGGGATTGGCGATGCGATTTTCTGGTTTACTTTGTTGCCGATTATGGCGGGAATTTGCTCATCATTTGCCAGCCAGGGAAACCTGC‑TGGGGCCGATTCTATTTTTCGCCGTTTACCTGCTTATCTTTTTCCTGCGCGt                                                                                                                                                                                          <  2:450217/239‑1 (MQ=255)
                                                          tGATGGGGTTATTAATTTCGATGGAAGAAAAAGGAGAAAACCGCGACACCATTAAAGGCCTCAAAGTGGCACTGTTTGGCCCAATCGCCGGGATTGGCGATGCGATTTTCTGGTTTACTTTGTTGCCGATTATGGCGGGAATTTGCTCATCATTTGCCAGCCAGGGAAACCTGC‑TGGGGCCGATTCTATTTTTCGCCGTTTACCTGCTTATCTTTTTCCTGCGCGTCGGCTGGACCCAc                                                                                                                                                                             <  1:572863/239‑1 (MQ=255)
                                                                    attaATTTCGATGGAAGAAAAAGGAGAAAACCGCGACACCATTAAAGGCCTCAAAGTGGCACTGTTTGGCCCAATCGCCGGGATTGGCGATGCGATTTTCTGGTTTACTTTGTTGCCGATTATGGCGGTGATTTTCTCATCACTTGCCCGCCCGGGGAACCTGCTTGGGGCCG‑TTCTATTTTTTCCCGTTTAACTCTCTATCTTTTTTCTGCGCGGGGGTGGGAACCACGTCCGTtttt                                                                                                                                                                   >  1:1084748/1‑239 (MQ=255)
                                                                          ttCGATGGAAGAAAAAGGAGAAAACCGCGACACCATTAAAGGCCTCAAAGTGGCACTGTTTGGCCCAATCGCCGGGATTGGCGATGCGATTTTCTGGTTTACTTTGTTGCCGATTATGGCGGGAATTTTCTCATCCTTTTCCAGCCACGGAAACCTGC‑TGGGGCCGGTTCTTTTTTTTGTCGCTTTCCTGCTTATTTTTTTTTTGCGCGCTGGCTGGACCCACGTCGGGTTTTCAGTCg                                                                                                                                                             >  1:865003/1‑239 (MQ=255)
                                                                            cGATGGAAGAAAAAGGAGAAAACCGCGACACCATTAAAGGCCTCAAAGTGGCACTGTTTGGCCCAATCGCCGGGATTGGCGATGCGATTTTCTGGTTTACTTTGTTGCCGATTATGGCGGGAAATTGCTCATTATTTTCCCCCCCGGGAAACCTGC‑TTGGGGCGGTTTTATTTTTTGCCGTTTTCCCTCTTATTTTTTTTCTGCGCGTCGGCTGGGCCCCCGCCCGTttttttctcgt                                                                                                                                                            >  1:289789/1‑231 (MQ=255)
                                                                               tGGAAGAAAAAGGAGAAAACCGCGACACCATTAAAGGCCTCAAAGTGGCACTGTTTGGCCCAATCGCCGGGATTGGCGATGCGATTTTCTGGTTTACTTTGTTGCCGATTATGGCGGGGAATTGTTCATCCTTTGCCCGCCATGGGAACCCGC‑TGGGGCCGAGTCTATTTTTTGCCGCTTTCCTGCCTTTTTTTTTTCTTCTCGTCGGCTGGGACCACGCCGGTTTTTCTCTTGGCGTg                                                                                                                                                        >  1:767665/1‑239 (MQ=255)
                                                                                 gaagaaAAAGGCGAAAACCGCGACACCATTAAAGGCCTCAAAGTGGCACTGTTTGGCCCAATCGCCGGGATTGGCGATGCGATTTTCTGGTTTACTTTGTTGCCGATTATGGCGGGAATTTGCTCATCATTTGCCAGCCAGGGAAACCTGC‑TGGGGCCGATTCTATTTTTCGCCGTTTACCTGCTTATCTTTTTCCTGCGCGTCGGCTGGACCCACGTCGGTTATTCAGTCGGCGTGa                                                                                                                                                       <  2:126848/238‑1 (MQ=255)
                                                                                      aaaaGGAGAAAACCGCGACACCATTAAAGGCCTCAAAGTGACACTGTTTGGCCCAATCGCCGGGATTGGCGATGCGATTTGCTGGTTTACTTTGTTGCCGATTATGGCGGGAATTTGCTCATCATTTGCCAGCCAGGGAAACCTGC‑TGGGGCCGATTCTATTTTTCGCCGTTTCCCTGCTTATCTTTTTCCTGCGCGTCGGCTGGACCCACGTCGGTTATTCAGTCGGCGTGAAGGc                                                                                                                                                   <  1:1054598/237‑1 (MQ=255)
                                                                                             gaAAACCGCGACACCATTAAAGGCCTCAAAGTGGCACTGTTTGGCCCAATCGCCGGGATTGGCGATGCGATTTTCTGGTTTACTTTGTTGCCGATTATGGCGGGAATTTGCTCATCATTTGCCAGCCAGGGAAACCTGC‑TGGGGCCGATTCTATTTTTCGCCGTTTACCTGCTTATCTTTTTCCTGCGCGTCGGCTGGACCCACGTCGGTTATTCAGTCGGCGTGAAGGCGATCGATa                                                                                                                                           <  1:668328/238‑1 (MQ=255)
                                                                                                              tAAAGGCCTCAAAGTGGCTCTGTTTCGCTCAATCGCCGGGATTGGCGATGCGATTTTCTGGTTTACTTTGTTTCCGATTATGGCGGGAATTTGCTCATCATTTGCCAGCCAGGGAAACCTGC‑TGGTGCCGATTCTATTTTTCGCCGTTTACCTGCTTATCTTTTTCCTGCGCGTCGGCTGGACCCACGTCGGTTATTCAGTCGGCGTGAAGGCGATCGATAAAGTGCGAGAGAACTCGc                                                                                                                         <  2:664892/239‑1 (MQ=255)
                                                                                                               aaaGGCCTCAAAGTGGCGCTGTTTGGCCCAATCGCCTGGATTGGCGATGCGATTTTCTGGTTTGCGTTGGTGCCGATTATGGCGGGAATTTGCTCATCATTTGCCAGCCAGGGAAACCTGC‑TGGGGCCGATTCTATTTTTCGCCGTTTACCTGCTTATCTTTTTCCTGCGCGTCGGCTGGACCCACGTCGGTTATTCAGTCGGCGTGAAGGCGATCGATAAAGTGCGAGAGAACTCGCa                                                                                                                        <  2:721168/239‑1 (MQ=255)
                                                                                                               aaaGGCCTCAAAGTGGCACTGTTTGGCCCAATCGCCGGGATTGGCGATGCGATTTTCTGGTTTACTTTGTTGCCGATTATGGCGGGAATTTGCTCATCATTTGCCAGCCAGGGAAACCTGC‑TGGGGCCGATTCTATTTTTCGCCGTTTACCTGCTTATCTTTTTCCTGCGCGTCGGCTGGACCCACGTCGGTTATTCAGTCGGCGTGAAGGCGATCGATAAAGTGCGAGAGAACTc                                                                                                                           <  1:609755/236‑1 (MQ=255)
                                                                                                                     cTCAAAGTGGCACTGTTTGGCCCAATCGCCGGGATTGGCGATGCGATTTTCTGGTTTACTTTGTTGCCGATTATGGCGGGAATTTTCTCATCATTTTTCCGCCCGGGGAACCTGC‑TTGGGCCGATTCTATTTTTTTCCGTTTTTCTGCTTTTCTTTTTTCTTCGCGTCGGCTGGACCCCCCTCGGTTATTCTGTCGGCGTGAAGGCGATCGATAAAATGCGGGAGAACTCGCAgatgat                                                                                                                  >  1:1023607/1‑239 (MQ=255)
                                                                                                                        aaaGTGGCACTGTTTGGCCCAATCGCCGGGATTGGCGATGCGATTTTCTGGTTTACTTTGTTGCCGATTATGGCGGGAATTTGCTCATCATTTGCCAGCCAGGGAAACCTGC‑TGGGGCCGATTCTATTTTTCGCCGTTTACCTGCTTATCTTTTTCCTGCGCGTCGGCTGGACCCACGTCGGTTATTCAGTCGGCGTGAAGGCGATCGATAAAGTGCGAGAGAACTCGCAGATGATTGc                                                                                                               <  1:756322/239‑1 (MQ=255)
                                                                                                                          aGTGGCACTGTTTGGCCCAATCGCCGGGATTGGCGATGCGATTTTCTGGTTTACTTTGTTGCCGATTATGGCGGGAATTTGCTCATCATTTGCCAGCCAGGGAAACCTGC‑TGGGGCCGATTCTATTTTTCGCCGTTTACCTGCTTATCTTTTTCCTGCGCGTCGGCTGGACCCAc                                                                                                                                                                             <  1:119103/175‑1 (MQ=255)
                                                                                                                                   gTTTGGCCCAATCGCCGGGATTGGCGATGCGATTTTCTGGTTTACTTTGTTGCCGATTATGGCGGGAATTTGCTCATCATTTGCCAGCCAGGGAAACCTGC‑TGGGGCCGATTCTATTTTTCGCCGTTTACCTGCTTATCTTTTTCCTGCGCGTCGGCTGGACCCACGTCGGTTATTCAGTCGGCGTGAAGGCGATCGATAAAGTGCGAGAGAACTCGCAGATGATTGCCCGTTCGGCaa                                                                                                    <  1:259864/239‑1 (MQ=255)
                                                                                                                                     ttGGCCCAATCGCCGGGATTGGCGATGCGATTTTCTGGTTTACTTTGTTGCCGATTATGGCGGGAATTTGCTCATCATTTGCCAGCCAGGGAAACCTGC‑TGGGGCCGATTCTATTTTTCGCCGTTTACCTGCTTATCTTTTTCCTGCGCGTCGGCTGGACCCACGTCGGTTATTCAGTCGGCGTGAAGGCGATCGATAAAGTGCGAGAGAACTCGCAGATGATTGCCCGTTCGGCAAcc                                                                                                  <  1:1104478/239‑1 (MQ=255)
                                                                                                                                                              tgcgatTTTCTGGTTTACTTTGTTGCCGATTATGGCGGGAATTTGCTCATCATTTGCCAGCCAGGGAAACCTGC‑TGGGGCCGATTCTATTTTTCGCCGTTTACCTGCTTATCTTTTTCCTGCGCGTCGGCTGGACCCACGGCGGTTATTCAGTCGGCGTGAAGGCGATCGATAAAGTGCGAGAGAACTCGCAGATGATTGCCCGTTCGGCAACCATCCTCGGGATCACGGTAATcggc                                                                          <  1:446358/238‑1 (MQ=255)
                                                                                                                                                                    tttCTGGTTTACTGTGTTGCCGATTATGGCGGGAATTTGCTCACCATTTGCCGGCCAGTGAAACCTGC‑TGGGGCAGATTCTATTTTTCGCCGTTTACCTGCTTATCTTTTTGCTGCGCGTCGGCTGGACCCACGTCGGTTATTCAGTCGGCGTGAAGGCGATCGATAAAGTGCGAGAGAACTCGCAGATGATTGCCCGTTCGGCAACCATCCTCGGGATCACGGTAATCGGCGGGCTg                                                                    <  2:453465/238‑1 (MQ=255)
                                                                                                                                                                     ttCTGCTTTACTTGGTTGCCGATTATGGCGGGAATGTGCTTATCATTTCCCAGGCAGGGAACCCTGC‑TGGTGTAGATTCGATTTTTCGCCGTTTACCTGCTTATCTCTTTCCTTCGCGTCGGCTGGACCCACGCCGGTTATTCAGTCGGCGTGAAGGCGATCGATTAAGTGCGAGAGAACTCGCAGATGATTGCGCGTTCGGCAACCATCCGCGAGAACACAGTAATCCGCGGGCTCAt                                                                  <  2:879544/239‑1 (MQ=255)
                                                                                                                                                                       gtgttttACTTTGTTGCCGATTATGGCGAGAATTTGCTCATCATTTGCCAGCCAGTGAATCCTGC‑TGGGGGCGATTCTATTTTTCGCCGTTTACCTGCGTATCTTTTTCCTGCGCTTCGGCTGGACCCACGTCGGTTATTCAGTCGGCGTGAAGGCGATCGATAAAGTGCGAGAGAACTCGCAGATGATTGCCCGTTCGGCAACCATCCTCGGGATCACGGTAATCGGCGGGCTGATCg                                                                <  2:430205/235‑1 (MQ=255)
                                                                                                                                                                         ggTTTACTTTGTTGCCGATTATGGCGGGAATTTGCTCATCATTTGCCAGCCAGGGAAACCTGC‑TGGGGCCGATTCTATTTTTCGCCGTTTACCTGCTTATCTTTTTCCTGCGCGTCGGCTGGACCCACGTCGGTTATTCAGTCGGCGTGAAGGCGATCGATAAAGTGCGAGAGAACTCGCAGATGATTGCCCGTTCGGCAACCATCCTCGGGATCACGGTAATCGGCGGGCTGATCGCt                                                              <  1:973120/239‑1 (MQ=255)
                                                                                                                                                                           tttACTTTGTTGCCGATTATGGCGGGAATTTGCTCATCATTTGCCAGCCAGGGAAACCTGC‑TGGGGCCGATTCTATTTTTCGCCGTTTACCTGCTTATCTTTTTCCTGCGCGTCGGCTGGACCCACGTCGGTTATTCAGTCGGCGTGAAGGCGATCGATAAAGTGCGAGAGAACTCGCAGATGATTGCCCGTTCGGCAACCATCCTCGGGATCACGGTAATCGGCGGGCTGATCGCtt                                                             <  1:141301/238‑1 (MQ=255)
                                                                                                                                                                             tACTTTGTTGCCGATTATGGCGGGAATTTGCTCATCATTTGCCAGCCAGGGAAACCTGC‑TGGGGCCGATTCTATTTTTCGCCGTTTACCTGCTTATCTTTTTCCTGCGCGTCGGCTGGACCCACGTCGGTTATTCAGTCGGCGTGAAGGCGATCGATAAAGTGCGAGAGAACTCGCAGATGATTGCCCGTTCGGCAACCATCCTCGGGATCACGGTAATCGGCGGGCTGATCGCTTCGt                                                          <  1:597869/239‑1 (MQ=255)
                                                                                                                                                                             gacttgtgtgCCGTTTATGGCGGAAATGTGCTCATCTTTGGCCAGCCAGGGGAACCTGC‑TTGGGCCGATTCTATTTTTCGCCGTTTTCCTGCTTATCTTTGTCCGGCCCGTCGGCTGGACCCACGTCGGTTATTCAGTCGGCGTGAAGGCGATCGATAAAGTGCGAGAGAACTCGCAGATGATTGCCCGTTCGGCAACCATCCTCGGGATCACGGTAAGCGGCGGGCTGAGCGCTGCGt                                                          <  2:1083752/231‑1 (MQ=255)
                                                                                                                                                                               cTTTGTTGCCGATTATGGCGGGAATTTGCTCATCATTTGCCAGCCAGGGAAACCTGC‑TGGGGCCGATTCTATTTTTCGCCGTTTACCTGCTTATCTTTTTCCTGCGCGTCGGCTGGACCCACGTCGGTTATTCAGTCGGCGTGAAGGCGATCGATAAAGTGCGAGAGAACTCGCAGATGATTGCCCGTTCGGCAACCATCCTCGGGATCACGGTAATCGGCGGGCTGATCGCTTCGt                                                          <  1:1024574/237‑1 (MQ=255)
                                                                                                                                                                                tttGTTGCCGATTATGGCGGGAATTTGCTCATCATTTGCCAGCCAGGGAAACCTGC‑TGGGGCCGATTCTATTTTTCGCCGTTTACCTGCTTATCTTTTTCCTGCGCGTCGGCTGGACCCACGTCGGTTATTCAGTCGGCGTGAAGGCGATCGATAAAGTGCGAGAGAACTCGCAGATGATTGCCCGTTCGGCAACCATCCTCGGGATCACGGTAATCGGCGGGCTGATCGCTTCGTAtg                                                       <  1:471160/239‑1 (MQ=255)
                                                                                                                                                                                 ttgttgCCGATTATGGCGGGAATTTGCTCATCATTTGCCAGCCAGGGAAACCTGC‑TGGGGCCGATTCTATTTTTCGCCGTTTACCTGCTTATCTTTTTCCTGCGCGTCGGCTGGACCCACGTCGGTTATTCAGTCGGCGTGAAGGCGATCGATAAAGTGCGAGAGAACTCGCAGATGATTGCCCGTTCGGCAACCATCCTCGGGATCACGGTAATCGGCGGGCTGATCGCTTCGTAtgt                                                      <  1:1015497/239‑1 (MQ=255)
                                                                                                                                                                                  tgttgCCGATTATGGCGGGAATTTGCTCATCATTTGCCAGCCAGGGAAACCTGC‑TGGGGCCGATTCTATTTTTCGCCGTTTACCTGCTTATCTTTTTCCTGCGCGTCGGCTGGACCCACGTCGGTTATTCAGTCGGCGTGAAGGCGATCGATAAAGTGCGAGAGAACTCGCAGATGATTGCCCGTTCGGCAACCATCCTCGGGATCACGGTAATCGGCGGGCTGATCGCTTCGTAtgtg                                                     <  1:1146835/239‑1 (MQ=255)
                                                                                                                                                                                     tgCCGATTATGGCGGGAATTTGCTCATCATTTGCCAGCCAGGGAAACCTGC‑TGGGGCCGATTCTATTTTTCGCCGTTTACCTGCTTATCTTTTTCCTGCGCGTCGGCTGGACCCACGTCGGTTATTCAGTCGGCGTGAAGGCGATCGATAAAGTGCGAGAGAACTCGCAGATGATTGCCCGTTCGGCAACCATCCTCGGGATCACGGTAATCGGCGGGCTGATCGCTTCGTATGTGCat                                                  <  1:1023347/239‑1 (MQ=255)
                                                                                                                                                                                      gCCGATTATGGCGGGAATTTGCTCATCATTTGCCAGCCAGGGAAACCTGC‑TGGGGCCGATTCTATTTTTCGCCGTTTGCCTGCTTATCTTTTTCCTGCGCGTCGGCTGGACCCACGTCGGTTATTCAGTCGGCGTGAAGGCGATCGATAAAGTGCGAGAGAACTCGCAGATGATTGCCCGTTCGGCAACCATCCTCGGGATCACGGTAATCGGCGGGCTGATCGCTTCGTATGTGCata                                                 <  1:992362/239‑1 (MQ=255)
                                                                                                                                                                                        cGATTATGGCGGGAATTTGCTCATCATTTGCCAGCCAGGGAAACCTGC‑TGGGGCCGATTCTATTTTTCGCCGTTTACCTGCTTATCTTTTTCCTGCGCGTCGGCTGGACCCACGTCGGTTATTCAGTCGGCGTGAAGGCGATCGATAAAGTGCGAGAGAACTCGCAGATGATTGCCCGTTCGGCAACCATCCTCGGGATCACGGTAATCGGCGGGCTGATCGCTTCGTATGTGCATAtt                                               <  1:819912/239‑1 (MQ=255)
                                                                                                                                                                                           ttATGGCGGGAATTTGGTCATCATTTGCCAGCCAGGGAACCCTGC‑TGGGGGCGATTCTATTTTTCGCCGTTTACCTGCTTATCTTTTTCCTGCGCGTCGGCTGGACCCACGTCGGTTATTCAGTCGGCGTGAAGGCGATCGATAAAGTGCGAGAGAACTCGCAGATGATTGCCCGTTCGGCAACCATCCTCGGGATCACGGTAATCGGCGGGCTGATCGCTTCGTATGTGCATATTaa                                             <  2:767665/238‑1 (MQ=255)
                                                                                                                                                                                                  gggAATTTGCTCATCATTTTCCAGCCAGGGAATCCTGC‑TGGTGCCGATTCTATTTTTCGCAGTTTACCTGCTTATCTTTTTCCTGCGCGTCGGCGGGACCCACGTCGGTTATTCAGTCGGCGTGAAGGCGATCGATAAAGTGCGAGAGAACTCGCAGATGATTGCCCGTTCGGCAACCATCCTCGGGATCACGGTAATCGGCGGGCTGATCGCTTCGTATGTGCATATTAACGTGGTGa                                     <  2:880709/239‑1 (MQ=255)
                                                                                                                                                                                                        ttGCTCATCATTTGCCAGCCAGGGAAACCTGC‑TGGGGCCGATTCTATTTTTCGCCGTTTACCTGCTTATCTTTTTCCTGCGCGTCGGCTGGACCCACGTCGGTTATTCAGTCGGCGTGAAGGCGATCGATAAAGTGCGAGAGAACTCGCAGATGATTGCCCGTTCGGCAACCATCCTCGGGATCACGGTAATCGGCGGGCTGATCGCTTCGTATGTGCATATTAACGTGGTGACATCGt                               <  1:335367/239‑1 (MQ=255)
                                                                                                                                                                                                             catcctTTGCCAGCCAGGGAAACCTGC‑TGGGGCCGATTCTATTTTTCGCCGTTTACCTGCTTATCTTTTTCCTGCGCGTCGGCTGGACCCACGTCGGTTATTCAGTCGTCGTGAAGGCGATCGATACAGTGCGAGAGAACTCGCAGATGATTGCCCGTTCGGCAACCATCCTCGGGATCACGTTAATCGGCGGGCTGATCGCTTCGTATGTGCATATTTACTGGGTGACATCTTGTGc                           >  2:720967/1‑238 (MQ=255)
                                                                                                                                                                                                               tcatTTGCCAGCCAGGGAAACCTGC‑TGGGGCCGATTCTATTTTTCGCCGTTTACCTGCTTATCTTTTTCCTGCGCGTCGGCTGGACCCACGTCGGTTATTCAGTCGGCGTGAAGGCGATCGATAAAGTGCGAGAGAACTCGCAGATGATTGCCCGTTCGGCAACCATCCTCGGGATCACGGTAATCGGCGGGCTGATCGCTTCGTATGTGCATATTAACGTGGTGACATCGTTTGCCAt                        <  1:445312/239‑1 (MQ=255)
                                                                                                                                                                                                               tcatTTGCCAGCCAGGGAAACCTGC‑TGGGGCCGATTCTATTTTTCGCCGTTTACCTGCTTATCTTTTTCCTGCGCGTCGGCTGGACCCACGTCGGTTATTCAGTCGGCGTGAAGGCGATCGATAAAGTGCGAGAGAACTCGCAGATGATTGCCCGTTCGGCAACCATCCTCGGGATCACGGTAATCGGCGGGCTGATCGCTTCGTATGTGCATATTAACGTGGTGACATCGTTTGCCAt                        <  1:949865/239‑1 (MQ=255)
                                                                                                                                                                                                                     gccagccaGGGAAACCTGC‑TGGGGCCGATTCTATTTTTCGCCGTTTACCTGCTTATCTTTTTCCTGCGCGTCGGCTGGACCCACGTCGGTTATTCAGTCGGCGTGAAGGCGATCGATAAAGTGCGAGAGAACTCGCAGATGATTGCCCGTTCGGCAACCATCCTCGGGATCACGGTAATCGGCGGGCTGAt                                                                  <  1:566656/191‑1 (MQ=255)
                                                                                                                                                                                                                     gccagccaGGGAAACCTGC‑TGGGGCCGATTCTATTTTTCGCCGTTTACCTGCTTATCTTTTTCCTGCGCGTCGGCTGGACCCACGTCGGTTATTCAGTCGGCGTGAAGGCGATCGATAAAGTGCGAGAGAACTCGCAGATGATTGCCCGTTCGGCAACCATCCTCGGGATCACGGTAATCGGCGGGCTGAt                                                                  >  2:566656/1‑191 (MQ=255)
                                                                                                                                                                                                                     gccagccaGGGAAACCTGC‑TGGGGCCGATTCTATTTTTCGCCGTTTACCTGCTTATCTTTTTCCTGCGCGTCGGCTGGACCCACGTCGGTTATTCAGTCGGCGTGAAGGCGATCGATAAAGTGCGAGAGAACTCGCAGATGATTGCCCGTTCGGCAACCATCCTCGGGATCACGGTAATCGGCGGGCTGATCGCTTCGTATGTGCCTATTAACGCGGTGACATCGTTTTCCATcaacaa                  >  1:701277/1‑239 (MQ=255)
                                                                                                                                                                                                                       cagccagGGAAACCTGC‑TGGGGCCGATTCTATTTTTCGCCGTTTACCTGCTTATCTTTTTCCTGCGCGTCGGCTGGACCCACGTCGGTTATTCAGTCGGCGTGAAGGCGATCGATAAAGTGCGAGAGAACTCGCAGATGATTGCCCGTTCGGCAACCATCCTCGGGATCACGGTAATCGGCGGGCTGATCGCTTCGTATGTGCATATTAACGTGGTGACATCGTTTGCCATCGACAAt                 <  1:319092/238‑1 (MQ=255)
                                                                                                                                                                                                                        agccagGGAAACCTGC‑TGGGGCCGATTCTATTTTTCGCCGTTTACCTGCTTATCTTTTTCCTGCGCGTCGGCTGGACCCACGTCGGTTATTCAGTCGGCGTGAAGGCGATCGATAAAGTGCGAGAGAACTCGCAGATGATTGCCCGTTCGGCAACCATCCTCGGGATCACGGTAATCGGCGGGCTGATCGCTTCGTATGTGCATATTAACGTGGTGACATCGTTTGCCATCGACAATa                >  1:333276/1‑238 (MQ=255)
                                                                                                                                                                                                                        agccagGGAAACCTGC‑TGGGGCCGATTCTATTTTTCGCCGTTTACCTGCTTATCTTTTTCCTGCGCGTCGGCTGGACCCACGTCGGTTATTCAGTCGGCGTGAAGGCGATCGATAAAGTGCGAGAGAACTCGCAGATGATTGCCCGTTCGGCAACCATCCTCGGGATCACGGTAATCGGCGGGCTGATCGCTTCGTATGTGCATATTAACGTGGTGACATCGTTTGCCATCGACAATa                >  1:374297/1‑238 (MQ=255)
                                                                                                                                                                                                                          ccagGGAAACCTGC‑TGGGGCCGATTCTATTTTTCGCCGTTTACCTGCTTATCTTTTTCCTGCGCGTCGGCTGGACCCACGTCGGTTATTCAGTCGGCGTGAAGGCGATCGATAAAGTGCGAGAGAACTCGCAGATGATTGCCCGTTCGGCAACCATCCTCGGGATCACGGTAATCGGCGGGCTGATCGCTTCGTATGTGCATATTAACGTGGTGACATCGTTTGCCATCGAc                    <  1:1147726/232‑1 (MQ=255)
                                                                                                                                                                                                                          ccagGGAAACCTGC‑TGGGGCCGATTCTATTTTTCGCCGTTTACCTGCTTATCTTTTTCCTGCGCGTCGGCTGGACCCACGTCGGTTATTCAGTCGGCGTGAAGGCGATCGATAAAGTGCGAGAGAACTCGCAGATGATTGCCCGTTCGGCAACCATCCTCGGGATCACGGTAATCGGCGGGCTGATCGCTTCGTATGTGCATATTAACGTGGTGACATCGTTTGCCATCGAc                    >  2:1147726/1‑232 (MQ=255)
                                                                                                                                                                                                                          ccagGGAAACCTGC‑TGGGGCCGATTCTATTTTTCGCCGTTTACCTGCTTATCTTTTTCCTGCGCGTCGGCTGGACCCACGTCGGTTATTCAGTCGGCGTGAAGGCGATCGATAAAGTGCGAGAGAACTCGCAGATGATTGCCCGTTCGGCAACCATCCTCGGGATCACGGTAATCGGCGGGCTGATCGCTTCGTATGTGCATATTAACGTGGTGACATCGTTTGCCATCGACAATAccc             >  1:694073/1‑239 (MQ=255)
                                                                                                                                                                                                                                   ccTGC‑TGGGGCCGCTTCTATTTTTCGCCGTTTACCTGCTTATCTTTTTCCTGCGCGTCGGCTGGACCCACGTCGGTTATTCAGTCGGCGTGAAGGCGATCGATAAAGTGCGAGAGAACTTGCAGATGATTGCCCGTTCGGCAACCATCCTCGGGATCAGGGTAATCGGCGGGATGATCGCTTCGCTTGTGCATATTATCGTGTGGACATCGGTTGTCCTCGACAATAccccacgccttg    >  2:902434/1‑230 (MQ=255)
                                                                                                                                                                                                                                     tgc‑tgGGGCCGATTCTATTTTTCGCCGTTTACCTGCTTATCTTTTTCCTGCGCGTCGGCTGGACCCACGTCGGTTATTCAGTCGGCGTGAAGGCGATCGATAAAGTGCGAGAGAACTCGCAGATGATTGCCCGTTCGGCAACCATCCTCGGGATCACGGTAATCGGCGGGCTGATCGCTTCGTATGTGCATATTAACGTGGTGACATCGTTTGCCATCGACAATACCCACAGCGTTgcg  <  1:613108/239‑1 (MQ=255)
                                                                                                                                                                                                                                      gc‑tgGGGCCGATTGTATTTATCGCCGTTTGCCTGCTTATCTTTTTCCTGCGCGTCCGCTGGACCCACGTCGGTTATTCAGTCGGCGTCAAGGCGATCGATAAAGTGCGAGCGAACTCGCAGATGATTGCCCGTTCGGCAACCATCCTAGGGATCTCGGCAATCGGCGGGCTGACTCCTTCGTATGAGAATCTTATCGAGCTGACATCGTGTGCCATCGACAATACCCACAGCGTTgc   >  2:366787/1‑237 (MQ=255)
                                                                                                                                                                                                                                              |                                                                                                                                                                                                                                      
GCGGCGATGAAAGATAACCTCGAATTTATTAATACCCACCCGAATCTGGTCGGATTCCTGATGGGGTTATTAATTTCGATGGAAGAAAAAGGAGAAAACCGCGACACCATTAAAGGCCTCAAAGTGGCACTGTTTGGCCCAATCGCCGGGATTGGCGATGCGATTTTCTGGTTTACTTTGTTGCCGATTATGGCGGGAATTTGCTCATCATTTGCCAGCCAGGGAAACCTGC‑TGGGGCCGATTCTATTTTTCGCCGTTTACCTGCTTATCTTTTTCCTGCGCGTCGGCTGGACCCACGTCGGTTATTCAGTCGGCGTGAAGGCGATCGATAAAGTGCGAGAGAACTCGCAGATGATTGCCCGTTCGGCAACCATCCTCGGGATCACGGTAATCGGCGGGCTGATCGCTTCGTATGTGCATATTAACGTGGTGACATCGTTTGCCATCGACAATACCCACAGCGTTGCG  >  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: ‑