Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,637,409 A→G K238K (AAA→AAG)  creC → sensory histidine kinase CreC

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,637,4090AG90.0% 127.3 / 5.1 50K238K (AAA→AAG) creCsensory histidine kinase CreC
Reads supporting (aligned to +/- strand):  ref base A (2/3);  new base G (17/28);  total (19/31)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.04e-01

CGGCTTATTGGCGTTTTGAGCGTAGGCAAACCGAACGCGGCGATGGCTCCGGTCATTAAGCGTAGCGAGCGGCGAATTTTATGGGCCAGCGCCATTTTGT‑TGGGGATTGCACTGGTGATTGGCGCAGGCATGGTTTGGTGGATCAACCGCTCTATTGCCCGGCTCACTCGCTATGCTGATTCCGTCACTGACAATAAGCCCGTTCCTCTCCCCGATCTCGGTAGTAGCGAGTTGCGTAAACTCGCGCAGGCGCTGGAAAGTATGCGCGTGAAGCTGGAAGGGAAAAACTATATTGAGCAGTATGTTTACGCATTAACTCATGAGCTAAAAAGCCCACTGGCGGCGATTCGTGGAGCGGCGGAAATTTTACGC‑GAAGGTCCGCCGCCGGAAGTGGTGGCTCGTTTTACTGACAACATTCTGACGCAAAATGCGCGTATGCAGGCATTGGTAGAAACGTTACTACGCCAGGCAAGACTGGAGAAT  >  NC_000913/4637170‑4637652
                                                                                                                                                                                                                                                |                                                                                                                                                                                                                                                    
