Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,604,517 C→T L120L (CTG→TTG)  bglJ → DNA‑binding transcriptional regulator BglJ

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,604,5170CT98.0% 162.3 / NA 51L120L (CTG→TTG) bglJDNA‑binding transcriptional regulator BglJ
Reads supporting (aligned to +/- strand):  ref base C (0/0);  major base T (30/20);  minor base G (1/0);  total (31/20)
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

CCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGACTGAACTGGCGATTAAGTTTCCGCGTACCCGGCGTTTAGTTATTGCGGATGATGATATTGAAGCTCGGCTGATTGGTTCATTGTCGCCATCACCGCTGGACGGTGTATTAAGTAAAGCGTCAAC‑GCTGGAGATTTTTCATCAGGAACTCTTTTTGTCATTAAATGGTGTACGTCAGGCGACCGACCGACTGAACAATCAGTGGTACATTAACCAAAGCCGGACGTTAAGCCCGACGGAGAGAGAAATATTGCGCTTTATGTCGCGTGGCTACTCAATGACACAAATTGCCGAGCAGCTTAAACGCAATATCAAAACGATCCGTGCACATAAATTTAATGTGATGTCGAAACTGGGCGTCAGTTCTGACGCA  >  NC_000913/4604273‑4604762
                                                                                                                                                                                                                                                     |                                                                                                                                                                                                                                                     
ccATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTGCCCTGAGCAGTGAGCGCAGAGAGGGATTATCTTGCCTGACTGAACTGGCGATTAAGTTTTCGCGTACCCGGCGTTTTGTTATTGCGGGTGGTTATATTGAAGCTCGGCCGATTGGTTCATTGTCGCCATCACCGCTGGGCGGGGGATTAAGGAAAGCGTCAAC‑GTTgg                                                                                                                                                                                                                                                    >  1:338317/1‑248 (MQ=255)
ccATACCTATACCAGGCATGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGACTGAACTGGCGATTAAGTTTCCGCGTACCCGGCGTTTAGTTATTGCGGATGATGATATTGAAGCTCGGCTGATTGGTTCATTGTCGCCATCACCGCTGGACGGTGTATGAAGTAAAGCTTCAAC‑GTTGgag                                                                                                                                                                                                                                                  <  1:443435/250‑1 (MQ=255)
 cATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCCTTTTTTTTTTTTTTTCCCTGAGAAGTGAGCGCAGAGAGGGGTTATCTTTCCTGACTTAACTGGCGATTAAATTTCCGCGTACCCCGCGTTTTGTTATTTCGGGTGATGATTTTGAAGCTCGGCTGGTTGGTTTATTTTCGCCCTCCCCGCTGGGCGGGGGATTAAGTAAAGCGGCAAC‑GTTGGGGAt                                                                                                                                                                                                                                                >  1:318271/1‑251 (MQ=255)
                   gAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGGCATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCCGGGAGGGGTTATCTTGCCTGGCTGAACTGGGGCTTTAGTTTCCGCGTACCCGGCGTTTTGTTTTTGCGGATGGTGATATTTAAGCTCGGGTGATTGGTTCATTTTCGCCCTCCCCGCTGGACGGTGTATTAAGTAAAACGTCAAC‑GTTGGAGGTTTTTTATCCGGAAATCtt                                                                                                                                                                                                                              >  1:129571/1‑251 (MQ=255)
                       cATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTTTTTTTTTTGCATGGGAAGGGAGCGCAGAGAGGGGTTATTTTGCCTTACTGAACTGGCGGTTAAGGTTCCGCGTAACCGGCGGTTTGTTTTTGCGGGTGATGATATTTAAGCTCGGCTGAGTGGTTCATTTTCGCCATCCCCGCTGGGCGGGGGAATAAATAAAGCGTCAAC‑GGTGGAGATTTTTTATCCGGAAATCtttttt                                                                                                                                                                                                                          >  2:94062/1‑250 (MQ=255)
                        aTTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCCTTTTTTTTTTTTCTTCCATGGGAAGTGAGCGCCGAGAGGGGTTATCTTTCCTGACCGAACTGGCGATTAAATTTCCGCGGACCCGGCGTTTTGTTATTTCGGGTGATGATATTGAAGCTCGGCTGGTTGGTTTCTTGTCGCCATCACCGCTGGGAGGGGTATTTAGTAAAACGTCAACAGTTGG‑GGTTTTTCATAAGGGACTCttttttt                                                                                                                                                                                                                         >  1:342140/1‑249 (MQ=255)
                          ttCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGACTGAACTGGCGATTAAGTTTCCGCGTACCCGGCGTTTAGTTATTGCGGATGATGATATTGAAGCTCGGCTGATTGGTTCATTGTCGCCATCACCGCTGGACGGTGTATTAAGTAAAGCGTCAAC‑GTTGGAGATTTTTCATCAGGAACTCTTTTTGTc                                                                                                                                                                                                                        <  2:512899/250‑1 (MQ=255)
                                  gCCATGGCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGACTGAACTGGCGATTAAGTTGCCGCGTACCCGGCGTTTAGTTATTGCGGATGATGATATTGAAGCTCGGCTGATTGGTTCATTGTCTCCATCACCGCTGTACGGTGTATTAATTAAAGCGTCAAC‑GTTGGAGATTTTTCATCAGGAACTCTTTTGGTCATTAAATg                                                                                                                                                                                                                <  1:87970/250‑1 (MQ=255)
                                                       gcggcgGTCATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGACTGAACTGGCGATTAAGTTTCCGCGTACCCGGCGTTTAGTTATTGCGGATGATGATATTGAAGCTCGGCTGATTGGTTCATTGTCGCCATCACCGCTGGACGGTGTATTAAGTAAAGCGTCAAC‑GTTGGAGATTTTTCATCAGGAACTCTTTTTGTCATTAAATGGTGTACGTCAGGcgaccgacc                                                                                                                                                                                           <  1:27812/250‑1 (MQ=255)
                                                         ggcggTCATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGGTTATCTTGCCTGACTGAACTGGCGATTAAGTTTCCGCGTACCCGGCGTTTAGTTATTGCGGATGATGATATTGAAGCTCGGCTGATTGGTTCATTGTCGCCATCACCGCTGGACGGTGTATTAAGTAAAGCGTCAAC‑GTTGGAGATTTTTCATCAGGAACTCTTTTTGTCATTAAATGGTGTACGTCAGGcgaccgaccgac                                                                                                                                                                                        >  1:335086/1‑251 (MQ=255)
                                                                aTTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGACTGAACTGGCGATTAAGTTTCCGCGTACCCGGCGTTTAGTTATTGCGGATGATGATATTGAAGCTCGGCTGATTGGTTCATTGTCGCCATCACCGCTGGACGGTGTATTAAGTAAAGCGTCAAC‑GTTGGAGATTTTTCATCAGTAACTCTTTTTGTCATTAAATGGTGTACGTCAGGCGACCGACCGACTGAACa                                                                                                                                                                                  >  1:411592/1‑250 (MQ=255)
                                                                                                             tCTTGCCTGACTGAACTGGCGATTAAGTTTCCGCGTACCCGGCGTTTAGTTATTGCGGATGATGATATTGAAGCTCGGCTGCTTGGTTCATTGTCGCCCTCACCGCTGGACGGTGTATTAAGTAAAGCGTCAAC‑GTTGGAGATTTTTCATCAGGAACTCTTTTTGTCATTAAATGGTGTACGTCAGGCGACCGACCGACTGAACAATCAGTGGTACATTAACCAAAGCCGGACGTTAAGCCCGACGgag                                                                                                                                      >  2:269830/1‑249 (MQ=255)
                                                                                                                tGCCTGACTGAACTGGCGATTAAGTTTCCGCGTACCCGGCGTTTAGTTATTGCGGATGATGATATTGAAGCTCGGCTGATTGGTTCATTGTCGCCATCACCGCTGGACGGTGTATTAAGTAAAGCGTCAAC‑GTTGGAGATTTTTCATCAGGAACTCTTTTTGTCATTAAATGGTGTACGTCAGGCGACCGACCGACTGAACAATCAGTGGTACATTAACCAAAGCCGGACGTTAAGCCCGACGgagaga                                                                                                                                   <  1:397002/249‑1 (MQ=255)
                                                                                                                         gaACTGGCGATTAAGTTTCCGCGTACCCGGCGTTTAGTTATTGCGGATGATGATATTGAAGCTCGGCTGATTGGTTCATTGTCGCCATCACCGCTGGACGGTGTATTAAGTAAAGCGTCAAC‑GTTGGAGATTTTTCATCAGGAACTCTTTTTGTCATTAAATGGTGTACGTCAGGCGACCGACCGACTGAACAATCAGTGGTACATTAACCAAAGCCGGACGTTAAGCCCGACGGAGAGAGAAATATTg                                                                                                                          <  1:164363/249‑1 (MQ=255)
                                                                                                                            cTGGCGATTAAGTTTCCGCGTACCCGGCGTTTAGTTATTGCGGATGATGATATTGAAGCTCGGCTGATTGGTTCATTGTCGCCATCACCGCTGGACGGTGTATTAAGTAAAGCGTCAAC‑GTTGGAGATTTTTCATCAGGAACTCTTTTTGTCATTAAATGGTGTACGTCAGGCGACCGACCGCCTGAACAATCAGTGGTACCTTAACCAAAGCCGGACGTTCAGCCCGACGGCGAGCGAAATATTGCGCtt                                                                                                                     >  2:175236/1‑251 (MQ=255)
                                                                                                                                   ttaggTTTCCGCGTACCCCGCGTTTCGTTCTTGCGGATGATGATATTGAAGCTCGGCTGATTGGTTCATTGTCGCCATCACCGCTGGACGGTGTATTAAGTAAAGCGTCAAC‑GTTGGAGATTTTTCATCAGGAACTCTTTTTGTCATTAAATGGAGTACGTCAGGCGACCGACCGACTGAACAATCAGTGGTACATTAACCAAAGCCGGACGTTAAGCCCGACGGAGCGAGAAATATTGCGCTTTATGTcg                                                                                                              >  2:101366/5‑251 (MQ=255)
                                                                                                                                   ttAAGTTTCCGCGTACCCGGCGTTTAGTTATTGCGGATGATGATATTGAAGCTCGGCTGATTGGTTCATTGTCGCCATCACCGCTGGACGGTGTATTAAGTAAAGCGTCAAC‑GTTGGAGATTTTTCATCAGGAACTCTTTTTGTCATTAAATGGTGTACGTCAGGCGACCGACCGACTGAACAATCAGTGGTACATTAACCAAAGCCGGACGTTAAGCCCGACGGAGAGCGAAATATTGCGCTTTATGTcg                                                                                                              >  1:345132/1‑251 (MQ=255)
                                                                                                                                             gcgTACCCGGCGTTTAGTTATTGCGGATGATGATATTGAAGCTCGGCTGATTGGTTCATTGTAGCCATCACCGCTGGACGGTGTATTAAGTAAAGCGTCAAC‑GTTGGAGATTTTTCATCAGGAACTCTTTTTGTCATTAAATGGTGTACGTCAGGCGACCGACCGACTGAACAATCAGTGGTACATTAACCAAAGCCGGACGTTAAGCCCGACGGCGAGAGAAATATTGCGCTTTATGTCGCGTGGCTACt                                                                                                    >  2:226516/1‑251 (MQ=255)
                                                                                                                                                  cccGGCGTTTAGTTATTGCGGATGATGATATTGAAGCTCGGCTGATTGGTTCATTGTCGCCATCACCGCTGGACGGTGTATTAAGTAAAGCGTCAAC‑GTTGGAGATTTTTCATCAGGAACTCTTTTTGTCATTAAATGGTGTACGTCAGGCGACCGACCGACTGAACAATCAGTGGTACATTAACCAAAGCCGGACGTTAAGCCCGACGGAGAGAGAAATATTGCGCTTTATGTCGCGTGGCTACTCAATg                                                                                               >  1:267108/1‑251 (MQ=255)
                                                                                                                                                  cccGGCGTTTAGTTATTGCGGATGATGATATTGAAGCTCGGCTGATTGGTTCATTGTCGCCATCACCGCTGGACGGTGTATTAAGTAAAGCGTCAAC‑GTTGGAGATTTTTCATCAGGAACTCTTTTTGTCATTAAATGGTGTACGTCAGGCGACCGACCGACTGAACAATCAGTGGTACATTAACCAAAGCCGGACGTTAAGCCCGACGGAGAGAGAAATATTGCGCTTTATGTCGCGTGGCTACTCAATg                                                                                               <  1:322605/251‑1 (MQ=255)
                                                                                                                                                      gCGTTTAGTGATTGCGGATGATGATATTGAAGCTCGGCTGATTGGTTCATTGTCGCCATCACCGCTGGACGGTGTATTAAGTAAATCGTCAAC‑GTTGGAGATTTTTCATCAGGAACTCTTTTTGTCATTAAATGGTGTACGTCAGGCGACCGACCGACTGAACAATCAGTGGTACATTAACCAAAGCCGGACGTTAAGCCCGACGGAGAGAGAAATATGGCGCTTTATGTCGCGTGTCTACTCAATGacac                                                                                           <  2:262495/251‑1 (MQ=255)
                                                                                                                                                         tttAGTTATTGCGGATGATGATATTGAAGCTCGGCTGATTGGTTCATTGTCGCCATCACCGCTGGACGGTGTATTAAGTAAAGCGTCAAC‑GTTGGAGATTTTTCATCAGGAACTCTTTTTGTCATTAAATGGTGTACGTCAGGCGACCGACCGACTGAACAATCAGTGGTACATTAACCAAAGCCGGACGTTAAGCCCGACGGAGAGAGAAATATTGCGCTTTATGTCGCGTGGCTACTCATTGACACaa                                                                                         >  2:318027/1‑250 (MQ=255)
                                                                                                                                                            aGTTATTGCGGATGATGATAGTGAAGCTCGGCTGATTGGTTCATTGTCGCCATCACCGCTGGACGGTGTATTAAGTAAAGCGTCAAC‑GTTGGAGATTTTTCATCAGGAACTCTTTTTGTCATTAAATGGTGTACGTCAGGCGACCGACCGACTGAACAATCAGTGGTACATTAACCAAAGCCGGACGTTAAGCCCGACGGAGAGAGAAATATTGCGCTTTATGTCGCGTGGCTAATCAATGACACAAAtt                                                                                      <  1:353531/250‑1 (MQ=255)
                                                                                                                                                             gTTATTGCGGATGATGATATTGAAGCTCGGCTGATTGGTTCATTGTCGCCATCACCGCTGGACGGTGTATTAAGTAAAGCGTCAAC‑GTTGGAGATTTTTCATCAGGAACTCTTTTTGTCATTAAATGGTGTACGTCAGGCGACCGACCGACTGAACAATCAGTGGTACATTAACCAAAGCCGGACGTTAAGCCCGACGGAGAGAGAAATATTGCGCTTTATGTCGCGTGGCTACTCAATGACACAAATTGc                                                                                    >  1:504776/1‑251 (MQ=255)
                                                                                                                                                                aTTGCGGATGATGATATTGAAGCTCGGCTGATTGGTTCATTGTCGCCATCACCGCTGGACGGTGTATTAAGTAAAGCGTCAAC‑GTTGGAGATTTTTCATCAGGAACTCTTTTTGTCATTAAATGGTGTACGTCAGGCGACCGACCGACTGAACAATCAGTGGTACATTAACCAAAGCCGGACGTTAAGCCCGACGGAGAGAGAAATATTGCGCTTTATGTCGCGTGGCTACTCAATGACACAAATTGCCg                                                                                  <  1:248523/250‑1 (MQ=255)
                                                                                                                                                                aTTGCGGATGATGAGATTGAAGCTCGGCTGATTGGTTCATTGTCGCCATCACCGATGGACGGTGTATTAAGTAAAGCGTCAAC‑GTTGGAGATTTTTCAGCAGGAACTCTTTTTGTCATTAAATGGTGTACGTCAGGCGACCGACCGACTGAACAATCTGTGGTACATTAACCAAAGCCGGACGTTAAGCCCGACGGAGAGAGAAATATTGCGCTTTATGTCGCGTGGCTACTCAATGACACAAATTGCCg                                                                                  <  1:243923/250‑1 (MQ=255)
                                                                                                                                                                  tGCGGATGATGATATTGAAGCTCGGCTGATTGGTTCATTGTCGCCATCACCGCTGGACGGTGTATTAAGTACAGCGTCCAC‑GTTGGAGATTTTTCATCAGGAACTCTTTTTGTCATTAAATGGTGTACGTCAGGCGACCGACCGACTGAACAATCAGTGGTACATTAACCAAAGCCGGAGGTTAAGCCCGACGGAGAGAGAAATATTGCGCTTTATGTCGCGTGGCGTCTCAATGAAACAAATTGCCgtgc                                                                               >  1:294107/1‑248 (MQ=255)
                                                                                                                                                                   gCGGATGATGATATTGAAGCTCGGCTGATTGGTTCATTGTCGCCATCACCGCTGGACGGTGTATTAAGTAAAGCGTCAAC‑GTTGGAGATTTTTCATCAGGAACTCTTTTTGTCATTAAATGGTGTACGTCAGGCGACCGACCGACTGAACAATCAGTGGTACATTAACCAAAGCCGGACGTTAAGCCCGACGGAGAGAGAAATATTGCGCTTTATGTCGCGTGGCTACTCAATGACACAAATTGCCGagc                                                                               >  1:225507/1‑250 (MQ=255)
                                                                                                                                                                       atgatgatATTGAAGCTCGGCTGATTGGTTCATTGTCGCCATCACCGCTGGACGGTGTATTAAGTAAAGCGTCAAC‑GTTGGAGATTTTTCATCAGGAACTCTTTTTGTCATTAAATGGTGTACGTCAGGCGACCGACCGACTGAACAATCAGTGGTACATTAACCAAAGTCGGACGTTAAGCCCGACGGAGAGAGAAATATTGCGCTTTATGTCGCGTGGCTACTCAATGACACAAATTGCCGAGCAGCtt                                                                          <  2:129571/251‑1 (MQ=255)
                                                                                                                                                                         gatgatATTGAAGCTCGGCTGATTGGTTCATTGTCGCCATCACCGCTGGACGGTGTATTAAGTAAAGCGTCAAC‑GTTGGAGATTTTTCATCAGGAACTCTTTTTGTCATTAAATGGTGTACGTCAGGCGACCGACCGACTGAACAATCAGTGGTACATTAACCAAAGCCGGACGTTAAGCCCGACGGAGAGAGAAATATTGCGCTTTATGTCGCGTGGCTACTCAATGACACAAATTGCCGAGCAGCTTa                                                                         >  2:243870/1‑250 (MQ=255)
                                                                                                                                                                              tatTGAAGCTCGGCTGATTGGTTCATTGTCGCCATCACCGCTGGACGGTGTATTAAGTAAAGCGTCAAC‑GTTGGAGATTTTTCATCAGGAACTCTTTTTGTCATTAAATGGTGTACGTCAGGCGACCGACCGACTGAACAATCAGTGGTACATTAACCAAAGCCGGACGTTAAGCCCGACGGAGAGAGAAATATTGCGCTTTATGTCGCGTGGCTACTCAATGACACAAATTGCCGagc                                                                               >  2:143435/1‑239 (MQ=255)
                                                                                                                                                                              tatTGAAGCTCGGCTGATTGGTTCATTGTCGCCATCACCGCTGGACGGTGTATTAAGTAAAGCGTCAAC‑GTTGGAGATTTTTCATCAGGAACTCTTTTTGTCATTAAATGGTGTACGTCAGGCGACCGACCGACTGAACAATCAGTGGTACATTAACCAAAGCCGGACGTTAAGCCCGACGGAGAGAGAAATATTGCGCTTTATGTCGCGTGGCTACTCAATGACACAAATTGCCGagc                                                                               <  1:143435/239‑1 (MQ=255)
                                                                                                                                                                               atTGAAGCTCGGCTGATTGGTTCATTGTCGCCATCACCGCTGGCCGGTGTATTAAGTAAAGCGTCAAC‑GTTGGAGATTTTTCATCAGGAACTCTTTTTGTCATTAAATGGTGTACGTCAGGCGACCGACCGACTGAACAATCAGTGGTACATTAACCAAAGCCGGACGTTAAGCCCGACGGCGAGAGACATATTGCGCTTTATGTCGCGTGGCTACTCAATGACACAAATTGCCGAGCAGCTTAAACGCa                                                                   >  1:315242/1‑250 (MQ=255)
                                                                                                                                                                               atTGAAGCTCGGCTGATTGGTTCATTGTCGCCATCACCGCTGGACGGTGTATTAAGTAAAGCGTCAAC‑GTTGGAGATTTTTCATCAGGAACTCTTTTTGTCATTAAATGGTGTACGTCAGGCGACCGACCGACTGAACAATCCGTGGTACATTACCCAAAGCCGGACGTTAAGCCCGACGGAGAGAGAAATATTGCGCTTTATGTCGCGTGGCTACTCAATGACACAAATTGCCGAGCAGCTTAAACGCa                                                                   >  2:179420/1‑250 (MQ=255)
                                                                                                                                                                                            tGATTGGTTCATTGTCGCCATCACCGCTGGACGGTGTATTAAGTAAAGCGTCAAC‑GTTGGAGATTTTTCATCAGGAACTCTTTTTGTCATTAAATGGTGTACGTCAGGCGACCGACCGACTGAACAATCAGTGGTACATTAACCAAAGCCGGACGTTAAGCCCGACGGAGAGAGAAATATTGCGCTTTATGTcg                                                                                                              >  2:199654/1‑194 (MQ=255)
                                                                                                                                                                                            tGATTGGTTCATTGTCGCCATCACCGCTGGACGGTGTATTAAGTAAAGCGTCAAC‑GTTGGAGATTTTTCATCAGGAACTCTTTTTGTCATTAAATGGTGTACGTCAGGCGACCGACCGACTGAACAATCAGTGGTACATTAACCAAAGCCGGACGTTAAGCCCGACGGAGAGAGAAATATTGCGCTTTATGTcg                                                                                                              <  1:199654/194‑1 (MQ=255)
                                                                                                                                                                                                           cGCCATCACCGCTGGACGGTGTATTAAGTAAAGCGTCAAC‑GTTGGAGATTTTTCATCAGGCACTCTTTTTGTCATTAAATGGTGTACGTCAGGCGACCGACCGACTGAACAATCAGTGGTACATTAACCAAAGCCGGACGTTAAGCCCGACGGAGAGAGAAATATTGCGCTTTATGTCGCGTGGCTACTCAATGACACAAATTGCCGAGCAGCTTAAACGCAATATCAAAACGATCCGTGCCCATAAAttt                                      >  2:434739/1‑251 (MQ=255)
                                                                                                                                                                                                           cGCCATCACCGCTGGACGGTGTATTAAGTAAAGCGTCAAC‑GTTGGAGATTTTTCATCAGGAACTCTTTTTGTCATTAAATGGTGTACGTCAGGCGACCGACCGACTGAACAATCAGTGGTACATTAACCAAAGCCGGACGTTAAGCCCGACGGAGAGAGAAATATTGCGCTTTATGTCGCGTGGCTACTCAATGACACAAATTGCCGAGCAGCTTAAACGCAATATCAAAACGATCCGAGCACATAAAttt                                      >  2:217417/1‑251 (MQ=255)
                                                                                                                                                                                                           cGCCATCACCGCTGGACGGTGTATTAAGTAAAGCGTCAAC‑GTTGGAGATTGTTCATCAGGAACTCTTTTTGTCATTAAATGGTGTACGTCAGGCGACCGACCGACTGAACAATCAGTGGTACATTAACCAAAGCCGGACGTTAAGCCCGACGGAGAGAGAAATATTGCGCTTTATGTCGCGTGGCTACTCAATGACACAAATTGCCGAGCAGCTTAAACGCAATATCAAAACGATCCGTGCACATAAAtt                                       <  1:226516/250‑1 (MQ=255)
                                                                                                                                                                                                               aTCAACGCTGGACGGTGTATTAAGTAAAGCGTCAAC‑GTTGGAGATTTTTCATCAGGAACTCTTTTTGTCATTAAATGGTGTACGTCAGGCGACCGACCGACTGACCAATCAGTGGTACATTAACCAAAGCCGGACGTTAAGCCCGACGGAGAGAGAAATATTGAGCTTTATGTCGCGTGGCTACTCAATGACACAAATTGCCGAGCAGCTTAAACGCAATATCAAAACACTCCGTGCACATAAATTTAAtg                                  >  2:498942/1‑251 (MQ=255)
                                                                                                                                                                                                                 cacagCTGGACGGAGTATTAAGTAAGGCGTCAAC‑GTTGGAGATTTTTCATCAGGAACTCTTTTTGTCATTCAATGGTGTACGTCAGGCGACCGACCGACTGAACAATCAGTGGTACATTAACCCAAGCCGGACGTTAAGCCCGACGGAGAGAGAAATATTGCGCTTTGTGTCGCGTGGCTACTCAATGACACAACTTGCCGAGCAGCTTAAACGCAATATCAAAACGATCCGTGCACATAAATTTACtgtg                                >  2:526644/5‑251 (MQ=255)
                                                                                                                                                                                                                     gCTGGACGGTGTATTAAGTAAAGCGTCAAC‑GTTGGAGATTTTTCATCAGGAACTCTTTTTGTCATTACATGGTGTACGTCAGGCGACCGACCGACTGAACAATCAGTGGTACATTAACCAAAGCCGGACGTTAAGCCCGACGGAGAGAGAAATATTGCGCTTTATGTCGCGTGGCTACTCAATGACACAAATTGCCGAGCAGCTTAAACGCAATATCAAAACGATCCGTGCACCTAAATTTAATGTGATGt                            >  1:249299/1‑251 (MQ=255)
                                                                                                                                                                                                                       tGGACGGTGTATTAAGTAAAGCGTCAAC‑GTTGGAGATTTTTCATCAGGAACTCTTTTTGTCATTAAATGGTGTACGTCAGGCGACCGACCGACTGAACAATCAGTGATACATTAACCAAAGCCGGACGTTAAGCCCGACGGAGAGAGAAATATTGCGCTTTATGTCGCGGGGCTACTCAATGACACAAATTGCCGAGCAGCTTAAACGCAATATCAAAACGATCCGTGCACATAAATTTAATGTGATGTCg                          >  2:156565/1‑251 (MQ=255)
                                                                                                                                                                                                                           cGGTGTATTAAGTAAAGCGTCAAC‑GTTGGAGATTTTTCATCAGGAACTCTTTTTGTCATTAAATGGTGTACGTCAGGCGACCGACTGACTGAACAATCAGTGGTACATTAACCAAAGCCGGACGTTAAGCCCGACGGAGAGAGAAATATTGCGCTTTATGTCGCGTGGCTACTCAATGACACAAATTGCCGAGCAGCTTAAACGCAATATCAAAACGATCCGTGCACATAAATTTAATGTGATGTCGAAAc                      >  1:552020/1‑251 (MQ=255)
                                                                                                                                                                                                                           cGGTGTATTAAGTAAAGCGTCAAC‑GTTGGAGATTTTTCATCAGGAACTCTTTTTGTCATTAAATGGTGTACGTCAGGCGACCGACCGACTGAACAATCAGTGGTACATTAACCAAAGCCGGACGTTAAGCCCGACGGAGAGAGAAATATTGCGCTTTATGTCGCGTGGCTACTCAATGACACAAATTGCCGAGCAGCTTAAACGCAATATCAAAACGATCCGTGCACATAAATTTAATGTGATGTCGAAAc                      <  2:225507/251‑1 (MQ=255)
                                                                                                                                                                                                                            ggTGTATTAAGTAAAGCGTCAAC‑GTTGGAGATTTTTCATCAGGAACTCTGTTTGTCATTAAATGGTGTACGTCAGGCGACCGACCGACTGAACAATCAGTGGTACATTAACCAAAGCCGGACGTTAAGCCCGACGGAGAGAGAAATATTGCGCTTTATTTCGCGTGGCTACTCAATGACACAAATTGCCGAGCAGCTTAAACGCAATATCAAAACGATCCGTGCACATAAATTGAATGTGAGGGCGAAACt                     <  2:294146/251‑1 (MQ=255)
                                                                                                                                                                                                                             gtgtATTAAGTAAAGCGTCAAC‑GTTGGAGATTTTTCATCAGGAACTCTTTTTGTCATTAAATGGTGTACGTCAGGCGACCGACCGACTGAACAATCAGTGGTACATTAACCAAAGcc                                                                                                                                                          <  1:326366/117‑1 (MQ=255)
                                                                                                                                                                                                                             gtgtATTAAGTAAAGCGTCAAC‑GTTGGAGATTTTTCATCAGGAACTCTTTTTGTCATTAAATGGTGTACGTCAGGCGACCGACCGACTGAACAATCAGTGGTACATTAACCAAAGcc                                                                                                                                                          >  2:326366/1‑117 (MQ=255)
                                                                                                                                                                                                                                           gCGTCAAC‑GTTGGAGATTTTTCATCAGGAACTCTTTTTGTCATTAAATGGTGTACGTCAGGCGACCGACCGACTGAACAATCAGTGGTACATTAACCAAAGCCGGACGTTAAGCCCGACGGAGAGAGAAATATTGCGCTTTATGTCGCGTGGCTACTCAATGACACAAATTGCCGAGCAGCTTAAACGCAATATCAAAACGATCCGTGCACATAAATTTAATGTGATGTCGAAACTGGGCGTCAGTTCTGa      <  1:243870/251‑1 (MQ=255)
                                                                                                                                                                                                                                              tCAAC‑GTTGGAGATTTTTCATCAGGAACTCTTTTTGTCATTAAATGGTGTACGTCAGGCGACCGACCGACTGAACAATCAGTGGTACATTAACCAAAGCCGGACGTTAAGCCCGACGGAGAGAGAAATATTGCGCTTTATGTCGCGTGGCTACTCAATGACACAAATTGCCGAGCAGCTTAAACGCACTATCAAAACGATCCGTGCACATAAATTTAATGTGATGTCGAAACTGGACGTCAGTTCTGACGc   >  1:359237/1‑251 (MQ=255)
                                                                                                                                                                                                                                               cAAC‑GTTGGAGATTTTTCATCAGGAACTCTTTTTGTCATTAAATGGTGTACGTCAGGCGACCGACCGACTGAACAATCAGTGGTACATTAACCAAAGCCGGACGTTAAGCCCGACGGAGAGAGAAATATTGCGCTTTATGTCGCGTGGCTACTCAATGACACAAATTGCCGAGCAGCTTAAACGCAATATCAAAACGATCCGTGCACATAAATTTAATGTGATGTCGAAACTGGGCGTCAGTTCTGACGCa  <  1:179420/251‑1 (MQ=255)
                                                                                                                                                                                                                                                     |                                                                                                                                                                                                                                                     
CCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGACTGAACTGGCGATTAAGTTTCCGCGTACCCGGCGTTTAGTTATTGCGGATGATGATATTGAAGCTCGGCTGATTGGTTCATTGTCGCCATCACCGCTGGACGGTGTATTAAGTAAAGCGTCAAC‑GCTGGAGATTTTTCATCAGGAACTCTTTTTGTCATTAAATGGTGTACGTCAGGCGACCGACCGACTGAACAATCAGTGGTACATTAACCAAAGCCGGACGTTAAGCCCGACGGAGAGAGAAATATTGCGCTTTATGTCGCGTGGCTACTCAATGACACAAATTGCCGAGCAGCTTAAACGCAATATCAAAACGATCCGTGCACATAAATTTAATGTGATGTCGAAACTGGGCGTCAGTTCTGACGCA  >  NC_000913/4604273‑4604762

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