Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 653,352 A→C Q373P (CAG→CCG)  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,3520AC98.6% 240.1 / ‑3.5 71Q373P (CAG→CCG) citAsensory histidine kinase in two‑component regulatory system with CitB
Reads supporting (aligned to +/- strand):  ref base A (1/0);  new base C (38/32);  total (39/32)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.79e-01

GTCTGAGCGTTATTGCCAACCGGGAAGCTATTCGTTCAGGTGATGATTTGCTGGGGGCCATTATCAGCTTTCGTAGTAAAGACGAAATTTCCACCCTCAATGCGCAACTGACGCAAATAAAACAATACGTTGAGAGCCTTCGTACATTGCGACACGAGCATCTCAATTGGATGTCGACGCTCAATGGTCTGTTGCAGATGAAAGAGTATGATCGCGTGCTGGCGATGGTGCAGGGGGAGTCTCAGGCCCAGCAACAGCTTATTGACAGCCTGCGCGAGGCGTTTGCCGATCGCCAGGTGGCGGGGCTGCTTTTTGGTAAAGTGCAGCGCGCCCGGGAACTGGGGCTAAAAATGATCATTGTCCCCGGTAGCCAGCTTTCGCAACTGCCGCCAGGACTGGATAGCACCGAGTTTGCAGCCATTGTGGGCAATTTACTTGATAACGCCTTCGAAGCCAGCCTGC  >  NC_000913/653121‑653582
                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                      
