Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,124,401 G→T A19E (GCA→GAA)  cytR ← DNA‑binding transcriptional repressor CytR

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,124,4010GT95.7% 125.5 / ‑1.1 47A19E (GCA→GAA) cytRDNA‑binding transcriptional repressor CytR
Reads supporting (aligned to +/- strand):  ref base G (1/1);  new base T (22/23);  total (23/24)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.95e-01

ACCAGATATCCGTGATTTGCCGCCGTAACTTCGATACCGCGAATAATTTCGCTAAAGAAGGGATCGCAGATATCCGGGACAATCACCAGAATGGTGCGGGATTCATTACGCTTGACGTTGCGCCCCATAGGCTGCGGTAAATAACCCACTTCCCGGGCCGCTTTTTCAACCCGATTACGGGTGGCCTGGGAGACTTTATCGGGATTCATTAATGCTCGGGAGACGGTCGCTGTAGAGACTTTTGCCTTGAGGGCAACGTCTTTCATGGTCGCGGCAGTTTCCTGCTTCTTCGCTTTCACACTCACTCCTCGCCTGGCACGTCAGGCGTACTACATCCATGTTTACTTCACATCGGCAACATTTTTAGCAGATAGCGCGTGAAAACGGTTACAGAATTTTCATGAAAAGTGTGATGAATATTGAATTTTTCGATCCGCCTCGCATCGTGAGCGGTTTAACCCTCA  >  NC_000913/4124158‑4124621
                                                                                                                                                                                                                                                   |                                                                                                                                                                                                                             
