Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 653,049 T→G 21.0% V272G (GTG→GGG)  citA → sensory histidine kinase in two‑component regulatory system with CitB

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913653,0490TG21.0% 122.3 / 2.6 67V272G (GTG→GGG) citAsensory histidine kinase in two‑component regulatory system with CitB
Reads supporting (aligned to +/- strand):  ref base T (32/18);  new base G (10/6);  total (43/24)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TCTGATGTTGCTGTCGTGGTTCCTGGCCGCGCATATCCGTCGGCAGATGATGGGCATGGAGCCAAAGCAAATCGCACGCGTGGTCCGTCAGCAAGAGGCGCTGTTTAGTTCGGTTTATGAAGGGCTGATTGCGGTGGATCCGCATGGTTACATTACCGCCATCAATCGTAACGCAAGAAAGATGCTGGGGCTGAGTTCCCCCGGACGGCAATGGTTGGGTAAACCCATTGTTGAAGTGGTCAGGCCCGCCGATTTCTTTACCGAACAGATTGATGAAAAACGTCAGGATGTGGTGGCGAACTTTAACGGTCTGAGCGTTATTGCCAACCGGGAAGCTATTCGTTCAGGTGATGATTTGCTGGGGGCCATTATCAGCTTTCGTAGTAAAGACGAAATTTCCACCCTCAATGCGCAACTGACGCAAATAAAACAATACGTTGAGAGCCTTCGTACATTGCGACACGAGCATCTCAA  >  NC_000913/652813‑653286
                                                                                                                                                                                                                                             |                                                                                                                                                                                                                                                 
tCTGATGTTGCTGTCGTGGTTCCTGGCCGCGCATATCCGTCGGCAGATGATGGGCATGGAGCCAAAGCAAATCGCACGCGTGGTCCGTCAGCAAGAGGCGCTGTTTAGTTCGGTTTATGAAGGGCTGATTGCGGTGGATCCGCATGGTTACATTACCGCCATCAATCGTAACGCAAGAAAGATGCTGGGGCTGAGTTCCCCCGGACGGCAATGGTTGGGTAAACCCATTGTTGAAgggg                                                                                                                                                                                                                                                 >  1:712772/1‑239 (MQ=255)
tCTGATGTTGCTGTCGTGGTTCCTGGCCGCGCATATCCGTCGGCAGATGATGGGCATGGAGCCAAAGCAAATCGCACGCGTGGTCCGTCAGCAAGAGGCGCTGTTTAGTTCGGTTTATGAAGGGCTGATTGCGGTGGATCCGCATGGTTACATTACCGCCATCAATCGTAACGCAAGAAAGATGCCGGGGCTGAGTTCACCCAGACGGCAATGGTTGGGAAAACCAATTGTGGAAgggg                                                                                                                                                                                                                                                 >  2:655455/1‑239 (MQ=255)
       ttGCTGTCGTGGTTCCTGGCCGCGCATATCCGTCGGCAGATGATGGGCATGGAGCCAAAGCAAATCGCACGCGTGGTCCGTCAGCAAGAGGCGCTGTTTAGTTCGGTTTATGAAGGGCTGATTGCGGTGGATCCGCATGGTTACATTACCGCCATCAATCGTAACGCAAGAAAGATGCTGGGGCTGAGTTCCCCCGGACGGCAATGGTTGGGTAAACCCATTGTTGAAGTGGTCAGGc                                                                                                                                                                                                                                          >  1:417537/1‑238 (MQ=255)
       ttGCTGTCGTGGTTCCTGGCCGCGCATATCCGTCGGCAGATGATGGGCATGGAGCCAAAGCAAATCGCACACGTGGTCCGTCAGCAAGAGGCGCTGTTTAGTTCGGTTTATGAAGGGCGGATTGCGGTGGATCCGCATGGTTACATTACCGCCATCAATCGTAACGCAAGAAAGATGCTGGGGCCGAGTTCCCCCGGCCGGCAATGGTTGGGTAAACCCATTGTTGAAAGGGTAAGGc                                                                                                                                                                                                                                          >  2:218625/1‑238 (MQ=255)
                       tGGCCGCGCATATCCGTCGGCAGATGATGGGCATGGAGCCACAGCAAATCGCACGCGTGGTCCGTCAGCAAGAGGCGCTGTTTAGTTCGGTTTATGAAGGGCTGATTGCGGTGGATCCGCATGGTTACATTACCGCCATCAATCGTAACGCAAGAAAGATGCTGGGGCTGAGTTCCCCCGGACGGCCATGGTTGGGTAAACCCCTTGTTGAGGGGTGACGGCCCGCCGCTTTTTTTaac                                                                                                                                                                                                                         >  2:44401/1‑237 (MQ=255)
                       tGGCCGCGCATATCCGTCGGCAGATGATGGGCATGGAGCCAAAGCAAATCGCACGCGTGGTCCGTCAGCAAGAGGCGCTGTTTAGTTCGGTTTATGAAGGGCTGATTGCGGTGGATCCGCATGGTTACATTACCGCCATCAATCGTAACGCAAGAAAGATGCGGGGGGTGAGTTCCCCCGGACGACAATGGTTGGGAAAAACCATTGTTGAAGTGGTTGGGCCCGCCGCTTCTTTTAcg                                                                                                                                                                                                                         >  2:222583/1‑238 (MQ=255)
                           cgcgCATACCCGTCGGAAGTTGAGGGGCAGGGAGCCAAAGCAAATCGCACGCGTGGTCCGTCAGCAAGAGGCGCTGTTTAGTTCGGTTTATGAAGGGCTGATTGCGGTGGATCCGCATGGTTACATTACCGCCATCAATCGTAACGCAAGAAAGATGCTGGGGCTGAGTTCCCCCGGACGGCAATGGTTGGGTAAACCCATTGTTGAAGTGGTCAGGCCCGCCGATTTCTTTACCGaa                                                                                                                                                                                                                      <  1:284784/238‑1 (MQ=255)
                            gcgcATATCCGTCGGCAGATGATGGGCATGGAGCCAAAGCAAATCGCACGCGTGGTCCGTCAGCAAGAGGCGCTGTTTAGTTCGGTTTATGAAGGGCTGATTGCGGTGGATCCGCATGGTTACATTACCGCCATCAATCGTAACGCAAGAAAGATGCTGGGGCTGAGTTCCCCCGGACGGCAATGGTTGGGTAAACCCATTGTTGAGGGGGTCAGGCCCGCCGATTTCTTTAGCggaaa                                                                                                                                                                                                                    >  2:530208/1‑235 (MQ=255)
                            gcgcATATCCGTCGGCAGATGATGGGCATGGAGCCAAAGCAAATCGCACGCGTGGTCCGTCAGCAAGAGGCGCTGTTTAGTTCGGTTTATGAAGGGCTGATTGCGGTGGATCCGCATGGTTACATTACCGCCATCAATCGTAACGCAAGAAAGATGCTGGGGCTGAGTTCCCCCGGACGGCAATGGTTGGGTAAACCCATTGTTGAAGTGGTCAGGCCCGCCGATTTCTTTACCGAACa                                                                                                                                                                                                                    >  1:547047/1‑239 (MQ=255)
                             cgcATATCCGTCGGCAGATGATGGGCATGGAGCCAAAGCAAATCGCACGCGTGGTCCGTCAGCAAGAGGCGCTGCTTAGTTGGGTTTATGAAGGGCTGATTGCGGGGGATCCGCACGTTTACATTACAGCCATCAATCGTAACGCAAGAAAGATGCTGGGGCTGAGTTCCCCCGGCAGCAAATGGATGGGGAAACCCATTTGTAAAGAGGTCAGGCCCCCGGCTTTTTTTACCGaaaaa                                                                                                                                                                                                                   >  2:349099/1‑236 (MQ=255)
                                              ctgatgGGCATGGAGCCAAAGCAAATCGCACGCGTGGTCCGTCAGCAAGAGGCGCTGGTTAGTTCGGTTTATGAAGGGCTGATTGCGGTGGATCCGCATGGTTACATTACCGCCATCAATCGTAACGCCAGAAAGATGCTGCGGCTGGGTTCCCCCGGACGGCACTGGTAGGGTAAAACCATTGTTGAAGGGGTAAGGCCCGCCGATTTCTTTACCGAAAAAATAGATGaaaaagataa                                                                                                                                                                                                  >  2:1001341/2‑234 (MQ=255)
                                              atgatgGGCATGGAGCCAAAGCAAATCGCACGCGTGGTCCGTCAGCAAGAGGCGCTGTTTAGTTCGGTTTATGAAGGGCTGATTGCGGAGGATCCGCATGGTTACATTACCGCCATCAATCGTAACGCAAGAAAGATGCTGGGGCGGAGTTCCCCCGGACGGCAATGGTTGGGTAAACCCATTGTTGAAGTGGTCAGGCCCGCCGATTTCTTTACCGACCGGATTGATGAAAAACGTCa                                                                                                                                                                                                  >  1:616327/1‑239 (MQ=255)
                                                   gggCATGGAGCCAAAGCAAATCGCACGCGTGGTCCGTCAGCAAGAGGCGCTGTTTAGTTCGGTTTATGAAGGGCTGATTGCGGTGGATCCGCATGGTTACATTACCGCCATCAATCGTAACGCAAGAAAGATGCTGGGGCTGAGTTCCCCCGGACGGCAATGGTTGGGTAAACCCATTGTTGAAGTGGTCAGGCCCGCCGCTTTCTTTACCGAACAGATTGTTGAAAAACGCCAGGCtg                                                                                                                                                                                             >  1:955750/1‑239 (MQ=255)
                                                            gcaaaaGCAAACCGCACGCGTGGGCCGTCAGCAAGGGGCGCTGTTTAGTTCGGTTTATGAAGGGCTGATGGGGGTGGTCCCGCATGGTTACATTACCGCCATCAATCGTAACGCAAGAAAGATGCTGGGGCTGAGTTCCCCCGGACGGCAATGGTTGGGTAAACCCATTGTTGAAGTGGTCAGGCCCGCCGATTTCTTTACCGAACAGATTGATGAAAAACGTCAGGATGTGGTGGCg                                                                                                                                                                                     <  1:900595/235‑1 (MQ=255)
                                                                    aaaCCGCACGGGGGGCCCGTCAGCAAGGGGCGCTGTTTAGTTCGGTTTATGAAGGGCTGATTGCGGTGGATCGGCATGGTTACATTACCGCCATCAATCGTAACGCAAGAAAGATGCTGGGGCTGAGTTCCCCCGGACGGCAATGGTTGGGTAAACCCATTGTTGAAGTGGTCAGGCCCGCCGATTTCTTTACCGAACAGATTGATGAAAAACGTCAGGATATGGTGGCGAACTTTa                                                                                                                                                                              <  1:218625/237‑1 (MQ=255)
                                                                     aaTCGCACGCGTGGTCCGTCAGCAAGAGGCGCTGTTTAGTTCGGTTTATGAAGGGCTGATTGCGGTGGATCCGCATGGTTACATTACCGCCATCAATCGTAACGCAAGAAAGATGCTGGGGCTGAGTTCCCCCGGACGGCAATGGTTGGGTAAACCCATTGTTGAAGTGGTCAGGCCCGCCGATTTCTTTACCGAACAGATTGATGAAAAACGTCAGGATGTGGTGGCGAACTTTAACg                                                                                                                                                                           >  1:387058/1‑239 (MQ=255)
                                                                       tgccgcgggTGGTCGGTCAGGAAGAGGGGCTGTTTGGTTCGTTTTAAGGAGGGCGGGTTGGGGGAGATCCGCTGGGTTAGTTTACCGCCGTTAGTCGTAGCGCAAGAAAGTTGGTGGGGTTGGGTCCCCCGGGAGGGCAATGGTTGGGTAAACCCATTGTTGAAGGGGTCAGGCCCGCCGATTTCTTTACCGAACAGATTGATGAAAAACGTCAGGATGTGGTGGCGAACTTTAACGGt                                                                                                                                                                         <  2:176181/231‑1 (MQ=255)
                                                                       tcccACGCGTGGTCGGTCAGCAAGAGGCGCTGTTTAGTCGGGTTTATAAAGGTCTGATTGCGGTGGATCCGCATGGTTACATGTCCGCCATCAATCGTAACGCAAGAAAGATGCTGGGGCTGAGTTCCCCCGGACGGCAATGGTTGGGTAACCCCATTGTTGAAGTGGTCAGGCCCGCCGATTTCTTTACCGAACAGATTGATGAAAAACGGCAGGATGTGGTGGCGAACTTTAACGGt                                                                                                                                                                         <  1:503516/236‑1 (MQ=255)
                                                                            cgcgTGGTCCGTCAGCAAGAGGCGCTGTTTAGTTCGGTTTATGAAGGGCTGATTGCGGTGGATCCGCATGGTTACATTACCGCCATCAATCGTAACGCAAGAAAGATGCTGGGGCTGAGTTCCCCCGGACGGCAATGGTTGGGTAAACCCATTGTTGAAGTGGTCAGGCCCGCCGATTTCTTTACCGAACAGATTAGTGAAAAACGTCAGGATGTGGGGGGGAACTTTAACGGTCTGAg                                                                                                                                                                    >  1:514066/1‑239 (MQ=255)
                                                                             gcgTGGTCCGTCAGCAAGAGGCGCTGTTTAGTTCGGTTTATGAAGGGCTGATTGCGGTGGATCCGCATGGTTACATTACCGCCATCAATCGTAACGCAAGAAAGATGCTGGGGCTGAGTTCCCCCGGACGGCAATGGTTGGGTAAACCCATTGTTGAAGGGGTCAGGCCCGCCGATTTCTTTTCCGACAAAATTGATAAAAAACGTCAGGAGGGGGTGGCGAACTTAAAGGGTcgaagg                                                                                                                                                                   >  2:503516/1‑234 (MQ=255)
                                                                                       tCAGCAAGAGGCGCTGTTTAGTTCGGTTTATGAAGGGCTGATTGCGGTGGATCCGCATGGTTACATTACCGCCATCAATCGTAACGCAAGAAAGATGCTGGGGCTGAGTTCCCCCGGACGGCAATGGTTGGGTAAACCCATTGTTGAAGTGGTCAGGCCCGCCGATTTCTTTACCGACCAGATTGATGAAAAACGTCAGGTGTTGGTGGCGAACTTTAACGGACTAAGCGTTATTGCCa                                                                                                                                                        >  1:417681/1‑239 (MQ=255)
                                                                                             agagGCGCTGTTTAGTTCGGTTTATGAAGGGCTGATTGCGGTGGATCCGCATGGTTACATTACCGCCATCAATCGTAACGCAAGAAAGATGCTGGGGCTGAGTTCCCCCGGACGGCAATGGTTGGGTAAACCCATTGTTGAAGTGGTCAGGCCCGCCGATTTCTTTACCGAACAGATTGATGAAAAACGTCCGGATGTGGTGGCGAACTTTAACGGTCTGAGCGTTTTTGCCAACCggg                                                                                                                                                  >  1:824523/1‑239 (MQ=255)
                                                                                                    cTGTTTAGTTCGGTTTATGAAGGGCTGATTGCGGTGGATCCGCATGGTTACATTACCGCCATCAATCGTAACGCAAGAAAGATGCTGGGGCTGAGTTCCCCCGGACGGCAATGGTTGGGTAAACCCATTGTTGAAGTGGTCAGGCCCGCCGATTTCTTTACCGAACAGATTGATGAAAAACGTCAGGATGTGGTGGCGAACTTTAACGGTCTGAGCGTTATTGCCAACCGGGAAGCTAt                                                                                                                                           >  1:492816/1‑239 (MQ=255)
                                                                                                     tttttttttggttttATAGAGGGCTGATGGGGGGGGTCCCGGAGGGTACAATTCCCGCCAACAAGCGAAACGCAAGAAAGATGGGGGGGTTGGGTTCCCCGGGCCGGAATTGGTGGGGTAAACCCATTGTGGAAGTGGTCAGGCCCGCCGATTTCTTTACCGAACAGATTGATGAAAAACGTCAGGATGTGGTGGCGAACTTTAACGGTCTGAGCGTTATTGCCAACCGGGAAGCTAt                                                                                                                                           <  2:387058/226‑1 (MQ=255)
                                                                                                     tGTTTAGTTCGGTTTATGAAGGGCTGATTGCGGTGGATCCGCATGGTTACATTACCGCCATCAATCGTAACGCAAGAAAGATGCTGGGGCTGAGTTCCCCCGGACGGCAATGGTTGGGTAAACCCATTGTTGAAGTGGTCAGGCCCGCCGATTTCTTTACCGAACAGATTGGTGAAAAACGTCAGGATGTGGTGGCCAACTTTAACGGCCTGAGCGTTATTTCCAACCGGGAAgacaa                                                                                                                                           >  2:83733/1‑234 (MQ=255)
                                                                                                         tAGTTCGGTTTAAGAAGGTCTGTTTGGGGTGGACCCGAATGGTTCAATTACCGCCATCAATCGAAACGCAAGAAAGATGCTGGGGCTGAGTTCCCCCGGACGGCAATGGTTGGGTAAACCCATTGTTGAAGTGGTCAGGCCCGCCGATTTCTTTACCGAACAGATTGATGAAAAACGTCAGGATGTGGTGGCGAACTTTAACGGTCTGAGCGTTATTGCCAACCGGGAAGCTATTCGtt                                                                                                                                      <  1:530208/239‑1 (MQ=255)
                                                                                                              cGGTTTATGAAGGGTTGAGTGGGGTGGTCCCGCATGGTTACATTACCGCCATCAAGCGTAACGCAAGAAAGATGCTGGGGCTGAGTTCCCCCGGACGGCATTGTGTGGGTAACCCCATTGTTGAAGTGGTCAGGCCCGCCGATTTCTTTACCGAACAGATTGATGAAAAACGTCAGGATGTGGTGGCGAACTTTAACGGTCTGAGCGTTATTGCCAACCGGGAAGCTATTCGTTCAGGt                                                                                                                                 <  1:218396/239‑1 (MQ=255)
                                                                                                                 tttATGAAGGGCTGATTGCGGTGGATCCGCATGGTTACATTACCGCCATCAATCGTAACGCAAGAAAGATGCTGGGGCTGAGTTCCCCCGGACGGCAATGGTTGGGTAAACCCATTGTTGAAGTGGTCAGGCCCGCCGATTTCTTTACCGAACAGATTGATGAAAAACGTCAGGATGTGGTGGCGAACTTTAACGGTCTGAGCGTTATTGCCAACCGGGAAGCTATTCGTTCAGGtgat                                                                                                                              >  1:606379/1‑239 (MQ=255)
                                                                                                                 tttATGAAGGGCTGATTGCGGTGGATCCGCATGGTTACATTACCGCCATCAATCGTAACGCAAGAAAGATGCTGGGGCTGAATTCCCCCGGACGGCAATGGTGGGGAAAAACCATTGTTGAGGGGGTCAGGCCCGCCGATTTCTTTACCGAACAGATTGATGACAAACGTCGGGCAGTGGTGGCGAACTTTAAAGGAGTAAGCGGTATTGCCAACCGGGAGGCAACTAGTTTAgatgat                                                                                                                              >  2:29588/1‑239 (MQ=255)
                                                                                                                 tttATGAAGGGCTGATTGCGGTGGATCCGCATGGTTAAATTACCGCCATCAATCGTACCGCAAGAAAAATGCTGGGGCTTAGTTCCCCAGGACGGCAATGGTTGGGTAAACCCATAGTTGCAGGGGCCAGGCCCGCCAATTCTTTTACCGAACAGCTTGTTGCAAAACGTACGGATGTGGTGTCGAACTTTAAGTGTCTGCACATTATTACCAACCCGGACGCAATTCGGTCGGGtggt                                                                                                                              >  2:507655/1‑237 (MQ=255)
                                                                                                                      gAAGGGCTGATTGCGGTGGATCCGCATGGTTACATTACCGCCATCAATCGTAACGCAAGAAAGATGCTGGGGCTGAGTTCCCCCGGACGGCAATGGTTGGGTAAACCCATTGTTGAAGTGGTCAGGCCCGCCGATTTCTTTACCGAACAGATTGATGAAAAACGTCAGGATGTGGTGGCGAACTTTAACGGTCTGAGCGTTATTGCCAACCGGGAAGCTATTCGGTAAGGGGATGAttt                                                                                                                        >  2:11572/1‑239 (MQ=255)
                                                                                                                          ggCTGATTGCGGTGGATCCGCATGGTTACATTACCGCCATCAATCGTAATGCAAGAAAGATGCTGGGGCTGAGTTCCCCCGGACGGCAATGTTTGGGTAAACCCATTGTTGAAGGGGTCAGGCCCGCCGATTTCTTTACCGAACAGATTGATGAAAAACGTCCGGTGTTGGGGGCGAAATTTAACGGTATGAGCGTTATTGCCAACCGGGAAGCTATTCGTTCAGGTGGTGGttggtgg                                                                                                                    >  2:656093/1‑234 (MQ=255)
                                                                                                                           gCTGATTGCGGTGGATCCGCATGGTTACATTACCGCCATCAGTCGTAACGCAAGAAAGGTGCTGGGGCTGAGTTCCCCCGGACGGCAATGGTTGGGGAAACCCATTTTTGAAGGGGTCAGGCCCGCCAATTTCTTTACCGAACAGATGTATGAAAAACGTCAGGATGTGGTGGCGAACTTTAACGGCTGGGGCGTTATTGCCAACCGTGAAGCCTTTCGTGCAGGTGATGATTTgtggg                                                                                                                   >  2:843488/1‑235 (MQ=255)
                                                                                                                           gCTGATTGCGGTGGATCCGCATGGTTACATTACCGCCATCAATCGTAACGCAAGAAAGATGCTGGGGCTGAGTTCCCCCGGACGGCAATGGTTGGGTAAACCCATTGTTGAAGTGGTCAGGCCCGCCGATTTCTTTACCGAACAGATTGATGAAAAACGTCAGGATGTGGTGGCGAACTTTAACGGTCTGGGCGTTATTGCCAACCGGGAAGCTATTCGTTTAGGGGATGATTTGCTgg                                                                                                                   >  2:716386/1‑239 (MQ=255)
                                                                                                                               aTTGCGGTGGATCCGCATGGTTACATTACCGCCATCAATCGTAACGCAAGAAAGATGCTGGGGCTGAGTTCCCCCGGACGGCAATGGTTGGGTAAACCCATTGTTGAAGTGGTCAGGCCCGCCGATTTCTTTACCGAACAGATTGATGAAAAACGTCAGGATGTGGTGGCGAACTTTAACGGTCTGAGCGTTATTGCCAACCGGGAAGCTATTCGTTCAGGGGATGATTTGCTggggg                                                                                                                >  1:96983/1‑238 (MQ=255)
                                                                                                                                tggtgggggATCCGGTTTGTTTCATTTCCCCCATGCAAGTGTAACCTAGGAATGGAGCTGGGGCCGGTTTCCCCCGGACGGAAATGGTTGGGTAAACCGATTGTGGAGGTGGTCAGGCCCGCCGATTTCTTTACCGAACAGATTGATGAAAAACGTCAGGATGTGGTGGCGAACTTTAACGGTCTGAGCGTTATTGCCAACCGGGAAGCTATTCGTTCAGGTGATGATTTGCTGGGGGc                                                                                                               <  2:616327/232‑1 (MQ=255)
                                                                                                                                         aaccGCATGGTTACATTACCGCCATCAATCGTAACGCAAGAAAGATGCTGGGGCTGAGTTCCCCCGGACGGCAATGGTTGGGTAAACCCATTGTTGAAGTGGTCAGGCCAGCCGATTTCTTTACCGAACAGATTGATGAAAAACGTCAGGCGGTGGTGGCGCCCTTTAACGGTCTGAGCGTTATTGCCCACCGGGGAGCACTTCGTTACGGTGATGATTTTCTGGGGGCCATTAtaa                                                                                                       >  2:765198/3‑235 (MQ=255)
                                                                                                                                                  gTTACATTACCGCCATCAATCGTAACGCAAGAAAGATGCTGGGGCTGAGTTCCCCCGGACGGCAATGGTTGGGTAAACCCATTGTTGAAGTGGTCAGGCCCGCCGATTTCTTTACCGACCAGATTGATGAAAAACGTCAGGATGTGGTGGCGAACTTTAACGGGCTGAGCGTTATTGCAAACCGGGAAGCTATTAGGTCAGGTGAGGTTTTGTTGGGGGCCATTATCAGCtttgtaggt                                                                                           >  2:691453/1‑233 (MQ=255)
                                                                                                                                                      cATTACCGCCATCAATCGTAACGCAAGAAAGATGCTGGGGCTGAGTTCCCCCGGACGGCAATGGTTGGGTAAACCCATTGTTGAAGTGGTCAGGCCCGCCGATTTCTTTACCGAACAGATTGATGAAAAACGTCAGGATGTGGTGGCGAACTTTAACGGTCTGAGCGTTATTGCCAACCGGGAAGCTATTCGTTCAGGTGATGATTTGCTGGGGGCCATTATCAGCTTTCGTAGTAAAg                                                                                       >  1:832949/1‑239 (MQ=255)
                                                                                                                                                         tACCGCCATCAATCGTAACGCAAGAAAGATGCTGGGGCTGAGTTCCCCCGGACGGCAATGGTTGGGTAAACCCATTGTTGAAGTGGTCAGGCCCGCCGATTTCTTTACCGAACAGATTGATGAAAAACGTCAGGATGTGGTGGCGAACTTTAACGGTCTGAGCGTTATTGCCAACCGGGAAGCTATTCGTTCAGGTGCTGATTTGCTGGGGGCCATTATCAGCTTTCGTAGTAAAGACg                                                                                    >  1:484024/1‑239 (MQ=255)
                                                                                                                                                            cGCCATCACCCGTAACGCAAGACAGATGCTGGGGCTGAGTTCCCCCGGACGGCAATGTGTGGGTAAACCCATTGTTGCAGGGGTCCGGCCCCCCGATTTCTTTACCGAACAGATTGTTGAAAAACACCGGGATGGGTGGTCGACTTTTAAGGACCGGGGGTTTATGCACAACCGGGAAGATATTCATCCGGGGATGGATTTGCGGGGGGCACTGTCAGCAGTTGGTAGAAAAGAGGaaa                                                                                 >  2:1029196/1‑239 (MQ=255)
                                                                                                                                                            cGCCATCAATCGTAACGCAAGAAAGATGCTGGGGCTGAGTTCCCCCGGACGGCAATGGTTGGGTAAACCCATTGTTGAAGTGGTCAGGCCCGCCGATTTCTTTACCGAACAGATTGATGAAAAACGTCAGGATGTGGTGGCGAACTTTAACGGTCTGAGCGTTATTGCCAACCGGGAAGCTATTCGTTCAGGTGATGATTTGCTGGGGGCCATTATCAGCTTTCGTAGTAAAGACGaaa                                                                                 <  1:44401/239‑1 (MQ=255)
                                                                                                                                                             gCCATCAATCGTAACGCAAGAAAGATGCTGGGGCTGAGTTCCCCCGGACGGCAATGGTTGGGTAAACCCATTGTTGAAGTGGTCAGGCCCGCCGATTTCTTTACCGAACAGATTGATGAAAAACGTCAGGATGTGGTGGCGAACTTTAACGGTCTGAGCGTTATTGCCAACCGGGGAGCTATTCGTTCAGGGGATGATTTGCTGGGGGCCATTATCAGCTTTCGTGGTAAAGCCGAAAt                                                                                >  1:549944/1‑239 (MQ=255)
                                                                                                                                                                aTCAATCGTAACGCAAGAAAGATGCTGGGGCTGAGTTCCCCCGGACGGCAATGGTTGGGTAAACCAATTGTTGAAGTGGTCAGGCCCGCCGATTTCTTTACCGAACAGATTGATGAAAAACGTCAGGATGTGGTGGCGAACTTTAACGGTCTGAGCGTTATTGCCAACCGGGAAGCTATTCGTTCAGGGGATGTTTTGCTGTGGGCCATTATCGTGTTTGCTATTAAAAACGAAAttt                                                                              >  2:227084/1‑238 (MQ=255)
                                                                                                                                                                 tCAATCGTAACGCAAGAAAGATGCTGGGGCTGAGTTCCCCCGGACGGCAATGGTTGGGTAAACCCATTGTTGAAGTGGTCAGGCCCGCCGATTTCTTTACCGAACAGATTGATGAAAAACGTCAGGATGTGGTGGCGAACTTTAACGGTCTGAGCGTTATTGCCAACCGGGAAGCTATTCGTTCAGGTGATGATTTGCTGGGGGCCATTATCAGCTTTCGTAGTAAAGCCGAAATTTc                                                                             >  1:835742/1‑238 (MQ=255)
                                                                                                                                                                       cGTAACGCAAGAAAGATGCTGGGGCTGAGTTCCCCCGGACGGCAATGGTTGGGTAAACCCATTGTTGAAGTGGTCAGGCCCGCCGATTTCTTTACCGAACAGATTGATGAAAAACGTCAGGCTGTGGTGGCGAACTTTAACGGTCTGAGCGTTATTGCCAACCGGGAAGCTATTCGTTCAGGTGATGATTTGCTGGGGGCCATTATCAGCTTTCGTAGTAAAGACGAAATTTCCACCCt                                                                       >  2:595992/1‑239 (MQ=255)
                                                                                                                                                                            cgaaagaaagaTGTTGGGGCTGAGTTCCCCCGGACGGCAATGGTTGGGTAACCCCATTGTTGAAGGGGTCAGGCCCGCCGATTTCTTTACCGAACAGATTGATGAAAAACGTCAGGATGTGGTGGCGAACTTTAACGGTCTGAGCGTTATTGCCAACCGGGAAGCTATTCGTTCAGGTGATGATTTGCTGGGGGCCATTATCAGCTTTCGTAGTAAAGACGAAATTTCCACCCTCAATg                                                                  <  1:929341/236‑1 (MQ=255)
                                                                                                                                                                                 gaaagaTGCTGGGGCTGAGTTCCCCCGGACGGCAATGGTTGGGTAAACCCATTGTTGAAGTGGTCAGGCCCGCCGATTTCTTTACCGAACAGATTGATGAAAAACGTCAGGATGTGGTGGCGAACTTTAACGGTCTGAGCGTTATTGCCAACCGGGAAGCTATTCGTTCAGGGGATGATTTGCTGGGGGCCATTATCAGCTTTCGTAGTAAAGACGAAATTTCCACCCTCAATGCGCa                                                              >  1:344203/1‑238 (MQ=255)
                                                                                                                                                                                       tgttGGGGCTTGTTTCCCCCGGGCGGAATTGGTTGGGTAAACCCGTTGTAGAAGGGGTAAGGCCCGCCGATTTCTTTACCGAACAGATTGATGAAAAACGTCAGGATGTGGTGGCGAACTTTAACGGTCTGAGCGTTATTGCCAACCGGGAAGCTATTCGTTCAGGTGATGATTTGCTGGGGGCCATTATCAGCTTTCGTAGTAAAGACGAAATTTCCACCCTCAATGCGCAACTGAc                                                        <  2:96983/235‑1 (MQ=255)
                                                                                                                                                                                              gtttgttccccccGGACGGCAATGGTTGGGTAAACCCATTGTTGAAGGGGTCAGGCCCGCCGATTTCTTTACCGAACAGATTGATGAAAAACGTCAGGATGTGGTGGCGAACTTTAACGGTCTGAGCGTTATTGCCAACCGGGAAGCTATTCGTTCAGGTGATGATTTGCTGGGGGCCATTATCAGCTTTCGTAGTAAAGACGAAATTTCCACCCTCAATGCGCAACTGACGCAAATa                                                 <  1:655455/230‑1 (MQ=255)
                                                                                                                                                                                                            gtcgtcattttttGGGTAACCCCTTTTTTGAGGTGGTCAGGCCCGTCGATTTCTTTACCGAACAGATTGATGAAAAACGTCAGGATGGGGTTGCTAACTTTAACGGTCTGAGCGTTATTGCCAACCGGGAAGCTATTCGTTCAGGTGATGATTTGCTGGGGGCCATTATCAGCTTTCGTAGTAAAGACGAAATTTCCACCCTCAATGCGCAACTGACGCAAATAAAACAATACGTTgag                                  <  2:514066/228‑1 (MQ=255)
                                                                                                                                                                                                            gACGGCAATGGTTGGGTAAACCCATTGTTGAAGTGGTCAGGCCCGCCGATTTCTTTACCGAACAGATTGATGAAAAACGTCAGGATGTGGTGGCGAACTTTAACGGTCTGAGCGTTATTGCCAACCGGGAAGCTATTCGTTCAGGTGATGATTTGCTGGGGGCCATTAt                                                                                                         <  1:554037/169‑1 (MQ=255)
                                                                                                                                                                                                            gACGGCAATGGTTGGGTAAACCCATTGTTGAAGTGGTCAGGCCCGCCGATTTCTTTACCGAACAGATTGATGAAAAACGTCAGGATGTGGTGGCGAACTTTAACGGTCTGAGCGTTATTGCCAACCGGGAAGCTATTCGTTCAGGTGATGATTTGATGGGGGCCATTAt                                                                                                         >  2:554037/1‑169 (MQ=255)
                                                                                                                                                                                                             aCGGCATTGGTTGGGTAACCCCATTTTTAATGTGGCCAGGCCCGCCGATTTCTTTACCGAACAGATTGATGAAAAACGGCAGGATGTGGTGGCGAACTTTAACGGTCTGAGCGTTATTGCCAACCGGGAAGCTATTCGTTCAGGTGATGATTTGCTGGGGGCCATTATCAGCTTTCGTAGTAAAGACGTAATTTCCACCCTCAATGCGCAACTGACGCAAATAAAACAATACGTTgag                                  <  1:545622/238‑1 (MQ=255)
                                                                                                                                                                                                             aCGGCAATGGTTGGGTAAACCCATTTTTGAAGGGGTCAGGCCCGCCGATTTCTTTACCGAACAGATTGATGAAAAACGTCAGGATGTGGTGGCGAACTTTAACGGTCTGAGCGTTATTGCCAACCGGGAAGCTATTCGTTCAGGTGATGATTTGCTGGGGGCCATTATCAGCTTTCGTAGTAAAGACGAAATTTCCACCCTCAATGCGCAACTGACGCAAATAAAACAATACGTTgag                                  <  1:691453/238‑1 (MQ=255)
                                                                                                                                                                                                               tgtattGGTTGAGTAAACCTTTTTTTGAATGGGCACGCCCCTCCGATTTCTTTACCGAACAGATTGATGAAAAACGTCAGGATGTGGTGGCGAACTTTAACGGTCTGAGCGTTATTGCCAACCGGGAAGCTATTCGTTCAGGTGATGATTTGCTGGGGGCCATTATCAGCTTTCGTAGTAAAGACGAAATTTCCACCCTCAATGCGCAACTGACGCAAATAAAACAATACGTTGagag                                <  2:835742/233‑1 (MQ=255)
                                                                                                                                                                                                               gtaaagggTTGGGTAAACCCATTGTTGAATGGGCAGGGCCCGCCGATTTCTTTACCGAACAGATTGATGAAAAACGTCAGGATGTGGTGGCGAACTTTAACGGTCTGAGCGTTATTGCCAACCGGGAAGCTATTCGTTCAGGTGATGATTTGCTGGGGGCCATTATCAGCTTTCGTAGTAAAGACGAAATTTCCACCCTCAATGCGCAACTGACGCAAATAAAACAATACGTTGagag                                <  1:1001341/232‑1 (MQ=255)
                                                                                                                                                                                                                                cccATTGTTGAAGTGGTCAGGCCCGCCGATTTCTTTACCGAACAGATTGATGAAAAACGTCAGGATGTGGTGGCGAACTTTAACGGTCTGAGCGTTATTGcc                                                                                                                                                         <  1:528258/102‑1 (MQ=255)
                                                                                                                                                                                                                                cccATTGTTGAAGTGGTCAGGCCCGCCGATTTCTTTACCGAACAGATTGATGAAAAACGTCAGGATGTGGTGGCGAACTTTAACGGTCTGAGCGTTATTGcc                                                                                                                                                         >  2:528258/1‑102 (MQ=255)
                                                                                                                                                                                                                                 ccATTGTTGAAGTGGTCAGGCCCGCCGATTTCTTTACCGAACAGATTGATGAAAAACGTCAGGATGTGGTGGCGAACTTTAACGGTCTGAGCGTTATTGCCAACCGGGAAGCTATTCGTTCAGGTGATGATTTGCTGGGGGCCATTATCAGCTTTCGTAGTAAAGACGAAATTTCCACCCTCAATGCGCAACTGACGCAAATAAAACAATACGTTGAGAGCCTTCGTACATTGCGacac             >  1:108207/1‑239 (MQ=255)
                                                                                                                                                                                                                                  ggttgttgAATTGGCCAGCCCCGCCGATTTCTTTACCGAACAGATTGATGAAAAACGTCAGGATGTGGTGGCGAACTTTAACGGTCTGAGCGTTATTGCCAACCGGGAAGCTATTCGTTCAGGTGATGATTTGCTGGGGGCCATTATCAGCTTTCGTAGTAAAGACGAAATTTCCGCCCTCAATGCGCAACTGACGCAAATAAAACAATACGTTGAGAGCCTTCGTACATTGCGACACg            <  1:519164/237‑1 (MQ=255)
                                                                                                                                                                                                                                    ttttagAAGGGGTCAGGCCCGCCGATTTCTTTACCGAACAGATTGATGAAAAACGTCAGGATGTGGTGGCGAACTTTAACGGTCTGAGCGTTATTGCCAACCGGGAAGCTATTCGTTCAGGTGATGATTTGCTGGGGGCCATTATCAGCTTTCGTAGTAAAGACGAAATTTCCACCCTCAATGCGCAACTGACGCAAATAAAACAATACGTTGAGAGCCTTCGTACATTGCGACACGAg          <  1:716386/234‑1 (MQ=255)
                                                                                                                                                                                                                                     tgttgAAGTGTTCAGGCCCGCCGATTTCTTTACCGAACAGATTGATGAAAAACGTCAGGATGTGGTGGCGAACTTTAACGGTCTGCGCGGTATTGCCAACCGGGAAGCTATTCGTTCAGGGGCTGATTTGCTAGGGTCCATTATCCGCTTTCGTAGTAAAAGCGAAATTTCCACCATCAATGCGTAACTTACGCACATAAAACAAAAAGATAAGACCCTTCGTACTTTGCGccaataca         >  2:945858/1‑234 (MQ=255)
                                                                                                                                                                                                                                     tgttgAAGTGGTCAGGCCCGCCGATTTCTTTACCGAACAGATTGATGAAAAACGTCAGGATGTGGTGGCGAACTTTAACGGTCTGAGCGTTATTGCCAACCGGGAAGCTATTCGTTCAGGTGATGATTTGCTGGGGGCCATTATCAGCTTTCGTAGTAAAGACGAAATTTCCACCCTCAATGCGCAACTGACGCAAATAAAACAATACGTTGCGAGCCTTCTTCCATTGCGACACGAGc         >  1:767929/1‑239 (MQ=255)
                                                                                                                                                                                                                                     tgttgAAGTGGTCAGGCCCGCCGATTTCTTTACCGAACAGATTGATGAAAAACGTCAGGATGTGGTGGCGAACTTTAACGGTCTGAGCGTTATTGCCAACCGGGAAGCTATTCGTTCAGGTGATGATTTGCTGGGGGCCATTATCAGCTTTCGTAGTAAAGACGAAATTTCCACCCTCAATGCGCAACTGACGCAAATAAAACAATACGTTGAGAGCCTTCGTACATTGCGACCCGAGc         >  1:293664/1‑239 (MQ=255)
                                                                                                                                                                                                                                     tgttgAAGTGGTCAGGCCCGCCGATTTCTTTACCGAACAGATTGATGAAAAACGTCAGGATGTGGTGGCGAACTTTAACGGTCTGAGCGTTATTGCCAACCGGGAAGCTATTCGTTCAGGTGATGATTTGCTGGGGGCCATTATCAGCTTTCGTAGTAAAGACGAAATTTCCACCCTCAATGCGCAACTGACGCAAATAAAACAATACGTTGAGACCCTTCGTCCATTTGGACACGAg          >  1:398278/1‑238 (MQ=255)
                                                                                                                                                                                                                                       ttgAAGTGGTCAGGCCCGCCGATTTCTTTACCGAACAGATTGATGAAAAACGTCAGGATGTGGTGGCGAACTTTAACGGTCTGAGCGTTATTGCCAACCGGGAAGCTATTCGTTCAGGTGATGATTTGCTGGGGGCCATTATCAGCTTTCGTAGTAAAGACGAAATTTCCACCCTCAATGCGCAACTGACGCAAATAAAACAATACGTTGAGAGCCTTCGTACATTGCGACACGAGCAt       <  1:11572/239‑1 (MQ=255)
                                                                                                                                                                                                                                       tgaaaGTGTTCAGGCCCGCCGATTTCTTTACCGAACAGATTGATGAAAAACGTCAGGATGTGGTGGCGAACTTTAACGGTCTGAGCGTTATTGCCAACCGGGAAGCTATTCGTTCAGGTGATGATTTGCTGGGGGCCATTATCAGCTTTCGTAGTAAAGACGAAATTTCCACCCTCAATGCGCAACTGACGCAAATAAAACAATACGTTGAGAGCCTTCGTACATTGCGACACGAGCAt       <  1:227084/236‑1 (MQ=255)
                                                                                                                                                                                                                                            gTGGTCAGGCCCGCCGATTTCTTTACCGAACAGATTGATGAAAAACGTCAGGATGTGGTGGCGAACTTTAACGGTCTGAGCGTTATTGCCAACCGGGAAGCTATTCGTTCAGGTGATGATTTGCTGGGGGCCATTATCAGCTTTCGTAGTAAAGCCGAAATTTCCACCCTCAATGCGCAACTGACGCAAATAAACAAAAACGTTGAGAGCCTTCGTACATAGCGACCCGAGccacacaa  >  2:267598/1‑232 (MQ=255)
                                                                                                                                                                                                                                             |                                                                                                                                                                                                                                                 
TCTGATGTTGCTGTCGTGGTTCCTGGCCGCGCATATCCGTCGGCAGATGATGGGCATGGAGCCAAAGCAAATCGCACGCGTGGTCCGTCAGCAAGAGGCGCTGTTTAGTTCGGTTTATGAAGGGCTGATTGCGGTGGATCCGCATGGTTACATTACCGCCATCAATCGTAACGCAAGAAAGATGCTGGGGCTGAGTTCCCCCGGACGGCAATGGTTGGGTAAACCCATTGTTGAAGTGGTCAGGCCCGCCGATTTCTTTACCGAACAGATTGATGAAAAACGTCAGGATGTGGTGGCGAACTTTAACGGTCTGAGCGTTATTGCCAACCGGGAAGCTATTCGTTCAGGTGATGATTTGCTGGGGGCCATTATCAGCTTTCGTAGTAAAGACGAAATTTCCACCCTCAATGCGCAACTGACGCAAATAAAACAATACGTTGAGAGCCTTCGTACATTGCGACACGAGCATCTCAA  >  NC_000913/652813‑653286

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