Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,637,346 T→C A217A (GCT→GCC creC → sensory histidine kinase CreC

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,637,3460TC93.9% 135.9 / ‑0.3 49A217A (GCT→GCCcreCsensory histidine kinase CreC
Reads supporting (aligned to +/- strand):  ref base T (2/1);  new base C (20/26);  total (22/27)
Fisher's exact test for biased strand distribution p-value = 5.81e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.63e-01

TGCAAAATCCTGCCGATCCCGAAAGTTCTGTGATGTATGTTGCCGCACCGATTATGGACGGCTCGCGGCTTATTGGCGTTTTGAGCGTAGGCAAACCGAACGCGGCGATGGCTCCGGTCATTAAGCGTAGCGAGCGGCGAATTTTATGGGCCAGCGCCATTTTGTTGGGGATTGCACTGGTGATTGGCGCAGGCATGGTTTGGTGGATCAACCGCTCTATTGCCCGGCTCACTCGCTATGCTGATTCCGTCACTGACAATAAGCCCGTTCCTCTCCCCGATCTCGGTAGTAGCGAGTTGCGTAAACTCGCGCAGGCGCTGGAAAGTATGCGCGTGAAGCTGGAAGGGAAAAACTATATTGAGCAGTATGTTTACGCATTAACTCATGAGCTAAAAAGCCCACTGGCGGCGATTCGTGGAGCGGCGGAAATTTTACGCGAAGGTCCGCCGCCGGAAGTGGTGGCTCGTTTTACTGACAACATTCTGACG  >  NC_000913/4637105‑4637592
                                                                                                                                                                                                                                                  |                                                                                                                                                                                                                                                      
tGCAAAATCCTGCCGATCCCGAAAGATCGGTGATGTATGTTGCCGCGCCGATTATGGACGGCTCGCGGCTTATTGGCGTTTTGAGCGTAGGCAAACCGAACGCGGCGATGGCTCCGGTCATTAAGCGTAGCGAGCGGCGAATTTTATGGGCCAGCGCCATTTTGCTGGGGATTGCACTGGTGATTGGCGCAGGCATGGTTTGGTGGATCAACCGCTCCATTGCCCGGCTCACTCGCTATGCCGATTCCGTc                                                                                                                                                                                                                                               >  1:299609/1‑251 (MQ=255)
   aaaaTCCTGCCGATCCCGAAAGTTCGGTGATGTATGTTGCCGCGCCGATTATGGACGGCTCGCGGCTTATTGGCGTTTTGAGCGTAGGCAAACCGAACGCGGCGATGGCTCCGGTCATTAAGCGTAGCGAGCGGCGAATTTTATGGGCCAGCGCCATTTTGCTGGGGATTGCACTGGTGATTGGCGCAGGCATGGTTTGGTGGATCAACCGCTCCATTGCCCGGCTCACTCGCTATGCCGATTCCGTCACt                                                                                                                                                                                                                                            <  1:179178/251‑1 (MQ=255)
     aaTCCTGCCGATCCCGAAAGGTCGGTGATGTATGTTGCCGCGCCGATTATGGACGGCTCGCGGCTTATTGGCGTTTTGAGCGTAGGCAAACCGAACGCGGCGATGGCTCCGGTCATTAAGCGTAGCGAGCGGCGAATTTTATGGGCCAGCGCCATTTTGCTGGGGATTGCACTGGTGATTGGCGCAGGCATGGTTTTGTGGATCAACCGCTCCATTGCCCGGCCCACTCGCCATGCCGATTCCGTCACTg                                                                                                                                                                                                                                           >  2:67100/1‑250 (MQ=255)
                 cccGAAAGTTCGGTGATGTATGTTGCCGCGCCGATTATGGACGGCTCGCGGCTTATTGGCGTTTTGAGCGTAGGCAAACCGAACGCGGCGATGGCTCCGGTCATTAAGCGTAGCGAGCGGCGAATTTTATGGGCCAGCGCCATTTTGCTGGGGATTGCACTGGTGATTGGCGCAGGCATGGTTTGGTGGATCAACCGCTCCATTGCCCGGCTCACTCGCTATGCCGATTCCGTCACTGACAATAAGCCCg                                                                                                                                                                                                                               <  1:171011/250‑1 (MQ=255)
                               gATGTATGTTGCCGCGCCGATTATGGACGGCTCGCGGCTTATTGGCGTTTTGAGCGTAGGCAAACCGAACGCGGCGATGGCTCCGGTCATTAAGCGTAGCGAGCGGCGAATTTTATGGGCCAGCGCCATTTTGCTGGGGATTGCACTGGTGATTGGCGCAGGCATGGTTTGGTGGATCAACCGCTCCATTGCCCGGCTCACTCGCTATGCCGATTCCGTCACTGACAATAAGCCCGTTCctcttcctgaa                                                                                                                                                                                                                 <  1:135463/250‑8 (MQ=255)
                                          ccGCGCCGTTTATGGACGGCTTGCGGCTTATTGGCGGTTTGAGCGTAGGCAAACCGAACGCGGCGATGGCTCCGGTCAATAAGCGTAGCGAGCGGCGAATTTTAGGGGCCAGTGCCATTTTGCTTGGGATTGCACTGGTGATTGGCGCAGGCATGGTTTGGTGGATCAACCGCTCCATTGCCCGGCTCACTCGCTATGCCGATTCCGTCACTGACAATAAGCCCGTTCCTCTTCCTGAACTCGGCAGTAGc                                                                                                                                                                                                     <  2:40187/251‑1 (MQ=255)
                                           cgcgccGATTATGGACGGCTCGCGGCTTATTGGCGTTTTGAGCGTAGGCAAACCGAACGCGGCGATGGCTCCGGTCATTAAGCGTAGCGAGCGGCGAATTTTATGGGCCAGCGCCATTTTGCTGGGGATTGCACTGGTGATTGGCGCAGGCATGGTTTGGTGGATCAACCGCGCCATTGCCCGGCTCACTCGCTATGCCGATTCCGTCACTGACAATAAGCCCGTTGCTATTCCTGAACTCGGCAGTAGCg                                                                                                                                                                                                    >  1:72919/5‑251 (MQ=255)
                                             cccagattaagggcggctttcgGCTTACTGGCGTTTTGAGCTTAGGGTAACCGAACTCGTGGATGGTTCCGGTAAGTAAGCTGAGCTGGCGGCGAATTTTATGTGCCAGCGCCATTTTGCTGGGGATTGCCCTGGTGATTGGCGCAGGCATGGTTTTGTGGATCAACCGCTCCATTGCCCGGCTCACTCGCTATGCCGATTCCGTCACTGACAATAAGCCCGTTCCTCTTCCTGAACTCGGCAGTAGctag                                                                                                                                                                                                  <  2:3618/231‑4 (MQ=255)
                                               ccGATTATGGACGGCTCGCGGCTTATTGGCGTTTTGAGCGTAGGCAAACCGAACGCGGCGATGGCTCCGGTCATTAAGCGTAGCGAGCGGCGAATTTTATGGGCCAGCGCCATTTTGCTGGGGATTGCACTGGTGATTGGCGCAGGCATGGTTTGGTGGATCAACCGCTCCATTGCCCGGCTCACTCGCTATGCCGATTCCGTCACTGACAATAAGCCCGTTCCTCTTCCTGAACTCGGCAGTAGCGAGtt                                                                                                                                                                                                >  1:63590/1‑251 (MQ=255)
                                                        gACGGCTCGCGGCTTATTGGCGTTTTGAGCGTAGGCAAACCGAACGCGGCGATGGCTCCGGTCATTAAGCGTAGCGAGCGGCGAATTTTATGGGCCAGCGCCATTTTGCTGGGGATTGCACTGGTGATTGGCGCAGGCATGGTTTGGTGGATCAACCGCTCCATTGCCCGGCTCACTCGCTATGCCGATTCCGTCACTGACAATAAGCCCGTTCCTCTTCCTGAACTCGGCAGTAGCGAGTTGCGTaagct                                                                                                                                                                                       <  1:147685/251‑4 (MQ=255)
                                                               cgcgGCTTATTGGCGTTTTGAGCGTAGGCAAACCGAACGCGGCGATGGCTCCGGTCATTAAGCGTAGCGAGCGGCGAATTTTATGGGCCAGCGCCATTTTGCTGGGGATTGCACTGGTGATTGGCGCAGGCATGGTTTGGTGGATCAACCGCTCCATTGCCCGGCTCACTCGCTATGCCGATTCCGTCACTGACAATAAGCCCGTTCCTCTTCCTGAACTCGGCAGTAGCGAGTTGCGTAAGCTTGCGCAg                                                                                                                                                                                <  1:341196/251‑1 (MQ=255)
                                                                  ggCTTATTGGCGTTTTGAGCGGAGGCAAGCCGAACGCGGCGATGGTTCCGGGCATTAAGCGTAGCGGGCGGCGAATTTTATGGGCCAGCGCCATTTTGCTGGGGATTGCACTGGTGATTGGCGCAGGCATGGTTTGGTGGATCAACCGCTCCATTGCCCGGCTCACTCGCTATGCCGATTCCGTCACTGACAATAAGCCCGTTCCTCTTCCTGAACTCGGCAGTAGCGAGTTGCGTAAGCTTGCGCAGgcg                                                                                                                                                                             <  2:131645/251‑1 (MQ=255)
                                                                      tATTGGCGTTTTGAGCGTAGGCAAACCGAACGCGGCGATGGCTCCGGTCATTAAGCGTAGCGAGCGGCGAATTTTATGGGCCAGCGCCATTTTGCTGGGGATTGCACTGGTGATTGGCGCAGGCATGGTTTGGTGGATCAACCGCTCCATTGCCCGGCTCACTCGCTATGCCGATTCCGTCACTGACAATAAGCCCGTTCCTCTTCCTGAACTCGGCAGTAGCGAGTTGCGTAAGCTTGCGCAGGCGCt                                                                                                                                                                           <  1:322398/249‑1 (MQ=255)
                                                                            cgttttggggtgtgtcAAACCGACCGTGGCGATGGCTTCGGTGATTAGGTCGAGTGAGCGGCGAAGTTTATGGGCCGGCGCCCTTTTGCTGGGGATTTCACTGGTGTTTGGCGCAGTCATGTTTTTGTGTATCTACCGCTCCATTGCCCGGCTCTCTCGCTATGCCGATTCCGTCACTGACAATAAGCCCGTTCCTCTTCCTGAACTCGGCAGTAGCGAGTTGCGTAAGATTGCGAAGGCGCTTGAAAGTa                                                                                                                                                                   <  2:109948/236‑1 (MQ=255)
                                                                             gTTTTGAGCGTAGGCAAACCGAACGCGGCGATGGCTCCGGTCATTAAGCGTAGCGAGCGGCGAATTTTATGGGCCAGCGCCATTTTGCTGGGGATTGCACTGGTGATTGGCGCAGGCATGGTTTGGTGGATCAACCGCTCCATTGCCCGGCTCACTCGCTATGCCGATTCCGTCACTGACAATAAGCCCGTTCCTCTTCCTGAACTCGGCAGTAGCGAGTTGCGTAAGCTTGCGCAGGCGCTGGAAAGTAt                                                                                                                                                                  <  1:67100/251‑1 (MQ=255)
                                                                                     cGTAGGCAAACCGAACGCGGCGATGGCTCCGGTCATTAAGCGTAGCGAGCGGCGAATTTTATGGGCCAGCGCCATTTTGCTGGGGATTGCACTGGTGATTGGCGCAGGCATGGTTTGGTGGATCAACCGCTCCATTGCCCGGCTCACTCGCTATGCCGATTCCGTCACTGACAATAAGCCCGTTCCTCTTCCTGAACTCGGCAGTAGCGAGTTGCGTAAGCTTGCGCAGGCGCTGGAAAGTATGCGCGTg                                                                                                                                                           >  1:173717/1‑250 (MQ=255)
                                                                                      gTAGGCAAACTGAACGCGGCGATGGCTCCGGTCATTAAGCGTAGCGAGCGGCGAATTTTATGGGCCAGCGCCATTTTGTTGGGGATTGCACTGGTGATTGGCGCAGGCATGGTTTGGTGGATCAACCGCTCTATTGCCCGGCTCACTCGCTATGCTGATTCCGTCACTGACAATAAGCCCGCTCCTCTCCCCGATCTCGGTAGTAGCGAGTTGCGTAAACTCGCGCAGGCGCTGGAAAGTATGAGCGTGa                                                                                                                                                          >  1:245673/1‑250 (MQ=255)
                                                                                            aaaCCGAACGCGGCGATGGCTCCGGTCATTAAGCGTAGCGAGCGGCGACTTTTATGGGCCAGCGCCATTTTGCTGGGGATTGCACTGGTGCTTGGCGCAGGCATGGTTTGGTGGATCAACCGCTCCATTGCCCGGCTCACTCGCTATGCCGATTCCGTCACTGAAAATAAGCCCGTTCCTCTTCCTGAACTCGGCAGTAGCGAGTTGCGTAAGCTTGCGCAGGCTCTGGTAAGTATGCGCGTGAAGATgg                                                                                                                                                    >  1:142258/1‑250 (MQ=255)
                                                                                                   aCGCGGCGATGGCTCCGGTCATTAAGCGTAGCGAGCGGCGAATTTTATGGGCCAGCGCCATTTTGCTGGGGATTGCACTGGTGATTGGCGCAGGCATGGTTTGGTGGATCAACCGCTCCATTGCCCGGCTCACTCGCTATGCCGATTCCGTCACTGACAATAAGCCCGTTCCTCTTCCTGAACTCGGCAGTAGCGAGTTGCGTAAGCTTGCGCAGGCGCTGGAAAGTATGCGCGTGAAGCTGGAAGGGaaa                                                                                                                                            <  1:336111/251‑1 (MQ=255)
                                                                                                                       aTTAAGCGTAGCGAGCGGCGAATTTTATGGGCCAGCGCCATTTTGCTGGGGCTTGCACTGGTGATTGGCGCAGGCATGGTTTGGTGGATCAACCGCTCCATTGCCCGGCTCACTCGCTATGCCGATTCCGTCACTGACAATAAGCCCGTTCCTCTTCCTGAACTCGGCAGTAGCGAGTTGCGTAAGCTTGCGCAGGCGCTGGAAAGTATGCGCGTGAAGCTGGAAGGGAAAAACTATATTGAGCAGTATGt                                                                                                                        >  1:306406/1‑251 (MQ=255)
                                                                                                                         gaaGCGGAGCGAGCGGCGAATTTTATGGGCAAGCGACATTTTGTTGGGGATTGCACTGGTGATTGGCGCAGGCATGGGTTGGTGGATCAACCGCTCTATTGCCCGGCTCACTCGCTATGCTGATTCCGTCACTGACAATAAGCCCGTTCCTCTCCCCGATCTCGGTAGTAGCGAGTTGCGTAAACTAGCGCAGGCGCTGGAAAGTATGCGCGTGAAGCTGGAAGGGAAAAACTATATTTAGCAGTATGtt                                                                                                                       <  2:9806/249‑1 (MQ=255)
                                                                                                                               gTAGCGAGCGGCGAATTTTATGGGCCAGCGCCATTTTGCTGGGGATTGCACTGGTGATTGGCGCAGGCATGGTTTGGTGGATCAACCGCTCCATTGCCCGGCTCACTCGCTATGCCGATTCCGTCACTGACAATAAGCCCGTTCCTCTTCCTGAACTCGGCAGTAGCGAGTTGCGTAAGCTTGCGCAGGCGCTGGAAAGTATGCGCGTGAAGCTGGAAGGGAAAAACTATATTGAGCAGTATGTTTatgcg                                                                                                                 >  1:37153/1‑247 (MQ=255)
                                                                                                                                  gcgagcgGCGAATTTTATGGGCCAGCGCCATTTTGCTGGGGATTGCACTGGTGATTGGCGCAGGCATGGTTTGGTGGATCAACCGCTCCATTGCCCGGCTCACTCGCTATGCCGATTCCGTCACTGACAATAAGCACGTTCCTCTTCCTGAACTCGGCAGTAGCGAGTTGCGTAAGCTTGCGCAGGCGCTGGAAAGTATGCGCGTGAAGCTGGAAGGGAAAAACTATATTGAGCAGTATGtttctgcgtt                                                                                                               >  1:54334/1‑243 (MQ=255)
                                                                                                                                    gagcgGCGAATTTTATGGGCCAGCGCCATTTGGCTGGGGATTGCACTGGTGATTGGCGCAGGCATGGTTTGGTGGATCAACCTCTCCATTGCCCGGCTCACTCGCTATGCCGATTCCGTCACTGACAATAAGCCCGTTCCTCTTCCTGAACTCGGCAGTAGCGAGTTGCGTAAGCTTGCTCAGGCGCTGGAAAGTATGCGCGTGAAGCTGGAAGGGAAAAACTATATTGAGCAGTATGTTTATGCGTTAAc                                                                                                            <  2:354692/251‑1 (MQ=255)
                                                                                                                                          cgAATTTTATGGGCCAGCGCCATTTTGTTGGGGATTGCACTGGTGATTGGCGCAGGCATGGTTTGGTGGATCAACCGCTCTATTGCCCGGCTCACTCGCTATGCTGATTCCGTCACTGACAATAAGCCCGTTCCTCTCCCCGATCTCGGTAGTAGCGAGTTGCGTAAACTCGCGCAGGCGCTGGAAAGTATGCGCGTGAAGCTGGAAGGGAAAAACTATATTGAGCAGTATGTTTACGCATTAACTCATg                                                                                                       >  1:198925/1‑250 (MQ=255)
                                                                                                                                          cgAATTTTATGGGCCAGCGCCATTTTGCTGGGGATTGCACTGGTGATTGGCGCAGGCATGGTTTGGTGGATCAACCGCTCCATTGCCCGGCTCACTCGCTATGCCGATTCCGTCACTGACAATAAGCCCGTTCCTCTTCCTGAACTCGGCAGTAGCGAGTTGCGTAAGCTTGCGCAGGCGCTGGAAAGTATGCGCGTGAAGCTGGAAGGGAAAAACTATATTGAGCAGTATGTTTATGCGTTAACCCATg                                                                                                       >  1:114885/1‑250 (MQ=255)
                                                                                                                                                         atcgcctgtttGCTGGGGATTGCACTGGTGATTGGCGCAGGCATGGTTTGGTGGATCAACCGCTCCATTGCCCGGCTCACGCGCTATGCCGATTCCGTCACTGACAATAAGCCCGTTCCTCTTCCTGAACTCGGCAGTAGCGAGTTGCGTAAGCTGGCGCAGGCGCTGGAAAGTATGCGCGTGAAGCTGGAAGGGAAAAACGATAGTGAGCAGTATGGTTATGCGTTAACCCATGAGCTAAAAAGCCCACt                                                                                       <  2:9928/243‑1 (MQ=255)
                                                                                                                                                          gcgcCATTTTGCTGGGGATTGCACTGGTGATTGGCGCAGGCATGGTTTGGTGGATCAACCGCTCCATTGCCCGGCTCACTCGCTATGCCGATTCCGTCACTGACAATAAGCCCGTTCCTCTTCCTGAACTCGGCAGTAGCGAGTTGCGTAAGCTTGCGCAGGCGCTGGAAAGTATGCGCGTGAAGCTGGAAGGGAAAAACTATATTGAGCAGTATGTTTATGCGTTAACCCATGAGCTAAAAAGCCCACTg                                                                                      <  1:205040/251‑1 (MQ=255)
                                                                                                                                                             ccATTTTGCTGGGGATTGCACTGGTGATTGGCGCAGGCATGGTTTGGTGGATCAACCGCTCCATTGCCCGGCTCACTCGCTATGCCGATTCCGTCACTGACAATAAGCCCGTTCCTCTTCCTGAACTCGGCAGTAGCGAGTTGCGTAAGCTTGCGCAGGCGCTGGAAAGTATGCGCGTGAAGCTGGAAGGGAAAAACTATATTGAGCAGTATGTTTATGCGTTAACCCATGAGCTAAAAAGCCCACTggcg                                                                                   >  1:355470/1‑251 (MQ=255)
                                                                                                                                                                          gATTGCACTGGTGATTGGCGCAGGCATGGTTTGGTGGATCAACCGCTCCATTGCCCGGCTCACTCGCTATGCCGATTCCGTCACTGACAATAAGCCCGTTCCTCTTCCTGAACTCGGCAGTAGCGAGTTGCGTAAGCTTGCGCAGGCGCTGGAAAGTATGCGCGTGAAGCTGGAAGGGAAAAACTATATTGAGCAGTATGTTTATGCGTTAACCCATGAGCTAAAAAGCCCACTGGCGGCGCTTCGTgg                                                                        >  1:168179/1‑249 (MQ=255)
                                                                                                                                                                                     tGATTGGCGCAGGCATGGTTTGGTGGATCAACCGCTCCATTGCCCGGCTCACTCGCTATGCCGATTCCGTCACTGACAATAAGCCCGTTCCTCTTCCTGAACTCGGCAGTAGCGAGTTGCGTAAGCTTGCGCAGGCGCTGGAAAGTATGCGCGTGAAGCTGGAAGGGAAAAACTATATTGAGCAGTATGTTTATGCGTTAACCCATGAGCTAAAAAGCCCACTGGCGGCGATTCGTGGCGCGGCGGAAAtt                                                           >  1:120724/1‑251 (MQ=255)
                                                                                                                                                                                       aTTGGCGCAGGCATGGTTTGGTGGATCAACCGCTCCATTGCCCGGCTCACTCGCTATGCCGATTCCGTCACTGACAATAAGCCCGTTCCTCTTCCTGAACTCGGCAGTAGCGAGTTGCGTAAGCTTGCGCAGGCGCTGGAAAGTATGCGCGTGAAGCTGGAAGGGAAAAACTATATTGAGCAGTATGTTTATGCGTTAACCCATGAGCTAAAAAGCCCACTGGCGGCGATTCGTGGCGCGGCGGAAAtttt                                                         <  1:173737/251‑1 (MQ=255)
                                                                                                                                                                                          ggCGCAGGCATGGTTTGGTGGATCAACCGCTCCATTGCCCGGCTCACTCGCTATGCCGATTCCGTCACTGACAATAAGCCCGTTCCTCTTCCTGAACTCGGCAGTAGCGAGTTGCGTAAGCTTGCGCAGGCGCTGGAAAGTATGCGCGTGAAGCTGGAAGGGAAAAACTATATTGAGCAGTATGTTTATGCGTTAACCCATGAGCTAAAAAGCCCACTGGCGGCGATTCGTGGCGCGGCGGAAATTTTAc                                                       <  1:152/250‑1 (MQ=255)
                                                                                                                                                                                            cgcAGGCATGGTTTGGTGGATCAACCGCTCCATTGCCCGGCTCACTCGCTATGCCGATTCCGTCACTGACAATAAGCCCGTTCCTCTTCCTGAACTCGGCAGTAGCGAGTTGCGTAAGCTTGCGCAGGCGCTGGAAAGTATGCGCGTGAAGCTGGAAGGGAAAAACTATATTGAGCAGTATGTTTATGCGTTAACCCATGAGCTAAAAAGCCCACTGGCGGCGATTCGTGGCGCGGCGGAAATTTTAcgcg                                                    >  1:188467/1‑251 (MQ=255)
                                                                                                                                                                                                ggCAGGGTTTGGTGGATCAACCGCTCCATTGCCCGGCTCACTCGCTAGGCCGAGTCCGTCACTGACAATAAGCCCGTTCCTCTTCCTGAACTCGGCAGTAGCGAGTTGCGTAAGCTTGCGCAGGCGCTGGAAAGTATGCGCGTGAAGCTGGAAGGGAAAAACTATATTGAGCAGTATGTTTATGCGTTAACCCATGAGCTAAAAAGCCCACTGGCGGCGATTCGTGGCGCGGCGGAAATTTTACGCGAAgg                                                <  2:188467/251‑1 (MQ=255)
                                                                                                                                                                                                      gTTTGGTGGATCAACCGCTCCATTGCCCGGCTCACTCGCTATGCCGATTCCGTCACTGACAATAAGCCCGTTCCTCTTCCTGAACTCGGCAGTAGCGAGTTGCGTAAGCTTGCGCAGGCGCTGGAAAGGATGCGCGTGAAGCTGGAAGGGAAAAACTATATTGAGCAGTATGTTTATGCGTTAACCCATGAGCTAAAAAGCCCACTGGCGGCGATTCGTGGCGCGGCGGAAATTTTACGCGAAGGTccgcc                                          >  1:315553/1‑251 (MQ=255)
                                                                                                                                                                                                        ttGGTGGATCAACCTCTCCATTGCGTGGCTCACTCGCTATGCCGATTCCGTCACTGTCAATAAGCCCGTTCCGCTTCCTGAACTCGGCAGTAGCGAGTTGCGTAAGCTTGCGCAGGCGCTGGAAAGTATGCGCGTGAAGCTGGAAGGGAAAAACTATATTGAGCAGTATGTTTATGCGTTAGCCCATGAGCTAAAAAGCCCCCTGGCGGCGATTCGTGGCGCGGCGGAAATTTTACGCGAAGGTccgcca                                         <  1:230411/250‑2 (MQ=255)
                                                                                                                                                                                                           gtggATCAACCGCTCCATTGCCCGGCTCACTCGCTATGCCGATTCCGTCACTGACAATAAGCCCGTTCCTCTTCCTGAACTCGGCAGTAGCGAGTTGCGTAAGCTTGCGCAGGCGCTGGAAAGTATGCGCGTGAAGCTGGAAGGGAAAAACTATATTGAGCAGTATGTTTATGCGTTAACCCATGAGCTAAAAAGCCCACTGGCGGCGATTCGTGGCGCGGCGGAAATTTTACGCGAAGGTCCGCCACCgg                                     <  1:366344/251‑1 (MQ=255)
                                                                                                                                                                                                              gATCAACCGCTCCATTGCCCGGCTCACTCGCTATGCCGATTCCGTCACTGACAATAAGCCCGTTCCTCTTCCTGAACTCGGCAGTAGCGAGTTGCGTAATCTTGCGCAGGCGCTGGAAAGTATGCGCGTGGAGCTGGAAGGGAAAAACTATATTGAGCAGTATGTTTatgcgttcccc                                                                                                           >  1:154379/1‑168 (MQ=255)
                                                                                                                                                                                                               aTCAACCGCTCCATTGCCCGGCTCACTCGCTATGCCGATTCCGTCACTGACAATAAGCCCGTTCCTCTTCCTGAACTCGGCAGTAGCGAGTTGCGTAAGCTTGCGCAGGCGCTGGAAAGTATGCGCGTGAAGCTGGAAGGGAAAAACTATATTGAGCAGTATGTTTATGCGTTAACCCATGAGCTAAAAAGCCCACTGGCGGCGATTCGTGGCGCGGCGGAAATTTTACGCGAAGGTCCGCCACCgggagg                                 >  1:228145/1‑247 (MQ=255)
                                                                                                                                                                                                                  aaCCGCTCCATTGCCCGGCTCACTCGCTATGCCGATTCCGTAACTGACAATAAGCCCGTTCCTCTTCCTGAACTCGGCAGTAGCGAGTTGCGTAAGCTTGCGCAGGCGCTGGAAAGTATGCGCGTGAAGCTGGAAGGGAAAAACTATATTGAGCAGTATGTTTATGCGTTAAACCATGAGCTAAAAAGCCCACTGGCGGCGATTCGTGGCGCGGCGGAAATTTTACGCGAAGGTCCGCCACCGGAAgtggt                              <  1:228749/251‑1 (MQ=255)
                                                                                                                                                                                                                     cGCTCCATTGCCCGGCTCACTCGCTATGCCGATTCCGTCACTGACAATAAGCCCGTTCCTCTTCCTGAACTCGGCAGTAGCGAGTTGCGTAAGCTTGCGCAGGCGCTGGAAAGAATGCGCGTGACGCTGGAAGGGAAAAACTATATTGAGAAGTATGTTTATGCGTTAACCCATGAGCTAAAAAGCCCACTGGCGGCGATTCGTGGCGCGGCGGAAATTTTACGCGAAGGTCCGCCACCGGAAGTGGTGGc                           >  1:370810/1‑251 (MQ=255)
                                                                                                                                                                                                                       cgccatggCCCGGCACCCTCGCTATGCCGTTTCCTTCACTGACAATAAGCCAGTTCCTCTTCCTGAACTCGGGAGGTGCGAGTTGCGGAAGCTTGCGCAGGCGCTGGAAAGTATGCGCGTGAAGCTGGAAGGGAAAAACTATATTGAGCAGTATGTGTATTCGTTAACCCATGAGCTAAAAAGCTCACTGGCGGCGATTCGTGGCGCGGCGGAAATTTTACGCGAAGGTCCGCCACCGGAAGTGGGGGCTc                         <  2:34605/244‑1 (MQ=255)
                                                                                                                                                                                                                           ggtGCCAGGCTCACTTGCTAGGCCGGTTCCCGCACTGACAATAAGCCCGTTCCTCTGCCTTAACTCGGCAGTAGCGAGTTGCGTAAGCTTACGCAGGCGCTGTAAAGTATGCGCGTGAAGCTGGAAGGGAAAAACTATATTGAGCAGTATGGTTATGCGTTAACCCATGAGCTAAAAAGCCCACTGGCGGCGATTCGTGGCGCGGCGGAAATTTTACGCGAAGGTCCGCCACCGGAAGTTGTGGCTCGttt                     <  2:167497/249‑1 (MQ=255)
                                                                                                                                                                                                                            aggctcGGCTCACGCGCTATGCCGATTCCGTCGCTGACAATAAGCCCGTTCCTTGTGCTGAACTCGGCAGTAGCGAGGTGCGTAAGCTTGCGCAGGCGCTGGAAAGTATGCGCGTGAAGCTGGAAGGGAAAAACTATATTGAGCAGTAGGTTTATGCGTTAACCCATGAGCTAAAAAGCCCACTGGCGGCGATTCTTGGCGCGGCGGAAATGTTACGCTAAGGTCCGCCACCTGAAGTGGTGGCTCGttt                     <  2:214204/245‑1 (MQ=255)
                                                                                                                                                                                                                              gCCCGGCTCACTCGCTATGCCGATTCCGTCACTGACAATAAGCCCGTTCCTCTTCCTGAACTCGGCAGTAGCGAGTTGCGTAAGCTTGCGCAGGCGCTGGAAAGTATGCGCGTGAAGCTGGAAGGGAAAAACTATATTGAGCAGTATGTTTATGCGTTAACCCATGAGCTAAAAAGCCCACTGGCGGCGATTCGTGGCGCGGCGGAAATTTTACGCGACGGTACGCCACCGGAAGTGGTGGCTCGttttcc                  >  2:122151/1‑249 (MQ=255)
                                                                                                                                                                                                                               cccggctgccccGCTATGCCGAGTCCGTCACTGACAATAAGCCCGTTCCTCTTCCTGAATTCGGCAGTAGCGAGTTGCGTAAGCTTGCGCAGGCGCGGGAAAGTATGCGCGTGAAGCTGGAAGGGAAAAACTATATTGAGCAGTATGTTTATGCGTTAACCCATGAGCTAAAAAGCCCACGGGCGGCGATTCGTGGCGCGGCGGAAATTTTACGCGAAGGGCCGCCACCTGAAGTGGTGGCTCGTTTTAcc                 <  2:345158/240‑2 (MQ=255)
                                                                                                                                                                                                                                      cACTCGCTATGCCGATTCCGTCACTGACAATAAGCCCGTTCCTCTTCCTGAACTCGGCAGTAGCGAGTTGCGTAAGCTTGCGCAGGCGCTGGAAAGTATGCGCGTGAAGCTGGAAGGGAAAAACTATATTGAGCAGTATGTTTATGCGTTAACCCATGAGCTAAAAAGCCCACTGGCGGCGATTCGTGGCGCGGCGGAAATTTTACGCGAAGGTCCGCCACCGGAAGTGGTGGCTCGTTTTACCGaccac           >  1:276643/1‑247 (MQ=255)
                                                                                                                                                                                                                                              aTGCCGATTCCGTCACTGACAATAAGCCCGTTCCTCTTCCTGAACTCGGCAGTAGCGAGTTGCGTAAGCTTGCGCAGGCGCTGGAAAGTATGCGCGTGAAGCTGGAAGGGAAAAACTATATTGAGCAGTATGTTTATGCGTTAACCCATGAGCTAAAAAGCCCACTGGCGGCGATTCGTGGCGCGGCGGAAATTTTACGCGAAGGTCCGCCACCGGAAGTGGTGGCTCGTTTTACCGACAACATTCTGACg  <  1:68967/251‑1 (MQ=255)
                                                                                                                                                                                                                                                  |                                                                                                                                                                                                                                                      
TGCAAAATCCTGCCGATCCCGAAAGTTCTGTGATGTATGTTGCCGCACCGATTATGGACGGCTCGCGGCTTATTGGCGTTTTGAGCGTAGGCAAACCGAACGCGGCGATGGCTCCGGTCATTAAGCGTAGCGAGCGGCGAATTTTATGGGCCAGCGCCATTTTGTTGGGGATTGCACTGGTGATTGGCGCAGGCATGGTTTGGTGGATCAACCGCTCTATTGCCCGGCTCACTCGCTATGCTGATTCCGTCACTGACAATAAGCCCGTTCCTCTCCCCGATCTCGGTAGTAGCGAGTTGCGTAAACTCGCGCAGGCGCTGGAAAGTATGCGCGTGAAGCTGGAAGGGAAAAACTATATTGAGCAGTATGTTTACGCATTAACTCATGAGCTAAAAAGCCCACTGGCGGCGATTCGTGGAGCGGCGGAAATTTTACGCGAAGGTCCGCCGCCGGAAGTGGTGGCTCGTTTTACTGACAACATTCTGACG  >  NC_000913/4637105‑4637592

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