Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,285,875 C→G 16.4% 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,8750CG16.4% 144.5 / 5.1 64P134A (CCG→GCG) agaDN‑acetylgalactosamine‑specific enzyme IID component of PTS
Reads supporting (aligned to +/- strand):  ref base C (14/38);  new base G (9/2);  total (23/41)
Fisher's exact test for biased strand distribution p-value = 1.11e-03
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GCGGCGATGAAAGATAACCTCGAATTTATTAATACCCACCCGAATCTGGTCGGATTCCTGATGGGGTTATTAATTTCGATGGAAGAAAAAGGAGAAAACCGCGACACCATTAAAGGCCTCAAAGTGGCACTGTTTGGCCCAATCGCCGGGATTGGCGATGCGATTTTCTGGTTTACTTTGTTGCCGATTATGGC‑GGGAATTTGCTCATCATTTGCCAGCCAGGGAAACCTGC‑TGGGGCCGATTCTATTTTTCGCCGTTTACCTGCTTATCTTTTTCCTGCGCGTCGGCTGGACCCACGTCGGTTATTCAGTCGGCGTGAAGGCGATCGATAAAGTGCGAGAGAACTCGCAGATGATTGCCCGTTCGGCAACCATCCTCGGGATCACGGTAATCGGCGGGCTGATCGCTTCGTATGTGCATATTAACGTGGTGACATCGTTTGCCATCGACAATACCCACAGCGTTGCG  >  NC_000913/3285638‑3286105
                                                                                                                                                                                                                                               |                                                                                                                                                                                                                                      
gcggcgATGAAAGATAACCTCGAATTTATTAATACCCACCCGAATCTGGTCGGATTCCTGATGGGGTTATTAATTTCGATGGAAGAAAAAGGAGAAAACCGCGACACCATTAAAGGCCTCAAAGTGGCACTGTTTGGCCCAATCGCCGGGATTGGCGATGCGATTTTCTGGTTTACTTTGTTGCCGATTATGGC‑GGGAATTTGCTCATCATTTGCCAGCCAGGGAAACCTGC‑TGGGGc                                                                                                                                                                                                                                        <  1:317352/238‑1 (MQ=255)
              tAACCTCGAATTTATTAATACCCACCCGAATCTGGTCGGATTCCTGATGGGGTTATTAATTTCGATGGAAGAAAAAGGAGAAAACCGCGACACCATTAAAGGCCTCAAAGTGGCACTGTTTGGCCCAATCGCCGGGATTGGCGATGCGATTTTCTGGTTTACTTTGTTGCCGATTATGGC‑GGGAATTTGCTCATCATTTGCCAGCCAGGGAAACCTGC‑TGGGGCCGATTCTAttttt                                                                                                                                                                                                                           <  1:416279/237‑1 (MQ=255)
              tAACCCCGAATTTATTAAAACCCACCCGAATCTGGTCGGATTCCTGATGGGGTTAGTAAGTTCGATGGAAGAAAAAGGAGAAAACCGCGACACCATTAAAGGCCTCAAAGTGGCACTGTTTGGCCCAATCGCCGGGATTGGCGATGCGATTTTCTGGTTTACTTTGTTGCCGATTATGGC‑GGGAATTTGCTCATCATTTGCCAGCCAGGGAAACCTGC‑TGGGGCCGATTCTAttttt                                                                                                                                                                                                                           <  1:984156/237‑1 (MQ=255)
                aCCTCGAATTTATTAATACCCACCCGAATCTGGTCGGATTCCTGATGGGGTTATTAATTTCGATGGAAGAAAAAGGAGAAAACCGCGACACCATTAAAGGCCTCAAAGTGGCACTGTTTGGCCCAATCGCCGGGATTGGCGATGCGATTTTCTGGTTTACTTTGTTGCCGATTATGGC‑GGGGATTTGGTCATTATTTTTCAGCCCGGGGAAACTTC‑TGGGGGCGGTTTTTTTTTTTGcc                                                                                                                                                                                                                       >  1:172346/1‑239 (MQ=255)
                    cGAATTTATTAATACCCACCCGAATCTGGTCGGATTCCTGATGGGGTTATTAATTTCGATGGAAGAAAAAGGAGAAAACCGCGACACCATTAAAGGCCTCAAAGTGGCACTGTTTGGCCCAATCGCCGGGATTGGCGATGCGATTTTCTGGTTTACTTTGTTGCCGATTATGGC‑GGGAATTTGCTCATCATTTGCCAGCCAGGGAAACCTGC‑TGGGGCCGATTCTATTTTTCGCCGttt                                                                                                                                                                                                                   <  1:157058/239‑1 (MQ=255)
                      gacagattacagccccgccGAATGTGGTCGGATTCCTGATGGGGTGATTAACTTCGAAGGAAGAAAAAGGAGAAAACCGCGACCCCATTAAAGGCCTCAAAGTAGCACTGTTTGGCCCAATCGCCGGGATTGGCGATGCGATTTTCTGGTTTAATTTGTTGCCGATTATGGC‑GAGAATTTGCTCATCATTTGCCAGCCAGGGAAACCTGC‑TGGTGCCGATTCTATTTTTCGCCGTTTAc                                                                                                                                                                                                                 <  2:238489/222‑1 (MQ=255)
                      aaTTTATTAATACCCACCCGAATCTGGTCGGATTCCTGATGGGGTTATTAATTTCGATGGAAGAAAAAGGAGAAAACCGCGACACCATTAAAGGCCTCAAAGTGGCACTGTTTGGCCCAATCGCCGGGATTGGCGATGCGATTTTCTGGTTTACTTTGTTGCCGATTATGGC‑GGGAAATTTCTCATCGTTTTGCATCCAAGGGAACACGC‑TGGGGGCGGTTCTtttttttttcccttt                                                                                                                                                                                                                  >  1:199853/1‑229 (MQ=255)
                                     acccGAATCTGGTCGGATTCCTGATGGGGTTATTAATTTCGATGGAAGAAAAAGGAGAAAACCGCGACACCATTAAAGGCCTCAAAGTGGCACTGTTTGGCCCAATCGCCGGGATTGGCGATGCGATTTTCTGGTTTACTTTGTTGCCGATTATGGC‑GGGGATTTGTTTTTCTTTTTCCCTCCCCGGGAAACTGC‑TTGGGGCGGGTTTTTTTTTTTCTGTTTTCCTTCTTTTCtttttt                                                                                                                                                                                                  >  2:572863/1‑238 (MQ=255)
                                     acccGAATCTGGTCGGATTCCTGATGGGGTTATTAATTTCGATGGAAGAAAAAGGAGAAAACCGCGACACCATTAAAGGCCTCAAAGTGGCACTGTTTGGCCCAATCGCCGGGATTGGCGATGCGATTTTCTGGTTTACTTTGTTGCCGATTATGGC‑GGGAGTTTGTTCATCCTTTGCCACCCAGGGAAACATGC‑TGGGGCCGCTTATCTTTTTTTCCGGTTATCTGCTTATCtttttt                                                                                                                                                                                                  >  2:1104478/1‑238 (MQ=255)
                                     acccGAATCTGGTCGGATTCCTGATGGGGTTATTAATTTCGATGGAAGAAAAAGGAGAAAACCGCGACACCATTAAAGGCCTCAAAGTGGCACTGTTTGGCCCAATCGCCGGGATTGGCGATGCGATTTTCTGGTTTACTTTGTTGCCGATTATGGC‑GGGAATATGCTCCTCCTTTTCCTACCACCGGAACATGC‑TGGGGGCGCTTCTATTTTTTCCCGTTTTGCCTCTTCTCtttttt                                                                                                                                                                                                  >  2:335367/1‑238 (MQ=255)
                                      cccGAATCTGGTCGGATTCCTGATGGGGTTATTAATTTCGATGGAAGAAAAAGGAGAAAACCGCGACACCATTAAAGGCCTCAAAGTGGCACTGTTTGGCCCAATCGCCGGGATTGGCGATGCGATTTTCTGGTTTACTTTGTTGCCGATTATGGC‑GGGGATTTGCTTATCTTTTTTCCCCCAGGGAAACCTTC‑TGGGGGCGTTTCTATTTTTTGCCGCTTTTCCGCtttttttttt                                                                                                                                                                                                   >  1:453465/1‑237 (MQ=255)
                                      cccGAATCTGGTCGGATTCCTGATGGGGTTATTAATTTCGATGGAAGAAAAAGGAGAAAACCGCGACACCATTAAAGGCCTCAAAGTGGCACTGTTTGGCCCAATCGCCGGGATTGGCGATGCGATTTTCTGGTTTACTTTGTTGCCGATTATGGC‑GGCGACTTTTTCAACTCTTTCCATTCCGGGGAAACTGC‑TAGGGGCGGTTCTTTTTTTTGTCCTTTATCCGCCTGTCttttttt                                                                                                                                                                                                 >  1:238489/1‑237 (MQ=255)
                                          aaTCTGGTCGGATTCCTGATGGGGTTATTAATTTCGATGGAAGAAAAAGGAGAAAACCGCGACACCATTAAAGGCCTCAAAGTGGCACTGTTTGGCCCAATCGCCGGGATTGGCGATGCGATTTTCTGGTTTACTTTGTTGCCGATTATGGC‑GGGGATTTGCTCATCATTTTCTAGCCCTGGGAGCCTGC‑TTGGGGCCCGTGTATTTTTTTCCCTTTTCCTTCTTATCTTTTTTCtccc                                                                                                                                                                                             >  2:597869/1‑236 (MQ=255)
                                           aTCTGGTCGGATTCCTGATGGGGTTATTAATTTCGATGGAAGAAAAAGGAGAAAACCGCGACACCATTAAAGGCCTCAAAGTGGCACTGTTTGGCCCAATCGCCGGGATTGGCGATGCGATTTTCTGGTTTACTTTGTTGCCGATTATGGC‑GGGAATTTGCTCATCATTTGCCAGCCAGGGAAACCTGC‑TGGGGCCGATTCTATTTTTCGCCGTTTACCTGCTTATCTTTTTCCTgcgc                                                                                                                                                                                            <  1:901008/239‑1 (MQ=255)
                                             cTGGGCGGATTCCGGATGGGGTTATTAATTGCGATGGAAGAAAAAGGAGAAATCCGCGACACCATTAAAGGCCGAAAAGTGGCACTGTTTGGCCCAATCGCCGGGATTGGCGATGCGATTTTCTGGTTTACTTTGTTGCCGATTATGGC‑GGGAATTTGCTCATCATTTGCCAGCCAGGGAAACCTGC‑TGGGGCCGATTCTATTTTTCGCCGTTTACCTGCTTATCTTTTTCCTGCGCGt                                                                                                                                                                                          <  2:450217/239‑1 (MQ=255)
                                                          tGATGGGGTTATTAATTTCGATGGAAGAAAAAGGAGAAAACCGCGACACCATTAAAGGCCTCAAAGTGGCACTGTTTGGCCCAATCGCCGGGATTGGCGATGCGATTTTCTGGTTTACTTTGTTGCCGATTATGGC‑GGGAATTTGCTCATCATTTGCCAGCCAGGGAAACCTGC‑TGGGGCCGATTCTATTTTTCGCCGTTTACCTGCTTATCTTTTTCCTGCGCGTCGGCTGGACCCAc                                                                                                                                                                             <  1:572863/239‑1 (MQ=255)
                                                                    attaATTTCGATGGAAGAAAAAGGAGAAAACCGCGACACCATTAAAGGCCTCAAAGTGGCACTGTTTGGCCCAATCGCCGGGATTGGCGATGCGATTTTCTGGTTTACTTTGTTGCCGATTATGGC‑GGTGATTTTCTCATCACTTGCCCGCCCGGGGAACCTGCTTGGGGCCG‑TTCTATTTTTTCCCGTTTAACTCTCTATCTTTTTTCTGCGCGGGGGTGGGAACCACGTCCGTtttt                                                                                                                                                                   >  1:1084748/1‑239 (MQ=255)
                                                                          ttCGATGGAAGAAAAAGGAGAAAACCGCGACACCATTAAAGGCCTCAAAGTGGCACTGTTTGGCCCAATCGCCGGGATTGGCGATGCGATTTTCTGGTTTACTTTGTTGCCGATTATGGC‑GGGAATTTTCTCATCCTTTTCCAGCCACGGAAACCTGC‑TGGGGCCGGTTCTTTTTTTTGTCGCTTTCCTGCTTATTTTTTTTTTGCGCGCTGGCTGGACCCACGTCGGGTTTTCAGTCg                                                                                                                                                             >  1:865003/1‑239 (MQ=255)
                                                                           tCGATGGAAGAAAAAGGAGAAAACCGCGACACCATTAAAGGCCTCAAAGTGGCACTGTTTGGCCCAATCGCCGGGATTGGCGATGCGATTTTCTGGTTTACTTTGTTGCCGATTATGGCGGGGAATTTG‑TCATCATTTGTCAACCCGGGGAAACTGC‑TGGGGGGGGTTATTTTTTTTGTCGTTTTCCTGCTTTTCTTTTTTCTGCTGCTCCGCCGGACTCACgcacgtgttttagttcc                                                                                                                                                            >  2:445312/1‑223 (MQ=255)
                                                                            cGATGGAAGAAAAAGGAGAAAACCGCGACACCATTAAAGGCCTCAAAGTGGCACTGTTTGGCCCAATCGCCGGGATTGGCGATGCGATTTTCTGGTTTACTTTGTTGCCGATTATGGC‑GGGAAATTGCTCATTATTTTCCCCCCCGGGAAACCTGC‑TTGGGGCGGTTTTATTTTTTGCCGTTTTCCCTCTTATTTTTTTTCTGCGCGTCGGCTGGGCCCCCGCCCGTttttttctcgt                                                                                                                                                            >  1:289789/1‑231 (MQ=255)
                                                                               tGGAAGAAAAAGGAGAAAACCGCGACACCATTAAAGGCCTCAAAGTGGCACTGTTTGGCCCAATCGCCGGGATTGGCGATGCGATTTTCTGGTTTACTTTGTTGCCGATTATGGC‑GGGGAATTGTTCATCCTTTGCCCGCCATGGGAACCCGC‑TGGGGCCGAGTCTATTTTTTGCCGCTTTCCTGCCTTTTTTTTTTCTTCTCGTCGGCTGGGACCACGCCGGTTTTTCTCTTGGCGTg                                                                                                                                                        >  1:767665/1‑239 (MQ=255)
                                                                                 gaagaaAAAGGCGAAAACCGCGACACCATTAAAGGCCTCAAAGTGGCACTGTTTGGCCCAATCGCCGGGATTGGCGATGCGATTTTCTGGTTTACTTTGTTGCCGATTATGGC‑GGGAATTTGCTCATCATTTGCCAGCCAGGGAAACCTGC‑TGGGGCCGATTCTATTTTTCGCCGTTTACCTGCTTATCTTTTTCCTGCGCGTCGGCTGGACCCACGTCGGTTATTCAGTCGGCGTGa                                                                                                                                                       <  2:126848/238‑1 (MQ=255)
                                                                                      aaaaGGAGAAAACCGCGACACCATTAAAGGCCTCAAAGTGACACTGTTTGGCCCAATCGCCGGGATTGGCGATGCGATTTGCTGGTTTACTTTGTTGCCGATTATGGC‑GGGAATTTGCTCATCATTTGCCAGCCAGGGAAACCTGC‑TGGGGCCGATTCTATTTTTCGCCGTTTCCCTGCTTATCTTTTTCCTGCGCGTCGGCTGGACCCACGTCGGTTATTCAGTCGGCGTGAAGGc                                                                                                                                                   <  1:1054598/237‑1 (MQ=255)
                                                                                             gaAAACCGCGACACCATTAAAGGCCTCAAAGTGGCACTGTTTGGCCCAATCGCCGGGATTGGCGATGCGATTTTCTGGTTTACTTTGTTGCCGATTATGGC‑GGGAATTTGCTCATCATTTGCCAGCCAGGGAAACCTGC‑TGGGGCCGATTCTATTTTTCGCCGTTTACCTGCTTATCTTTTTCCTGCGCGTCGGCTGGACCCACGTCGGTTATTCAGTCGGCGTGAAGGCGATCGATa                                                                                                                                           <  1:668328/238‑1 (MQ=255)
                                                                                                              tAAAGGCCTCAAAGTGGCTCTGTTTCGCTCAATCGCCGGGATTGGCGATGCGATTTTCTGGTTTACTTTGTTTCCGATTATGGC‑GGGAATTTGCTCATCATTTGCCAGCCAGGGAAACCTGC‑TGGTGCCGATTCTATTTTTCGCCGTTTACCTGCTTATCTTTTTCCTGCGCGTCGGCTGGACCCACGTCGGTTATTCAGTCGGCGTGAAGGCGATCGATAAAGTGCGAGAGAACTCGc                                                                                                                         <  2:664892/239‑1 (MQ=255)
                                                                                                               aaaGGCCTCAAAGTGGCGCTGTTTGGCCCAATCGCCTGGATTGGCGATGCGATTTTCTGGTTTGCGTTGGTGCCGATTATGGC‑GGGAATTTGCTCATCATTTGCCAGCCAGGGAAACCTGC‑TGGGGCCGATTCTATTTTTCGCCGTTTACCTGCTTATCTTTTTCCTGCGCGTCGGCTGGACCCACGTCGGTTATTCAGTCGGCGTGAAGGCGATCGATAAAGTGCGAGAGAACTCGCa                                                                                                                        <  2:721168/239‑1 (MQ=255)
                                                                                                               aaaGGCCTCAAAGTGGCACTGTTTGGCCCAATCGCCGGGATTGGCGATGCGATTTTCTGGTTTACTTTGTTGCCGATTATGGC‑GGGAATTTGCTCATCATTTGCCAGCCAGGGAAACCTGC‑TGGGGCCGATTCTATTTTTCGCCGTTTACCTGCTTATCTTTTTCCTGCGCGTCGGCTGGACCCACGTCGGTTATTCAGTCGGCGTGAAGGCGATCGATAAAGTGCGAGAGAACTc                                                                                                                           <  1:609755/236‑1 (MQ=255)
                                                                                                                     cTCAAAGTGGCACTGTTTGGCCCAATCGCCGGGATTGGCGATGCGATTTTCTGGTTTACTTTGTTGCCGATTATGGC‑GGGAATTTTCTCATCATTTTTCCGCCCGGGGAACCTGC‑TTGGGCCGATTCTATTTTTTTCCGTTTTTCTGCTTTTCTTTTTTCTTCGCGTCGGCTGGACCCCCCTCGGTTATTCTGTCGGCGTGAAGGCGATCGATAAAATGCGGGAGAACTCGCAgatgat                                                                                                                  >  1:1023607/1‑239 (MQ=255)
                                                                                                                        aaaGTGGCACTGTTTGGCCCAATCGCCGGGATTGGCGATGCGATTTTCTGGTTTACTTTGTTGCCGATTATGGC‑GGGAATTTGCTCATCATTTGCCAGCCAGGGAAACCTGC‑TGGGGCCGATTCTATTTTTCGCCGTTTACCTGCTTATCTTTTTCCTGCGCGTCGGCTGGACCCACGTCGGTTATTCAGTCGGCGTGAAGGCGATCGATAAAGTGCGAGAGAACTCGCAGATGATTGc                                                                                                               <  1:756322/239‑1 (MQ=255)
                                                                                                                          aGTGGCACTGTTTGGCCCAATCGCCGGGATTGGCGATGCGATTTTCTGGTTTACTTTGTTGCCGATTATGGC‑GGGAATTTGCTCATCATTTGCCAGCCAGGGAAACCTGC‑TGGGGCCGATTCTATTTTTCGCCGTTTACCTGCTTATCTTTTTCCTGCGCGTCGGCTGGACCCAc                                                                                                                                                                             <  1:119103/175‑1 (MQ=255)
                                                                                                                                   gTTTGGCCCAATCGCCGGGATTGGCGATGCGATTTTCTGGTTTACTTTGTTGCCGATTATGGC‑GGGAATTTGCTCATCATTTGCCAGCCAGGGAAACCTGC‑TGGGGCCGATTCTATTTTTCGCCGTTTACCTGCTTATCTTTTTCCTGCGCGTCGGCTGGACCCACGTCGGTTATTCAGTCGGCGTGAAGGCGATCGATAAAGTGCGAGAGAACTCGCAGATGATTGCCCGTTCGGCaa                                                                                                    <  1:259864/239‑1 (MQ=255)
                                                                                                                                     ttGGCCCAATCGCCGGGATTGGCGATGCGATTTTCTGGTTTACTTTGTTGCCGATTATGGC‑GGGAATTTGCTCATCATTTGCCAGCCAGGGAAACCTGC‑TGGGGCCGATTCTATTTTTCGCCGTTTACCTGCTTATCTTTTTCCTGCGCGTCGGCTGGACCCACGTCGGTTATTCAGTCGGCGTGAAGGCGATCGATAAAGTGCGAGAGAACTCGCAGATGATTGCCCGTTCGGCAAcc                                                                                                  <  1:1104478/239‑1 (MQ=255)
                                                                                                                                                              tgcgatTTTCTGGTTTACTTTGTTGCCGATTATGGC‑GGGAATTTGCTCATCATTTGCCAGCCAGGGAAACCTGC‑TGGGGCCGATTCTATTTTTCGCCGTTTACCTGCTTATCTTTTTCCTGCGCGTCGGCTGGACCCACGGCGGTTATTCAGTCGGCGTGAAGGCGATCGATAAAGTGCGAGAGAACTCGCAGATGATTGCCCGTTCGGCAACCATCCTCGGGATCACGGTAATcggc                                                                          <  1:446358/238‑1 (MQ=255)
                                                                                                                                                                    tttCTGGTTTACTGTGTTGCCGATTATGGC‑GGGAATTTGCTCACCATTTGCCGGCCAGTGAAACCTGC‑TGGGGCAGATTCTATTTTTCGCCGTTTACCTGCTTATCTTTTTGCTGCGCGTCGGCTGGACCCACGTCGGTTATTCAGTCGGCGTGAAGGCGATCGATAAAGTGCGAGAGAACTCGCAGATGATTGCCCGTTCGGCAACCATCCTCGGGATCACGGTAATCGGCGGGCTg                                                                    <  2:453465/238‑1 (MQ=255)
                                                                                                                                                                     ttCTGCTTTACTTGGTTGCCGATTATGGC‑GGGAATGTGCTTATCATTTCCCAGGCAGGGAACCCTGC‑TGGTGTAGATTCGATTTTTCGCCGTTTACCTGCTTATCTCTTTCCTTCGCGTCGGCTGGACCCACGCCGGTTATTCAGTCGGCGTGAAGGCGATCGATTAAGTGCGAGAGAACTCGCAGATGATTGCGCGTTCGGCAACCATCCGCGAGAACACAGTAATCCGCGGGCTCAt                                                                  <  2:879544/239‑1 (MQ=255)
                                                                                                                                                                       gtgttttACTTTGTTGCCGATTATGGC‑GAGAATTTGCTCATCATTTGCCAGCCAGTGAATCCTGC‑TGGGGGCGATTCTATTTTTCGCCGTTTACCTGCGTATCTTTTTCCTGCGCTTCGGCTGGACCCACGTCGGTTATTCAGTCGGCGTGAAGGCGATCGATAAAGTGCGAGAGAACTCGCAGATGATTGCCCGTTCGGCAACCATCCTCGGGATCACGGTAATCGGCGGGCTGATCg                                                                <  2:430205/235‑1 (MQ=255)
                                                                                                                                                                         ggTTTACTTTGTTGCCGATTATGGC‑GGGAATTTGCTCATCATTTGCCAGCCAGGGAAACCTGC‑TGGGGCCGATTCTATTTTTCGCCGTTTACCTGCTTATCTTTTTCCTGCGCGTCGGCTGGACCCACGTCGGTTATTCAGTCGGCGTGAAGGCGATCGATAAAGTGCGAGAGAACTCGCAGATGATTGCCCGTTCGGCAACCATCCTCGGGATCACGGTAATCGGCGGGCTGATCGCt                                                              <  1:973120/239‑1 (MQ=255)
                                                                                                                                                                           tttACTTTGTTGCCGATTATGGC‑GGGAATTTGCTCATCATTTGCCAGCCAGGGAAACCTGC‑TGGGGCCGATTCTATTTTTCGCCGTTTACCTGCTTATCTTTTTCCTGCGCGTCGGCTGGACCCACGTCGGTTATTCAGTCGGCGTGAAGGCGATCGATAAAGTGCGAGAGAACTCGCAGATGATTGCCCGTTCGGCAACCATCCTCGGGATCACGGTAATCGGCGGGCTGATCGCtt                                                             <  1:141301/238‑1 (MQ=255)
                                                                                                                                                                             tACTTTGTTGCCGATTATGGC‑GGGAATTTGCTCATCATTTGCCAGCCAGGGAAACCTGC‑TGGGGCCGATTCTATTTTTCGCCGTTTACCTGCTTATCTTTTTCCTGCGCGTCGGCTGGACCCACGTCGGTTATTCAGTCGGCGTGAAGGCGATCGATAAAGTGCGAGAGAACTCGCAGATGATTGCCCGTTCGGCAACCATCCTCGGGATCACGGTAATCGGCGGGCTGATCGCTTCGt                                                          <  1:597869/239‑1 (MQ=255)
                                                                                                                                                                             gacttgtgtgCCGTTTATGGC‑GGAAATGTGCTCATCTTTGGCCAGCCAGGGGAACCTGC‑TTGGGCCGATTCTATTTTTCGCCGTTTTCCTGCTTATCTTTGTCCGGCCCGTCGGCTGGACCCACGTCGGTTATTCAGTCGGCGTGAAGGCGATCGATAAAGTGCGAGAGAACTCGCAGATGATTGCCCGTTCGGCAACCATCCTCGGGATCACGGTAAGCGGCGGGCTGAGCGCTGCGt                                                          <  2:1083752/231‑1 (MQ=255)
                                                                                                                                                                               cTTTGTTGCCGATTATGGC‑GGGAATTTGCTCATCATTTGCCAGCCAGGGAAACCTGC‑TGGGGCCGATTCTATTTTTCGCCGTTTACCTGCTTATCTTTTTCCTGCGCGTCGGCTGGACCCACGTCGGTTATTCAGTCGGCGTGAAGGCGATCGATAAAGTGCGAGAGAACTCGCAGATGATTGCCCGTTCGGCAACCATCCTCGGGATCACGGTAATCGGCGGGCTGATCGCTTCGt                                                          <  1:1024574/237‑1 (MQ=255)
                                                                                                                                                                                tttGTTGCCGATTATGGC‑GGGAATTTGCTCATCATTTGCCAGCCAGGGAAACCTGC‑TGGGGCCGATTCTATTTTTCGCCGTTTACCTGCTTATCTTTTTCCTGCGCGTCGGCTGGACCCACGTCGGTTATTCAGTCGGCGTGAAGGCGATCGATAAAGTGCGAGAGAACTCGCAGATGATTGCCCGTTCGGCAACCATCCTCGGGATCACGGTAATCGGCGGGCTGATCGCTTCGTAtg                                                       <  1:471160/239‑1 (MQ=255)
                                                                                                                                                                                 ttgttgCCGATTATGGC‑GGGAATTTGCTCATCATTTGCCAGCCAGGGAAACCTGC‑TGGGGCCGATTCTATTTTTCGCCGTTTACCTGCTTATCTTTTTCCTGCGCGTCGGCTGGACCCACGTCGGTTATTCAGTCGGCGTGAAGGCGATCGATAAAGTGCGAGAGAACTCGCAGATGATTGCCCGTTCGGCAACCATCCTCGGGATCACGGTAATCGGCGGGCTGATCGCTTCGTAtgt                                                      <  1:1015497/239‑1 (MQ=255)
                                                                                                                                                                                  tgttgCCGATTATGGC‑GGGAATTTGCTCATCATTTGCCAGCCAGGGAAACCTGC‑TGGGGCCGATTCTATTTTTCGCCGTTTACCTGCTTATCTTTTTCCTGCGCGTCGGCTGGACCCACGTCGGTTATTCAGTCGGCGTGAAGGCGATCGATAAAGTGCGAGAGAACTCGCAGATGATTGCCCGTTCGGCAACCATCCTCGGGATCACGGTAATCGGCGGGCTGATCGCTTCGTAtgtg                                                     <  1:1146835/239‑1 (MQ=255)
                                                                                                                                                                                     tgCCGATTATGGC‑GGGAATTTGCTCATCATTTGCCAGCCAGGGAAACCTGC‑TGGGGCCGATTCTATTTTTCGCCGTTTACCTGCTTATCTTTTTCCTGCGCGTCGGCTGGACCCACGTCGGTTATTCAGTCGGCGTGAAGGCGATCGATAAAGTGCGAGAGAACTCGCAGATGATTGCCCGTTCGGCAACCATCCTCGGGATCACGGTAATCGGCGGGCTGATCGCTTCGTATGTGCat                                                  <  1:1023347/239‑1 (MQ=255)
                                                                                                                                                                                      gCCGATTATGGC‑GGGAATTTGCTCATCATTTGCCAGCCAGGGAAACCTGC‑TGGGGCCGATTCTATTTTTCGCCGTTTGCCTGCTTATCTTTTTCCTGCGCGTCGGCTGGACCCACGTCGGTTATTCAGTCGGCGTGAAGGCGATCGATAAAGTGCGAGAGAACTCGCAGATGATTGCCCGTTCGGCAACCATCCTCGGGATCACGGTAATCGGCGGGCTGATCGCTTCGTATGTGCata                                                 <  1:992362/239‑1 (MQ=255)
                                                                                                                                                                                        cGATTATGGC‑GGGAATTTGCTCATCATTTGCCAGCCAGGGAAACCTGC‑TGGGGCCGATTCTATTTTTCGCCGTTTACCTGCTTATCTTTTTCCTGCGCGTCGGCTGGACCCACGTCGGTTATTCAGTCGGCGTGAAGGCGATCGATAAAGTGCGAGAGAACTCGCAGATGATTGCCCGTTCGGCAACCATCCTCGGGATCACGGTAATCGGCGGGCTGATCGCTTCGTATGTGCATAtt                                               <  1:819912/239‑1 (MQ=255)
                                                                                                                                                                                           ttATGGC‑GGGAATTTGGTCATCATTTGCCAGCCAGGGAACCCTGC‑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:949865/239‑1 (MQ=255)
                                                                                                                                                                                                                tcatTTGCCAGCCAGGGAAACCTGC‑TGGGGCCGATTCTATTTTTCGCCGTTTACCTGCTTATCTTTTTCCTGCGCGTCGGCTGGACCCACGTCGGTTATTCAGTCGGCGTGAAGGCGATCGATAAAGTGCGAGAGAACTCGCAGATGATTGCCCGTTCGGCAACCATCCTCGGGATCACGGTAATCGGCGGGCTGATCGCTTCGTATGTGCATATTAACGTGGTGACATCGTTTGCCAt                        <  1:445312/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:374297/1‑238 (MQ=255)
                                                                                                                                                                                                                         agccagGGAAACCTGC‑TGGGGCCGATTCTATTTTTCGCCGTTTACCTGCTTATCTTTTTCCTGCGCGTCGGCTGGACCCACGTCGGTTATTCAGTCGGCGTGAAGGCGATCGATAAAGTGCGAGAGAACTCGCAGATGATTGCCCGTTCGGCAACCATCCTCGGGATCACGGTAATCGGCGGGCTGATCGCTTCGTATGTGCATATTAACGTGGTGACATCGTTTGCCATCGACAATa                >  1:333276/1‑238 (MQ=255)
                                                                                                                                                                                                                           ccagGGAAACCTGC‑TGGGGCCGATTCTATTTTTCGCCGTTTACCTGCTTATCTTTTTCCTGCGCGTCGGCTGGACCCACGTCGGTTATTCAGTCGGCGTGAAGGCGATCGATAAAGTGCGAGAGAACTCGCAGATGATTGCCCGTTCGGCAACCATCCTCGGGATCACGGTAATCGGCGGGCTGATCGCTTCGTATGTGCATATTAACGTGGTGACATCGTTTGCCATCGAc                    >  2:1147726/1‑232 (MQ=255)
                                                                                                                                                                                                                           ccagGGAAACCTGC‑TGGGGCCGATTCTATTTTTCGCCGTTTACCTGCTTATCTTTTTCCTGCGCGTCGGCTGGACCCACGTCGGTTATTCAGTCGGCGTGAAGGCGATCGATAAAGTGCGAGAGAACTCGCAGATGATTGCCCGTTCGGCAACCATCCTCGGGATCACGGTAATCGGCGGGCTGATCGCTTCGTATGTGCATATTAACGTGGTGACATCGTTTGCCATCGAc                    <  1:1147726/232‑1 (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)
                                                                                                                                                                                                                                               |                                                                                                                                                                                                                                      
GCGGCGATGAAAGATAACCTCGAATTTATTAATACCCACCCGAATCTGGTCGGATTCCTGATGGGGTTATTAATTTCGATGGAAGAAAAAGGAGAAAACCGCGACACCATTAAAGGCCTCAAAGTGGCACTGTTTGGCCCAATCGCCGGGATTGGCGATGCGATTTTCTGGTTTACTTTGTTGCCGATTATGGC‑GGGAATTTGCTCATCATTTGCCAGCCAGGGAAACCTGC‑TGGGGCCGATTCTATTTTTCGCCGTTTACCTGCTTATCTTTTTCCTGCGCGTCGGCTGGACCCACGTCGGTTATTCAGTCGGCGTGAAGGCGATCGATAAAGTGCGAGAGAACTCGCAGATGATTGCCCGTTCGGCAACCATCCTCGGGATCACGGTAATCGGCGGGCTGATCGCTTCGTATGTGCATATTAACGTGGTGACATCGTTTGCCATCGACAATACCCACAGCGTTGCG  >  NC_000913/3285638‑3286105

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