cccAGATATCCGTGATTTGCCGCCGTAACTTCGATACCGCGAATAATTTCGCTAAAGAAGGGATCGCAGATATCCGGGACAAGCACCAGAATGGTGCGGGATTCATTACGCTTGACGTTGCGCCCCATAGGCTGCGGTAAATAACCCACTTCCCGGGCCGCTTTTTCAACCCGATTACGGGTGGCCTGGGAGACTTTATCGGGATTCATTAATGCTCGGGAGACGGTCGCTGTAGAGACTTTTTCCTTGAg                                                                                                                                                                                                                        <  1:416150/250‑1 (MQ=255)
                   ccgccgTAACTTCGATACCGCGAATAATTTCGCTAAAGAAGGGATCGCAGATATCCGGGACAATCACCAGCATGATGCGGGATTCATTACGCTTGCCGTTGCGCCCCATAGGCTGCGGTAAATAACCCACTTCCCGGGCCGCTTTTTCAACCCGCTTACGGGTGGCCTGGGAGACTTTATCGGGATTCATTAATGCTCGGGAGACGGTCGCTGTAGAGACTTTTTCCTTGCGGGCAACGTCTTTCATGGTc                                                                                                                                                                                                     >  2:416010/1‑251 (MQ=255)
                    cgccgTACCTTCGATACCGCGAATAATTTCGCTCAAGAAGGGATCGCAGATATACGGGACAATCACCAGAATGGTGCGGGATTCATTACGCTTGACGTTGCGCACCATAGGCTGCGGTAAATAACCCACTTCCCGGGCCGCTTTTTCAACCCGATTACGGGTGGCCTGGGAGACTTTATCGGGCTTCATTAATGCTCGGGAGACGGTCGCTGTAGAGACTTTTTCCTTGAGGGCACCGTCTTTCATGGTcg                                                                                                                                                                                                    >  2:277223/1‑251 (MQ=255)
                      ccgTAACTTCGATACCGCGAATAATTTCGCTAAAGAAGGGATCGCAGATATCCGGGACAATCACCAGAATGGTGCGGGATTCATTACGCTTGACGTTGCGCCCCATAGGCTGCGGTAAATAACCCACTTCCCGGGCCGCTTTTTCAACCCGATTACGGGTGGCCTGGGAGACTTTATCGGGATTCATTAATGCTCGGGAGACGGTCGCTGTAGAGACTTTTTCCTTGAGGGCAACGTCTTTCATGGTcgcg                                                                                                                                                                                                  <  1:354602/251‑1 (MQ=255)
                       cgTAACTTCGATACCGCGAATAATTTCGCTAAAGAAGGGATCGCAGATATCCGGGACAATCACCAGAATGGTGCGGGATTCATTACGCTTGACGTTGCGCCCCATAGGCTGCGGTAAATAACCCACTTCCCGGGCCGCTTTTTCAACCCGATTACGGGTGGCCTGGGCGACTTTATCGGGTTTCATTAATGCTCGGGAGCCGGTCGCTGTAGAGACTTTTTCCTTGAGGGCAACGTCTTTCATGGTCCCgg                                                                                                                                                                                                 >  2:166420/1‑251 (MQ=255)
                        gTAACTTCGATACCGCGAATAATTTCGCTAAAGAAGGGATCGCAGATATCCGGGACAATCACCAGAATGGTGCGGGATTCATTACGCTTGACGTTGCGCCCCATAGGCTGCGGTAAATAACACACTTCCCGGGCCGCTTTTTCAACCCGATTACGGGTGGCCTGGGAGACTTTATCGGGCTTCATTAATGCTCGGGCGACGGTCGCTGTAGAGACTTTTTCCTTGAGGGCAACGTCTTTCATGGTCGCGGc                                                                                                                                                                                                >  1:443993/1‑251 (MQ=255)
                           aCTTCGATACCGCGAATAATTTCGCTAAAGAAGGGATCGCAGATATCCGGGACAATCACCAGAATGGTGCGGGATTCATTACGCTTGACGTTGCGCCCCATAGGCTGCGGTAAATAACCCACTTCCCGGGCCGCTTTTTCAACCCGATTACGGGTGGCCTGGGAGACTTTATCGGGATTCATTAATGCTCGGGAGACGGTCGCTGTAGAGACTTTTTCCTTGAGGGCAACGTCTTTCATGGTCGCGGCAGt                                                                                                                                                                                             >  1:50670/1‑251 (MQ=255)
                              tCGATACCGCGAATAATTTCGCTAAAGAAGGGATCGCAGATATCCGGGACAATCACCAGAATGGTGCGGGATTCATTACGCTTGACGTTGCGCCCCATAGGCTGCGGTAAATAACCCACTTCCCGGGCCGCTTTTTCAACCCGATTACGGGTGGCCTGGGAGACTTTATCGGGATTCATTAATGCTCGGGAGACGGTCGCTGTAGAGACTTTTTCCTTGAGGGCAACGTCTTTCATGGTCGCGGCAGttt                                                                                                                                                                                           <  1:242738/250‑1 (MQ=255)
                               cGATACCGCGAATAATTTCGCTAAAGAAGGGATCGCAGATATCCGGGACAATCACCAGAATGGTGCGGGATTCATTACGCTTGACGTTGCGCCCCATAGGCTGCGGTAAATAACCCACTTCCCGGGCCGCTTTTTCAACCCGATTACGGGTGGCCTGGGAGACTTTATCGGGATTCATTAATGCTCGGGAGACGGTCGCTGTAGAGACTTTTTCCTTGAGGGCAACGTCTTTCATGGTCGCGGCAGTTTcc                                                                                                                                                                                         >  1:347903/1‑251 (MQ=255)
                                      gcgAATAATTTCGCTAAAGAAGGGATCGCAGATATCCGGGACAATCACCAGAATGGTGCGGGATTCATTACGCTTGACGTTGCGCCCCATAGGCTGCGGTAAATAACCCACTTCCCGGGCCGCTTTTTCAACCCGATTACGGGTGGCCTGGGAGACTTTATCGGGATTCATTAATGCTCGGGAGACGGTCGCTGTAGAGACTTTTTCCTTGAGGGCAACGTCTTTCATGGTCGCGGCAGTTTCCTGcttct                                                                                                                                                                                  <  1:10041/251‑1 (MQ=255)
                                                             gATCGCAGATATCCGGGACAATCACCAGAATGGTGCGGGATTCATTACGCTTGACGTTGCGCCCCATAGGCTGCGGTAAATAACCCACTTCCCGGGCCGCTTTTTCAACCCGATTACGGGTGGCCTGGGAGACTTTATCGGGATTCATTAATGCTCGGGAGACTGTCGCTGTAGAGACTTTTTCCTTGAGGGCAACGTCTTTCATGGTCGCGGCAGTTTCCTGCTTCTTCGCTTTCACACTCACTCCTCGc                                                                                                                                                           <  1:277223/251‑1 (MQ=255)
                                                                    gATAGCCGGGACAATCACCAGAATGGGGCGGGATTCATTACGCTTGACGTTGCGCCCCATAGGCTGCGGTAAATAACCCACTTCCCGGGCCGCTTTTTCAACCCGATTACGGGTGGCCTGGGAGACTTTATCGGGATTCATTAATGCTCGGGAGACGGTCGCTGTAGAGACTTTTTCCTTGAGGGCAACGTCTTTCATGGTCGCGGCAGTTTCCTGCTTCTTCGCTTTCACACTCACTCCTCGCCTGGCAc                                                                                                                                                    <  1:105446/251‑1 (MQ=255)
                                                                       gtCCGGGAGAATCACCAGATTGGTGCGTGATTCATTACGCTTGACGTTGCGCCCCATAGGCTGCGGTAAATAGCCCACTGTCCGGGCCGCTGTGTCAACCAGATTACGGGTGGCCTGGGAGACTTTATCGGGATTCATTAATGCTCGGGAGACGGTCGCTTTAGAGACTTTTTCCTTGAGGGCAACGTCTTTCATGGGCGCGGCAGTTGCCTGCTTCTTCGCTTTCACACTCACTCCTCGTCTGGCTCGTc                                                                                                                                                 <  2:79683/250‑1 (MQ=255)
                                                                             gACAATCACCAGAATGGTGCGGGCTTCATTACGCTTGACGTTGCGCCCCATAGGCTGCGGTAAATAACCCACTTCCCGGGCCGCTTTTTCAACCCGATTACGGGTGGCCTGGGAGACTTTATCGGGATTCATTAATGCTCGGGAGACGGTCGCTGTAGAGACTTTTTCCTTGAGGGCAACGTCTTTCATGGTCGCGGCAGTTTCCTGCTTCTTCGCTTTCACACTCACTCCTCGCCTGGCACGTCAGGCGt                                                                                                                                           >  1:370906/1‑251 (MQ=255)
                                                                                         aaTGGTGCGGGATTCATCACTCTTGACGTTGCGCCCCATAGGCTGCGGTTAATAACCCACGTCCCGGGCCGCTCTTTCAACCCGAGTACGGGGTGCCTGGGAGACTTTATCGGTATTCAGTAATGCTCGTGAGACGGTCGCTGTAGAGACTTGTTCCTTGAGGTCAACGTCGTTCATGGTCGCGGCAGTTTCCGGCTGCTTCGCTTTCACACTCACTCCTCGCCTGGCACGTCAGGCCTACTACATCCATg                                                                                                                               <  2:427302/251‑1 (MQ=255)
                                                                                            ggTGCGGGATTCATTACGCTTGACGTTGCGCCCCATAGGCTGCGGTAAATAACCCACTTCCCGGGCCGCTTTTTCAACCCGATTACGGGTGGCCTGGGAGACTTTATCGGGATTCATTAATGCTCGGGAGACGGTCGCTGTAGAGACTTTTTCCTTGAGGGCAACGTCTTTCATGGTCGCGGCAGTTTCCTGCTTCTTCGCTTTCACACTCACTCCTCGCCTGGCACGTCAGGCGTACTACATCCATGttt                                                                                                                            >  1:220462/1‑251 (MQ=255)
                                                                                               gCGGGATTCATTACGCTTGACGTTGCCCCCCATAGGCTGCGGTAAATAACCCACTTCCCGGGCCGCTTTTTCAACCCGAGTACGGGTGGCCTGGGAGACTTTATCGGGATTCATTAATGCTCGGGAGACGGTCGCTGTAGAGACTTTTTCCTTGAGGGCAACGTCTTTCATGGTCGCGGCAGTTTCCTGCTTCTTCGCTTTCACACTCACTCCTCGCCTGGCACGTCAGGCGTACTACATCCATGTTTACt                                                                                                                         <  1:462713/251‑1 (MQ=255)
                                                                                                    aTTCATTACGCTTGACGTTGCGCCCCATAGGCTGCGGTAAATAACCCACTTCCCGGGCCGCTTTTTCAACCCGATTACGGGTGGCCTGGGAGACTTTATCGGGATTCATTAATGCTCGGGAGACGGTCGCTGTAGAGACTTTTGCCTTGAGGGCAACGTCTTTCATGGTCGCGGCAGTTTCCTGCTTCTTCGCTTTCACACTCACTCCTCGCCTGGCACGTCAGGCGTACTACATCCATGTTTACTTcaca                                                                                                                    <  1:189111/251‑1 (MQ=255)
                                                                                                             gCTTGACGTTGAGCCACATAGGCTGCGGTAACTAACCCACTTCCCGGGCCGCTTTTTCAACCCGATTACGGGTGGCCTGGGAGACTTTATCGGGATTCATTAATGCTCGGGGGACGGTCGCTGTAGAGACTTTTTCCTTGAGGGCAACGTCTTTCATGGTCGCGGCAGTTTCCTGCTTCTTCGCTTTCACACTCACTCCTCGCCTGGCACGTCAGGCGTACTACATCCATGTTTACTTCACATAGGCAAc                                                                                                            >  2:67937/1‑250 (MQ=255)
                                                                                                                  aCGTTGCGCCCCATAGGCTGCGGTAAATAACCCACTTCCCGGGCCGCTTTTTCAACCCGATTACGGGTGGCCTGGGAGACTTTATCGGGATTCATTAATGCTCGGGAGACGGTCGCTGTAGAGACTTTTTCCTTGAGGGCAACGTCTTTCATGGTCGCGGCAGTTTCCTGCTTCTTCGCTTTCACACTCACTCCTCGCCTGGCCCGTCAGGCGGACTACATCCATGTTTACTTCACATCGGCAACAttttt                                                                                                      >  2:420785/1‑251 (MQ=255)
                                                                                                                   cGTTGCGCCCCATAGGCTGCGGTAAATAACCCACTTCCCGGGCCGCTTTTTCAACCCGATTACGGGTGGCCTGGGAGACTTTATCGGGATTCATTAATGCTCGGGAGACGGTCGCTGTAGAGACTTTTTCCTTGAGGGCAACGTCTTTCATGGTCGCGGCAGTTTCCTGCTTCTTCGCTTTCACACTCACTCCTCGCCTGGCACGTCAGGCGTACTACATCCATGTTTACTTCACATCGGCAACAttttt                                                                                                      >  1:150678/1‑250 (MQ=255)
                                                                                                                    gTTGCGCCCCATAGGCTGCGGTAAATAACCCACTTCCCGGGCCGCTTTTTCAACCCGCTTACGGGTGGCCTGGGAGACTTTATCGGGATTCATTAATGCTCGGGAGACGGTCGCTGTAGAGACTTTTTCCTTGAGGGCAACGTCTTTCATGGTCGCGGCAGTTTCCTGCTTCTTCGCTTTCACACTCACTCCTCGCCTGGCACGTCAGGCGTACTACATCCATGTTTACTTCACATCGGCAACCttttt                                                                                                      >  1:240256/1‑249 (MQ=255)
                                                                                                                     ttGCGCCCCATAGGCTGCGGTAAATAACCCACTTCCCGGGCCGCTTTTTCAACCCGATTACGGGTGGCCTGGGAGACTTTATCGGGATTCATTAATGCTCGGGAGACGGTCGCTGTAGAGACTTTTTCCTTGAGGGCAACGTCTTTCATGGTCGCGGCAGTTTCCTGCTTCTTCGCTTTCACACTCACTCCTCGCCTGGCACGTCAGGCGTACTACATCCATGTTTACTTCACATCGGCAACATTTTTAGc                                                                                                   >  1:36403/1‑251 (MQ=255)
                                                                                                                                      cGGTAAATAACCCACTTCCCGGGCCTCTTTGTCAACCCGATTACGGGTGGCTTGGGAGACGTTATCGGGATTCAGTAATGCTCGGGAGCCGGTCGCTGTAGAGACTTTTTCCTTGAGGGCAACGTCTTTCATGGTCGCGGCAGTTTCCTGCTTCTTCGCTTTCACACTCACTCCTCGCCCGGCACGTCAGGCGTACTACATCCATGTTTACTTCACATCGGCAACATTTTTAGCAGATAGCGCGTGAAAAc                                                                                  <  2:347903/251‑1 (MQ=255)
                                                                                                                                       ggTAAATAAACCCCTTCCCGGGCCCTTTTTTCAACCCCATTCCGGGTGGCCTGGGAGACTTTATCGGGATTCATTTATGCTCGGGAGACGGTCGCTGTAGAGACTTTTTCCTTGAGGGCAACGTCTTTCATGGTCGCGGCAGTTTCCTGCTTCTTCGCTTTCACACTCACTCCTCGCCTGGCACGTCAGGCGTACTACATCCATGTTTACTTCACATCGGCAACATTTTTAGCAGATAGCGCGTGAAAACg                                                                                 <  1:412525/251‑1 (MQ=255)
                                                                                                                                            aTAACCCACTTCCCGGGCCGCTTTTTCAACCCGATTACGGGTGGCCTGGGAGACTTTATCGGGATTCATTAATGCTCGGGAGACGGTCGCTGTAGAGACTTTTTCCTTGAGGGCAACGTCTTTCATGGTCGCGGCAGTTTCCTGCTTCTTCGCTTTCACACTCACTCCTCGCCTGGCACGTCAGGCGTACTACATCCATGTTTACTTCACATCGGCAACATTTTTAGCAGATAGCGCGTGAAAACGGTTAc                                                                            >  1:315698/1‑251 (MQ=255)
                                                                                                                                             tAACCCACTTCCCGGGCCGCTTTTTCAACCCGATTACGGGTGGCCTGGGAGACATTATCGGGATTCATTAATGCTCGGGAGACGGTCGCAGTAGAGACTTTCTCCTTGAGGGCAACGTCTTTCATGGTCGCGGAAGTTTCCTGCTTCTTCGCTTTCACACTCACTCCTCGCCTGGCACGTCAGGCGTACTACATCCATGTTTACTTCACATCGGCAACCTTTTTAGCAGATAGCGCGTGAAAACGGTTAc                                                                            >  1:236356/1‑250 (MQ=255)
                                                                                                                                                 ccACTTCCCGGGCCGCGTTTTCAACCCGATTACGGGTGGCCTGGGAGACTTTATCGGGATTCATTAATGCTCGGGAGACGGTCGCTGTAGAGACTTTTTCCTTGAGGGCAACGTCTTTCATGGTCGCGGCAGTTTCCTGCTTCTTCGCTTTCACACTCACTCCTCGCCTGGCACGTCAGGCGTACTACATCCATGTTTACTTCACATCGGCAACATTTTTAGCAGATAGCGCGTGAAAACGGTTACAGAAt                                                                      <  1:430032/251‑1 (MQ=255)
                                                                                                                                                     ttaccGGGCCGATTTTTCACCCCGATTACGGGTGGCCTGGGAGACTTTATCGGGATTCATTAATGCTCGGGAGACGGTCGCCGTCGAGACTTTTTCCTTGAGGGCAACGTCTTTCATGGTCGCGGACGTATACTGCTTCTTCGCTTTCACACTCACTCCTCGCCTGGCACGTCAGGCGTACTACATCCAAGTTTACTTCACATCGGCAACATATTTAGAAGATAGCGCGTGCAAACGTTACCAGAAttatc                                                                  >  2:455229/4‑248 (MQ=255)
                                                                                                                                                          gggCCGCTTTTTCAACCCGATTACGGGTGGCCTTGGAGACTTTATCGGGATTCATTAATGCTCGGGAGACGGTCGCTGTAGAGACTTTTTCCTTGAGGGCAACGTCTTTCATGGTCGCGGCAGTTTCCTGCTTCTTCGCTTTCACACTCACTCCTCGCCTGGCACGTCAGGCGTACTACATCCATGTTTACTTCACATCGGCAACATTTTTAGCAGATAGCGCGTGAAAACGGTTACAGAATTTTCATGaa                                                             <  1:234202/251‑1 (MQ=255)
                                                                                                                                                                  ttgtCAACCCGATTACGGGTGGCCTGGGAGACTTTATCGGGATTCATTAATGCTCGAGAGACGGTCGCTGAAGAGACTTTTGCCTTGGGGGCAACGTCTTTCATGGTCGCGGCAGTTTCCTTCTTGTTCGCTTTCACACTCAATCCTCGCCTGGCACGGCCGGCGGACTACTTCCATGTTTAATTCACATAGGCAACATTTTTGGCTGCTAGGGCGTAAAAAAGGTTAAAACATATTCATGAAAAgtgtat                                                     >  2:17891/4‑249 (MQ=255)
                                                                                                                                                                            gATTACGGGTGGCCTGGGAGACTTTATCGGGATTCATTAATGCTCGGGAGACGGTCGCTGTAGAGACTTTTTCCTTGAGGGCAACGTCTTTCATGGTCGCGGCAGTTTCCTGCTTCTTCGCTTTCACACTCACTCCTCGCCTGGCACGTCAGGCGTACTACATCCATGTTTACTTCACATc                                                                                                                  <  2:255815/181‑1 (MQ=255)
                                                                                                                                                                            gATTACGGGTGGCCTGGGAGACTTTATCGGGATTCATTAATGCTCGGGAGACGGTCGCTGTAGAGACTTTTTCCTTGAGGGCAACGTCTTTCATGGTCGCGGCAGTTTCCTGCTTCTTCGCTTTCACACTCACTCCTCGCCTGGCACGTCAGGCGTACTACATCCATGTTTACTTCACATc                                                                                                                  >  1:255815/1‑181 (MQ=255)
                                                                                                                                                                                 cGGGTGGCCTGGGAGTCTTTATCGGTATTCATTAATGCTCGGGAGACGGTCGCTGTAGAGACTTTTTCCTTGAGGGCAACGTCTTTCATTGTCGCGTCAGTTTCCTGCTTCTTCGCTTTCACACTCACTCCTCGCCTGGCACGTCAGGCGTACTACATCCATGTTTACTTCTCATCGGCAACATTTTTAGCAGATAGCGCGTGAAAACTGTTACAGAATTTTCATGAAAAGTGTGATGAATATTGAAtttt                                      <  2:396756/251‑1 (MQ=255)
                                                                                                                                                                                 cGGGTGGCCTGGGAGACTTTATCTGGATTCATTAATGCTCGGGAGACGGTCGCTGTAGAGACTTTTTCCTTGAGGGCAACGTCTTTCATGGTCGCGGCAGTTTCCTGCTTCTTCGCTTTCACACTCACTCCTCGCCTGGCACGTCAGGCGTACTACATCCATGTTTACTTCACATCGGCAACATTTTTAGCAGATAGCGCGTGAAAACGGTTACAGAATTTTCATGAAAAGTGTGATGAATATTGAAtttt                                      <  1:319822/251‑1 (MQ=255)
                                                                                                                                                                                  gggTGGCCTGGGAGACTTTATCGGGATTCATTAATGCTCGGGAGACGGTCGCTGTAGAGACTTTTTCCTTGAGGGCAACGTCTTTCATGGTCGCGGCAGTTTCCTGCTTCTTCGCTTTCACACTCACTCCTCGCCTGGCACGTCAGGCGTACTACATCCATGTTTACTTCACATCGGCAACATTTTTAGCAGATAGCGCGTGAAAACGGTTACAGAATTTTCATGAAAAGTGTGATGAATATTGAAttttt                                     <  1:198156/251‑1 (MQ=255)
                                                                                                                                                                                           gggAGTCTTTATCGGGATTCATTAATGCTCGGGAGACGGTCGCTGTAGAGACTTTTTCCTTGAGGGCAACGTCTTTCATGGTCGCGGCAGTTTCCTGCTTCTTCGCTTTCACACTCACTCCTCGCCTGGCACGTCAGGCGTACTACATCCATGTTTACTTCACATCGGCAACATTTTTAGCAGATAGCGCGTGAAAACGGTTACAGAATTTTCATGAAAAGTGTGATGAATATTGAATTTTTCGATCCg                              <  2:370906/249‑1 (MQ=255)
                                                                                                                                                                                            ggAGACTTTATCGGGATTCATTAATGCTCGGGAGACGGTCGCTGTAGAGACTTTTTCCTTGAGGGCAACGTCTTTCATGGTCGCGGCAGTTTCCTGCTTCTTCGCTTTCACACTCACTCCTCGCCTGGCACGTCAGGCGTACTACATCCATGTTTACTTCACATCGGCAACATTTTTAGCAGATAGCGCGTGAAAACGGTTACAGAATTTTCATGAAAAGTGTGATGAATATTGAATTTTTCGATCCGCCt                           <  2:97238/251‑1 (MQ=255)
                                                                                                                                                                                              gtaCTTTATCGGGATTCATTAATGCTCGGTAGTCGGTCGCTGTAGAGTCTTTTTCATTGAGGGCAACGTCTTTCATGGTTGCGGCAGTGTCCTGCTTCTTCGCTTTCACACTCACTCCTCGCCTGGCACGTCAGGCGTACTACAACCATGTTTACTTCACATCGGCAACATTTTTAGCAGATAGCGCGTGAAACCGGTTACAGAATTTTCATGAAAAGTGTGATGAATATTGAATTTTTCGATCCGCCTCg                         <  2:86092/249‑1 (MQ=255)
                                                                                                                                                                                                       cGGGATTCATTAATGCTCGGGAGACGGTCGCTGTAGAGACTTTTTCCTTGAGGGCAACGTCTTTCATGGTCGCGGCAGTTTCCTGCTTCTTCGCTTTCACACTCACTCCTCGCCTGGCACGTCAGGCGTACTACATCCATGTTTACTTCACATCGGCAACATTTTTAGCAGATAGCGCGTGAAAACGGTTACAGAATTTTCATGAAAAGTGTGATGAATATTGAATTTTTCGATCCGCCTCGCATCGTg                  >  1:203617/1‑249 (MQ=255)
                                                                                                                                                                                                             tCATTAATGCTCGGGAGACGGTCGCTGTAGAGACTTTTTCCTTGAGGGCAACGTCTTTCATGGTCGCGGCAGTTTCCTGCTTCTTCGCTTTCACACTCACTCCTCGCCTGGCACGTCAGGCGtactac                                                                                                                                      >  1:97734/1‑128 (MQ=255)
                                                                                                                                                                                                             tCATTAATGCTCGGGAGACGGTCGCTGTAGAGACTTTTTCCTTGAGGGCAACGTCTTTCATGGTCGCGGCAGTTTCCTGCTTCTTCGCTTTCACACTCACTCCTCGCCTGGCACGTCAGGCGtactac                                                                                                                                      <  2:97734/128‑1 (MQ=255)
                                                                                                                                                                                                              cATTAATGCTCGGGAGACGGTCGCTGTAGAGACTTTTTCCTTGAGGGCAACGTCTTTCATGGTCGCGGCAGTTTCCTGCTTCTTCGCTTTCACACTCACTCCTCGCCTGGCACGTCAGGCGTACTACATCCATGTTTACTTCACATCGGCAACATTTTTAGCAGATAgcgc                                                                                          >  1:439908/1‑171 (MQ=255)
                                                                                                                                                                                                              cATTAATGCTCGGGAGACGGTCGCTGTAGAGACTTTTTCCTTGAGGGCAACGTCTTTCATGGTCGCGGCAGTTTCCTGCTTCTTCGCTTTCACACTCACTCCTCGCCTGGCACGTCAGGCGTACTACATCCATGTTTACTTCACATCGGCAACATTTTTAGCAGATAgcgc                                                                                          <  2:439908/171‑1 (MQ=255)
                                                                                                                                                                                                              cATTAATGCTCGGGAGACGGTCGCTGTAGAGACTTTTTCCTTGAGGGCAACGTCTTTCATGGTCGCGGCAGTTTCCTGCTTCTTCGCTTTCACACTCACTCCTCGCCTGGCACGTCAGGCGTACTACATCCATGTTTACTTCACATCGGCAACATTTTTAGCAGATAGCGCGTGAAAACGGTTACAGAATTTTCATGAAAAGTGTGATGAATATTGAATTTTTCGATCCGCCTCGCATCGTGAGCGGttt          >  1:306498/1‑250 (MQ=255)
                                                                                                                                                                                                              cATTAATGCTCGGGAGACGGTCGCTGTAGAGACTTTTTCCTTGAGGGCAACGTCTTTCATGGTCGCGGCAGTTTCCTGCTTCTTCGCTTTCACACTCACTCCTCGCCTGGCACGTCAGGCGTACTACATCCATGTTTACTTCACATCGGCAACATTTTTAGCAGATAGCGCGTGAAAACGGTTACAGAATTTTCATGAAAAGTGTGATGAATATTGAATTTTTCGATCCGCCTCGCATCGTGAGCGGttt          >  1:421257/1‑250 (MQ=255)
                                                                                                                                                                                                                     gCTCGGGAGACGGTCGCTGTAGAGACTTCTTCCTTGATGGCAACGTCTTTCATTGTCGCGGCAGTTTCCTGCTTCTTCGCTTTCACACTCACTCCTCGCCTGGCACGGCAGGCGTACTACATCCATGTTTACTTCACATCGGCAACATTTTTAGCAGATAGCGCGTGAAAACGGTTACAGAATTTTCATGAAAAGTGTGATGAATATTGAATTTTTCGATCCGCCTCGCATCGTGAGCGGTTTATCCCTCa  <  2:224490/251‑1 (MQ=255)
                                                                                                                                                                                                                                                   |                                                                                                                                                                                                                             
ACCAGATATCCGTGATTTGCCGCCGTAACTTCGATACCGCGAATAATTTCGCTAAAGAAGGGATCGCAGATATCCGGGACAATCACCAGAATGGTGCGGGATTCATTACGCTTGACGTTGCGCCCCATAGGCTGCGGTAAATAACCCACTTCCCGGGCCGCTTTTTCAACCCGATTACGGGTGGCCTGGGAGACTTTATCGGGATTCATTAATGCTCGGGAGACGGTCGCTGTAGAGACTTTTGCCTTGAGGGCAACGTCTTTCATGGTCGCGGCAGTTTCCTGCTTCTTCGCTTTCACACTCACTCCTCGCCTGGCACGTCAGGCGTACTACATCCATGTTTACTTCACATCGGCAACATTTTTAGCAGATAGCGCGTGAAAACGGTTACAGAATTTTCATGAAAAGTGTGATGAATATTGAATTTTTCGATCCGCCTCGCATCGTGAGCGGTTTAACCCTCA  >  NC_000913/4124158‑4124621

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: