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,5170CT97.7% 152.5 / NA 45L120L (CTG→TTG) bglJDNA‑binding transcriptional regulator BglJ
Reads supporting (aligned to +/- strand):  ref base C (1/0);  major base T (20/23);  minor base G (1/0);  total (22/23)
Fisher's exact test for biased strand distribution p-value = 4.77e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGACTGAACTGGCGATTAAGTTTCCGCGTACCCGGCGTTTAGTTATTGCGGATGATGATATTGAAGCTCGGCTGATTGGTTCATTGTCGCCATC‑ACCGCTGGACGGTGTATTAAGTAAAGCGTCAACGCTGGAGATTTTTCATCAGGAACTCTTTTTGTCATTAAATGGTGTACGTCAGGCGACCGACCGACTGAACAATCAGTGGTACATTAACCAAAGCCGGACGTTAAGCCCGACGGAGAGAGAAATATTGCGCTTTATGTCGCGTGGCTACTCAATGACACAAATTGCCGAGCAGCTTAAACGCAATATCAAAACGATCCGTGCACATAAATTTAATGTGATGTCGAAACTGGGCGTCAGTTCTGACGC  >  NC_000913/4604273‑4604761
                                                                                                                                                                                                                                                     |                                                                                                                                                                                                                                                    
ccATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTGCCATGAGAAGGGAGCGCAGAGAGGGATTATCTTGCCTGACTGAACTGGCGATTAAGTTTCCGCGTACCCGGCGTTTAGTTATTGCGGATGATGAGATTGAAGCTCGGCTGATTGGTTCATTGTCGCCATC‑ACCGCTGGACGGTGTATTAAGTAAAGCGTCAACGTTGgag                                                                                                                                                                                                                                                 <  1:348787/250‑1 (MQ=255)
        ataCCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTTTGCCAAGGGGAGTGAGCGCCGAGAGGGATTATCTTGCCTGACTGGACTTGCGGTTAAGTTTCCGCGTACCCGGCGGTTAGTTATTGCGGATGATGATATTTAAGCTCGGCTGATTTGTTCCTTGTCGCCATC‑ACCGCCGGGCGGGTTATTAAGTAAAACGTCAACGGTGGAGAttttttat                                                                                                                                                                                                                                        >  1:401307/1‑248 (MQ=255)
          aCCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGACTGAACTGGCGATTAAGTTTCCGCGTACCCGGCGTTTAGTTATTGCGGATGATGATATTGAAGCTCGGCTGATTGGTTCATTGTCGCCATC‑ACCGCTGGACGGTGTATTAAGTAAAGCGTCAACGTTGGAGATTTTtcatc                                                                                                                                                                                                                                       <  1:531701/250‑1 (MQ=255)
           ccAGTCAGGAGTCATTTCAGGATGCCATGGCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGACTGAACTGGCGATTAAGTTTCCGCGTACCCGGCGTTTAGTTATTTCGGATGATGATATTGAAGCTCGGCTGATTGGTTCATTGTCGCCATC‑ACCGCTGGACTGTGTATTAAGTAAAGCGTCAACGTTGGAGATTTTTCATCAg                                                                                                                                                                                                                                     <  2:566220/251‑1 (MQ=255)
                     gttaTTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTGTTCTGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGACTGAACTGGCGATTAAGTTTCCGCGTACCCGGCGTTTAGTTATTGCGGATGATGATATTGAAGCTCGGCTGATTGGTTCATTGTCGCCATC‑ACCGCTGGACGGTGTATTAAGTAAAGCGTCAACGTTGGAGATTTTTCATCAGGAACTCtttt                                                                                                                                                                                                                           <  2:146249/248‑1 (MQ=255)
                        aTTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGACTGAACTGGCGATTAAGTTTCCGCGTACCCGGCGTTTAGTTATTGCGGATGATGATATTGAAGCTCGGCTGATTGGTTCATTGTCGCCATC‑ACCGCTGGACGGTGTATTAAGTAAAGCGTCAACGTTGGAGATTTTTCATCAGGAACTCTTTTTGt                                                                                                                                                                                                                        <  1:319942/251‑1 (MQ=255)
                                    cATGTCGCGGATCTCGTTTGCGGCGGTCAATTTTTCTTTTTTTACCACTAGAAGGGCGTGGAGGGAGGGGGTATCTTTCCTGACTGACCTGGCGAATAAGTTTCCGCGTAACCCGCGGTTTGTTATTTCGGGTGGTGATATTGAAGCGCGGCTGATTGGGTCCTTGTCGCCATCTCCCGC‑GGGCGGTGTATTAAATAAAACGGTAAAGCTGGGGATTTTTTCTCAGGAACTCTTTTTTTCATTAAATgggg                                                                                                                                                                                                            >  2:422502/1‑249 (MQ=255)
                                                       gcggcgGTCATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGACTGAACGGGCGATTAAGTTTCCGCGTCCCCGGCGTTTAGTTATTGCGGATGATGATATTGAAGCTCGGCTGATTGGTTCATTGTCGCCATC‑ACCGCTGGACGGTGTATTAAGTAAAGCGTCAACGTTGGAGATTTTTCATCAGGAACTCTTTTTGTCATTAAATGGTGTACGTCAGGcgaccgacc                                                                                                                                                                                          <  2:205902/250‑1 (MQ=255)
                                                           cggTCATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGACTGAACTGGCGATTAAGTTTCCGCGTACCCGGCGTTTAGTTATTGCGGATGATGATATTGAAGCTCGGCTGATTGGTTCATTGTCGCCATC‑ACCGCTGGACGGTGTATTAAGTAAAGCGTCAACGTTGGAGATTTTTCATCAGGAACTCTTTTTGTCATTAAATGGTGTACTTCAGGCGACCGACCGACTg                                                                                                                                                                                     <  1:451047/251‑1 (MQ=255)
                                                                          tttctGCAATGAGAAGTGAGCGAAGAGAGGGATTATCTTGCCTGACTGAACTGGCGATTAAGTTTCCGCGTACCCGGCGTTTAGTTATTGCGGATGATGATATTGAAGCTCGGCTGATTGGTTCATTGTAGCCATC‑ACCGCTGGACGGTGTATTAAGTAAAGCGTCAACGTTGGAGATTTTTCATCAGGAACTATTTTTGTCATTAAATGGTGTACGTCAGGCGACCGACCGACTGAACAATCAATGGTa                                                                                                                                                                       >  2:426457/1‑250 (MQ=255)
                                                                            tctGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGACTGAACTGGCGATTAAGTTTCCGCGTACCCGGCGTTTAGTTATTGCGGATGATGATATTGAAGCTCGGCTGATTGGTTCATTGTCGCCATC‑ACCGCTGGACGGTGTATTAAGTAAAGCGTCAACGTTGGAGATTTTTCATCAGGAACTCTTTTTGTCATTAAATGGTGTACGTCAGGCGACCGACCGACTGAACAATCAGTGGTACAt                                                                                                                                                                    >  2:103997/1‑251 (MQ=255)
                                                                              tGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGACTGAACTGGCGATTAAGTTTCCGCGTACCCGGCGTTTAGTTATTGCGGATGATGATATTGAAGCTCGGCTGATTGGTTCATTGTCGCCATC‑ACCGCTGGACGGTGTATTAAGTAAAGCGTCAACGTTGGAGATTTTTCATCAGGAACTCTTTTTGTCATTAAATGGTGTACGTCAGGCGACCGACCGACTGAACAATCAGTGGTACAtt                                                                                                                                                                   >  2:17868/1‑250 (MQ=255)
                                                                               gCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGACTGAACTGGCGATTAAGTTTCCGCGTACCCGGCGTTTAGTTATTGCGGATGATGATATTGAAGCTCGGCTGATTGGTTCATTGTCGCCATC‑ACCGCTGGACGGTGTATTAAGTAAAGCGTCAACGTTGGAGATTTTTCATCAGGAACTCTTTTTGTCATTAAATGGTGTACGTCAGGCGACCGACCGACTGAACAATCAGTGGTACATTaa                                                                                                                                                                 <  1:156299/251‑1 (MQ=255)
                                                                                gcATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGAGTGAACTGGCGATTAAGTTTCCGCGTGCCGGGTGGTTAGTTATTTCGGAGTATGTTGTTGAGGCGGTGTGGATGGGTTCATTGTGGACATC‑ACCAGTGGACGGTGGATGAAGTAAAGCTTCAACTTTGGAGATGTTTCATCAGGAACTGTGTTGGTCAGTAAATGTTGTACGTCGGGTGACGTACCGACTAAGCAATCAGTGGTagtttaac                                                                                                                                                                <  2:536073/250‑8 (MQ=255)
                                                                                                     agGGATTATCTTGCCTGACTGAACTGGCGATTAAGTTTCCGCGTACCCGGCGTTTAGTTATTGCGGATGATGATATTGAAGCTCGGCTGATTGGTTCATTGTCGCCATC‑ACCGCTGGACGGTGTATTAAGTAAAGCGTCAACGTTGGAGATTTTTCATCAGGAACTCTTTTTGTCATTAAATGGTGTACGTCAGGCGACCGACCGACTGAACAATCAGTGGTACATTAACCAAAGCCGGACGTTAAGCCCg                                                                                                                                           >  1:425434/1‑251 (MQ=255)
                                                                                                       ggATTATCTTGCCTGACTGAACTGGCGATTAAGTTTCCGCGTACCCGGCGTTTAGTTATTGCGGATGATGATATTGAAGCTCGGCTGATTGGTTCATTGTCGCCATC‑ACCGCTGGACGGTGTATTAAGTAAAGCGTCAACGTTGGCGATTTTTCATCAGGACCTCTTTTTGTCATTAAATGGTGTACGTCAGGCGACCGACCGACTGAACAATCAGTGGTACATTAACCAAAGCCGGACGTTAAGCCCGAc                                                                                                                                         >  1:87987/1‑251 (MQ=255)
                                                                                                                tGCCTGACTGAACTGGCGATTAAGTTTCCGCGTACCCGGCGTTTAGTTATTGCGGATGATGATATTGAAGCTCGGCTGATTGGTTCATTGTCGACATC‑ACCGCTGGACGGTGTATTAAGTAAAGCGTCAACGTTGGAGATTTTTCATCAGGAACTCTTTTTGTCATTAAATGGTGTACGTCAGGCGACCGACCGACTGAACAATCAGTGGTACATTAACCAAAGCCGGACGTTAAGCCCGACGgagcg                                                                                                                                   >  1:93386/1‑246 (MQ=255)
                                                                                                                          aaCTGGCGATTAAGTTTCCGCGTACCCGGCGTTTAGTTATTGCGGATGATGATATTGAAGCTCGGCTGATTGGTTCATTGTCGCCATC‑ACCGCTGGACGGTGTATTAAGTAAAGCGTCAACGTTGGAGATTTTTCATCAGGAACTCTTTTTGTCATTAAATGGTGTACGTCAGGCGACCGACCGACTGAACAATCAGTGGTACATTAACCAAAGCCGGACGTTAAGCCCGACGGAGAGAGAAATATTgcgc                                                                                                                      <  1:450638/251‑1 (MQ=255)
                                                                                                                                     aaGTTTCCGCGTACCCGGCGTTTAGTTATTGCGGATGATGATATTGAAGCTCGGCTGATTGGTTCATTGTCGCCATC‑ACCGCTGGACGGTGTATTAAGTAAAGCGTCAACGTTGGAGATTTTTCATCAGGAACTCTTTTTGTCATTAAATGGTGTACGTCAGGCGACCGACCGACTGAACAATCAGTGGTACATTAACCACAGCCGGACGTTAAGCCCGACGGAGAGAGACATATTGCGCTTTATGTcgcg                                                                                                           >  2:6057/1‑251 (MQ=255)
                                                                                                                                              cgTACCCGGCGTTTAGTTATTGCGGATGATGATATTGAAGCTCGGCTGATTGGTTCATTGTCGCCATC‑ACCGCTGGACGGTGTATTAAGTAAAGCGTCAACGTTGGAGATTTTTCATCAGGAACTCTTTTTGTCATTAAATGGTGTACGTCAGGCGACCGACCGACTGAACAATCAGTGGTACATTAACCAAAGCCGGACGTTAAGCCCGACGGAGAGAGAAATATTGCGCTTTATGTCGCGTGGCTACTc                                                                                                  >  1:371772/1‑251 (MQ=255)
                                                                                                                                                gtcccGGCGTTTAGTTAGTGCGGATGATGATATTGAAGCTCGGCTGATTGGTTCATTGTCGCCATC‑ACCGCTGGACGGTGTATTAAGTAAAGCGTCATCGTTGTAGATTTTTCATCAGGAACTGTGTTTGTCATTAAATTGTGGACGTCAGGCGACCGACCGACTGAACAATCAGTTGTACATTAACCAAAGCCGTACGTTAAGCCCGACGGAGAGAGAAATATTGCGCTTTATGTCGCGTGGTTACTCaa                                                                                                <  2:245465/249‑1 (MQ=255)
                                                                                                                                                    cGGCGTTTAGTGATTGCGGATGAGGATATTGAAGCTCGGCTGATTGGTTCATTGTCGCCATC‑ACCGCTGGACGGTGTATTAAGTAAAGCGTCAACGTTGGAGATTTTTCATCAGGAACTCTTTTTGTCATTAAATGGTGTACGTCAGGCGACCGACCGACTGAACAATCAGTGGTACATTAACCAAAGCCGGACGTTAAGCCCGACGGAGAGAGAAATATTGCGCTTTATGTCGCGTGGCTACTCAATGac                                                                                            <  1:574648/251‑1 (MQ=255)
                                                                                                                                                     ggCGTTTAGTTATTGCGGATGATGATATTGAAGCTCGGCTGATTGGTTCATTGTCGCCATC‑ACCGCTGGACGGTGTATTAAGTAAAGCGTCAACGTTGGAGATTTTTCATCAGGAACTCTTTTTGTCATTAAATGGTGTACGTCAGGCGACCGACCGACTGAACAATCAGTGGTACATTAACCAAAGCCGGACGTTAAGCCCGACGGAGAGAGAAATATTGCGCTTTATGTCGCGTGGCTACTCAATGac                                                                                            >  1:181766/1‑250 (MQ=255)
                                                                                                                                                      gCGTTTAGTTATTGCGGATGATGATATTGAAGCTCGGCTGATTGGTTCATTGTCGCCATC‑ACCGCTGGACGGTGTATTAAGTAAAGCGTCAACGTTGGAGATTTTTCATCAGGAACTCTTTTTGTCATTAAATGGTGTACGTCAGGCGACCGACCGACTGATCAATCAGTGGTACATTAACCAGAGCCGGACGTTAAGCCCGACGGAGAGAGAAATATTGCGCTTTATGTCGCGTGGCTACTCAATGacac                                                                                          <  1:422502/251‑1 (MQ=255)
                                                                                                                                                       cGTTTAGTTATTGCGGATGATGATATTGAAGCTCGGCTGATTGGTTCATTGTCGCCATC‑ACCGCTGGACGGTGTATTAAGTAAAGCGTCAACGTTGGAGATTTTTCATCAGGAACTCTTTTTGTCATTAAATGGTGTACGTCAGGCGACCGACCGACTGAACAATCAGTGGTACATTAACCAAAGCCGGACGTTAAGCCCGACGGAGAGAGAAATATTGCGCTTTATGTCGCGTGGCTACTCAATGacac                                                                                          >  2:447850/1‑250 (MQ=255)
                                                                                                                                                                       atgatgatATTGAAGCTCGGCTGATTGGTTCATTGTCGCCATC‑ACCGCTGGACGGTGTATTAAGTAAAGCGTCAACGTTGGAGATTTTTCATCAGGAACTCTTTTTGTCATTAAATGGTGTACGTCAGGCGACCGACCGACTGAACAATCAGTGGTACATTAACCAAAGCCGGACGTTAAGCCCGACGGAGATAGAAATATTGCGCTTTATGTCGCGTGGCTACTCAATGACACAAATTGCCGAGCAGCtt                                                                         <  2:529022/251‑1 (MQ=255)
                                                                                                                                                                           tgatATTGAAGCTCGGCTGATTGGTTCATTGTCGCCATC‑ACCGCTGGACGGTGTATTAAGTAAAGCGTCAACGTTGGAGATTTTTCATCAGGAACTCTTTTTGTCATTAAATGGTGTACGTCAGGCGACCGACCGACTGAACAATCAGTGGTACATTAACCAAAGCCGGACGTTAAGCCCGACGGAGAGAGAAATATTGCGCTTTATGTCGCGTGGCTAc                                                                                                    <  2:329061/220‑1 (MQ=255)
                                                                                                                                                                           tgatATTGAAGCTCGGCTGATTGGTTCATTGTCGCCATC‑ACCGCTGGACGGTGTATTAAGTAAAGCGTCAACGTTGGAGATTTTTCATCAGGAACTCTTTTTGTCATTAAATGGTGTACGTCAGGCGACCGACCGACTGAACAATCAGTGGTACATTAACCAAAGCCGGACGTTAAGCCCGACGGAGAGAGAAATATTGCGCTTTATGTCGCGTGGCTAc                                                                                                    >  1:329061/1‑220 (MQ=255)
                                                                                                                                                                               atTGAAGCTCGGCTGATTGGTTCATTGTCGCCATC‑ACCGCTGGACGGTGTATTAAGTAAAGCGTCAACGTTGGAGATTTTTCATCAGGAACTCTTTTTGTCATTAAATGGTGTACGTCAGGCGACCGACCGACTGAACAATCAGTGGTACATTAACCAAAGCCGGACGTTAAGCCCGACGGAGAGAGAAATATTGCGCTTTATGTCGCGTGGCTACTCAATGACACAAATTGCCGAGCAGCTTAAACGCa                                                                  <  1:17868/250‑1 (MQ=255)
                                                                                                                                                                                        cGGCTTATTGGTTCATTGTCGCCATC‑ACCGCTGGACGGTGTATTAAGTAAAGCGTCAACGTTGGAGATTTTTCATCAGGAACTCTTTTTGTCATTAAATGGTGTACGTCAGGCGACCGACCGACTGAACAATCAGTGGTACATTAACCAAAGCCGTACGTTAAGCCCGACGGAGAGAGAAATATTGCGCTTTATGTCTCGTGGCTACTCAATGACACAAATTGCCGAGCAGCTTAAACGCAATATCAAAAc                                                        <  2:238533/251‑1 (MQ=255)
                                                                                                                                                                                                 ggTTCATTGTCGCCATC‑ACCGCTGGACGGTGTATTAAGTAAAGCGTCAACGTTGGAGATTTTTCATCAGGAACTCTTTTTGTCATTAAATGGTGTACGTCAGGCGACCGACCGACTGAACAATCAGTGGTACATTAACCAAAGCCGGACGTTAAGCCCGACGGAGAGAGAAATATTGCGCTTTATGTCGCGTGGCTACTCAATGACACAAATTGCCGAGCAGCTTAAACGCAATATCAAAACGATCCGTg                                                <  1:585295/250‑1 (MQ=255)
                                                                                                                                                                                                     cATTGTCGCCATC‑ACCGCTGGACGGTGTATTAAGTAAAGCGTCAACGTTGGAGATTTTTCATCAGGAACTCTTTTTGTCATTAAATGGTGTACGTCAGGCGACCGACCGACTGAACAATCAGTGGTACATTAACCAAAGCCGGACGTTAAGCCCGACGGAGAGAGAAATATTGCGCTTTATGTCGCGTGGCTACTCAATGACACAAATTGCCGAGCAGCTTAAACGCAATATCAAAACGATCCGTGCACAt                                           >  1:332614/1‑251 (MQ=255)
                                                                                                                                                                                                                   aCCGCTGGACGGTGTATTAAGTAAAGCGTCAACGTTGGAGATTTTTCATCAGGAACTCTTTTTGTCATTAAATGGTGTACGTCAGGCGACCGACCGACTGAACAATCAGTGGTACATTAACCAAAGCCGGACGTTAAGCCCGACGGAGAGAGAAATATTGCGCTTTATGTCGCGTGGCTACTCAATGACACAAATTGCCGAGCAGCTTAAACGCAATATCAAAACGATCCGTGCACATAAATTTAATGTGa                              <  1:453344/251‑1 (MQ=255)
                                                                                                                                                                                                                    ccGCTGGACGGTGTATTAAGTAAAGCGTCAACGTTGGAGATTTTTCATCAGGAACTCTTTTTGTCATTAAATGGTGTACGTCAGGCGACCGACCGACTGAACAATCAGTGGTACATTAACCAAAGCCGGACGTTAAGCCCGACGGAGAGAGAAATATTGCGCTTTATGTCGCGTGGCTACTCAATGACACAAATTGCCGAGCAGCTTAAACGCAATATCAAAACGATCCGTGCACATAAATTTAATGTGAt                             <  1:441671/251‑1 (MQ=255)
                                                                                                                                                                                                                       cTGGACGGTGTATTAAGTAAAGCGTCAACGTTGGAGATTTTTCATCAGGAACTCTTTTTGTCATTAAATGGTGTACGTCAGGCGACCGACCGACTGAACAATCAGAGGTACATTAACCAAAGCCGGACGTTAAGCCCGACGGAGAGAGAAATATTGCGCTTTATGTCGCGTGCCTACTAAATGACACAAATTACCGAGCAGCTTAAACGCAATATCAAAACGATCCGTGCACATAAATTTAATGTGATGTc                          >  1:20883/1‑251 (MQ=255)
                                                                                                                                                                                                                         ggACGGTGTATTAAGTAAAGCGTCAACGTTGGAGATTTTTCATCAGGAACTCTTTTTGTCATTAACTGGTGTACGTCAGGCGACCGACCGACTGAACATTCAGTGGTACATTAACCAAAGCCGGACGTTAAGCCCGACGGAGAGAGAAATATTGCGCTTTATGTCGCGAGGCTACTCAATGACACAAATTGCCGAGCAGCTTAAACGCAATATCAAAACGTTCCGCGCACATAAATTTAACGTGATGTCg                         >  2:167796/1‑250 (MQ=255)
                                                                                                                                                                                                                          gACGGTGTATTAAGTAAAGCGTCAACGTTGGAGATTTTTCATCAGGAACTCTTTTTGTCATTAAATGGTGTACGTCAGGCGACCGACCGACTGAACAATCAGTGGTACATTAACCAAAGCCGGACGTTAAGCCCGACGGAGAGAGAAATATTGCGCTTTATGTCGCGTGGCTACTCAATGACAAAAATTGCCGAGCAGCTTAAACGCAATTTCAAAAAGATCCGTGCACATAAATTTAATGTGATGTCGa                        >  2:192272/1‑250 (MQ=255)
                                                                                                                                                                                                                            cGGTGTATTAAGTAAAGCGTCAACGTTGGAGATTTTTCATCAGGAACTCTTTTTGTCATTAAATGGTGTACGTCAGGCGACCGACCGACTGAACAATCAGTGGTACATTAACCAAAGCCGGACGTTAAGCCCGACGGAGAGAGAAATATTGCGCTTTATGTCGCGTGGCTACTCAATGACACAAATTGCCGAGCAGCTTAAACGCAATATCAAAACGATCCGTGCACATAAATTTAATGTGATGTCGAAAc                     >  1:580053/1‑251 (MQ=255)
                                                                                                                                                                                                                                  aTTAAGTACAGCGTCAACGTTGGAGATTTTTCATCCGGCACTCTTTTTGTCATTAAATGATGTACGTCAGGCGAACGACCGACTGAACAATCAGTGGTACATTAACCACAGCCGGACGTTAAGCCCGACGGAGAGAGAAAAATGGAGCTTTATGTCGCGTGGCTAATCAATGCCACAAATTGCCGAGCAGCTTAAACGCAATATCCAAACGATCCGGGCCAATAAATTTAATGAGATATCGTAACTAGGCg               >  2:355219/1‑251 (MQ=255)
                                                                                                                                                                                                                                       gTAAAGCGTCAACGTTGGAGATTTTTCATCAGGAACTCTTTTTGTCATTAAATGGTGTACGTCAGGCGACCGACCGACTGAACAATCAGTGGTACATTAACCAAAGCCGGACGTTAAGCCCGACGGAGAGAGAAATATTGCGCTTTATGTCGCGTGGCTACTCAATGACACAAATTGCCGAGCAGCTTAAACGCAATATCAAAACGATCCGTGCACATAAATTTAATGTGATGTCGAAACTGGGCGTCAGt          <  1:273365/251‑1 (MQ=255)
                                                                                                                                                                                                                                          aaGCGTCAACGTTGGAGATTTTTCATCAGGAACTCTTTTTGTCATTAAATGGTGTACGTCAGGCGACCGACCGACTGAACAATCAGTGGTACATTAACCAAAGCCGGACGTTAAGCCCGACGGAGAGAGAAATATTGCGCTTTATGTCGCGTGGCTACTCAATGACACAAATTGCCGAGCAGCTTAAACGCAATATCAAAACGATCCGTGCACATAAATTTAATGTGATGTCGAAACTGGGCGTCAGTTCt       <  1:196216/251‑1 (MQ=255)
                                                                                                                                                                                                                                           aGCGTCAACGTTGGAGATTTTTCATCAGGAACTCTTTTTGTCATTAAATGGTGTACGTCAGGCGACCGACCGACTGAACAATCAGTGGTACATTAACCAAAGCCGGACGTTCAGCCCGACGGAGAGAGAAATATTGCGCTTTATGTCGCGTGGCTACTCAATGACACAAATTGCCGAGCAGCTTAAACGCAAAATCAAAACGATCCGTGCACATAAATTTAATGTGATGTCGAAACTGGGCGTCAGTTCTg      >  1:215174/1‑251 (MQ=255)
                                                                                                                                                                                                                                            gCGTCAACGTTGGAGATTCTTCATCAGGAACTCTTTTTGTCATTAAATGGTGTACGTCAGGCGACCGACCGACTGAACAATCAGTGGTACATTAACCAAAGCCGGACGTTAAGCCCGACGGAGAGAGAAATATTGCGCTTTATGTCGCGTGGCTACTCAATGACACAAATTGCCGAGCAGCTTAAACGCAATATCAAAACGATCCGTGCACATAAATTTAATGTGATGTCGAAACTGGGCGTCAGTTCTGa     <  1:302747/251‑1 (MQ=255)
                                                                                                                                                                                                                                               tCAACGTTGGAGATTTTTCATCGGGAACTCTTTTTGTCATTAAATGGTGTACGTCAGGCGACCGACCGACTGCACAATCAGTGGTACATTAACCAAAGCCGGACGTTAAGCCCGACGGAGAGAGAAATATTGCGCTTTATGTCGCGTGGCTACTCAATGACCCAAATTGCCGAGCAGCTTAAACGCAATATCAAAACGATCCGTGCACATAAATTTAATGTGATGTCGAAACTGGGCGTCAGTTCAGACGc  >  2:285473/1‑251 (MQ=255)
                                                                                                                                                                                                                                                 aaCGTTGGAGATTTTTCATCAGGAACTCTTTTTGTCATTAAATGGTGTACGTCAGGCGACCGACCGACTGAACAATCAGTGGTACATTAACCAAAGCCGGACGTTAAGCCCGACGGAGAGAGAAATATTGCGCTTTATGTCGCGTGGCTACTCAATGACACAAATTGCCGAGCAGCTTAAACGCAATATCAAAACGATCCGTGCACATAAATTTAATGTGATGTCGAAACTGGGCGTCAGTTCTGACGc  >  1:521032/1‑249 (MQ=255)
                                                                                                                                                                                                                                                     |                                                                                                                                                                                                                                                    
CCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGACTGAACTGGCGATTAAGTTTCCGCGTACCCGGCGTTTAGTTATTGCGGATGATGATATTGAAGCTCGGCTGATTGGTTCATTGTCGCCATC‑ACCGCTGGACGGTGTATTAAGTAAAGCGTCAACGCTGGAGATTTTTCATCAGGAACTCTTTTTGTCATTAAATGGTGTACGTCAGGCGACCGACCGACTGAACAATCAGTGGTACATTAACCAAAGCCGGACGTTAAGCCCGACGGAGAGAGAAATATTGCGCTTTATGTCGCGTGGCTACTCAATGACACAAATTGCCGAGCAGCTTAAACGCAATATCAAAACGATCCGTGCACATAAATTTAATGTGATGTCGAAACTGGGCGTCAGTTCTGACGC  >  NC_000913/4604273‑4604761

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: ‑