cgGCTTATTGTCGTTTTGAGCGTAGGCAAACCGAACGCGGCGATGGCTCCGGTCATTAAGCGTAGCGAGCGGCGAATTTTATGGGCCAGCGCCATTTTGC‑TGGGGATTGCACTGGTGATTGGCGCAGGCATGGTTTGGTGGATCAACCGCTCCATTGCCCGGCTCACTCGCTATGCCGATTCCGTCACTGACAATAAGCCCGTTCCTCTTCCTGAACTCGTCAGTAGCGAGTTGCGTAAGCTTGCGCAGgc                                                                                                                                                                                                                                           <  1:124615/251‑1 (MQ=255)
cgGCTTATTGGCGTTTTGAGCGTAGGCAAACCGAACGCGGCGATGGCTCCGGTCATTAAGCGTAGCGAGCGGCGAATTTTATGGGCCAGCGCCATTTTGC‑TGGGGATTGCACTGGTGATTGGCGCAGGCATGGTTTGGTGGATCAACCGCTCCATTGCCCGGCTCACTCGCTATGCCGATTCCGTCACTGACAATAAGCCCGTTCCTCTTCCTGAACTCGGCAGTAGCGAGTTGCGTAAGCTTGCGCAGgc                                                                                                                                                                                                                                           >  1:267975/1‑251 (MQ=255)
           cGTTTTGAGCGTAGGCAAACCGAACGCGGCGATGGCTCCGGTCATTAAGCGTAGCGAGCGGCGAATTTTATGGGCCAGCGCCATTTTGC‑TGGGGATTGCACTGGGGATTGGCGCAGGCATGGTTTGGTGGATCAACCGCTCCCTTGCCCGGCTCCCTCGCTATGCCGATTCCGTCACTGACAATAAGCCCGTTCCTCTTCCTGAACTCGGCAGTAGCGAGTTGCGTAAGCTTGCGCAGGCGCTggcaggg                                                                                                                                                                                                                                 >  2:188987/1‑245 (MQ=255)
              tttGATCGTAGGCAAACCGAACGCGGCGATGGCTCCGGTCATTAAGCGTAGCGAGCTGCGAATTTTATGGGCCAGCGCCATTTTGC‑TGGGGATTGCACTGGTGATTGGCGCAGGCATGGTTTGGTGGATCAACCGCTCCATTGCCCGGCTCACTCGCTATGCCGATTCCGTCACTGACAATAAGCCCGTTCCTCTTCCTGAACTCGGCAGTAGCGAGTTGCGTAAGCTTGCGCAGGCGCTGGAAAGTATg                                                                                                                                                                                                                              <  1:501976/250‑1 (MQ=255)
                   gCGTAGGCAAACCGAACGCGGCGATGGCTCCGGTCATTAAGCGTAGCGAGCGGCGAATTTTATGGGCCAGCGCCATTTTGC‑TGGGGATTGCACTGGTGATTGGCGCAGGCATGGTTTGGTGGATCAACCGCTCCATTGCCCGGCTCACTCGCTATGCCGATTCCGTCACTGACAATAAGCCCGTTCCTCTTCCTGAACTCGGCAGTAGCGAGTTGCGTAAGCTTGCGCAGGCGCTGGAAAGTATGCGCGTg                                                                                                                                                                                                                        <  1:358633/251‑1 (MQ=255)
                     gTAGGCAAACCGAACGCGGCGATGGCTCCGGTCATTAAGCGTAGCGAGCGGCGAATTTTATGGGCCAGCGCCATTTTGC‑TGGGGATTGCACTGGTGATTGGCGCAGGCATGGTTTGGTGGATCAACCGCTCCATTGCCCGGCTCACTCGCTATGCCGATTCCGTCACTGACAATAAGCCCGTTCCTCTTCCTGAACTCGGCAGTAGCGAGTTGCGTAAGCTTGCGCAGGCGCTGGAAAGTATGCGCGTGa                                                                                                                                                                                                                       >  1:52485/1‑250 (MQ=255)
                                 aaCGCGGCGATGGCCCCGGTCATTAAGCGTAGCGAGCGGCGAATTTTATGGGCCAGCGCCATTTTGC‑TGGGGATTGCACTGGTGATTGGCGCAGGCATGGTTTGGTGGATCAACCGCTCCATTGCCCGGCTCACTCGCTATGCCGATTCCGTCACTGACAATAAGCCCGTTCCTCTTCCTGAACTCGGCAGTAGCGAGTTGCGTAAGCTTGCGCAGGCGCTGGAAAGTATGCGCGTGAAGCTGGAAGGGaa                                                                                                                                                                                                          <  2:52485/251‑1 (MQ=255)
                                         gATGGCTGCGGTCATTAAGCGTAGCGAGCGGCGAATTTTATGGGCCAGCGCCATTTTGT‑TGGGGATTGCACTGGTGATTGGCGCAGGCATGGTTTGGTGGATCAACCGCTCTATTGCCCGGCTCACTCGCTATGCTGATTCCGTCACTGACAATAAGCCCGTTCCTCTCCCAGATCTCGGTAGTAGCGAGTTGCGTAAACTCGCGCAGGCGCTGGAAAGTATGCGCGTGAAGCTGGAAGGGAAAACCtat                                                                                                                                                                                                   >  2:37518/1‑250 (MQ=255)
                                         gATGGCTCCGGTCATTAAGCGTAGCGAGCGGCGAATTTTATGGGCCAGCGCCATTTTGC‑TGGGGATTGCACTGGTGATTGGCGCAGGCATGGTTTGGTGGATCAACCGCTCCATTGCCCGGCTCACTCGCTATGCCGATTCCGTCACTGACCATAAGCCCGTTCCTCTTCCTGCACTCGGCAGTAGCGAGTTGCGTAAGCTTGCGCAGGCGCTGGAAAGTATGCGCGTGAAGCTGGAAGGGAAAAACtat                                                                                                                                                                                                   >  2:312706/1‑250 (MQ=255)
                                          aTGGCTCCGGTCATTAAGCGTAGCGAGCGGCGAATTTTATGGGCCAGCGCCATTTTGC‑TGGGGATTGCACTGGTGATTGGCGCAGGCATGGTTTGGTGGATCAACCGCTCCATTGCCCGGCTCACTCGCTATGCCGATTCCGTCACTGACAATAAGCCCGTTCCTCTTCCTGAACTCGGCAGTAGCGAGTTGCGTAAGCTTGCGCAGGCGCTGGAAAGTATGCGCGTGAAGCTGGAAGGGAAAAACTatat                                                                                                                                                                                                 >  1:502227/1‑251 (MQ=255)
                                               gccGGTCATTAAGCGTAGCGAGCGGCGAATTTGTTGGGCCAGCGCCATTT‑GTCTGGGGAGTGCACTGGTGATGGGCGCAGGCATGGTTTGGTGGATCAACCGTACCATTGCCCGGCTCACTCGCTATGCCGATTCCGTCACTGACAATAAGCCCGTTCCTCTTCCTGAACTCGGCAGTAGCGAGTTGCGTAAGCTTGCGCAGGCGCTGGAAAGTATGCGCGTGAAGCTGGAAGGGAAAAACTATATTGAGc                                                                                                                                                                                            <  2:384630/250‑1 (MQ=255)
                                                 cGGTCATTAAGCGTAGCGAGCGGCGAATTTTATGGGCCAGCGCCATTTTGC‑TGGGGATTGCACTGGTGATTGGCGCAGGCATGGTTTGGTGGATCAACCGCTCCATTGCCCGGCTCACTCGCTATGCCGATTCCGTCACTGACAATAAGCCCGTTCCTCGTCCTGAACTCGGCAGTAGCGAGTTGCGTAAGCTTGCGCAGGCGCTGGAAAGTATGCGCGTGAAGCTGGAAGGGAAAAACTATATTGAGCAg                                                                                                                                                                                          <  2:385606/251‑1 (MQ=255)
                                                   gTCATTAAGCGTAGCGAGCGGCGAATTTTATGGGCCAGCGCCATTTTGC‑TGGGGATTGCACTGGTGATTGGCGCAGGCATGGTTTGGTGGATCAACCGCTCCATTGCCCGGCTCACTCGCTATGCCGATTCCGTCACTGACAATAAGCCCGTTCCTCTTCCTGAACTCGGCAGTAGCGAGTTGCGTAAGCTTGCGCAGGCGCTGGAAAGTATGCGCGTGAAGCTGGAAGGGAAAAACTATATTGAGCAGt                                                                                                                                                                                         >  1:245368/1‑250 (MQ=255)
                                                          aGCGTAGCGAGCGGCGAATTTTATGGGCCAGCGCCATTTTGC‑TGGGGATTGCACTGGTGATTGGCGCAGGCATGGTTTGGTGGATCATCCGCTCCATTGCCCGGCTCACTCGCTATGCCGATTCCGTCACTGCCAATAAGCCCGTTCCTCTTCCTGAACTCGGCAGTAGCGAGTTGCGTAAGCTTGCGCAGGCGCTGGAAAGTATGAGCGGGAAGCTGGAAGGGAGAAACTATATTGAGCAGTATGTTTat                                                                                                                                                                                 >  2:400921/1‑250 (MQ=255)
                                                                gcgagcgGCGAATTTTATGGGCCAGCGCCATTTTGT‑TGGGGATTGCACTGGTGATTGGCGCAGGCATGGTTTGGTGGATCAACCGCTCTATTGCCCGGCTCACTCGCTATGCTGATTCCGTCACTGACAATAAGCCCGTTCCTCTCCCCGATCTCGGTAGTAGCGAGTTGCGTAAACTCGCGCAGGCGCTGGAAAGTATGCGCGTGAAGCTGGAAGGGAAAAACTATATTGAGCAGTATGTTTACGCATTa                                                                                                                                                                           <  1:37518/251‑1 (MQ=255)
                                                                 cgagcgGCGAATTTTATGGGCCAGCGCCATTTTGC‑TGGGGATTGCACTGGTGATTGGCGCAGGCATGGTTTGGTGGATCAACCGCTCCATTGCCCGGCTCACTCGCTATGCCGATTCCGTCACTGACAATAAGCCCGTTCCTCTTCCTGAACTCGGCAGTAGCGAGTTGCGTAAGCTTGCGCAGGCGCTGGAAAGTATGCGCGTGAAGCTGGAAGGGAAAAACTATATTGAGCAGTATGTTTATGCGTTaa                                                                                                                                                                          <  2:268062/251‑1 (MQ=255)
                                                                           aTTTTATGGGCCAGCGCCATTTTGC‑TGGGGATTGCACTGGTGATTGGCGCAGGCATGGTTTGGTGGATCAACCGCTCCATTGCCCGGCTCACTCGCTATGCCGATTCCGTCACTGACAATAAGCCCGTTCCTCTTCCTGAACTCGGCAGTAGCGAGTTGCGTAAGCTTGCGCAGGCGCTGGAAAGTATGCGCGTGAAGCTGGAAGGGAAAAACTATATTGAGCAGTATGTTTATGCGTTAACCCATGAGCt                                                                                                                                                                >  1:360273/1‑251 (MQ=255)
                                                                                aTGGGCCAGCGCCATTTTGC‑TGGGGATTGCACTGGTGATTGGCGCAGGCATGGTTTGGTGGATCAACCGCTCCATTGCCCGGCTCACTCGCTATGCCGATTCCGTCACTGACAATAAGCCCGTTCCTCTTCCTGAACTCGGCAGTAGCGAGTTGCGTAAGCTTGCGCAGGCGCTGGAAAGTATGCGCGTGAAGCTGGAAGGGAAAAACTATATTGAGCAGTATGTTTATGCGTTAACCCATGAGCTaaaaa                                                                                                                                                           <  2:374366/251‑1 (MQ=255)
                                                                                  gggCCAGCGCCATTTTGC‑TGGGGATTGCACTGGTGATTGGCGCAGGCATGGTTTGGTGGATCAACCGCTCCATTGCCCGGCTCACTCGCTATGCCGATTCCGTCACTGACAATAAGCCCGTTCCTCTTCCTGAACTCGGCAGTAGCGAGTTGCGTAAGCTTGCGCAGGCGCTGGAAAGTATGCGCGTGAAGCTGGAAGGGAAAAACTATATTGAGCAGTATGTTTATGCGTTAACCCATGAGCTAAAAAg                                                                                                                                                          <  1:213748/250‑1 (MQ=255)
                                                                                   ggCCAGCGCCATTTTGC‑TGGGGATTGCACTGGTGATTGGCGCAGGCATGGTTTGGTGGATCAACCGCTCCATTGCCCGGCTCACTCGCTATGCCGATTCCGTCACTGACAATAAGCCCGTTCCTCTTCCTGAACTCGGCAGTAGCGAGTTGCGTAAGCTTGCGCAGGCGCTGGAAAGTATGCGCGTGAAGCTGGAAGGGAAAAACTATATTGAGCAGTATGTTTATGCGTTAACCCATGAGCTAAAAAGcc                                                                                                                                                        <  1:82642/251‑1 (MQ=255)
                                                                                            cATTTTGC‑TGGGGATTGCACTGGTGATTGGCGCAGGCATGGTTTGGTGGATCAACCGCTCCATTGCCCGGCTCACTCGCTATGCCGATTCCGTCACTGACAATAAGCCCGTTCCTCTTCCTGAACTCGGCAGTAGCGAGTTGCGTAAGCTTGCGCAGGCGCTGGAAAGTATGCGCGTGACGCTGGAAGGGAAAAACTATATTGAGCAGTATGTTTATGCGTTAACCCATGAGCTAAAAAGCCCACTggcgg                                                                                                                                               >  1:326377/1‑251 (MQ=255)
                                                                                                   gctgGGGATTTCACTGGTGATTGGCGCAGGCATGGTGTGGTGGATCACCCGCTCCATTGCCCGGCTCACTCGCTATGCCGATTCCGTCACTGACAATAAGCCCGTTCCTCTTCCTGAACTCGGCAGTAGCGAGTTGCGTAAGCGTGCGCAGGCGCTGGAAAGTATGCGCGTGAAGCTGGAAGGGAAAAACTATATTGAGCAGTATGTTTATGCGTTAACCCATGAGCTAAAAAGCCCTCTGGCGGCGAtt                                                                                                                                          <  2:124021/248‑1 (MQ=255)
                                                                                                   gctgGGGATTGCACTGGTGATTGGCGCAGGCATGGTTTGGTGGATCAACCGCTCCATTGCCCGGCTCACTCGCTATGCCGATTCCGTCACTGACAATAAGCCCGTTCCTCTTCCTGAACTCGGCAGTAGCGAGTTGCGTAAGCTTGCGCAGGCGCTGGAAAGTATGCGCGTGAAGCTGGAAGGGAAAAACTATATTGAGCAGTATGTTTATGCGTTAACCCATGAGCTAAAAAGACCACTGGCGGCGGTTc                                                                                                                                         >  2:378929/3‑251 (MQ=255)
                                                                                                              cACTGGTGATTGGCGCAGGCATGGTTTGGTGGATCAACCGCTCCATTGCCCGGCTCACTCGCTATGCCGATTCCGTCACTGACAATAAGCCCGTTCCTCTTCCTGAACTCGGCAGTAGCGAGTTGCGTAAGCTTGCGCAGGCGCTGGAAAGTATGCGCGTGAAGCTGGAAGGGAAAAACTATATTGAGCAGTATGTTTATGCGTTAACCCATGAGCTAAAAAGCCCACTGGCGGCGATTCGTGGCgcggcg                                                                                                                              >  1:250549/1‑251 (MQ=255)
                                                                                                                     gATTGGCGCAGGCATGGTTTGGTGGATCAACCGCTCCATTGCCCGGCTCACTCGCTATGCCGATTCCGTCACTGACAATAAGCCCGTTCCTCTTCCTGAACTCGGCAGTAGCGAGTTGCGTAAGCTTGCGCAGGCGCTGGAAAGTATGCGCGTGAAGCTGGaa                                                                                                                                                                                                               <  1:402559/163‑1 (MQ=255)
                                                                                                                     gATTGGCGCAGGCATGGTTTGGTGGATCAACCGCTCCATTGCCCGGCTCACTCGCTATGCCGATTCCGTCACTGACAATAAGCCCGTTCCTCTTCCTGAACTCGGCAGTAGCGAGTTGCGTAAGCTTGCGCAGGCGCTGGAAAGTATGCGCGTGAAGCTGGaa                                                                                                                                                                                                               >  2:402559/1‑163 (MQ=255)
                                                                                                                              aGGCATGGTTTGGTGGATCAACCGCTCCATTGCCCGGCTCACTCGCTATGCCGATTCCGTCACTGACAATAAGCCCGTTCCTCTTCCTGAACTCGGCAGTAGCGAGTTGCGTAAGCTTGCGCAGGCGCTGGAAAGTATGCGCGTGAAGCTGGAAGGGAAAAACTATATTGAGCAGTATGTTTATGCGTTAACCCATGAGCTAAAAAGCCCACTGGCGGCGATTCGTGGCGCGGCGGAAATTTTACGC‑Ga                                                                                                               >  1:238339/1‑249 (MQ=255)
                                                                                                                                 gaTGGTTTGGTGGATCAACCGCTCCATTGCCCGGCTCACTCGCTATGCCGATTCCGTCACTGACAATAAGCCCGTTCCTCTTCCTGAACTCGGCAGTAGCGAGTTGCGTAAGCTTGCGCAGGCGCTTGAAAGGATGCGCGTGAAGCTGGAGGGGAAAAACTATATTGCGCGGTATGTTTATGCGTTAACCCATGAGCTAAAAAGCCCCCTGCAGGCGGTTCGTGGCGCGGGGGTACTTTTAcgcacaggtc                                                                                                           >  2:90991/2‑244 (MQ=255)
                                                                                                                                        tggtggATCAACCGCTCCATTGCCCGGCTCACTCTCTATGCCGATTCCGTCACTGACAATAAGCCCGTTCCTCTTCCTGAACTCGGCAGTAGCGAGTTGCGTAAGCTTGCGCAGGCGCTGGAAAGTATGCGCGTGAAGCTGGAAGGGACAAAATATATTGAGCAGTATGTTTATGCGTTAACCCATGAGCTAAAAAGCCCACTTGCGGCGATTCGTGGCGCGGCGGAAATTTTTCGCAGAAggcccccccc                                                                                                    >  2:336735/1‑243 (MQ=255)
                                                                                                                                                   ccGCTCCATTGCCCGGCTCACTCGCTATGCCTATTCCGTCACTGACAATAAGCCCGTTCCTCTTCCTGAACTCGGCAGTAGCGAGTTGCGTAAGCTTGCGCAGGCGCTGGAAAGTATGCGCGTGAAGCTGGAAGGGAAAAACTATATTGAGCAGTATGTTTATGCGTTAACCCATGAGCTAAAAAGCCCACTGGCGGCGATTCTTGGCGCGTCGGAAATTTTACGC‑GAAGGTCCGCTTCCGGAAGtggtgg                                                                                        <  1:423398/251‑1 (MQ=255)
                                                                                                                                                               ccGGCTCACTCGCTATGCCGATTCCGTCACTGACAATAAGCCCGTTCCTCTTCCTGAACTCGGCAGTAGCGAGTTGCGTAAGCTTGCGCAGGCGCTGGAAAGTATGCGCGTGAAGCTGGAAGGGAAAAACTATATTGAGCAGTATGTTTATGCGTTAACCCATGAGCTAAAAAGCCCACTGGCGGCGATTCGTGGCGCGGCGGAAATTTTACGC‑GAAGGTCCGCCACCGGAAGTGGTGGCTCGTTTTAccg                                                                            <  1:254724/251‑3 (MQ=255)
                                                                                                                                                                  gCTCACTCGCTATGCTGATTCCGTCACTGACAATAAGCCCGTTCCTCTCCCCGATCTCGGTAGTAGCGAGTTGCGTAAACTCACGCAGGCGCTGGAAAGTATGCGCGTGAAGCTGGAAGGGAAAAACTATATTGAGCAGTATGTTTACGCATTAACTCATGAGCTAAAAAGCCCACTGGCGGCGATTCGTGGAGCGGCGGAAATTTTACGC‑GAAGGTCCGCCGCCGGAAGTGGTGGCTCGTTTTACTGaca                                                                         <  1:5563/251‑1 (MQ=255)
                                                                                                                                                                     cACTCGCTATGCCGATTCCGTCACTGACAATAAGCCCGTTCCTCTTCCTGAACTCGGCAGTAGCGAGTTGCGTAAGCTTGCGCAGGCGCTGGAAAGTATGCGCGTGAAGCTGGAAGGGAAAAACTATATTGAGCAGTATGTTTATGCGTTAACCCATGAGCTAAAAAGCCCACTGGCGGCGATTCGTGGCGCGGCGGAAATTTTACGC‑GAAGGTCCGCCGCCGGAAGTGGTGGCTCGTTTTACCGacaaca                                                                      <  1:132725/251‑1 (MQ=255)
                                                                                                                                                                         cGCTATGCCGATTCCGGCACAGACAATAAGCCCGGTCCTCGTCCTGATCTCGGCAGTAGCGAGTTGCGTAAGCTTGCGCAGGCGCTGGAAAGTATGCGCGTGAAGCTGGAAGGGAAAAACTATATTGAGCAGGATGTTTATGCGTTAACCCATGAGCTAAAAAGCCCACTGGCGGCGATTCGTGGCGCGGCGGAAATTTTACGC‑GAAGGTCCGCCACCGGAAGTGGTGGCTCGTTTTACCGACAACAtt                                                                    <  2:135564/249‑1 (MQ=255)
                                                                                                                                                                              tgccgATTCCGTCACTGACAATAAGCCCGTTCCTCTTCCTGAACTCGGCAGTAGCGAGTTGCGTAAGCTTGCGCAGGCGCTGGAAAGTATGCGCGTGAAGCTGGAAGGGAAAAACTATATTGAGCAGTATGTTTATGCGTTAACCCATGAGCTAAAAAGCCCACTGGCGGCGATTCGTGGCGCGGCGGAAATTTTACGC‑GAAGGTCCGCCACCGGAAgtggt                                                                                          <  2:296264/218‑1 (MQ=255)
                                                                                                                                                                              tgccgATTCCGTCACTGACAATAAGCCCGTTCCTCTTCCTGAACTCGGCAGTAGCGAGTTGCGTAAGCTTGCGCAGGCGCTGGAAAGTATGCGCGTGAAGCTGGAAGGGAAAAACTATATTGAGCAGTATGTTTATGCGTTAACCCATGAGCTAAAAAGCCCACTGGCGGCGATTCGTGGCGCGGCGGAAATTTTACGC‑GAAGGTCCGCCACCGGAAggggt                                                                                          >  1:296264/5‑222 (MQ=255)
                                                                                                                                                                                   aTTCCGTCACTGACAATAAGCCCGTTCCTCTTCCTGAACGCGGCAGTAGCGAGTTGCGTAAGCTTGCGCAGGCGCTGGAAAGTATGCGCGTGAAGCTGGAAGGGAAAAACTATATTGAGCAGTATGTTTATGCGTTAACCCATGAGCTAAAAAGCCCACTGGCGGCGATTCGTGGCGCGGCGGAAATTTTACGC‑GAAGGTCCGCCACCGGAAGTGGTGGCTCGTTTTACCGACAACATTCTGACGCAAAAt                                                        <  1:335246/251‑1 (MQ=255)
                                                                                                                                                                                      ccGTCACTGACAATAAGCCCGTTCCTCTTCCTGAACTCGGCAGTAGCGAGTTGCGTAAGCTTGCGCAGGCGCTGGAAAGTATGCGCGTGAAGCTGGAAGGGAAAAACTATATTGAGCAGTATGTTTATGCGTTAACCCATGAGCTAAAAAGCCCACTGGCGGCGATTCGTGGCGCGGCGGAAATTTTACGC‑GAAGGTCCGCCACCGGAAGTGGTGGCTCGTTTTACCGACAACATTCTGACGCAAAATgcg                                                     <  1:176683/251‑1 (MQ=255)
                                                                                                                                                                                      ccGTCACTGACAATAAGCCCGTTCCTCTTCCTGAACTCGGCAGTAGCGAGTTGCGTAAGCTTGCGCAGGCGCTGGAAAGTATGCGCGTGAAGCTGGAAGGGAAAAACTATATTGAGCAGTATGTTTATGCGTTAACCCATGAGCTAAAAAGCCCACTGGCGGCGATTCGTGGCGCGGCGGAAATTTTACGC‑GAAGGTCCGCCACCGGAAGTGGTGGCTCGTTTTACCGACAACATTCTGACGCAAAATgc                                                      <  1:56542/250‑1 (MQ=255)
                                                                                                                                                                                              gACAATAATCCCGTTCCTCTTCCTGAACTCGGCAGTAGCGAGTTGCGTAAGCTTGCGCAGGCGCTGGAAAGTATGCGCGTGAAGCTGGAAGGGAAAAACTATATTGAGCAGTATGTTTATGCGTTAACCCATGAGCTAAAAAGCCCACTGGCGGCGATTCGTGGCGCGGCGGAAATTTTACGC‑GAAGGTCCGCCACCGGAAGTGGTGGCTCGTTTTACCGACAACATTCTGACGCAAAATGCGCGTATg                                               <  2:502227/249‑1 (MQ=255)
                                                                                                                                                                                              gACAATAAGCCCGTTCCTCTTCCTGAACTCGGCAGTAGCGAGTTGCGTAAGCTTGCGCAGGCGCTGGAAAGTATGCGCGTGAAGCTGGAAGGGAAAAACTATATTGAGCAGTATGTTTATGCGTTAACCCATGAGCTAAAAAGCCCACTGGCGGCGATTCGTGGCGCGGCGGAAATTTTACGC‑GAAGGTCCGCCACCGGAAGTGGTGGCTCGTTTTACCGACAACATTCTGACGCAAAATGCGCGTATGCa                                             <  1:312706/251‑1 (MQ=255)
                                                                                                                                                                                                      gCCCGTTCCTCTTCCTGAACTCGGCAGTAGCGAGTTGCGTAAGCTTGCGCAGGCGCTGGAAAGTATGCGCGTGAAGCTGGAAGGGAAAAACTATATTGAGCAGTATGTTTATGCGTTAACCCATGAGCTAAAAAGCCCACTGGCGGCGATTCGTGGCGCGGCGGAAATTTTACGC‑GAAGGTCCGCCAACGGAAGTGGTGGCGCGGTTTACCGACCACAATCTGACGCCAAATGCGCGTATGCAGGCATTgg                                     >  1:364678/1‑251 (MQ=255)
                                                                                                                                                                                                       cccGTTCCTCTTCCTGAACTCGGCAGTAGCGAGTTGCGTAAGCTTGCGCAGGCGCTGGAAAGTATGCGCGTGAAGCTGGAAGGGAAAAACTATATTGAGCAGTATGTTTATGCGTTAACCCATGAGCTAAAAAGCCCACTGGCGGCGATTCGTGGCGCGGCGGAAATTTTACGC‑GAAGGTCCGCCACCGGAAGTGGTGGCTCGTTTTACCGACAACATTCTGACGCAAAATGCGCGTATGCAGGCATTGGt                                    <  2:245368/251‑1 (MQ=255)
                                                                                                                                                                                                                cttcctgaactcGGCAGTAGCGAGTTGCGTAAGCTTGCGCAGGCGCTGGAAAGTATGCGCGTGAAGCTGGAAGGGAAAAACTATATTGAGCAGTATGTTTATGCGTTAACCCATGAGCTAAAAAGCCCACTGGCGGCGATTCGTGGCGCGGCGGAAATTTTACGC‑GAAGGTCCGCCACCGGAAGTGGTGGCTCGTTTTACCGACAACATTCTGACGCAAAATGCGCGTATTCAGGCATTGGTAGAAACg                             <  1:1841/240‑1 (MQ=255)
                                                                                                                                                                                                                   cctgaactcGGCAGTACCGAGTTGCGTAAGCTTGCGCAGGCGCTGGAAAGTATGCGCGTGAATCTTGAAGGGAAAAACTATATTGAGCAGTATGGTTAGGCTTTAACCCATGAGCTAAAAAGCCCACTGGCGGCGATTCGTGGCGCGGCGGAAATTTTACGC‑GAAGGTCCGCCACCGGAAGTGGTGGCTCGTTTTACCGACAACATTCGGACGCAAAATGCGCGTATGCAGGCCTTGGTAGAAACGTtact                        <  2:447644/245‑1 (MQ=255)
                                                                                                                                                                                                                        tctcGGTAGTAGCGAGTTGCGTAAACTCGCGCAGGCGCTGGAAAGTATGCGCGTGAAGCTGGAAGGGAAAAACTATATTGAGCAGTATGTTTACGCATTAACTCATGAGCTAAAAAGCCCACTGGCGGCGATTCGTGGAGCGGCGGAAATTTTACGC‑GAAGGTccgccg                                                                                                     >  1:36069/1‑169 (MQ=255)
                                                                                                                                                                                                                        tctcGGTAGTAGCGAGTTGCGTAAACTCGCGCAGGCGCTGGAAAGTATGCGCGTGAAGCTGGAAGGGAAAAACTATATTGAGCAGTATGTTTACGCATTAACTCATGAGCTAAAAAGCCCACTGGCGGCGATTCGTGGAGCGGCGGAAATTTTACGC‑GAAGGTccgccg                                                                                                     <  2:36069/169‑1 (MQ=255)
                                                                                                                                                                                                                                  agCGAGTTGCGTAAGCTTGCGCAGGCGCTGGAAAGTATGCGCGTGAAGCTGGAAGGGAAAAACTATATTGAGCAGTATGTTTATGCGTTAACCCATGAGCTAAAAAGCCCACTGGCGGCGATTCGTGGCGCGGCGGAAATTTTACGC‑GAAGGTCCGCCACCGGAAGTGGTGGCTCGTTTTACCGACAACATTCTGACGCAAAATGCGCGTATGCAGGCATTGGTAGAAACGTTACTACGCCAGGCAAGAc          <  1:205046/250‑1 (MQ=255)
                                                                                                                                                                                                                                         tGCGTAAGCTTGCGCAGGCGCTGGAAAGTATGCGCGTGAAGCTGGAAGGGAAAAACTATATTGAGCAGTATGTTTATGCGTTAACCCATGAGCTAAAAAGCCCACTGGCGGCGATTCGTGGCGCGGCGGAAATTTTACGC‑GAAGGTCCGCCACCGGAAGTGGTGGCTCGTTTTACCGACAACATTCTGACGCAAAATGCGCGTATGCAGGCATTGGTAGAAACGTTACTACGCCAGGCAAGACTGGAGaa   <  1:447304/250‑1 (MQ=255)
                                                                                                                                                                                                                                          gCGTAAGCTTGCGCAGGCGCTGGAAAGTATGCGCGTGAAGCTGGAAGGGAAAAACTATATTGAGCAGTATGTTTATGCGTTAACCCATGAGCTAAAAAGCCCACTGGCGGCGATTCGTGGCGCGGCGGAAATTTTACGC‑GAAGGTCCGCCACCGGAAGTGGTGGCTCGTTTTACCGACAACATTCTGACGCAAAATGCGCGTATGCAGGCATTGGTAGAAACGTTACTACGCCAGGCAAGACTGGAGAAt  <  2:360273/250‑1 (MQ=255)
                                                                                                                                                                                                                                                |                                                                                                                                                                                                                                                    
CGGCTTATTGGCGTTTTGAGCGTAGGCAAACCGAACGCGGCGATGGCTCCGGTCATTAAGCGTAGCGAGCGGCGAATTTTATGGGCCAGCGCCATTTTGT‑TGGGGATTGCACTGGTGATTGGCGCAGGCATGGTTTGGTGGATCAACCGCTCTATTGCCCGGCTCACTCGCTATGCTGATTCCGTCACTGACAATAAGCCCGTTCCTCTCCCCGATCTCGGTAGTAGCGAGTTGCGTAAACTCGCGCAGGCGCTGGAAAGTATGCGCGTGAAGCTGGAAGGGAAAAACTATATTGAGCAGTATGTTTACGCATTAACTCATGAGCTAAAAAGCCCACTGGCGGCGATTCGTGGAGCGGCGGAAATTTTACGC‑GAAGGTCCGCCGCCGGAAGTGGTGGCTCGTTTTACTGACAACATTCTGACGCAAAATGCGCGTATGCAGGCATTGGTAGAAACGTTACTACGCCAGGCAAGACTGGAGAAT  >  NC_000913/4637170‑4637652

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