gTCTGAGCGTTATTGCCAACCGGGAAGCTATTCGTTCAGGTGATGATTTGCTGGGGGCCATTATCAGCTTTCGTAGTAAAGACGAAATTTCCACCCTCAATGCGCAACTGACGCAAATAAAACAATACGTTGAGAGCCTTCGTACATTGCGACACGAGCATCTCAATTGGATGTCGACGCTCAATGGTCTGTTGCAGATGAAAGAGTATGATCGCGTGCTGGCGATGGTGCCgggg                                                                                                                                                                                                                                    <  2:524/236‑1 (MQ=255)
gTCTGAGCGTTATTGCCAACCGGGAAGCTATTCGTTCAGGTGATGATTTGCTGGGGGCCATTATCAGCTTTCGTAGTAAAGACGAAATTTCCACCCTCAATGCGCAACTGACGCAAATAAAACAATACGTTGAGAGCCTTCGTACATTGCGACACGAGCATCTCAATTGGATGTCGACGCTCAATGGTCTGTTGCAGATGAAAGAGTATGATCGCGTGCTGGCGATGGTGCCgggg                                                                                                                                                                                                                                    <  1:290717/236‑1 (MQ=255)
     aGCGTTATTGCCAACCGGGAAGCTATTCGTTCAGGTGATGATTTGCTGGGGGCCATTATCAGCTTTCGTAGTAAAGACGAAATTTCCACCCTCAATGCGCAACTGACGCAAATAAAACAATACGTTGAGAGCCTTCGTACATTGCGACACGAGCATCTCAATTGGATGTCGACGCTCAATGGTCTGTTGCAGATGAAAGAGTATGATCGCGTGCTGGCGATGGTGCCGGGTGAGtc                                                                                                                                                                                                                               >  1:158400/1‑236 (MQ=255)
     aGCGTTATTGCCAACCGGGAAGCTATTCGTTCAGGTGATGATTTGCTGGGGGCCATTATCAGCTTTCGTAGTAAAGACGAAATTTCCACCCTCAATGCGCAACTGACGCAAATAAAACAATACGTTGAGAGCCTTCGTACATTGCGACACGAGCATCTCAATTGGATGTCGACGCTCAATGGTCTGTTGCAGATGAAAGAGTATGATCGCGTGCTGGCGATGGTGCCCGGGGATtc                                                                                                                                                                                                                               >  1:732330/1‑236 (MQ=255)
                 aaCCGGGAAGCTATTCGTTCAGGTGATGATTTGCTGGGGGCCATTATCAGCTTTCGTAGTAAAGACGAAATTTCCACCCTCAATGCGCAACTGACGCAAATAAAACAATACGTTGAGAGCCTTCGTACATTGCGACACGAGCATCTCAATTGGATGTCGACGCTCAATGGTCTGTTGCAGATGAAAGAGTATGATCGCGTGCTGGCGATGGTGCCGGGGGAGTCTCAGGCCCAGCa                                                                                                                                                                                                                   <  1:541875/236‑1 (MQ=255)
                 aaCCGGGAAGCTATTCGTTCAGGTGATGATTTGCTGGGGGCCATTATCAGCTTTCGTAGTAAAGACGAAATTTCCACCCTCAATGCGCAACTGACGCAAATAAAACAATACGTTGAGAGCCTTCGTACATTGCGACACGAGCATCTCAATTGGATGTCGACGCTCAATGGTCTGTTGCAGATGAAAGAGTATGATCGCGTGCTGGCGATGGTGCCGGGGGAGTCTCAGGCCCAGCa                                                                                                                                                                                                                   >  1:263054/1‑236 (MQ=255)
                  aCCGGGAAGCTATTCGTTCAGGTGATGATTTGCTGGGGGCCATTATCAGCTTTCGTAGTAAAGACGAAATTTCCACCCTCAATGCGCAACTGACGCAAATAAAACAATACGTTGAGAGCCTTCGTACATTGCGACACGAGCATCTCAATTGGATGTCGACGCTCAATGGTCTGTTGCAGATGAAAGAGTATGATCGCGTGCTGGCGATGGTGCCGGGGGAGTCTCAGGCCCa                                                                                                                                                                                                                      >  1:744673/1‑232 (MQ=255)
                   ccGGGAAGCTATTCGTTCAGGTGATGATTTGCTGGGGGCCATTATCAGCTTTCGTAGTAAAGACGAAATTTCCACCCTCAATGCGCAACTGACGCAAATAAAACAATACGTTGAGAGCCTTCGTACATTGCGACACGAGCATCTCAATTGGATGTCGACGCTCAATGGTCTGTTGCAGATGAAAGAGTATGATCGCGTGCTGGCGATGGTGCCGGGGGAGTCTCAGGCCCAGCAAc                                                                                                                                                                                                                 <  2:230975/236‑1 (MQ=255)
                       gAAGCTATTCGTTCAGGTGATGATTTGCTGGGGGCCATTATCAGCTTTCGTAGTAAAGACGAAATTTCCACCCTCAATGCGCAACTGACGCAAATAAAACAATACGTTGAGAGCCTTCGTACATTGCGACACGAGCATCTCAATTGGATGTCGACGCTCAATGGTCTGTTGCAGATGAAAGAGTATGATCGCGTGCTGGCGATGGTGCCGGGGGAGTCCCAGGCCCAGCAACAGCt                                                                                                                                                                                                             >  1:254967/1‑236 (MQ=255)
                              ttCGTTCAGGTGATGATTTGCTGGGGGCCATTATCAGCTTTCGTAGTAAAGACGAAATTTCCACCCTCAATGCGCAACTGACGCAAATAAAACAATACGTTGAGAGCCTTCGTACATTGCGACACGAGCATCTCAATTGGATGTCGACGCTCAATGGTCTGTTGCAGATGAAAGAGTATGATCGCGTGCTGGCGATGGTGCCGGGGGAGTCTCAGGCCCAGCAACAGCTTATTGa                                                                                                                                                                                                       >  1:182943/1‑235 (MQ=255)
                               tCGTTCAGGTGATGATTTGCTGGGGGCCATTATCAGCTTTCGTAGTAAAGACGAAATTTCCACCCTCAATGCGCAACTGACGCAAATAAAACAATACGTTGAGAGCCTTCGTACATTGCGACACGAGCATCTCAATTGGATGTCGACGCTCAATGGTCTGTTGCAGATGAAAGAGTATGATCGCGTGCTGGCGATGGTGCCGGGGGAGTCTCAGGCCCAGCAACAGCTTATTGa                                                                                                                                                                                                       <  1:582231/234‑1 (MQ=255)
                                  tgcAGGTGATGATTAGCTGGGGGCCATTATCAGCTTTCGTAGTAAAGACGAAATGTCCACCCTCAATGCGCAACTGACGCAAATAAAACAATACGTTGAGAGCATTCGTACATTGCGACACGAGCATCTCAATTGGATGTCGACGCTCAATGGTCTGTTGCAGATGAAAGAGTATGATCGCGTGCTGGCGATGGTGCCGGGGGAGTCTCAGGCCCAGCAACAGCTTATTGACAGcc                                                                                                                                                                                                  <  2:107966/234‑1 (MQ=255)
                                      ggTGATGATTTGCTGGGGGCCATTATCAGCTTTCGTAGTAAAGACGAAATTTCCACCCTCAATGCGCAACTGACGCAAATAAAACAATACGTTGAGAGCCTTCGTACATTGCGACACGAGCATCTCAATTGGATGTCGACGCTCAATGGTCTGTTGCAGATGAAAGAGTATGATCGCGTGCTGGCGATGGTGCCGGGGGAGTCTCAGGCCCAGCAACAGCTTATTGACAGCCTgcg                                                                                                                                                                                              >  1:727224/1‑236 (MQ=255)
                                        tgatgaTTTGCTGGGGGCCATTATCAGCTTTCGTAGTAAAGACGAAATTTCCACCCTCAATGCGCAACTGACGCAAATAAAACAATACGTTGAGAGCCTTCGTACATTGCGACACGAGCATCTCAATTGGATGTCGACGCTCAATGGTCTGTTGCAGATGAAAGAGTATGATCGCGTGCTGGCGATGGTGCCGGGGGAGTCTCAGGCCCAGCAACAGCTTATTGACAGCCTGcgcg                                                                                                                                                                                            >  2:78853/1‑236 (MQ=255)
                                            ggtttGCTGGGAGCCATTAACAGCTTTCGTAGTAAAGACGAAATTTCCACCCTCAATGCGCAACTGACGCAAATAAAACAATACGTTGAGAGCCTTCGTACATTGCGACACGAGCATCTCAATTGGATGTCGACGCTCAATGGTCTGTTGCAGATGAAAGAGTATGATCGCGTGCTGGCGATGGTGCCGGGGGAGTCTCAGGCCCAGCAACAGCTTATTGACAGCCTGCGCGAgg                                                                                                                                                                                         <  2:119267/233‑1 (MQ=255)
                                            gatTTGCTGGGGGCCATTATCAGCTTTCGTAGTAAAGACGAAATTTCCACCCTCAATGCGCAACTGACGCAAATAAAACAATACGTTGAGAGCCTTCGTACATTGCGACACGAGCATCTCAATTGGATGTCGACGCTCAATGGTCTGTTGCAGATGAAAGAGTATGATCGCGTGCTGGCGATGGTGCCgggg                                                                                                                                                                                                                                    <  2:579000/192‑1 (MQ=255)
                                            gatTTGCTGGGGGCCATTATCAGCTTTCGTAGTAAAGACGAAATTTCCACCCTCAATGCGCAACTGACGCAAATAAAACAATACGTTGAGAGCCTTCGTACATTGCGACACGAGCATCTCAATTGGATGTCGACGCTCAATGGTCTGTTGCAGATGAAAGAGTATGATCGCGTGCTGGCGATGGTGCCgggg                                                                                                                                                                                                                                    >  1:579000/1‑192 (MQ=255)
                                                 gCTGGGGGCCATTATCAGCTTTCGTAGTAAAGACGAAATTTCCACCCTCAATGCGCAACTGACGCAAATAAAACAATACGTTGAGAGCCTTCGTACATTGCGACACGAGCATCTCAATTGGATGTCGACGCTCAATGGTCTGTTGCAGATGAAAGAGTATGATCGCGTGCTGGCGATGGTGCCGGGGGAGTCTCAGGCCCAGCAACAGCTTATTGACAGCCTGCGCGAGGCGTTTg                                                                                                                                                                                   >  2:105829/1‑236 (MQ=255)
                                                   tGGGGGCCCTTAGCAGCTTTAGTAGTAAAGACGAAATTTCCACCCTCAATGCGCAACTGACGCAAATAAACCAATACGTTGAGAGCCTTCGTACATTGCGACACGAGCATCTCAATTGGATGTCGACGCTCAATGGTCTGTTGCAGATGAAAGAGTATGATCGCGTGCTGGCGATGGTGCCGGGGGAGTCTCAGGCCCAGCAACAGCTTATTGACAGCCTGCGCGAGGCGTTTGcc                                                                                                                                                                                 <  2:48028/236‑1 (MQ=255)
                                                    gggggCCATTATCAGCTTTCGTAGTAAAGACGAAATTTCCACCCTCAATGCGCAACTGACGCAAATAAAACAATACGTTGAGAGCCTTCGTACATTGCGACACGAGCATCTCAATTGGATGTCGACGCTCAATGGTCTGTTGCAGATGAAAGAGTATGATCGCGTGCTGGCGATGGTGCCGGGGGAGTCTCAGGCCCAGCAACAGCTTATTGACAGCCTGCGCGAGGCGTTTGCCg                                                                                                                                                                                <  2:263054/236‑1 (MQ=255)
                                                      gggCCATTATCAGCTTTCGTAGTAAAGACGAAATTTCCACCCTCAATGCGCAACTGACGCAAATAAAACAATACGTTGAGAGCCTTCGTACATTGCGACACGAGCATCTCAATTGGATGTCGACGCTCAATGGTCTGTTGCAGATGAAAGAGTATGATCGCGTGCTGGCGATGGTGCCGGGGGAGTCTCAGGCCCAGCAACAGCTTATTGACAGCCTGCGCGAGGCGTTTGCCGa                                                                                                                                                                               >  1:423172/1‑235 (MQ=255)
                                                       ggCCATTATCAGCTTTCGTAGTAAAGACGAAATTTCCACCCTCAATGCGCAACTGACGCAAATAAAACAATACGTTGAGAGCCTTCGTACATTGCGACACGAGCATCTCAATTGGATGTCGACGCTCAATGGTCTGTTGCAGATGAAAGAGTATGATCGCGTGCTGGCGATGGTGCCgggg                                                                                                                                                                                                                                    >  2:647741/1‑181 (MQ=255)
                                                       ggCCATTATCAGCTTTCGTAGTAAAGACGAAATTTCCACCCTCAATGCGCAACTGACGCAAATAAAACAATACGTTGAGAGCCTTCGTACATTGCGACACGAGCATCTCAATTGGATGTCGACGCTCAATGGTCTGTTGCAGATGAAAGAGTATGATCGCGTGCTGGCGATGGTGCCgggg                                                                                                                                                                                                                                    <  1:647741/181‑1 (MQ=255)
                                                                 aGCTTTCGTAGTAAAGACGAAATTTCCACCCTCAATGCGCAACTGACGCAAATAAAACAATACGTTGAGAGCCTTCGTACATTGCGACACGAGCATCTCAATTGGATGTCGACGCTCAATGGTCTGTTGCAGATGAAAGAGTATGATCGCGTGCTGGCGATGGTGCCGGGGGAGTCTCAGGCCCAGCAACAGCTTATTGACAGCCTGCGCGAGGCGTTTGCCGATCGCCAGGTGGc                                                                                                                                                                   <  1:105829/236‑1 (MQ=255)
                                                                      tCGTAGTAAAGACGAAATTTCCACCCTCAATGCGCAACTGACGCAAATAAAACAATACGTTGAGAGCCTTCGTACATTGCGACACGAGCATCTCAATTGGATGTCGACGCTCAATGGTCTGTTGCAGATGAAAGAGTATGATCGCGTGCTGGCGATGGTGCCGGGGGAGTCTCAGGCCCAGCAACAGCTTATTGa                                                                                                                                                                                                       >  2:747951/1‑195 (MQ=255)
                                                                      tCGTAGTAAAGACGAAATTTCCACCCTCAATGCGCAACTGACGCAAATAAAACAATACGTTGAGAGCCTTCGTACATTGCGACACGAGCATCTCAATTGGATGTCGACGCTCAATGGTCTGTTGCAGATGAAAGAGTATGATCGCGTGCTGGCGATGGTGCCGGGGGAGTCTCAGGCCCAGCAACAGCTTATTGa                                                                                                                                                                                                       <  1:747951/195‑1 (MQ=255)
                                                                       cGTAGTAAAGACGAAATTTCCACCCTCAATGCGCAACTGACGCAAATAAAACAATACGTTGAGAGCCTTCGTACATTGCGACACGAGCATCTCAATTGGATGTCGACGCTCAATGGTCTGTTGCAGATGAAAGAGTATGATCGCGTGCTGGCGATGGTGCCGGGGGAGTCTCAGGCCCAGCAACAg                                                                                                                                                                                                               <  2:441954/186‑1 (MQ=255)
                                                                       cGTAGTAAAGACGAAATTTCCACCCTCAATGCGCAACTGACGCAAATAAAACAATACGTTGAGAGCCTTCGTACATTGCGACACGAGCATCTCAATTGGATGTCGACGCTCAATGGTCTGTTGCAGATGAAAGAGTATGATCGCGTGCTGGCGATGGTGCCGGGGGAGTCTCAGGCCCAGCAACAg                                                                                                                                                                                                               >  1:441954/1‑186 (MQ=255)
                                                                        gtagtaAAGACGAAATTTCCACCCTCAATGCGCAACTGACGCAAATAAAACAATACGTTGAGAGCCTTCGTACATTGCGACACGAGCATCTCAATTGGATGTCGACGCTCAATGGTCTGTTGCAGATGAAAGAGTATGATCGCGTGCTGGCGATGGTGCCGGGGGAGTCTCAGGCCCAGCAACAGCTTATTGACAGCCTGCGCGAGGCGTTTGCCGATCGCCAGGTGGCGGGgctg                                                                                                                                                            <  1:78853/236‑1 (MQ=255)
                                                                                  cGAAATTTCCACCCTCAATGCGCAACTGACGCAAATAAAACAATACGTTGAGAGCCTTCGTACATTGCGACACGAGCATCTCAATTGGATGTCGACGCTCAATGGTCTGTTGCAGATGAAAGAGTATGATCGCGTGCTGGCGATGGTGCCGGGGGAGTCTCAGGCCCAGCAACAGCTTATTGACAGCCTGCGCGAGGCGTTTGCCGATCGCCAGGTGGCGGGGCTGCTTTTTGGTa                                                                                                                                                  <  1:365173/236‑1 (MQ=255)
                                                                                            aCCCTAAATGCGCAACTGACGAAAATAAAGCAATACGTTGAGAGCCTTCGTACATTGCGACACGAGCATCTCAATTGGATGTCGACGCTCAATGGTCTGTTGCAGATGAAAGAGTATGATCGCGTGCTGGCGATGGTGCCGGGGGAGTCTCAGGCCCAGCAACAGCTTATTGACAGCCTGCGCGAGGCGTTTGCCGATCGCCAGGTGGCGGGGCTGCTTTTTGGTAAAGTGCAgcg                                                                                                                                        <  2:744673/236‑1 (MQ=255)
                                                                                                        cAACTGACGCAAATAAAACAATACGTTGAGAGCCTTCGTACATTGCGACACGAGCATCTCAATTGGATGTCGACGCTCAATGGTCTGTTGCAGATGAAAGAGTATGATCGCGTGCTGGCGATGGTGCCGGGGGAGTCTCAGGCCCAGCAACAGCTTATTGACAGCCTGCGCGAGGCGTTTGCCGATCGCCAGGGGGCGGGGCTGCCTTTTGGTAACTGGCAGCGGGCCCGGGAACt                                                                                                                            >  2:175701/1‑236 (MQ=255)
                                                                                                          aCTGACGCAAATAAAACAATACGTTGAGAGCCTTCGTACATTGCGACACGAGCATCTCAATTGGATGTCGACGCTCAATGGTCTGTTGCAGATGAAAGAGTATGATCGCGTGCTGGCGATGGTGCCGGGGGAGTCTCAGGCCCAGCAACAGCTTATTGTCAGCCTGCGCGAGGCGTTGGCGGACTGCCAGGCAGCGGGGCGGCTTTTTGGTAAAGAGCAGAGCGCACgggcaat                                                                                                                            >  2:398712/1‑230 (MQ=255)
                                                                                                            tGACGCAAATAAAACAATACGTTGAGAGCCTTCGTACATTGCGACACGAGCATCTCAATTGGATGTCGACGCTCAATGGTCTGTTGCAGATGAAAGAGTATGATCGCGTGCTGGCGATGGTGCCGGGGGAGTCTCAGGCCCAGCAACAGCTTATTGACAGCCTGCGCGAGGCGTTTGCCGATCGCCAGGTGGCGGGGCTGCTTTTTGGTAAAGTGCAGCGCGCCCGGGAACTgggg                                                                                                                        >  1:465578/1‑236 (MQ=255)
                                                                                                                         gcAATTCGTTGAGAGCCTTCGTACATTGCGACACGAGCATCGCAATTGGATGTCGACGCTCAATGGTCTGTTGCAGATGAAAGAGTATGATCGCGTGCTGGCGATGGTGCCGGGGGAGTCTCAGGCCCAGCAACAGCTTATTGACAGCCTGCGCGAGGCGTTTGCCGATCGCCAGGTGGCGGGGCTGCTTTTTGGTAAAGTGCAGCGCGCCCGGGAACTGGGGCTAAAAATGATCa                                                                                                           <  1:747254/235‑1 (MQ=255)
                                                                                                                           aaTACGTTGAGAGCCTTCATACATTGCGACACGAGCATCGCAATTGGATGTCGACGCTCAATGGTCTGTTGCAGATGAAAGAGTATGATCGCGTGCTGGCGATGGTGCCGGGGGAGTCTCAGGCCCAGCAACAGCTTATTGACAGCCTGCGCGAGGCGTTTGCCGATCGCCAGGTGGCGGGGCTGCTTTTTGGTAAAGTGCAGCGCGCCCGGGAACTGGGGCTAAAAATGATCAtt                                                                                                         <  1:333588/236‑1 (MQ=255)
                                                                                                                             tACGTTGAGAGCCTTCGTACATTGCGACACGAGCATCTCAATTGGATGTCGACGCTCAATGGTCTGTTGCAGATGAAAGAGTATGATCGCGTGCTGGCGATGGTGCCGGGGGAGTCTCAGGCCCAGCAACAGCTTATTGACAGCCTGCGCGAGGCGTTTGCCGATCGCCAGGGGGCGGGGCGGCTTTTTGGTAAAGTGCAGCGCGCCCGGGAACTGGGGCTAAAAATGATCATTg                                                                                                        >  1:686044/1‑235 (MQ=255)
                                                                                                                               cGTTGAGAGCCTTCGTACATTGCGACACGAGCATCTCAATTGGATGTCGACGCTCAATGGTCTGTTGCAGATGAAAGAGTATGATCGCGTGCTGGCGATGGTGCAGGGGGAGTCTCAGGCCCAGCAACAGCTTATTGACAGCCTGCGCGAGGCGTTTGCCGATCGCCAGGTGGCGGGGCTGCTTTTTGGTCACGTGCAGCGCGCTCGGGAACTGGTGCGACAACTGATCAtt                                                                                                         >  2:178888/1‑232 (MQ=255)
                                                                                                                                  tGAGAGCCTTCGTACATTGCGACACGAGCATCTCAATTGGATGTCGACGCTCAATGGTCTGTTGCAGATGAAAGAGTATGATCGCGTGCTGGCGATGGTGCCGGGGGAGTCTCAGGCCCAGCAACAGCTTATTGACAGCCTGCGCGAGGCGTTTGCCGATCGCCAGGTGGCGGGGCTGCTTTTTGGTAAAGTGCAGCGCGCCCGGGAACTGGGGCTAAAAATGATCATTGTccc                                                                                                    >  1:622909/1‑234 (MQ=255)
                                                                                                                                  tGAGAGCCTTCGTACATTGCGACACGAGCATCTCAATTGGATGTCGACGCTCAATGGTCTGTTGCAGATGAAAGAGTATGATCGCGTGCTGGCGATGGTGCCGGGGGAGTCTCAGGCCCAGCAACAGCTTATTGACAGCCTGCGCGAGGCGTTTGCCGATCGCCAGGTGGCGGGGCTGCTTTTTGGTAAAGTGCAGCGCGCCCGGGAACTGGGGCTAAAAATGATCATTGTccc                                                                                                    <  2:622909/234‑1 (MQ=255)
                                                                                                                                  tGAGAGCCTTCGTACATTGCGACACGAGCATCTCAATTGGATGTCGACGCTCAATGGTCTGTTGCAGATGAAAGAGTATGATCGCGTGCTGGCGATGGTGCCGGGGGAGTCTCAGGCCCAGCAACAGCTTATTGACAGCCTGCGCGAGGCGTTTGCCGATCGCCAGGTGGCGGGGCTGCTTTTTGGTAAAGTGCAGCGCGCCCGGGAACTGGGGCTAAAAATGATCATTGTCCCCg                                                                                                  >  1:671965/1‑236 (MQ=255)
                                                                                                                                             atACATTGCGACAAGTGCATGTCAATTGGATGTCGACGCTCAATGGTCTGGTGAAGATGAAAGAGTATGATCGCGTGCTGGCGATGGTGCCGGGGGAGTCTCAGGCCCAGCAACAGCTTATTGACAGCCTGCGCGAGGCGTTTGCCGATCGCCAGGTGGCGGGGCTGCTTTTTGGTAAAGTGCAGCGCGCCCGGGAACTGGGGCTAAAAATGATCATTGTCCCCAGTTGccgggcg                                                                                       <  2:254967/235‑6 (MQ=255)
                                                                                                                                                        cacGAGCATCTCAATTGGATGTCGACGCTCAATGGTCTGTTGCAGATGAAAGAGTATGATCGCGTGCTGGCGATGGTGCCGGGGGAGTCTCAGGCCCAGCAACAGCTTATTGACAGCCTGCGCGAGGCGTTTGCCGATCGCCAGGTGGCGGGGCTGCTTTTTGGTAAAGTGCAGCGCGCCCGGGAACTGGGGCTAAAAATGATCATTGTCCCCGGTAGCCAGCTTTCGCAACTgcc                                                                            >  1:126459/1‑236 (MQ=255)
                                                                                                                                                           gAGCATCTCAATTGGATGTCGACGCTCAATGGTCTGTTGCAGATGAAAGAGTATGATCGCGTGCTGGCGATGGTGCCGGGGGAGTCTCAGGCCCAGCAACAGCTTATTGACAGCCTGCGCGAGGCGTTTGCCGATCGCCAGGTGGCGGGGCTGCTTTTTGGTAAAGTGCAGCGCGCCCGGGAACTGGGGCTAAAAATGATCATTGTCCCCGGTAGCCAGCTTTCGCAACTgccgc                                                                          <  1:175701/235‑1 (MQ=255)
                                                                                                                                                                tctcAATTGGATGTCGACGCTCAATGGTCTGTTGCAGATGAAAGAGTATGATCGCGTGCTGGCGATGGTGCCGGGGGAGTCTCAGGCCCAGCAACAGCTTATTGACAGCCTGCGCGAGGCGTTTGCCGATCGCCAGGTGGCGGGGCTGCTTTTTGGTAAAGTGCAGCGCGCCCGGGAACTGGGGCTAAAAATGATCATTGTCCCCGGTAGCCAGCTTTCGCAACTGCCGCCAGGAc                                                                    >  2:654079/1‑236 (MQ=255)
                                                                                                                                                                tcacAATTGGAAGTCGACGGTCAATGGGCTGTTGCAGATGAAAGAGTTTGATCGCGGGCTGGCGATGGTGCCGGGGGAGTCTCAGGCCCAGCAACAGCTTATTGACAGCCTGCGCGAGGCGTTTGCCGATCGCCAGGTGGCGGGGCTGCTTTTTGGTAAAGTGCAGCGCGCCCGGGAACTGGGGCTAAAAATGATCATTGTCCCCGGTAGCCAGCTTTCGCAACTGCCGCCAGGa                                                                     <  2:727224/232‑1 (MQ=255)
                                                                                                                                                                  tcAATTGGATGTCGACGCTCAATGGTCTGTTGCAGATGAAAGAGTATGATCGCGTGCTGGCGATGGTGCCGGGGGAGTCTCAGGCCCAGCAACAGCTTATTGACAGCCTGCGCGAGGCGTTTGCCGATCGCCAGGTGGCGGGGCTGCTTTTTGGTAAAGTGCAGCGCGCCCGGGAACTGGGGCTAAAAATGATCATTGTCCCCGGTAGCCAGCTTTCGCAACTGCCGCCACGACTg                                                                  >  1:258952/1‑236 (MQ=255)
                                                                                                                                                                    aaTTGGATGTCGACGCTCAATGGTCTGTTGCAGATGAAAGAGTATGATCGCGTGCTGGCGATGGTGCCGGGGGAGTCTCAGGCCCAGCAACAGCTTATTGACAGCCTGCGCGAGGCGTTTGCCGATCGCCAGGTGGCGGGGCTGCTTTTTGGTAAAGTGCAGCGCGCCCGGGAACTGGGGCTAAAAATGATCATTGTCCCCGGTAGCCAGCTTTCGCAACTGCCGCCAGGACTGGa                                                                <  2:126459/236‑1 (MQ=255)
                                                                                                                                                                           tGTCGACGCTCAATGGTCTGTTGCAGATGAAAGAGTATGATCGCGTGCTGGCGATGGTGCCGGGGGAGTCTCAGGCCCAGCAACAGCTTATTGACAGCCTGCGCGAGGCGTTTGCCGATCGCCAGGTGGCGGGGCTGCTTTTTGGTAAAGTGCAGCGCGCCCGGGAACTGGGGCTAAAAATGATCATTGTCCCCGGTAGCCAGCTTTCGCAACTGCCGCCAGGACTGGATAGCAcc                                                         >  1:76444/1‑236 (MQ=255)
                                                                                                                                                                           tGTCGACGCTCAATGGTCTGTTGCAGATGAAAGAGTATGATCGCGTGCTGGCGATGGTGCCGGGGGAGTCTCAGGCCCAGCAACAGCTTATTGACAGCCTGCGCGAGGCGTTTGCCGATCGCCAGGTGGCGGGGCTGCTTTTTGGTAAAGTGCAGCGCGCCCGGGAACTGGGGCTAAAAATGATCATTGTCCCCGGTAGCCAGCTTTCGCAACTGCCGCCAGGACTGGATAGCAcc                                                         >  1:483452/1‑236 (MQ=255)
                                                                                                                                                                           tGTCGACGCTCAATGGTCTGTTGCAGATGAAAGAGTATGATCGCGTGCTGGCGATGGTGCCGGGGGAGTCTCAGGCCCAGCAACAGCTTATTGACAGCCTGCGCGAGGCGTTTGCCGATCGCCAGGTGGCGGGGCTGCTTTTTGGTAAAGTGCAGCGCGCCCGGGAACTGGGGCTAAAAATGATCATTGTCCCCGGTAGCCAGCTTTCGCAACTGCCGCCAGGACTGGATAGCAcc                                                         >  1:63759/1‑236 (MQ=255)
                                                                                                                                                                              caaCGCTCAATGGTCTGTTGCAGATGAAAGAGTATGATCGCGTGCTGGCGATGGTGCCGGGGGAGTCTCAGGCCCAGCAACAGCTTATTGACAGCCTGCGCGAGGCGTTTGCCGATCGCCAGGTGGCGGGGCTGCTTTTTGGTAAAGTGCAGCGCGCCCGGGAACTGGGGCTAAAAATGATCATTGTCCCCGGTAGCCAGCTTTCGCAACTGCCGCCAGGACTGGATAGCACCGa                                                       <  2:182943/233‑1 (MQ=255)
                                                                                                                                                                              cGACGCTCAATGGTCTGTTGCAGATGAAAGAGTATGATCGCGTGCTGGCGATGGTGCCGGGGGAGTCTCAGGCCCAGCAACAGCTTATTGACAGCCTGCGCGAGGCGTTTGCCGATCGCCAGGTGGCGGGGCTGCTTTTTGGTAAAGTGCAGCGCGCCCGGGAACTGGGGCTAAAAATGATCATTGTCCCCGGTGGCCAGCTTTCGCAACGGCCGCCAGGACTgggtgtgaacagg                                                      >  2:430927/1‑225 (MQ=255)
                                                                                                                                                                              cGACGCTCAATGGTCTGTTGCAGATGAAAGAGTATGATCGCGTGCTGGCGATGGTGCCGGGGGAGTCTCAGGCCCAGCAACAGCTTATTGACAGCCTGCGCGAGGCGTTTGCCGATCGCCAGGTGGCGGGGCTGCTTTTTGGTAAAGTGCAGCGCGCCCGGGAACTGGGGCTAAAAATGATCATTGTCCCCGGTAGCCACCTTTCGCAACTGACGCCCGGACTGGATAGCAccagg                                                      >  2:723489/1‑233 (MQ=255)
                                                                                                                                                                                                    gtagAAAGAGTATGATCGCGTGCTGGCGATGGTGCCGGGGGATTCTCAGGCCCAGCAACAGCTTATTGACAGCCTGCGCGAGGCGTTTGCCGATCGCCAGGTGGCGGGGCTGCTTTTTGGTAAAGTGCAGCGCGCCCGGGAACTGGGGCTAAAAATGATCATTGTCCCCGGTAGCCAGCTTTCGCAACTGCCGCCAGGACTGGATAGCACCGAGTTTGCAGCCATTGTGGGCAAtt                                <  2:76444/233‑1 (MQ=255)
                                                                                                                                                                                                    gATGAAAGAGTATGATCGCGTGCTGGCGATGGTGCCGGGGGAGTCTCAGGCCCAGCAACAGCTTATTGACAGCCTGCGCGAGGCGTTTGCCGATCGCCAGGTGGCGTGGCTGCTTTTTGGTAAAGTGCAGCGCGCCCGGGAACTGGGGCTAAAAATGATCATTGTCCCCGGTAGCCAGCTTTCGCAACTGCCGCCAGGACTGGATAGCACCGAGTTTGCAGCCATTGTGGGCaa                                  >  1:744830/1‑234 (MQ=255)
                                                                                                                                                                                                    gATGAAAGAGTATGATCGCGTGCTGGCGATGGTGCCGGGGGAGTCTCAGGCCCAGCAACAGCTTATTGACAGCCTGCGCGAGGCGTTTGCCGATCGCCAGGTGGCGGGGCTGCTTTTTGGTAAAGTGCAGCGCGCCCGGGAACTGGGGCTAAAAATGATCATTGTCCCCGGTAGCCAGCTTTCGCAACTGCCGCCAGGACTGGATAGCACCGAGTTTGCAGCCATTGTGGGCAAtt                                <  2:483452/236‑1 (MQ=255)
                                                                                                                                                                                                     aTGAAAGAGTATGATCGCGTGCTGGCGATGGTGCCGGGGGAGTCTCAGGCCCAGCAACAGCTTATTGACAGCCTGCGCGAGGCGTTTGCCGATCGCCAGGTGGCGGGGCTGCTTTTTGGTAAAGTGCAGCGCGCCCGGGAACTGGGGCTAAAAATGATCATTGTCCCCGGTAGCCAGCTTTCGCAACTGCCGCCAGGa                                                                     <  1:496668/198‑1 (MQ=255)
                                                                                                                                                                                                     aTGAAAGAGTATGATCGCGTGCTGGCGATGGTGCCGGGGGAGTCTCAGGCCCAGCAACAGCTTATTGACAGCCTGCGCGAGGCGTTTGCCGATCGCCAGGTGGCGGGGCTGCTTTTGGTAAAGTGCAGCGCGCCCGGGAACTGGGGCTAAAAATGATCATTGTCCCCGGTAGCCAGCTTTCGCAACTGCCGCCAGGACTGGATAGCACCGAGGTTGCAGCCATTGTGGGCCATTTa                              >  2:338201/1‑236 (MQ=255)
                                                                                                                                                                                                     aTGAAAGAGTATGATCGCGTGCTGGCGATGGTGCCGGGGGAGTCTCAGGCCCAGCAACAGCTTATTGACAGCCTGCGCGAGGCGTTTGCCGAACGCCAGGTGGCGGGGCTGCTTTTTGGTAAAGTGCAGCGCGCCCGGGAACTGGGGCTAAAAATGATCATTGTCCCCGGTAGCCAGCTTTCGCAACTGCCGCCAGGa                                                                     >  2:496668/1‑198 (MQ=255)
                                                                                                                                                                                                         aaGAGTATGATCGCGTGCTGGCGATGGTGCCGGGGGAGTCTCAGGCCCAGCAACAGCTTATTGACAGCCTGCGCGAGGCGTTTGCCGATCGCCAGGTGGCGGGGCTGCTTTTTGGTAAAGTGCAGCGCGCCCGGGAACTGGGGCTAAAAATGATCATTGTCCCCGGTAGCCAGCTTTCGCAACTGCCGCCAGGACTGGATAGCACCGAGTTTGCAGCCATTGTGGGCAATTTACtt                           >  1:248776/1‑236 (MQ=255)
                                                                                                                                                                                                            agTATGATCGCGTGCTGGCGATGGTGCCGGGGGAGTCTCAGGCCCAGCAACAGCTTATTGACAGCCTGCGCGAGGCGTTTGCCGATCGCCAGGTGGCGGGGCTGCTTTTTGGTAAAGTGCAGCGCGCCCGGGAACTGGGGCTAAAAATGATCATTGTCCCCGGTAGCCAGCTTTCGCAACTGCCGCCAGGACTGGATAGCACCGAGTTTGCAGCCATTGTGGGCAATTTACTTGa                         >  1:137064/1‑235 (MQ=255)
                                                                                                                                                                                                                tGATCGCGTGCTGGCGATGGTGCCGGGGGAGTCTCAGGCCCAGCAACAGCTTATTGACAGCCTGCGCGAGGCGTTTGCCGATCGCCAGGTGGCGGGGCTGCTTTTTGGTAAAGTGCAGCGCGCCCGGGAACTGGGGCTAAAAATGATCATTGTCCCCGGTAGCTAGCTTTCGCAACTGCCGCCAGGACTGGATATCACCGAGTTTTCAGCCATTGTGGTCAATTTACTTGCTAACg                    >  2:229159/1‑236 (MQ=255)
                                                                                                                                                                                                                  aTCGGGTGCGGGGGATGGTGCCGGGGGAGTCTCAGGCCTGGCAACAGCTTATTGACAGCCTGCGCGAGGCGTTTGCCGATCGACAGGTGGCGGGGCTGCTTTTTGGTAAAGTGCAGCGCGCCCGGGAACTGGGGATAAAAATGATCATTGTCACCGGTAGCCAGCTTTCGAAACTGCCGCCAGGACTGGATAGCACCGAGTTTGCAGCCATTGTGGGCAATTTACTTGATAACg                    <  1:229159/234‑1 (MQ=255)
                                                                                                                                                                                                                        tGCTGGCGATGGTGCCGGGGGAGTCTCAGGCCCAGCAACAGCTTATTGACAGCCTGCGCGAGGCGTTTGCCGATCGCCAGGTGGCGGGGCTGCTTTTTGGTAAAGTGCAgcg                                                                                                                                        >  1:379305/1‑112 (MQ=255)
                                                                                                                                                                                                                        tGCTGGCGATGGTGCCGGGGGAGTCTCAGGCCCAGCAACAGCTTATTGACAGCCTGCGCGAGGCGTTTGCCGATCGCCAGGTGGCGGGGCTGCTTTTTGGTAAAGTGCAgcg                                                                                                                                        <  2:379305/112‑1 (MQ=255)
                                                                                                                                                                                                                        tGCTGGCGATGGTGCCGGGGGAGTCTCAGGCCCAGCAACAGCTTATTGACAGCCTGCGCGAGGCGTTTGCCGATCGCCAGGTGGCGGGGCTGCTTTTTGGTAAAGTGCAGCGCGCCCGGGAACTGGGGCTAAAAATGATCATTGTCCCCGGTAGCCAGCTTTCGCAACTGCCGCCAGGACTGGATAGCACCGAGTTTGCAGCCATTGTGGGCAATTTACTTGATAACGCCTTCGaa            >  2:536883/1‑236 (MQ=255)
                                                                                                                                                                                                                               gagtgTGCCGGGGGAGTCTCGGGCCCAGCAACAGCTTAGTGACAGCCTGCGCGAGGCGTTTGCCGATCGCCAGGTGGCGGGGCTGCTTTTTGGTAAAGTGCAGCGCGCCCGGGAACTGGGGCTAAAAATGATCATTGTCCCCGGTAGCCAGCTTTCGCAACTGCCGCCAGGACTGGATAGCACCGAGTTTGCAGCCATTGTGGGCAATTTACTTGATAACGCCTTCGAagccagc      <  2:258952/231‑1 (MQ=255)
                                                                                                                                                                                                                                  ggTGCCGGGGGAGTCTCAGGCCCAGCAACAGCTTATTGACAGCCTGCGCGAGGCGTTTGCCGATCGCCAGGTGGCGGGGCTGCTTTTTGGTAAAGTGCAGCGCGCCCGGGAACTGGGGCTAAAAATGATCATTGTCCCCGGTAGCCAGCTTTCGCAACTGCCGCCAGGa                                                                     >  1:215079/1‑169 (MQ=255)
                                                                                                                                                                                                                                  ggTGCCGGGGGAGTCTCAGGCCCAGCAACAGCTTATTGACAGCCTGCGCGAGGCGTTTGCCGATCGCCAGGTGGCGGGGCTGCTTTTTGGTAAAGTGCAGCGCGCCCGGGAACTGGGGCTAAAAATGATCATTGTCCCCGGTAGCCAGCTTTCGCAACTGCCGCCAGGa                                                                     <  2:215079/169‑1 (MQ=255)
                                                                                                                                                                                                                                  ggTGCCGGGGGAGTCTCAGGCCCAGCAACAGCTTATTGACAGCCTGCGCGAGGCGTTTGCCGATCGCCAGGTGGCGGGGCTGCTTTTTGGTAAAGTGCAGCGCGCCCGGGAACTGGGGCTAAAAATGATCATTGTCCCCGGTAGCCAGCTTTCGCAACTGCCGCCAGGACTGGATAGCACCGAGTTTGCAGCCATTGTGGGCCATTTACTTGGTAACGCCTTCCAAGCCAGGCTGc  >  2:341129/1‑236 (MQ=255)
                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                      
GTCTGAGCGTTATTGCCAACCGGGAAGCTATTCGTTCAGGTGATGATTTGCTGGGGGCCATTATCAGCTTTCGTAGTAAAGACGAAATTTCCACCCTCAATGCGCAACTGACGCAAATAAAACAATACGTTGAGAGCCTTCGTACATTGCGACACGAGCATCTCAATTGGATGTCGACGCTCAATGGTCTGTTGCAGATGAAAGAGTATGATCGCGTGCTGGCGATGGTGCAGGGGGAGTCTCAGGCCCAGCAACAGCTTATTGACAGCCTGCGCGAGGCGTTTGCCGATCGCCAGGTGGCGGGGCTGCTTTTTGGTAAAGTGCAGCGCGCCCGGGAACTGGGGCTAAAAATGATCATTGTCCCCGGTAGCCAGCTTTCGCAACTGCCGCCAGGACTGGATAGCACCGAGTTTGCAGCCATTGTGGGCAATTTACTTGATAACGCCTTCGAAGCCAGCCTGC  >  NC_000913/653121‑653582

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