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

CCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGACTGAACTGGCGATTAAGTTTCCGCGTACCCGGCGTTTAGTTATTGCGGATGATGATATTGAAGCTCGGCTGATTGGTTCATTGTCGCCATCACCGCTGGACGGTGTATTAAGTAAAGCGTCAACGCTGGAGATTTTTCATCAGGAACTCTTTTTGTCATTAAATGGTGTACGTCAGGCGACCGACCGACTGAACAATCAGTGGTACATTAACCAAAGCCGGACGTTAAGCCCGACGGAGAGAGAAATATTGCGCTTTATGTCGCGTGGCTACTCAATGACACAAATTGCCGAGCAGCTTAAACGCAATATCAAAACGATCCGTGCACATAAATTTAATGTGATGTCGAAACTGGGCGTCAGTTCTGACGCA  >  NC_000913/4604273‑4604762
                                                                                                                                                                                                                                                     |                                                                                                                                                                                                                                                     
ccATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTGCCCTGAGCAGTGAGCGCAGAGAGGGATTATCTTGCCTGACTGAACTGGCGATTAAGTTTTCGCGTACCCGGCGTTTTGTTATTGCGGGTGGTTATATTGAAGCTCGGCCGATTGGTTCATTGTCGCCATCACCGCTGGGCGGGGGATTAAGGAAAGCGTCAACGTTgg                                                                                                                                                                                                                                                    >  1:338317/1‑248 (MQ=255)
ccATACCTATACCAGGCATGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGACTGAACTGGCGATTAAGTTTCCGCGTACCCGGCGTTTAGTTATTGCGGATGATGATATTGAAGCTCGGCTGATTGGTTCATTGTCGCCATCACCGCTGGACGGTGTATGAAGTAAAGCTTCAACGTTGgag                                                                                                                                                                                                                                                  <  1:443435/250‑1 (MQ=255)
 cATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCCTTTTTTTTTTTTTTTCCCTGAGAAGTGAGCGCAGAGAGGGGTTATCTTTCCTGACTTAACTGGCGATTAAATTTCCGCGTACCCCGCGTTTTGTTATTTCGGGTGATGATTTTGAAGCTCGGCTGGTTGGTTTATTTTCGCCCTCCCCGCTGGGCGGGGGATTAAGTAAAGCGGCAACGTTGGGGAt                                                                                                                                                                                                                                                >  1:318271/1‑251 (MQ=255)
                   gAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGGCATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCCGGGAGGGGTTATCTTGCCTGGCTGAACTGGGGCTTTAGTTTCCGCGTACCCGGCGTTTTGTTTTTGCGGATGGTGATATTTAAGCTCGGGTGATTGGTTCATTTTCGCCCTCCCCGCTGGACGGTGTATTAAGTAAAACGTCAACGTTGGAGGTTTTTTATCCGGAAATCtt                                                                                                                                                                                                                              >  1:129571/1‑251 (MQ=255)
                       cATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTTTTTTTTTTGCATGGGAAGGGAGCGCAGAGAGGGGTTATTTTGCCTTACTGAACTGGCGGTTAAGGTTCCGCGTAACCGGCGGTTTGTTTTTGCGGGTGATGATATTTAAGCTCGGCTGAGTGGTTCATTTTCGCCATCCCCGCTGGGCGGGGGAATAAATAAAGCGTCAACGGTGGAGATTTTTTATCCGGAAATCtttttt                                                                                                                                                                                                                          >  2:94062/1‑250 (MQ=255)
                        aTTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCCTTTTTTTTTTTTCTTCCATGGGAAGTGAGCGCCGAGAGGGGTTATCTTTCCTGACCGAACTGGCGATTAAATTTCCGCGGACCCGGCGTTTTGTTATTTCGGGTGATGATATTGAAGCTCGGCTGGTTGGTTTCTTGTCGCCATCACCGCTGGGAGGGGTATTTAGTAAAACGTCAACAGTTGGGGTTTTTCATAAGGGACTCttttttt                                                                                                                                                                                                                         >  1:342140/1‑249 (MQ=255)
                          ttCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGACTGAACTGGCGATTAAGTTTCCGCGTACCCGGCGTTTAGTTATTGCGGATGATGATATTGAAGCTCGGCTGATTGGTTCATTGTCGCCATCACCGCTGGACGGTGTATTAAGTAAAGCGTCAACGTTGGAGATTTTTCATCAGGAACTCTTTTTGTc                                                                                                                                                                                                                        <  2:512899/250‑1 (MQ=255)
                                  gCCATGGCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGACTGAACTGGCGATTAAGTTGCCGCGTACCCGGCGTTTAGTTATTGCGGATGATGATATTGAAGCTCGGCTGATTGGTTCATTGTCTCCATCACCGCTGTACGGTGTATTAATTAAAGCGTCAACGTTGGAGATTTTTCATCAGGAACTCTTTTGGTCATTAAATg                                                                                                                                                                                                                <  1:87970/250‑1 (MQ=255)
                                                       gcggcgGTCATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGACTGAACTGGCGATTAAGTTTCCGCGTACCCGGCGTTTAGTTATTGCGGATGATGATATTGAAGCTCGGCTGATTGGTTCATTGTCGCCATCACCGCTGGACGGTGTATTAAGTAAAGCGTCAACGTTGGAGATTTTTCATCAGGAACTCTTTTTGTCATTAAATGGTGTACGTCAGGcgaccgacc                                                                                                                                                                                           <  1:27812/250‑1 (MQ=255)
                                                         ggcggTCATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGGTTATCTTGCCTGACTGAACTGGCGATTAAGTTTCCGCGTACCCGGCGTTTAGTTATTGCGGATGATGATATTGAAGCTCGGCTGATTGGTTCATTGTCGCCATCACCGCTGGACGGTGTATTAAGTAAAGCGTCAACGTTGGAGATTTTTCATCAGGAACTCTTTTTGTCATTAAATGGTGTACGTCAGGcgaccgaccgac                                                                                                                                                                                        >  1:335086/1‑251 (MQ=255)
                                                                aTTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGACTGAACTGGCGATTAAGTTTCCGCGTACCCGGCGTTTAGTTATTGCGGATGATGATATTGAAGCTCGGCTGATTGGTTCATTGTCGCCATCACCGCTGGACGGTGTATTAAGTAAAGCGTCAACGTTGGAGATTTTTCATCAGTAACTCTTTTTGTCATTAAATGGTGTACGTCAGGCGACCGACCGACTGAACa                                                                                                                                                                                  >  1:411592/1‑250 (MQ=255)
                                                                                                             tCTTGCCTGACTGAACTGGCGATTAAGTTTCCGCGTACCCGGCGTTTAGTTATTGCGGATGATGATATTGAAGCTCGGCTGCTTGGTTCATTGTCGCCCTCACCGCTGGACGGTGTATTAAGTAAAGCGTCAACGTTGGAGATTTTTCATCAGGAACTCTTTTTGTCATTAAATGGTGTACGTCAGGCGACCGACCGACTGAACAATCAGTGGTACATTAACCAAAGCCGGACGTTAAGCCCGACGgag                                                                                                                                      >  2:269830/1‑249 (MQ=255)
                                                                                                                tGCCTGACTGAACTGGCGATTAAGTTTCCGCGTACCCGGCGTTTAGTTATTGCGGATGATGATATTGAAGCTCGGCTGATTGGTTCATTGTCGCCATCACCGCTGGACGGTGTATTAAGTAAAGCGTCAACGTTGGAGATTTTTCATCAGGAACTCTTTTTGTCATTAAATGGTGTACGTCAGGCGACCGACCGACTGAACAATCAGTGGTACATTAACCAAAGCCGGACGTTAAGCCCGACGgagaga                                                                                                                                   <  1:397002/249‑1 (MQ=255)
                                                                                                                         gaACTGGCGATTAAGTTTCCGCGTACCCGGCGTTTAGTTATTGCGGATGATGATATTGAAGCTCGGCTGATTGGTTCATTGTCGCCATCACCGCTGGACGGTGTATTAAGTAAAGCGTCAACGTTGGAGATTTTTCATCAGGAACTCTTTTTGTCATTAAATGGTGTACGTCAGGCGACCGACCGACTGAACAATCAGTGGTACATTAACCAAAGCCGGACGTTAAGCCCGACGGAGAGAGAAATATTg                                                                                                                          <  1:164363/249‑1 (MQ=255)
                                                                                                                            cTGGCGATTAAGTTTCCGCGTACCCGGCGTTTAGTTATTGCGGATGATGATATTGAAGCTCGGCTGATTGGTTCATTGTCGCCATCACCGCTGGACGGTGTATTAAGTAAAGCGTCAACGTTGGAGATTTTTCATCAGGAACTCTTTTTGTCATTAAATGGTGTACGTCAGGCGACCGACCGCCTGAACAATCAGTGGTACCTTAACCAAAGCCGGACGTTCAGCCCGACGGCGAGCGAAATATTGCGCtt                                                                                                                     >  2:175236/1‑251 (MQ=255)
                                                                                                                                   ttaggTTTCCGCGTACCCCGCGTTTCGTTCTTGCGGATGATGATATTGAAGCTCGGCTGATTGGTTCATTGTCGCCATCACCGCTGGACGGTGTATTAAGTAAAGCGTCAACGTTGGAGATTTTTCATCAGGAACTCTTTTTGTCATTAAATGGAGTACGTCAGGCGACCGACCGACTGAACAATCAGTGGTACATTAACCAAAGCCGGACGTTAAGCCCGACGGAGCGAGAAATATTGCGCTTTATGTcg                                                                                                              >  2:101366/5‑251 (MQ=255)
                                                                                                                                   ttAAGTTTCCGCGTACCCGGCGTTTAGTTATTGCGGATGATGATATTGAAGCTCGGCTGATTGGTTCATTGTCGCCATCACCGCTGGACGGTGTATTAAGTAAAGCGTCAACGTTGGAGATTTTTCATCAGGAACTCTTTTTGTCATTAAATGGTGTACGTCAGGCGACCGACCGACTGAACAATCAGTGGTACATTAACCAAAGCCGGACGTTAAGCCCGACGGAGAGCGAAATATTGCGCTTTATGTcg                                                                                                              >  1:345132/1‑251 (MQ=255)
                                                                                                                                             gcgTACCCGGCGTTTAGTTATTGCGGATGATGATATTGAAGCTCGGCTGATTGGTTCATTGTAGCCATCACCGCTGGACGGTGTATTAAGTAAAGCGTCAACGTTGGAGATTTTTCATCAGGAACTCTTTTTGTCATTAAATGGTGTACGTCAGGCGACCGACCGACTGAACAATCAGTGGTACATTAACCAAAGCCGGACGTTAAGCCCGACGGCGAGAGAAATATTGCGCTTTATGTCGCGTGGCTACt                                                                                                    >  2:226516/1‑251 (MQ=255)
                                                                                                                                                  cccGGCGTTTAGTTATTGCGGATGATGATATTGAAGCTCGGCTGATTGGTTCATTGTCGCCATCACCGCTGGACGGTGTATTAAGTAAAGCGTCAACGTTGGAGATTTTTCATCAGGAACTCTTTTTGTCATTAAATGGTGTACGTCAGGCGACCGACCGACTGAACAATCAGTGGTACATTAACCAAAGCCGGACGTTAAGCCCGACGGAGAGAGAAATATTGCGCTTTATGTCGCGTGGCTACTCAATg                                                                                               >  1:267108/1‑251 (MQ=255)
                                                                                                                                                  cccGGCGTTTAGTTATTGCGGATGATGATATTGAAGCTCGGCTGATTGGTTCATTGTCGCCATCACCGCTGGACGGTGTATTAAGTAAAGCGTCAACGTTGGAGATTTTTCATCAGGAACTCTTTTTGTCATTAAATGGTGTACGTCAGGCGACCGACCGACTGAACAATCAGTGGTACATTAACCAAAGCCGGACGTTAAGCCCGACGGAGAGAGAAATATTGCGCTTTATGTCGCGTGGCTACTCAATg                                                                                               <  1:322605/251‑1 (MQ=255)
                                                                                                                                                      gCGTTTAGTGATTGCGGATGATGATATTGAAGCTCGGCTGATTGGTTCATTGTCGCCATCACCGCTGGACGGTGTATTAAGTAAATCGTCAACGTTGGAGATTTTTCATCAGGAACTCTTTTTGTCATTAAATGGTGTACGTCAGGCGACCGACCGACTGAACAATCAGTGGTACATTAACCAAAGCCGGACGTTAAGCCCGACGGAGAGAGAAATATGGCGCTTTATGTCGCGTGTCTACTCAATGacac                                                                                           <  2:262495/251‑1 (MQ=255)
                                                                                                                                                         tttAGTTATTGCGGATGATGATATTGAAGCTCGGCTGATTGGTTCATTGTCGCCATCACCGCTGGACGGTGTATTAAGTAAAGCGTCAACGTTGGAGATTTTTCATCAGGAACTCTTTTTGTCATTAAATGGTGTACGTCAGGCGACCGACCGACTGAACAATCAGTGGTACATTAACCAAAGCCGGACGTTAAGCCCGACGGAGAGAGAAATATTGCGCTTTATGTCGCGTGGCTACTCATTGACACaa                                                                                         >  2:318027/1‑250 (MQ=255)
                                                                                                                                                            aGTTATTGCGGATGATGATAGTGAAGCTCGGCTGATTGGTTCATTGTCGCCATCACCGCTGGACGGTGTATTAAGTAAAGCGTCAACGTTGGAGATTTTTCATCAGGAACTCTTTTTGTCATTAAATGGTGTACGTCAGGCGACCGACCGACTGAACAATCAGTGGTACATTAACCAAAGCCGGACGTTAAGCCCGACGGAGAGAGAAATATTGCGCTTTATGTCGCGTGGCTAATCAATGACACAAAtt                                                                                      <  1:353531/250‑1 (MQ=255)
                                                                                                                                                             gTTATTGCGGATGATGATATTGAAGCTCGGCTGATTGGTTCATTGTCGCCATCACCGCTGGACGGTGTATTAAGTAAAGCGTCAACGTTGGAGATTTTTCATCAGGAACTCTTTTTGTCATTAAATGGTGTACGTCAGGCGACCGACCGACTGAACAATCAGTGGTACATTAACCAAAGCCGGACGTTAAGCCCGACGGAGAGAGAAATATTGCGCTTTATGTCGCGTGGCTACTCAATGACACAAATTGc                                                                                    >  1:504776/1‑251 (MQ=255)
                                                                                                                                                                aTTGCGGATGATGATATTGAAGCTCGGCTGATTGGTTCATTGTCGCCATCACCGCTGGACGGTGTATTAAGTAAAGCGTCAACGTTGGAGATTTTTCATCAGGAACTCTTTTTGTCATTAAATGGTGTACGTCAGGCGACCGACCGACTGAACAATCAGTGGTACATTAACCAAAGCCGGACGTTAAGCCCGACGGAGAGAGAAATATTGCGCTTTATGTCGCGTGGCTACTCAATGACACAAATTGCCg                                                                                  <  1:248523/250‑1 (MQ=255)
                                                                                                                                                                aTTGCGGATGATGAGATTGAAGCTCGGCTGATTGGTTCATTGTCGCCATCACCGATGGACGGTGTATTAAGTAAAGCGTCAACGTTGGAGATTTTTCAGCAGGAACTCTTTTTGTCATTAAATGGTGTACGTCAGGCGACCGACCGACTGAACAATCTGTGGTACATTAACCAAAGCCGGACGTTAAGCCCGACGGAGAGAGAAATATTGCGCTTTATGTCGCGTGGCTACTCAATGACACAAATTGCCg                                                                                  <  1:243923/250‑1 (MQ=255)
                                                                                                                                                                  tGCGGATGATGATATTGAAGCTCGGCTGATTGGTTCATTGTCGCCATCACCGCTGGACGGTGTATTAAGTACAGCGTCCACGTTGGAGATTTTTCATCAGGAACTCTTTTTGTCATTAAATGGTGTACGTCAGGCGACCGACCGACTGAACAATCAGTGGTACATTAACCAAAGCCGGAGGTTAAGCCCGACGGAGAGAGAAATATTGCGCTTTATGTCGCGTGGCGTCTCAATGAAACAAATTGCCgtgc                                                                               >  1:294107/1‑248 (MQ=255)
                                                                                                                                                                   gCGGATGATGATATTGAAGCTCGGCTGATTGGTTCATTGTCGCCATCACCGCTGGACGGTGTATTAAGTAAAGCGTCAACGTTGGAGATTTTTCATCAGGAACTCTTTTTGTCATTAAATGGTGTACGTCAGGCGACCGACCGACTGAACAATCAGTGGTACATTAACCAAAGCCGGACGTTAAGCCCGACGGAGAGAGAAATATTGCGCTTTATGTCGCGTGGCTACTCAATGACACAAATTGCCGagc                                                                               >  1:225507/1‑250 (MQ=255)
                                                                                                                                                                       atgatgatATTGAAGCTCGGCTGATTGGTTCATTGTCGCCATCACCGCTGGACGGTGTATTAAGTAAAGCGTCAACGTTGGAGATTTTTCATCAGGAACTCTTTTTGTCATTAAATGGTGTACGTCAGGCGACCGACCGACTGAACAATCAGTGGTACATTAACCAAAGTCGGACGTTAAGCCCGACGGAGAGAGAAATATTGCGCTTTATGTCGCGTGGCTACTCAATGACACAAATTGCCGAGCAGCtt                                                                          <  2:129571/251‑1 (MQ=255)
                                                                                                                                                                         gatgatATTGAAGCTCGGCTGATTGGTTCATTGTCGCCATCACCGCTGGACGGTGTATTAAGTAAAGCGTCAACGTTGGAGATTTTTCATCAGGAACTCTTTTTGTCATTAAATGGTGTACGTCAGGCGACCGACCGACTGAACAATCAGTGGTACATTAACCAAAGCCGGACGTTAAGCCCGACGGAGAGAGAAATATTGCGCTTTATGTCGCGTGGCTACTCAATGACACAAATTGCCGAGCAGCTTa                                                                         >  2:243870/1‑250 (MQ=255)
                                                                                                                                                                              tatTGAAGCTCGGCTGATTGGTTCATTGTCGCCATCACCGCTGGACGGTGTATTAAGTAAAGCGTCAACGTTGGAGATTTTTCATCAGGAACTCTTTTTGTCATTAAATGGTGTACGTCAGGCGACCGACCGACTGAACAATCAGTGGTACATTAACCAAAGCCGGACGTTAAGCCCGACGGAGAGAGAAATATTGCGCTTTATGTCGCGTGGCTACTCAATGACACAAATTGCCGagc                                                                               >  2:143435/1‑239 (MQ=255)
                                                                                                                                                                              tatTGAAGCTCGGCTGATTGGTTCATTGTCGCCATCACCGCTGGACGGTGTATTAAGTAAAGCGTCAACGTTGGAGATTTTTCATCAGGAACTCTTTTTGTCATTAAATGGTGTACGTCAGGCGACCGACCGACTGAACAATCAGTGGTACATTAACCAAAGCCGGACGTTAAGCCCGACGGAGAGAGAAATATTGCGCTTTATGTCGCGTGGCTACTCAATGACACAAATTGCCGagc                                                                               <  1:143435/239‑1 (MQ=255)
                                                                                                                                                                               atTGAAGCTCGGCTGATTGGTTCATTGTCGCCATCACCGCTGGCCGGTGTATTAAGTAAAGCGTCAACGTTGGAGATTTTTCATCAGGAACTCTTTTTGTCATTAAATGGTGTACGTCAGGCGACCGACCGACTGAACAATCAGTGGTACATTAACCAAAGCCGGACGTTAAGCCCGACGGCGAGAGACATATTGCGCTTTATGTCGCGTGGCTACTCAATGACACAAATTGCCGAGCAGCTTAAACGCa                                                                   >  1:315242/1‑250 (MQ=255)
                                                                                                                                                                               atTGAAGCTCGGCTGATTGGTTCATTGTCGCCATCACCGCTGGACGGTGTATTAAGTAAAGCGTCAACGTTGGAGATTTTTCATCAGGAACTCTTTTTGTCATTAAATGGTGTACGTCAGGCGACCGACCGACTGAACAATCCGTGGTACATTACCCAAAGCCGGACGTTAAGCCCGACGGAGAGAGAAATATTGCGCTTTATGTCGCGTGGCTACTCAATGACACAAATTGCCGAGCAGCTTAAACGCa                                                                   >  2:179420/1‑250 (MQ=255)
                                                                                                                                                                                            tGATTGGTTCATTGTCGCCATCACCGCTGGACGGTGTATTAAGTAAAGCGTCAACGTTGGAGATTTTTCATCAGGAACTCTTTTTGTCATTAAATGGTGTACGTCAGGCGACCGACCGACTGAACAATCAGTGGTACATTAACCAAAGCCGGACGTTAAGCCCGACGGAGAGAGAAATATTGCGCTTTATGTcg                                                                                                              >  2:199654/1‑194 (MQ=255)
                                                                                                                                                                                            tGATTGGTTCATTGTCGCCATCACCGCTGGACGGTGTATTAAGTAAAGCGTCAACGTTGGAGATTTTTCATCAGGAACTCTTTTTGTCATTAAATGGTGTACGTCAGGCGACCGACCGACTGAACAATCAGTGGTACATTAACCAAAGCCGGACGTTAAGCCCGACGGAGAGAGAAATATTGCGCTTTATGTcg                                                                                                              <  1:199654/194‑1 (MQ=255)
                                                                                                                                                                                                           cGCCATCACCGCTGGACGGTGTATTAAGTAAAGCGTCAACGTTGGAGATTTTTCATCAGGCACTCTTTTTGTCATTAAATGGTGTACGTCAGGCGACCGACCGACTGAACAATCAGTGGTACATTAACCAAAGCCGGACGTTAAGCCCGACGGAGAGAGAAATATTGCGCTTTATGTCGCGTGGCTACTCAATGACACAAATTGCCGAGCAGCTTAAACGCAATATCAAAACGATCCGTGCCCATAAAttt                                      >  2:434739/1‑251 (MQ=255)
                                                                                                                                                                                                           cGCCATCACCGCTGGACGGTGTATTAAGTAAAGCGTCAACGTTGGAGATTTTTCATCAGGAACTCTTTTTGTCATTAAATGGTGTACGTCAGGCGACCGACCGACTGAACAATCAGTGGTACATTAACCAAAGCCGGACGTTAAGCCCGACGGAGAGAGAAATATTGCGCTTTATGTCGCGTGGCTACTCAATGACACAAATTGCCGAGCAGCTTAAACGCAATATCAAAACGATCCGAGCACATAAAttt                                      >  2:217417/1‑251 (MQ=255)
                                                                                                                                                                                                           cGCCATCACCGCTGGACGGTGTATTAAGTAAAGCGTCAACGTTGGAGATTGTTCATCAGGAACTCTTTTTGTCATTAAATGGTGTACGTCAGGCGACCGACCGACTGAACAATCAGTGGTACATTAACCAAAGCCGGACGTTAAGCCCGACGGAGAGAGAAATATTGCGCTTTATGTCGCGTGGCTACTCAATGACACAAATTGCCGAGCAGCTTAAACGCAATATCAAAACGATCCGTGCACATAAAtt                                       <  1:226516/250‑1 (MQ=255)
                                                                                                                                                                                                               aTCAACGCTGGACGGTGTATTAAGTAAAGCGTCAACGTTGGAGATTTTTCATCAGGAACTCTTTTTGTCATTAAATGGTGTACGTCAGGCGACCGACCGACTGACCAATCAGTGGTACATTAACCAAAGCCGGACGTTAAGCCCGACGGAGAGAGAAATATTGAGCTTTATGTCGCGTGGCTACTCAATGACACAAATTGCCGAGCAGCTTAAACGCAATATCAAAACACTCCGTGCACATAAATTTAAtg                                  >  2:498942/1‑251 (MQ=255)
                                                                                                                                                                                                                 cacagCTGGACGGAGTATTAAGTAAGGCGTCAACGTTGGAGATTTTTCATCAGGAACTCTTTTTGTCATTCAATGGTGTACGTCAGGCGACCGACCGACTGAACAATCAGTGGTACATTAACCCAAGCCGGACGTTAAGCCCGACGGAGAGAGAAATATTGCGCTTTGTGTCGCGTGGCTACTCAATGACACAACTTGCCGAGCAGCTTAAACGCAATATCAAAACGATCCGTGCACATAAATTTACtgtg                                >  2:526644/5‑251 (MQ=255)
                                                                                                                                                                                                                     gCTGGACGGTGTATTAAGTAAAGCGTCAACGTTGGAGATTTTTCATCAGGAACTCTTTTTGTCATTACATGGTGTACGTCAGGCGACCGACCGACTGAACAATCAGTGGTACATTAACCAAAGCCGGACGTTAAGCCCGACGGAGAGAGAAATATTGCGCTTTATGTCGCGTGGCTACTCAATGACACAAATTGCCGAGCAGCTTAAACGCAATATCAAAACGATCCGTGCACCTAAATTTAATGTGATGt                            >  1:249299/1‑251 (MQ=255)
                                                                                                                                                                                                                       tGGACGGTGTATTAAGTAAAGCGTCAACGTTGGAGATTTTTCATCAGGAACTCTTTTTGTCATTAAATGGTGTACGTCAGGCGACCGACCGACTGAACAATCAGTGATACATTAACCAAAGCCGGACGTTAAGCCCGACGGAGAGAGAAATATTGCGCTTTATGTCGCGGGGCTACTCAATGACACAAATTGCCGAGCAGCTTAAACGCAATATCAAAACGATCCGTGCACATAAATTTAATGTGATGTCg                          >  2:156565/1‑251 (MQ=255)
                                                                                                                                                                                                                           cGGTGTATTAAGTAAAGCGTCAACGTTGGAGATTTTTCATCAGGAACTCTTTTTGTCATTAAATGGTGTACGTCAGGCGACCGACTGACTGAACAATCAGTGGTACATTAACCAAAGCCGGACGTTAAGCCCGACGGAGAGAGAAATATTGCGCTTTATGTCGCGTGGCTACTCAATGACACAAATTGCCGAGCAGCTTAAACGCAATATCAAAACGATCCGTGCACATAAATTTAATGTGATGTCGAAAc                      >  1:552020/1‑251 (MQ=255)
                                                                                                                                                                                                                           cGGTGTATTAAGTAAAGCGTCAACGTTGGAGATTTTTCATCAGGAACTCTTTTTGTCATTAAATGGTGTACGTCAGGCGACCGACCGACTGAACAATCAGTGGTACATTAACCAAAGCCGGACGTTAAGCCCGACGGAGAGAGAAATATTGCGCTTTATGTCGCGTGGCTACTCAATGACACAAATTGCCGAGCAGCTTAAACGCAATATCAAAACGATCCGTGCACATAAATTTAATGTGATGTCGAAAc                      <  2:225507/251‑1 (MQ=255)
                                                                                                                                                                                                                            ggTGTATTAAGTAAAGCGTCAACGTTGGAGATTTTTCATCAGGAACTCTGTTTGTCATTAAATGGTGTACGTCAGGCGACCGACCGACTGAACAATCAGTGGTACATTAACCAAAGCCGGACGTTAAGCCCGACGGAGAGAGAAATATTGCGCTTTATTTCGCGTGGCTACTCAATGACACAAATTGCCGAGCAGCTTAAACGCAATATCAAAACGATCCGTGCACATAAATTGAATGTGAGGGCGAAACt                     <  2:294146/251‑1 (MQ=255)
                                                                                                                                                                                                                             gtgtATTAAGTAAAGCGTCAACGTTGGAGATTTTTCATCAGGAACTCTTTTTGTCATTAAATGGTGTACGTCAGGCGACCGACCGACTGAACAATCAGTGGTACATTAACCAAAGcc                                                                                                                                                          <  1:326366/117‑1 (MQ=255)
                                                                                                                                                                                                                             gtgtATTAAGTAAAGCGTCAACGTTGGAGATTTTTCATCAGGAACTCTTTTTGTCATTAAATGGTGTACGTCAGGCGACCGACCGACTGAACAATCAGTGGTACATTAACCAAAGcc                                                                                                                                                          >  2:326366/1‑117 (MQ=255)
                                                                                                                                                                                                                                           gCGTCAACGTTGGAGATTTTTCATCAGGAACTCTTTTTGTCATTAAATGGTGTACGTCAGGCGACCGACCGACTGAACAATCAGTGGTACATTAACCAAAGCCGGACGTTAAGCCCGACGGAGAGAGAAATATTGCGCTTTATGTCGCGTGGCTACTCAATGACACAAATTGCCGAGCAGCTTAAACGCAATATCAAAACGATCCGTGCACATAAATTTAATGTGATGTCGAAACTGGGCGTCAGTTCTGa      <  1:243870/251‑1 (MQ=255)
                                                                                                                                                                                                                                              tCAACGTTGGAGATTTTTCATCAGGAACTCTTTTTGTCATTAAATGGTGTACGTCAGGCGACCGACCGACTGAACAATCAGTGGTACATTAACCAAAGCCGGACGTTAAGCCCGACGGAGAGAGAAATATTGCGCTTTATGTCGCGTGGCTACTCAATGACACAAATTGCCGAGCAGCTTAAACGCACTATCAAAACGATCCGTGCACATAAATTTAATGTGATGTCGAAACTGGACGTCAGTTCTGACGc   >  1:359237/1‑251 (MQ=255)
                                                                                                                                                                                                                                               cAACGTTGGAGATTTTTCATCAGGAACTCTTTTTGTCATTAAATGGTGTACGTCAGGCGACCGACCGACTGAACAATCAGTGGTACATTAACCAAAGCCGGACGTTAAGCCCGACGGAGAGAGAAATATTGCGCTTTATGTCGCGTGGCTACTCAATGACACAAATTGCCGAGCAGCTTAAACGCAATATCAAAACGATCCGTGCACATAAATTTAATGTGATGTCGAAACTGGGCGTCAGTTCTGACGCa  <  1:179420/251‑1 (MQ=255)
                                                                                                                                                                                                                                                     |                                                                                                                                                                                                                                                     
CCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGACTGAACTGGCGATTAAGTTTCCGCGTACCCGGCGTTTAGTTATTGCGGATGATGATATTGAAGCTCGGCTGATTGGTTCATTGTCGCCATCACCGCTGGACGGTGTATTAAGTAAAGCGTCAACGCTGGAGATTTTTCATCAGGAACTCTTTTTGTCATTAAATGGTGTACGTCAGGCGACCGACCGACTGAACAATCAGTGGTACATTAACCAAAGCCGGACGTTAAGCCCGACGGAGAGAGAAATATTGCGCTTTATGTCGCGTGGCTACTCAATGACACAAATTGCCGAGCAGCTTAAACGCAATATCAAAACGATCCGTGCACATAAATTTAATGTGATGTCGAAACTGGGCGTCAGTTCTGACGCA  >  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: