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

TCTGATGTTGCTGTCGTGGTTCCTGGCCGCGCATATCCGTCGGCAGATGATGGGCATGGAGCCAAAGCAAATCGCACGCGTGGTCCGTCAGCAAGAGGCGCTGTTTAGTTCGGTTTATGAAGGGCTGATTGCGGTGGATCCGCATGGTTACATTACCGCCAT‑CAATCGTAACGCAAGAAAGATGCTGGGGCTGAGTTCCCCCGGACGGCAATGGTTGGGTAAACCCATTGTTGAAGTGGTCA‑GGCCCGCCGATTTCTTTACCGAACAGATTGATGAAAAACGTCAGGATGTGGTGGCGAACTTTAACGGTCTGAGCGTTATTGC‑CAACCGGGAAGCTATTCGTTCAGGTGAT‑GATTTGCTGGGGGCCATTATCAGC‑TTTCGTAGTAAAGACGAAATTTCCACCCTCAATGCGCAACTGACGCAAATAAAACAATACGTTGAGAGCCTTCGTACATTGCGACACGAGCATCTCAA  >  NC_000913/652813‑653286
                                                                                                                                                                                                                                             |                                                                                                                                                                                                                                                 
tCTGATGTTGCTGTCGTGGTTCCTGGCCGCGCATATCCGTCGGCAGATGATGGGCATGGAGCCAAAGCAAATCGCACGCGTGGTCCGTCAGCAAGAGGCGCTGTTTAGTTCGGTTTATGAAGGGCTGATTGCGGTGGATCCGCATGGTTACATTACCGCCAT‑CAATCGTAACGCAAGAAAGATGCTGGGGCTGAGTTCCCCCGGACGGCAATGGTTGGGTAAACCCATTGTTGAAgggg                                                                                                                                                                                                                                                 >  1:712772/1‑239 (MQ=255)
tCTGATGTTGCTGTCGTGGTTCCTGGCCGCGCATATCCGTCGGCAGATGATGGGCATGGAGCCAAAGCAAATCGCACGCGTGGTCCGTCAGCAAGAGGCGCTGTTTAGTTCGGTTTATGAAGGGCTGATTGCGGTGGATCCGCATGGTTACATTACCGCCAT‑CAATCGTAACGCAAGAAAGATGCCGGGGCTGAGTTCACCCAGACGGCAATGGTTGGGAAAACCAATTGTGGAAgggg                                                                                                                                                                                                                                                 >  2:655455/1‑239 (MQ=255)
       ttGCTGTCGTGGTTCCTGGCCGCGCATATCCGTCGGCAGATGATGGGCATGGAGCCAAAGCAAATCGCACGCGTGGTCCGTCAGCAAGAGGCGCTGTTTAGTTCGGTTTATGAAGGGCTGATTGCGGTGGATCCGCATGGTTACATTACCGCCAT‑CAATCGTAACGCAAGAAAGATGCTGGGGCTGAGTTCCCCCGGACGGCAATGGTTGGGTAAACCCATTGTTGAAGTGGTCA‑GGc                                                                                                                                                                                                                                          >  1:417537/1‑238 (MQ=255)
       ttGCTGTCGTGGTTCCTGGCCGCGCATATCCGTCGGCAGATGATGGGCATGGAGCCAAAGCAAATCGCACACGTGGTCCGTCAGCAAGAGGCGCTGTTTAGTTCGGTTTATGAAGGGCGGATTGCGGTGGATCCGCATGGTTACATTACCGCCAT‑CAATCGTAACGCAAGAAAGATGCTGGGGCCGAGTTCCCCCGGCCGGCAATGGTTGGGTAAACCCATTGTTGAAAGGGTAA‑GGc                                                                                                                                                                                                                                          >  2:218625/1‑238 (MQ=255)
                       tGGCCGCGCATATCCGTCGGCAGATGATGGGCATGGAGCCACAGCAAATCGCACGCGTGGTCCGTCAGCAAGAGGCGCTGTTTAGTTCGGTTTATGAAGGGCTGATTGCGGTGGATCCGCATGGTTACATTACCGCCAT‑CAATCGTAACGCAAGAAAGATGCTGGGGCTGAGTTCCCCCGGACGGCCATGGTTGGGTAAACCCCTTGTTG‑AGGGGTGACGGCCCGCCGCTTTTTTTaac                                                                                                                                                                                                                         >  2:44401/1‑237 (MQ=255)
                       tGGCCGCGCATATCCGTCGGCAGATGATGGGCATGGAGCCAAAGCAAATCGCACGCGTGGTCCGTCAGCAAGAGGCGCTGTTTAGTTCGGTTTATGAAGGGCTGATTGCGGTGGATCCGCATGGTTACATTACCGCCAT‑CAATCGTAACGCAAGAAAGATGCGGGGGGTGAGTTCCCCCGGACGACAATGGTTGGGAAAAACCATTGTTGAAGTGGTTG‑GGCCCGCCGCTTCTTTTAcg                                                                                                                                                                                                                         >  2:222583/1‑238 (MQ=255)
                           cgcgCATACCCGTCGGAAGTTGAGGGGCAGGGAGCCAAAGCAAATCGCACGCGTGGTCCGTCAGCAAGAGGCGCTGTTTAGTTCGGTTTATGAAGGGCTGATTGCGGTGGATCCGCATGGTTACATTACCGCCAT‑CAATCGTAACGCAAGAAAGATGCTGGGGCTGAGTTCCCCCGGACGGCAATGGTTGGGTAAACCCATTGTTGAAGTGGTCA‑GGCCCGCCGATTTCTTTACCGaa                                                                                                                                                                                                                      <  1:284784/238‑1 (MQ=255)
                            gcgcATATCCGTCGGCAGATGATGGGCATGGAGCCAAAGCAAATCGCACGCGTGGTCCGTCAGCAAGAGGCGCTGTTTAGTTCGGTTTATGAAGGGCTGATTGCGGTGGATCCGCATGGTTACATTACCGCCAT‑CAATCGTAACGCAAGAAAGATGCTGGGGCTGAGTTCCCCCGGACGGCAATGGTTGGGTAAACCCATTGTTGAGGGGGTCA‑GGCCCGCCGATTTCTTTAGCggaaa                                                                                                                                                                                                                    >  2:530208/1‑235 (MQ=255)
                            gcgcATATCCGTCGGCAGATGATGGGCATGGAGCCAAAGCAAATCGCACGCGTGGTCCGTCAGCAAGAGGCGCTGTTTAGTTCGGTTTATGAAGGGCTGATTGCGGTGGATCCGCATGGTTACATTACCGCCAT‑CAATCGTAACGCAAGAAAGATGCTGGGGCTGAGTTCCCCCGGACGGCAATGGTTGGGTAAACCCATTGTTGAAGTGGTCA‑GGCCCGCCGATTTCTTTACCGAACa                                                                                                                                                                                                                    >  1:547047/1‑239 (MQ=255)
                             cgcATATCCGTCGGCAGATGATGGGCATGGAGCCAAAGCAAATCGCACGCGTGGTCCGTCAGCAAGAGGCGCTGCTTAGTTGGGTTTATGAAGGGCTGATTGCGGGGGATCCGCACGTTTACATTACAGCCAT‑CAATCGTAACGCAAGAAAGATGCTGGGGCTGAGTTCCCCCGGCAGCAAATGGATGGGGAAACCCATTTGTAAAGAGGTCA‑GGCCCCCGGCTTTTTTTACCGaaaaa                                                                                                                                                                                                                   >  2:349099/1‑236 (MQ=255)
                                              ctgatgGGCATGGAGCCAAAGCAAATCGCACGCGTGGTCCGTCAGCAAGAGGCGCTGGTTAGTTCGGTTTATGAAGGGCTGATTGCGGTGGATCCGCATGGTTACATTACCGCCAT‑CAATCGTAACGCCAGAAAGATGCTGCGGCTGGGTTCCCCCGGACGGCACTGGTAGGGTAAAACCATTGTTGAAGGGGTAA‑GGCCCGCCGATTTCTTTACCGAAAAAATAGATGaaaaagataa                                                                                                                                                                                                  >  2:1001341/2‑234 (MQ=255)
                                              atgatgGGCATGGAGCCAAAGCAAATCGCACGCGTGGTCCGTCAGCAAGAGGCGCTGTTTAGTTCGGTTTATGAAGGGCTGATTGCGGAGGATCCGCATGGTTACATTACCGCCAT‑CAATCGTAACGCAAGAAAGATGCTGGGGCGGAGTTCCCCCGGACGGCAATGGTTGGGTAAACCCATTGTTGAAGTGGTCA‑GGCCCGCCGATTTCTTTACCGACCGGATTGATGAAAAACGTCa                                                                                                                                                                                                  >  1:616327/1‑239 (MQ=255)
                                                   gggCATGGAGCCAAAGCAAATCGCACGCGTGGTCCGTCAGCAAGAGGCGCTGTTTAGTTCGGTTTATGAAGGGCTGATTGCGGTGGATCCGCATGGTTACATTACCGCCAT‑CAATCGTAACGCAAGAAAGATGCTGGGGCTGAGTTCCCCCGGACGGCAATGGTTGGGTAAACCCATTGTTGAAGTGGTCA‑GGCCCGCCGCTTTCTTTACCGAACAGATTGTTGAAAAACGCCAGGCtg                                                                                                                                                                                             >  1:955750/1‑239 (MQ=255)
                                                            gcaaaaGCAAACCGCACGCGTGGGCCGTCAGCAAGGGGCGCTGTTTAGTTCGGTTTATGAAGGGCTGATGGGGGTGGTCCCGCATGGTTACATTACCGCCAT‑CAATCGTAACGCAAGAAAGATGCTGGGGCTGAGTTCCCCCGGACGGCAATGGTTGGGTAAACCCATTGTTGAAGTGGTCA‑GGCCCGCCGATTTCTTTACCGAACAGATTGATGAAAAACGTCAGGATGTGGTGGCg                                                                                                                                                                                     <  1:900595/235‑1 (MQ=255)
                                                                    aaaCCGCACGGGGGGCCCGTCAGCAAGGGGCGCTGTTTAGTTCGGTTTATGAAGGGCTGATTGCGGTGGATCGGCATGGTTACATTACCGCCAT‑CAATCGTAACGCAAGAAAGATGCTGGGGCTGAGTTCCCCCGGACGGCAATGGTTGGGTAAACCCATTGTTGAAGTGGTCA‑GGCCCGCCGATTTCTTTACCGAACAGATTGATGAAAAACGTCAGGATATGGTGGCGAACTTTa                                                                                                                                                                              <  1:218625/237‑1 (MQ=255)
                                                                     aaTCGCACGCGTGGTCCGTCAGCAAGAGGCGCTGTTTAGTTCGGTTTATGAAGGGCTGATTGCGGTGGATCCGCATGGTTACATTACCGCCAT‑CAATCGTAACGCAAGAAAGATGCTGGGGCTGAGTTCCCCCGGACGGCAATGGTTGGGTAAACCCATTGTTGAAGTGGTCA‑GGCCCGCCGATTTCTTTACCGAACAGATTGATGAAAAACGTCAGGATGTGGTGGCGAACTTTAACg                                                                                                                                                                           >  1:387058/1‑239 (MQ=255)
                                                                       tgccgcgggTGGTCGGTCAGGAAGAGGGGCTGTTTGGTTCGTTTTAAGGAGGGCGGGTTGGGGGAGATCCGCTGGGTTAGTTTACCGCCGT‑TAGTCGTAGCGCAAGAAAGTTGGTGGGGTTGGGTCCCCCGGGAGGGCAATGGTTGGGTAAACCCATTGTTGAAGGGGTCA‑GGCCCGCCGATTTCTTTACCGAACAGATTGATGAAAAACGTCAGGATGTGGTGGCGAACTTTAACGGt                                                                                                                                                                         <  2:176181/231‑1 (MQ=255)
                                                                       tcccACGCGTGGTCGGTCAGCAAGAGGCGCTGTTTAGTCGGGTTTATAAAGGTCTGATTGCGGTGGATCCGCATGGTTACATGTCCGCCAT‑CAATCGTAACGCAAGAAAGATGCTGGGGCTGAGTTCCCCCGGACGGCAATGGTTGGGTAACCCCATTGTTGAAGTGGTCA‑GGCCCGCCGATTTCTTTACCGAACAGATTGATGAAAAACGGCAGGATGTGGTGGCGAACTTTAACGGt                                                                                                                                                                         <  1:503516/236‑1 (MQ=255)
                                                                            cgcgTGGTCCGTCAGCAAGAGGCGCTGTTTAGTTCGGTTTATGAAGGGCTGATTGCGGTGGATCCGCATGGTTACATTACCGCCAT‑CAATCGTAACGCAAGAAAGATGCTGGGGCTGAGTTCCCCCGGACGGCAATGGTTGGGTAAACCCATTGTTGAAGTGGTCA‑GGCCCGCCGATTTCTTTACCGAACAGATTAGTGAAAAACGTCAGGATGTGGGGGGGAACTTTAACGGTCTGAg                                                                                                                                                                    >  1:514066/1‑239 (MQ=255)
                                                                             gcgTGGTCCGTCAGCAAGAGGCGCTGTTTAGTTCGGTTTATGAAGGGCTGATTGCGGTGGATCCGCATGGTTACATTACCGCCAT‑CAATCGTAACGCAAGAAAGATGCTGGGGCTGAGTTCCCCCGGACGGCAATGGTTGGGTAAACCCATTGTTGAAGGGGTCA‑GGCCCGCCGATTTCTTTTCCGACAAAATTGATAAAAAACGTCAGGAGGGGGTGGCGAACTTAAAGGGTcgaagg                                                                                                                                                                   >  2:503516/1‑234 (MQ=255)
                                                                                       tCAGCAAGAGGCGCTGTTTAGTTCGGTTTATGAAGGGCTGATTGCGGTGGATCCGCATGGTTACATTACCGCCAT‑CAATCGTAACGCAAGAAAGATGCTGGGGCTGAGTTCCCCCGGACGGCAATGGTTGGGTAAACCCATTGTTGAAGTGGTCA‑GGCCCGCCGATTTCTTTACCGACCAGATTGATGAAAAACGTCAGGTGTTGGTGGCGAACTTTAACGGACTAAGCGTTATTGC‑Ca                                                                                                                                                        >  1:417681/1‑239 (MQ=255)
                                                                                             agagGCGCTGTTTAGTTCGGTTTATGAAGGGCTGATTGCGGTGGATCCGCATGGTTACATTACCGCCAT‑CAATCGTAACGCAAGAAAGATGCTGGGGCTGAGTTCCCCCGGACGGCAATGGTTGGGTAAACCCATTGTTGAAGTGGTCA‑GGCCCGCCGATTTCTTTACCGAACAGATTGATGAAAAACGTCCGGATGTGGTGGCGAACTTTAACGGTCTGAGCGTTTTTGC‑CAACCggg                                                                                                                                                  >  1:824523/1‑239 (MQ=255)
                                                                                                    cTGTTTAGTTCGGTTTATGAAGGGCTGATTGCGGTGGATCCGCATGGTTACATTACCGCCAT‑CAATCGTAACGCAAGAAAGATGCTGGGGCTGAGTTCCCCCGGACGGCAATGGTTGGGTAAACCCATTGTTGAAGTGGTCA‑GGCCCGCCGATTTCTTTACCGAACAGATTGATGAAAAACGTCAGGATGTGGTGGCGAACTTTAACGGTCTGAGCGTTATTGC‑CAACCGGGAAGCTAt                                                                                                                                           >  1:492816/1‑239 (MQ=255)
                                                                                                     tttttttttggttttATAGAGGGCTGATGGGGGGGGTCCCGGAGGGTACAATTCCCGCCAA‑CAAGCGAAACGCAAGAAAGATGGGGGGGTTGGGTTCCCCGGGCCGGAATTGGTGGGGTAAACCCATTGTGGAAGTGGTCA‑GGCCCGCCGATTTCTTTACCGAACAGATTGATGAAAAACGTCAGGATGTGGTGGCGAACTTTAACGGTCTGAGCGTTATTGC‑CAACCGGGAAGCTAt                                                                                                                                           <  2:387058/226‑1 (MQ=255)
                                                                                                     tGTTTAGTTCGGTTTATGAAGGGCTGATTGCGGTGGATCCGCATGGTTACATTACCGCCAT‑CAATCGTAACGCAAGAAAGATGCTGGGGCTGAGTTCCCCCGGACGGCAATGGTTGGGTAAACCCATTGTTGAAGTGGTCA‑GGCCCGCCGATTTCTTTACCGAACAGATTGGTGAAAAACGTCAGGATGTGGTGGCCAACTTTAACGGCCTGAGCGTTATTTC‑CAACCGGGAAgacaa                                                                                                                                           >  2:83733/1‑234 (MQ=255)
                                                                                                         tAGTTCGGTTTAAGAAGGTCTGTTTGGGGTGGACCCGAATGGTTCAATTACCGCCAT‑CAATCGAAACGCAAGAAAGATGCTGGGGCTGAGTTCCCCCGGACGGCAATGGTTGGGTAAACCCATTGTTGAAGTGGTCA‑GGCCCGCCGATTTCTTTACCGAACAGATTGATGAAAAACGTCAGGATGTGGTGGCGAACTTTAACGGTCTGAGCGTTATTGC‑CAACCGGGAAGCTATTCGtt                                                                                                                                      <  1:530208/239‑1 (MQ=255)
                                                                                                              cGGTTTATGAAGGGTTGAGTGGGGTGGTCCCGCATGGTTACATTACCGCCAT‑CAAGCGTAACGCAAGAAAGATGCTGGGGCTGAGTTCCCCCGGACGGCATTGTGTGGGTAACCCCATTGTTGAAGTGGTCA‑GGCCCGCCGATTTCTTTACCGAACAGATTGATGAAAAACGTCAGGATGTGGTGGCGAACTTTAACGGTCTGAGCGTTATTGC‑CAACCGGGAAGCTATTCGTTCAGGt                                                                                                                                 <  1:218396/239‑1 (MQ=255)
                                                                                                                 tttATGAAGGGCTGATTGCGGTGGATCCGCATGGTTACATTACCGCCAT‑CAATCGTAACGCAAGAAAGATGCTGGGGCTGAGTTCCCCCGGACGGCAATGGTTGGGTAAACCCATTGTTGAAGTGGTCA‑GGCCCGCCGATTTCTTTACCGAACAGATTGATGAAAAACGTCAGGATGTGGTGGCGAACTTTAACGGTCTGAGCGTTATTGC‑CAACCGGGAAGCTATTCGTTCAGGtgat                                                                                                                              >  1:606379/1‑239 (MQ=255)
                                                                                                                 tttATGAAGGGCTGATTGCGGTGGATCCGCATGGTTACATTACCGCCAT‑CAATCGTAACGCAAGAAAGATGCTGGGGCTGAATTCCCCCGGACGGCAATGGTGGGGAAAAACCATTGTTGAGGGGGTCA‑GGCCCGCCGATTTCTTTACCGAACAGATTGATGACAAACGTCGGGCAGTGGTGGCGAACTTTAAAGGAGTAAGCGGTATTGC‑CAACCGGGAGGCAACTAGTTTAgatgat                                                                                                                              >  2:29588/1‑239 (MQ=255)
                                                                                                                 tttATGAAGGGCTGATTGCGGTGGATCCGCATGGTTAAATTACCGCCAT‑CAATCGTACCGCAAGAAAAATGCTGGGGCTTAGTTCCCCAGGACGGCAATGGTTGGGTAAACCCATAGTTGCAGGGGCCA‑GGCCCGCCAATTCTTTTACCGAACAGCTTGTTGCAAAACGTACGGATGTGGTGTCGAACTTTAAGTGTCTGCACATTATTAC‑CAACCCGGACGCAATTCGGTCGGGtggt                                                                                                                              >  2:507655/1‑237 (MQ=255)
                                                                                                                      gAAGGGCTGATTGCGGTGGATCCGCATGGTTACATTACCGCCAT‑CAATCGTAACGCAAGAAAGATGCTGGGGCTGAGTTCCCCCGGACGGCAATGGTTGGGTAAACCCATTGTTGAAGTGGTCA‑GGCCCGCCGATTTCTTTACCGAACAGATTGATGAAAAACGTCAGGATGTGGTGGCGAACTTTAACGGTCTGAGCGTTATTGC‑CAACCGGGAAGCTATTCGGTAAGGGGAT‑GAttt                                                                                                                        >  2:11572/1‑239 (MQ=255)
                                                                                                                          ggCTGATTGCGGTGGATCCGCATGGTTACATTACCGCCAT‑CAATCGTAATGCAAGAAAGATGCTGGGGCTGAGTTCCCCCGGACGGCAATGTTTGGGTAAACCCATTGTTGAAGGGGTCA‑GGCCCGCCGATTTCTTTACCGAACAGATTGATGAAAAACGTCCGGTGTTGGGGGCGAAATTTAACGGTATGAGCGTTATTGC‑CAACCGGGAAGCTATTCGTTCAGGTGGT‑GGttggtgg                                                                                                                    >  2:656093/1‑234 (MQ=255)
                                                                                                                           gCTGATTGCGGTGGATCCGCATGGTTACATTACCGCCAT‑CAGTCGTAACGCAAGAAAGGTGCTGGGGCTGAGTTCCCCCGGACGGCAATGGTTGGGGAAACCCATTTTTGAAGGGGTCA‑GGCCCGCCAATTTCTTTACCGAACAGATGTATGAAAAACGTCAGGATGTGGTGGCGAACTTTAACGGCTGGGGCGTTATTGC‑CAACCGTGAAGCCTTTCGTGCAGGTGAT‑GATTTgtggg                                                                                                                   >  2:843488/1‑235 (MQ=255)
                                                                                                                           gCTGATTGCGGTGGATCCGCATGGTTACATTACCGCCAT‑CAATCGTAACGCAAGAAAGATGCTGGGGCTGAGTTCCCCCGGACGGCAATGGTTGGGTAAACCCATTGTTGAAGTGGTCA‑GGCCCGCCGATTTCTTTACCGAACAGATTGATGAAAAACGTCAGGATGTGGTGGCGAACTTTAACGGTCTGGGCGTTATTGC‑CAACCGGGAAGCTATTCGTTTAGGGGAT‑GATTTGCTgg                                                                                                                   >  2:716386/1‑239 (MQ=255)
                                                                                                                               aTTGCGGTGGATCCGCATGGTTACATTACCGCCAT‑CAATCGTAACGCAAGAAAGATGCTGGGGCTGAGTTCCCCCGGACGGCAATGGTTGGGTAAACCCATTGTTGAAGTGGTCA‑GGCCCGCCGATTTCTTTACCGAACAGATTGATGAAAAACGTCAGGATGTGGTGGCGAACTTTAACGGTCTGAGCGTTATTGC‑CAACCGGGAAGCTATTCGTTCAGGGGAT‑GATTTGCTggggg                                                                                                                >  1:96983/1‑238 (MQ=255)
                                                                                                                                tggtgggggATCCGGTTTGTTTCATTTCCCCCATGCAAGTGTAACCTAGGAATGGAGCTGGGGCCGGTTTCCCCCGGACGGAAATGGTTGGGTAAACCGATTGTGGAGGTGGTCA‑GGCCCGCCGATTTCTTTACCGAACAGATTGATGAAAAACGTCAGGATGTGGTGGCGAACTTTAACGGTCTGAGCGTTATTGC‑CAACCGGGAAGCTATTCGTTCAGGTGAT‑GATTTGCTGGGGGc                                                                                                               <  2:616327/232‑1 (MQ=255)
                                                                                                                                         aaccGCATGGTTACATTACCGCCAT‑CAATCGTAACGCAAGAAAGATGCTGGGGCTGAGTTCCCCCGGACGGCAATGGTTGGGTAAACCCATTGTTGAAGTGGTCA‑GGCCAGCCGATTTCTTTACCGAACAGATTGATGAAAAACGTCAGGCGGTGGTGGCGCCCTTTAACGGTCTGAGCGTTATTGC‑CCACCGGGGAGCACTTCGTTACGGTGAT‑GATTTTCTGGGGGCCATTAtaa                                                                                                       >  2:765198/3‑235 (MQ=255)
                                                                                                                                                  gTTACATTACCGCCAT‑CAATCGTAACGCAAGAAAGATGCTGGGGCTGAGTTCCCCCGGACGGCAATGGTTGGGTAAACCCATTGTTGAAGTGGTCA‑GGCCCGCCGATTTCTTTACCGACCAGATTGATGAAAAACGTCAGGATGTGGTGGCGAACTTTAACGGGCTGAGCGTTATTGC‑AAACCGGGAAGCTATTAGGTCAGGTGAG‑GTTTTGTTGGGGGCCATTATCAGC‑tttgtaggt                                                                                           >  2:691453/1‑233 (MQ=255)
                                                                                                                                                      cATTACCGCCAT‑CAATCGTAACGCAAGAAAGATGCTGGGGCTGAGTTCCCCCGGACGGCAATGGTTGGGTAAACCCATTGTTGAAGTGGTCA‑GGCCCGCCGATTTCTTTACCGAACAGATTGATGAAAAACGTCAGGATGTGGTGGCGAACTTTAACGGTCTGAGCGTTATTGC‑CAACCGGGAAGCTATTCGTTCAGGTGAT‑GATTTGCTGGGGGCCATTATCAGC‑TTTCGTAGTAAAg                                                                                       >  1:832949/1‑239 (MQ=255)
                                                                                                                                                         tACCGCCAT‑CAATCGTAACGCAAGAAAGATGCTGGGGCTGAGTTCCCCCGGACGGCAATGGTTGGGTAAACCCATTGTTGAAGTGGTCA‑GGCCCGCCGATTTCTTTACCGAACAGATTGATGAAAAACGTCAGGATGTGGTGGCGAACTTTAACGGTCTGAGCGTTATTGC‑CAACCGGGAAGCTATTCGTTCAGGTGCT‑GATTTGCTGGGGGCCATTATCAGC‑TTTCGTAGTAAAGACg                                                                                    >  1:484024/1‑239 (MQ=255)
                                                                                                                                                            cGCCAT‑CACCCGTAACGCAAGACAGATGCTGGGGCTGAGTTCCCCCGGACGGCAATGTGTGGGTAAACCCATTGTTGCAGGGGTCC‑GGCCCCCCGATTTCTTTACCGAACAGATTGTTGAAAAACACCGGGATG‑GGTGGTCGACTTTTAAGGACCGGGGGTTTATGCACAACCGGGAAGATATTC‑ATCCGGGGATGGATTTGC‑GGGGGGCACTGTCAGCAGTTGGTAGAAAAGAGGaaa                                                                                 >  2:1029196/1‑239 (MQ=255)
                                                                                                                                                            cGCCAT‑CAATCGTAACGCAAGAAAGATGCTGGGGCTGAGTTCCCCCGGACGGCAATGGTTGGGTAAACCCATTGTTGAAGTGGTCA‑GGCCCGCCGATTTCTTTACCGAACAGATTGATGAAAAACGTCAGGATGTGGTGGCGAACTTTAACGGTCTGAGCGTTATTGC‑CAACCGGGAAGCTATTCGTTCAGGTGAT‑GATTTGCTGGGGGCCATTATCAGC‑TTTCGTAGTAAAGACGaaa                                                                                 <  1:44401/239‑1 (MQ=255)
                                                                                                                                                             gCCAT‑CAATCGTAACGCAAGAAAGATGCTGGGGCTGAGTTCCCCCGGACGGCAATGGTTGGGTAAACCCATTGTTGAAGTGGTCA‑GGCCCGCCGATTTCTTTACCGAACAGATTGATGAAAAACGTCAGGATGTGGTGGCGAACTTTAACGGTCTGAGCGTTATTGC‑CAACCGGGGAGCTATTCGTTCAGGGGAT‑GATTTGCTGGGGGCCATTATCAGC‑TTTCGTGGTAAAGCCGAAAt                                                                                >  1:549944/1‑239 (MQ=255)
                                                                                                                                                                aT‑CAATCGTAACGCAAGAAAGATGCTGGGGCTGAGTTCCCCCGGACGGCAATGGTTGGGTAAACCAATTGTTGAAGTGGTCA‑GGCCCGCCGATTTCTTTACCGAACAGATTGATGAAAAACGTCAGGATGTGGTGGCGAACTTTAACGGTCTGAGCGTTATTGC‑CAACCGGGAAGCTATTCGTTCAGGGGAT‑GTTTTGCTGTGGGCCATTATCGTG‑TTTGCTATTAAAAACGAAAttt                                                                              >  2:227084/1‑238 (MQ=255)
                                                                                                                                                                 t‑CAATCGTAACGCAAGAAAGATGCTGGGGCTGAGTTCCCCCGGACGGCAATGGTTGGGTAAACCCATTGTTGAAGTGGTCA‑GGCCCGCCGATTTCTTTACCGAACAGATTGATGAAAAACGTCAGGATGTGGTGGCGAACTTTAACGGTCTGAGCGTTATTGC‑CAACCGGGAAGCTATTCGTTCAGGTGAT‑GATTTGCTGGGGGCCATTATCAGC‑TTTCGTAGTAAAGCCGAAATTTc                                                                             >  1:835742/1‑238 (MQ=255)
                                                                                                                                                                       cGTAACGCAAGAAAGATGCTGGGGCTGAGTTCCCCCGGACGGCAATGGTTGGGTAAACCCATTGTTGAAGTGGTCA‑GGCCCGCCGATTTCTTTACCGAACAGATTGATGAAAAACGTCAGGCTGTGGTGGCGAACTTTAACGGTCTGAGCGTTATTGC‑CAACCGGGAAGCTATTCGTTCAGGTGAT‑GATTTGCTGGGGGCCATTATCAGC‑TTTCGTAGTAAAGACGAAATTTCCACCCt                                                                       >  2:595992/1‑239 (MQ=255)
                                                                                                                                                                            cgaaagaaagaTGTTGGGGCTGAGTTCCCCCGGACGGCAATGGTTGGGTAACCCCATTGTTGAAGGGGTCA‑GGCCCGCCGATTTCTTTACCGAACAGATTGATGAAAAACGTCAGGATGTGGTGGCGAACTTTAACGGTCTGAGCGTTATTGC‑CAACCGGGAAGCTATTCGTTCAGGTGAT‑GATTTGCTGGGGGCCATTATCAGC‑TTTCGTAGTAAAGACGAAATTTCCACCCTCAATg                                                                  <  1:929341/236‑1 (MQ=255)
                                                                                                                                                                                 gaaagaTGCTGGGGCTGAGTTCCCCCGGACGGCAATGGTTGGGTAAACCCATTGTTGAAGTGGTCA‑GGCCCGCCGATTTCTTTACCGAACAGATTGATGAAAAACGTCAGGATGTGGTGGCGAACTTTAACGGTCTGAGCGTTATTGC‑CAACCGGGAAGCTATTCGTTCAGGGGAT‑GATTTGCTGGGGGCCATTATCAGC‑TTTCGTAGTAAAGACGAAATTTCCACCCTCAATGCGCa                                                              >  1:344203/1‑238 (MQ=255)
                                                                                                                                                                                       tgttGGGGCTTGTTTCCCCCGGGCGGAATTGGTTGGGTAAACCCGTTGTAGAAGGGGTAA‑GGCCCGCCGATTTCTTTACCGAACAGATTGATGAAAAACGTCAGGATGTGGTGGCGAACTTTAACGGTCTGAGCGTTATTGC‑CAACCGGGAAGCTATTCGTTCAGGTGAT‑GATTTGCTGGGGGCCATTATCAGC‑TTTCGTAGTAAAGACGAAATTTCCACCCTCAATGCGCAACTGAc                                                        <  2:96983/235‑1 (MQ=255)
                                                                                                                                                                                              gtttgttccccccGGACGGCAATGGTTGGGTAAACCCATTGTTGAAGGGGTCA‑GGCCCGCCGATTTCTTTACCGAACAGATTGATGAAAAACGTCAGGATGTGGTGGCGAACTTTAACGGTCTGAGCGTTATTGC‑CAACCGGGAAGCTATTCGTTCAGGTGAT‑GATTTGCTGGGGGCCATTATCAGC‑TTTCGTAGTAAAGACGAAATTTCCACCCTCAATGCGCAACTGACGCAAATa                                                 <  1:655455/230‑1 (MQ=255)
                                                                                                                                                                                                            gtcgtcattttttGGGTAACCCCTTTTTTGAGGTGGTCA‑GGCCCGTCGATTTCTTTACCGAACAGATTGATGAAAAACGTCAGGATGGGGTTGCTAACTTTAACGGTCTGAGCGTTATTGC‑CAACCGGGAAGCTATTCGTTCAGGTGAT‑GATTTGCTGGGGGCCATTATCAGC‑TTTCGTAGTAAAGACGAAATTTCCACCCTCAATGCGCAACTGACGCAAATAAAACAATACGTTgag                                  <  2:514066/228‑1 (MQ=255)
                                                                                                                                                                                                            gACGGCAATGGTTGGGTAAACCCATTGTTGAAGTGGTCA‑GGCCCGCCGATTTCTTTACCGAACAGATTGATGAAAAACGTCAGGATGTGGTGGCGAACTTTAACGGTCTGAGCGTTATTGC‑CAACCGGGAAGCTATTCGTTCAGGTGAT‑GATTTGCTGGGGGCCATTAt                                                                                                         <  1:554037/169‑1 (MQ=255)
                                                                                                                                                                                                            gACGGCAATGGTTGGGTAAACCCATTGTTGAAGTGGTCA‑GGCCCGCCGATTTCTTTACCGAACAGATTGATGAAAAACGTCAGGATGTGGTGGCGAACTTTAACGGTCTGAGCGTTATTGC‑CAACCGGGAAGCTATTCGTTCAGGTGAT‑GATTTGATGGGGGCCATTAt                                                                                                         >  2:554037/1‑169 (MQ=255)
                                                                                                                                                                                                             aCGGCATTGGTTGGGTAACCCCATTTTTAATGTGGCCA‑GGCCCGCCGATTTCTTTACCGAACAGATTGATGAAAAACGGCAGGATGTGGTGGCGAACTTTAACGGTCTGAGCGTTATTGC‑CAACCGGGAAGCTATTCGTTCAGGTGAT‑GATTTGCTGGGGGCCATTATCAGC‑TTTCGTAGTAAAGACGTAATTTCCACCCTCAATGCGCAACTGACGCAAATAAAACAATACGTTgag                                  <  1:545622/238‑1 (MQ=255)
                                                                                                                                                                                                             aCGGCAATGGTTGGGTAAACCCATTTTTGAAGGGGTCA‑GGCCCGCCGATTTCTTTACCGAACAGATTGATGAAAAACGTCAGGATGTGGTGGCGAACTTTAACGGTCTGAGCGTTATTGC‑CAACCGGGAAGCTATTCGTTCAGGTGAT‑GATTTGCTGGGGGCCATTATCAGC‑TTTCGTAGTAAAGACGAAATTTCCACCCTCAATGCGCAACTGACGCAAATAAAACAATACGTTgag                                  <  1:691453/238‑1 (MQ=255)
                                                                                                                                                                                                               tgtattGGTTGAGTAAACCTTTTTTTGAA‑TGGGCACGCCCCTCCGATTTCTTTACCGAACAGATTGATGAAAAACGTCAGGATGTGGTGGCGAACTTTAACGGTCTGAGCGTTATTGC‑CAACCGGGAAGCTATTCGTTCAGGTGAT‑GATTTGCTGGGGGCCATTATCAGC‑TTTCGTAGTAAAGACGAAATTTCCACCCTCAATGCGCAACTGACGCAAATAAAACAATACGTTGagag                                <  2:835742/233‑1 (MQ=255)
                                                                                                                                                                                                               gtaaagggTTGGGTAAACCCATTGTTGAA‑TGGGCAGGGCCCGCCGATTTCTTTACCGAACAGATTGATGAAAAACGTCAGGATGTGGTGGCGAACTTTAACGGTCTGAGCGTTATTGC‑CAACCGGGAAGCTATTCGTTCAGGTGAT‑GATTTGCTGGGGGCCATTATCAGC‑TTTCGTAGTAAAGACGAAATTTCCACCCTCAATGCGCAACTGACGCAAATAAAACAATACGTTGagag                                <  1:1001341/232‑1 (MQ=255)
                                                                                                                                                                                                                                cccATTGTTGAAGTGGTCA‑GGCCCGCCGATTTCTTTACCGAACAGATTGATGAAAAACGTCAGGATGTGGTGGCGAACTTTAACGGTCTGAGCGTTATTGc‑c                                                                                                                                                         <  1:528258/102‑1 (MQ=255)
                                                                                                                                                                                                                                cccATTGTTGAAGTGGTCA‑GGCCCGCCGATTTCTTTACCGAACAGATTGATGAAAAACGTCAGGATGTGGTGGCGAACTTTAACGGTCTGAGCGTTATTGc‑c                                                                                                                                                         >  2:528258/1‑102 (MQ=255)
                                                                                                                                                                                                                                 ccATTGTTGAAGTGGTCA‑GGCCCGCCGATTTCTTTACCGAACAGATTGATGAAAAACGTCAGGATGTGGTGGCGAACTTTAACGGTCTGAGCGTTATTGC‑CAACCGGGAAGCTATTCGTTCAGGTGAT‑GATTTGCTGGGGGCCATTATCAGC‑TTTCGTAGTAAAGACGAAATTTCCACCCTCAATGCGCAACTGACGCAAATAAAACAATACGTTGAGAGCCTTCGTACATTGCGacac             >  1:108207/1‑239 (MQ=255)
                                                                                                                                                                                                                                  ggttgttgAATTGGCCA‑GCCCCGCCGATTTCTTTACCGAACAGATTGATGAAAAACGTCAGGATGTGGTGGCGAACTTTAACGGTCTGAGCGTTATTGC‑CAACCGGGAAGCTATTCGTTCAGGTGAT‑GATTTGCTGGGGGCCATTATCAGC‑TTTCGTAGTAAAGACGAAATTTCCGCCCTCAATGCGCAACTGACGCAAATAAAACAATACGTTGAGAGCCTTCGTACATTGCGACACg            <  1:519164/237‑1 (MQ=255)
                                                                                                                                                                                                                                    ttttagAAGGGGTCA‑GGCCCGCCGATTTCTTTACCGAACAGATTGATGAAAAACGTCAGGATGTGGTGGCGAACTTTAACGGTCTGAGCGTTATTGC‑CAACCGGGAAGCTATTCGTTCAGGTGAT‑GATTTGCTGGGGGCCATTATCAGC‑TTTCGTAGTAAAGACGAAATTTCCACCCTCAATGCGCAACTGACGCAAATAAAACAATACGTTGAGAGCCTTCGTACATTGCGACACGAg          <  1:716386/234‑1 (MQ=255)
                                                                                                                                                                                                                                     tgttgAAGTGTTCA‑GGCCCGCCGATTTCTTTACCGAACAGATTGATGAAAAACGTCAGGATGTGGTGGCGAACTTTAACGGTCTGCGCGGTATTGC‑CAACCGGGAAGCTATTCGTTCAGGGGCT‑GATTTGCTAGGGTCCATTATCCGC‑TTTCGTAGTAAAAGCGAAATTTCCACCATCAATGCGTAACTTACGCACATAAAACAAAAAGATAAGACCCTTCGTACTTTGCGccaataca         >  2:945858/1‑234 (MQ=255)
                                                                                                                                                                                                                                     tgttgAAGTGGTCA‑GGCCCGCCGATTTCTTTACCGAACAGATTGATGAAAAACGTCAGGATGTGGTGGCGAACTTTAACGGTCTGAGCGTTATTGC‑CAACCGGGAAGCTATTCGTTCAGGTGAT‑GATTTGCTGGGGGCCATTATCAGC‑TTTCGTAGTAAAGACGAAATTTCCACCCTCAATGCGCAACTGACGCAAATAAAACAATACGTTGCGAGCCTTCTTCCATTGCGACACGAGc         >  1:767929/1‑239 (MQ=255)
                                                                                                                                                                                                                                     tgttgAAGTGGTCA‑GGCCCGCCGATTTCTTTACCGAACAGATTGATGAAAAACGTCAGGATGTGGTGGCGAACTTTAACGGTCTGAGCGTTATTGC‑CAACCGGGAAGCTATTCGTTCAGGTGAT‑GATTTGCTGGGGGCCATTATCAGC‑TTTCGTAGTAAAGACGAAATTTCCACCCTCAATGCGCAACTGACGCAAATAAAACAATACGTTGAGAGCCTTCGTACATTGCGACCCGAGc         >  1:293664/1‑239 (MQ=255)
                                                                                                                                                                                                                                     tgttgAAGTGGTCA‑GGCCCGCCGATTTCTTTACCGAACAGATTGATGAAAAACGTCAGGATGTGGTGGCGAACTTTAACGGTCTGAGCGTTATTGC‑CAACCGGGAAGCTATTCGTTCAGGTGAT‑GATTTGCTGGGGGCCATTATCAGC‑TTTCGTAGTAAAGACGAAATTTCCACCCTCAATGCGCAACTGACGCAAATAAAACAATACGTTGAGACCCTTCGTCCATTTGGACACGAg          >  1:398278/1‑238 (MQ=255)
                                                                                                                                                                                                                                       ttgAAGTGGTCA‑GGCCCGCCGATTTCTTTACCGAACAGATTGATGAAAAACGTCAGGATGTGGTGGCGAACTTTAACGGTCTGAGCGTTATTGC‑CAACCGGGAAGCTATTCGTTCAGGTGAT‑GATTTGCTGGGGGCCATTATCAGC‑TTTCGTAGTAAAGACGAAATTTCCACCCTCAATGCGCAACTGACGCAAATAAAACAATACGTTGAGAGCCTTCGTACATTGCGACACGAGCAt       <  1:11572/239‑1 (MQ=255)
                                                                                                                                                                                                                                       tgaaaGTGTTCA‑GGCCCGCCGATTTCTTTACCGAACAGATTGATGAAAAACGTCAGGATGTGGTGGCGAACTTTAACGGTCTGAGCGTTATTGC‑CAACCGGGAAGCTATTCGTTCAGGTGAT‑GATTTGCTGGGGGCCATTATCAGC‑TTTCGTAGTAAAGACGAAATTTCCACCCTCAATGCGCAACTGACGCAAATAAAACAATACGTTGAGAGCCTTCGTACATTGCGACACGAGCAt       <  1:227084/236‑1 (MQ=255)
                                                                                                                                                                                                                                            gTGGTCA‑GGCCCGCCGATTTCTTTACCGAACAGATTGATGAAAAACGTCAGGATGTGGTGGCGAACTTTAACGGTCTGAGCGTTATTGC‑CAACCGGGAAGCTATTCGTTCAGGTGAT‑GATTTGCTGGGGGCCATTATCAGC‑TTTCGTAGTAAAGCCGAAATTTCCACCCTCAATGCGCAACTGACGCAAATAAACAAAAACGTTGAGAGCCTTCGTACATAGCGACCCGAGccacacaa  >  2:267598/1‑232 (MQ=255)
                                                                                                                                                                                                                                             |                                                                                                                                                                                                                                                 
TCTGATGTTGCTGTCGTGGTTCCTGGCCGCGCATATCCGTCGGCAGATGATGGGCATGGAGCCAAAGCAAATCGCACGCGTGGTCCGTCAGCAAGAGGCGCTGTTTAGTTCGGTTTATGAAGGGCTGATTGCGGTGGATCCGCATGGTTACATTACCGCCAT‑CAATCGTAACGCAAGAAAGATGCTGGGGCTGAGTTCCCCCGGACGGCAATGGTTGGGTAAACCCATTGTTGAAGTGGTCA‑GGCCCGCCGATTTCTTTACCGAACAGATTGATGAAAAACGTCAGGATGTGGTGGCGAACTTTAACGGTCTGAGCGTTATTGC‑CAACCGGGAAGCTATTCGTTCAGGTGAT‑GATTTGCTGGGGGCCATTATCAGC‑TTTCGTAGTAAAGACGAAATTTCCACCCTCAATGCGCAACTGACGCAAATAAAACAATACGTTGAGAGCCTTCGTACATTGCGACACGAGCATCTCAA  >  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: ‑