Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,333,520 C→T G144D (GGC→GAC)  basR ← DNA‑binding transcriptional dual regulator BasR

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,333,5200CT91.5% 134.2 / 2.2 47G144D (GGC→GAC) basRDNA‑binding transcriptional dual regulator BasR
Reads supporting (aligned to +/- strand):  ref base C (1/3);  new base T (27/16);  total (28/19)
Fisher's exact test for biased strand distribution p-value = 2.89e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.57e-01

CAGATTCAATTAGTTTTCCTCATTCGCGACCAGCATATAGCCAAAGCCGCGCACGGTGCGGATACGGGCTTTGCCCACTTTGTCGCGCAGATTGTGGATATGCACTTCCAGGGTGTTGGTCGAGGGTTCATTGTCCCAGTTATAGATGTCGTTGTAGAGAATTTCCCGATGCACCGGACTGCCTGCTTTGAGCATTAACCGTGACAGCAGAGCATATTCTTTGGGCGTCAGAATCAACTCTTCACCGCCCATCCATACCTGACGGCGACCCATGTTCAGCGTCAGATTGCCAACAATCAGCTCACTTTCGCCCTGATTATTATGGCGTCGTAGCAGGGCGCGGATACGGGCATGTAACTCTTCCAGCGCAAAAGGCTTCACCAGATAGTCGTCGGCACCGACATCCAGCCCGGCGATTTTGTCGGTCAGCGTATCGCGAGCGGTGAGGATCAGTACCGGCAGGGTATATTTTTTCTGCCGGATACGGGC  >  NC_000913/4333273‑4333761
                                                                                                                                                                                                                                                        |                                                                                                                                                                                                                                                  
cAGATTCAATTAGTTTTCCTCATTCGCGACCAGCATATAGCCAAAGCCGCGCACGGTGCGGATACGGGCTTTGCCCACTTTGTCGCGCAGATTGTGGATATGCACTTCCAGGGTGTTGGTCGAGGGTTCATTGTCCCAGTTATAGATGTCGTTGTAGAGAATTTCCCGATGCACCGGACTGCCTGCTTTGAGCATTAACCGTGACAGCAGAGCATATTCTTTGGGCGTCAGAATCAACTCTTCACCGTcca                                                                                                                                                                                                                                                 <  1:222856/251‑1 (MQ=255)
  gATTCAATTAGTTTTCCTCATTCGCGACCAGCATATAGCCAAAGCCGCGCACGGTGCGGATACGGGCTTTGCCCACTTTGTCGCGCAGATTGTGGATATGCACTTCCAGGGTGTTGGTCGAGGGTTCATTGTCCCAGTTATAGATGTCGTTGTAGAGAATTTCCCGATGCACCGGACTGCCTGCTTTGAGCATTAACCGTGACAGCAGAGCATATTCTTTGGGCGTCAGAATCAACTCTTCACCGTccatc                                                                                                                                                                                                                                               >  2:105594/1‑251 (MQ=255)
    ttCAATTAGTTTTCCTCATTCGCGACCAGCATATAGCCAAAGCCGCGCACGGTGCGGATACGGGCTTTGCCCACTTTGTCGCGCAGATTGTGGATATGCACTTCCAGGGTGTTGGTCGAGGGTTCATTGTCCCAGTTATAGATGTCGTTGTAGAGAATTTCCCGATGCACCGGACTGCCTGCTTTGAGCATTAACCGTGACAGCAGAGCATATTCTTTGGGCGTCAGAATCAACTCTTCACCGTccatcc                                                                                                                                                                                                                                              >  1:165942/1‑250 (MQ=255)
       aaTTAGTTTTCCTCATTCGCGACCAGCATATAGCCAAAGCCGCGCACGTTGCGGATACGGGCTTTGCCCACTTTGTCGCGCAGATTGTGGATATGCACTTCCAGGGTGTTGGTCGAGGGTTCATTGTCCCAGTTATAGATGTCGTTGTAGAGAATTTCCCGATGCACCGGACTGCCTGCTTTGAGCATTAACCGTGACAGCCGAGCATATTCTTTGGGCGTCAGAATCAACTCTTCACCGTCCATCCATAc                                                                                                                                                                                                                                          >  2:110020/1‑251 (MQ=255)
       aaTTAGTTTTCCTCATTCGCGACCAGCATATAGCCAAAGCCGCGCACGGTGCGGATACGGGCTTTGCCCACTTTGTCGCGCAGATTGTGGATATGCACTTCCAGGGTGTTGGTCGAGGGTTCATTGTCCCAGTTATAGATGTCGTTGTAGAGAATTTCCCGATGCACCGGACTGCCTGCTTTGAGCATTAACCGTGACAGCAGAGCATATTCTTTGGGCGTCAGAATCAACTCTTCACCGTCCATCCATAc                                                                                                                                                                                                                                          >  1:305027/1‑251 (MQ=255)
          tAGTTTTCCTCATTCGCGACCAGCATATAGCCAAAGCCGCGCACGGTGCGGATACGGGCTTTGCCCACTTTGTCGCGCAGATTGTGGATATGCACTTCCAGGGTGTTGGTCGAGGGTTCATTGTCCCAGTTATAGATGTCGTTGTAGAGAATTTCCCGATGCACCGGACTGCCTGCTTTGAGCATTAACCGTGACAGCAGAGCATATTCTTTGGGCGTCAGAATCAACTCTTCACCGTCCATCCATACCt                                                                                                                                                                                                                                        <  1:206198/250‑1 (MQ=255)
                 ccTCATTCGCGACCAGCATATAGCCAAAGCCGCGCACGGTGCGGATACGGGCTTTGCCCACTTTGTCGCGCAGATTGTGGATATGCACTTCCAGGGTGTTGGTCGAGGGTTCATTGTCCCAGTTATAGATGTCGTTGTAGAGAATTTCCCGATGCACCGGACTGCCTGCTTTGAGCATTAACCGTGACAGCAGAGCATATTCTTTGGGCGTCAGAATCAACTCTTCACCGCCCATCCATACCTGACGGCg                                                                                                                                                                                                                                 <  2:186789/250‑1 (MQ=255)
                            aCCAGCATATAGCCAAAGCCGCGCACGGTGCGGATACGGGCTTTGCCCACTTTGTCGCGCAGATTGTGGATATGCACTTCCAGGGTGTTGGTCGAGGGTTCATTGTCCCAGTTATAGATGTCGTTGTAGAGAATTTCCCGATGCACCGGACTGCCTGCTTTGAGCATTAACCGTGACAGCAGAGCATATTCTTTGGGCGTCAGAATCAACTCTTCACCGTCCATCCATACCTGACGGCGACCCATGTTc                                                                                                                                                                                                                       >  1:67404/1‑249 (MQ=255)
                                               cgcgCACGGTGCGGATACGGGCTTTGCCCACTTTGTCGCGCAGATTGTGGATATGCACTTCCAGGGTGTTGGTCGAGGGTTCATTGTCCCAGTTATAGATGTCGTTGTAGAGAATTTCCCGATGCACCGGACTGCCTGCTTTGAGCATTAACCGTGACAGCAGAGCATATTCTTTGGGCGTCAGAATCAACTCTTCACCGTCCATCCATACCTGACGGCGACCCATGTTCAGCGTCAGATTGCCAACAATc                                                                                                                                                                                                 <  1:147560/251‑1 (MQ=255)
                                                      ggTGCGGATACGGGCTTGGCCCACTTTGTCGCGCAGATTGGGGATATGCACTTCCAGGGTGTTGGTCGAGGGTTCATTGTCCCAGTTATAGATGTCGTTGTAGAGAATTTCCCGATGCACCGGACTGCCTGCTTTGAGCATTAACCGTGACAGCAGAGCATATTCTTTGGGCGTCAGAATCAACTCTTCACCGCCCATCCATCCCTGACGGCGACCCATGTTCAGCGTCAGATTGCCAACAATCAGCTCAc                                                                                                                                                                                          <  2:415263/251‑1 (MQ=255)
                                                        tGCGGATACGGGCTTTGCCCACTTTGTCGCGCAGATTGTGGATATGCACTTCCAGGGTGTTGGTCGAGGGTTCATTGTCCCAGTTATAGATGTCGTTGTAGAGAATTTCCCGATGCACCGGACTGCCTGCTTTGAGCATTAACCGTGACAGCAGAGCATATTCTTTGGGCGTCAGAATCAACTCTTCACCGTCCATCCATACCTGACGGCGACCCATGTTCAGCGTCAGATTGCCAACAATCAGCTCACtt                                                                                                                                                                                        >  1:211237/1‑251 (MQ=255)
                                                                        gCCCACTTTGTCGCGCAGATTGTGGATATGCACTTCCAGGGTGTTGGTCGAGGGTTCATTGTCCCAGTTATAGATGGCGTTGTAGAGAATTTCCCGATGCACCGGACTGCCTGCTTTGAGCATTAACCGTGACAGCAGAGCATATTCTTTGGGCGTCAGAATCAACTCTTCACCGTCCATCCATACCTGACGGCGACCCATGTTCAGCGTCAGATTGCCAACAATCAGCTCACTTTCGCCCTGAttatta                                                                                                                                                                         <  1:221936/250‑1 (MQ=255)
                                                                                tGTCGCGCAGATTGTGGATATGCACTTCCAGGGTGTTGGTCGAGGGTTCATTGTACCAGTTATAGATGACGTTGTAGAGAATTTCCCGCTGCACCGGACTGCCTGCTTTGCGCATTAACCGTGACAGCAGAGCATATTCTTTGGGCGTCAGAATCAACTCTTCACCGCCCATCCATACCTGCCGGCGACCCATGTTCAGCGTCAGATTGCCAACAATCAGCTCACTTTCGCCCTGATTATTATGGcgtcgt                                                                                                                                                                >  1:391269/1‑251 (MQ=255)
                                                                                     cgcAGATTGTGGATATGCACTTCCAGGGTGTTGGTCGAGGGTTCATTGTCCCAGTTATAGATGTCGTTGTAGAGAATTCCCCGCTGCACCGGACTGCCTGCTTTGAGCATTAACCGTGACAGCAGAGCATATTCTGTGGGCGTCAGAATCAACTCTTCACCGTCCATCCATACCTGACGGCGACCCATGTTCAGCGTCAGATTGCCAACAATCAGCTCACTTTCGCACTGATTATTATGGCGTCGTAGCAg                                                                                                                                                           <  2:145431/251‑1 (MQ=255)
                                                                                      gcAGATTGTGGATATGCACTTCCAGGGTGTTGGTCGAGGGTTCATTGTCCCAGTTATAGATGTCGTTGTAGAGAATTTGCCGATGCACCGGACTGCCTGCTTTGTGCATTAATCGAGACAGCAGAGCATATTCTTTGGGCGTCAGAATCAAATCTTAACAGTCCCTCCATACCTAACGGCGTCCCATGTTCAGCGTAAGATTGACACAATTCAGATCACTTTCGTCCTGTTTATTATAGCGTCGTAGccgg                                                                                                                                                          >  2:294916/1‑248 (MQ=255)
                                                                                        aGATTGTGGATATGCACTTCCAGGGTGTCGGTAGAGGGTTCATTGTCCCAGTTATAGATGTCGTTGTAGAGAATTTCCCGATGCACCGGACTGCCTGCTTTGAGCATTAACCGTGACAGCAGAGCATATTCTTTGGGCGTCAGAATCAACTCTTCACCGCCCATCCATACCTGACGGCGACCCATGTTCAGCGTCAGATTGCCAACAATCAGCTCACTTTCGCCCTGATTATTATGGCGTCGTAGCAGGgc                                                                                                                                                        <  1:208630/251‑1 (MQ=255)
                                                                                             gtgGATATGCACTTCCAGGGTGTTGGTCGAGGGTTCATTGTCCCAGTTATAGATGTCGTTGTAGAGAATTTCCCGATGCACCGGACTGCCTGCTTTGAGCATTAACCGTGACAGCAGAGCATATTCTTTGGGCGTCAGAATCAACTCTTCACCGTCCATCCATACCTGACGGCGACCCATGTTCAGCGTCAGATTGCCAACAATCAGCTCACTTTCGCCCTGATTATTATGGCGTCGTAGCAGGGCGCgg                                                                                                                                                    >  1:409463/1‑250 (MQ=255)
                                                                                                 atatGCACTTCCAGGGTGTTGGTCGAGGGTTCATTGTCCCAGTTGTAGATGTCGTTGTAGAGAATTTCGCGATGCACCGGACTGCCTGCTTTGAGCATTAAACGTGACAGCAGAGCATATTCTTTGGGCGTCAGAATCAACTCTTCACCGTCCATCCATACCTGACGGCGACCCATGTTCAGCGTCAGATTGCCAACAATCAGCTCACTTTCGCCCTGATTATTATGGCGTCGTAGCAGGGCGCGGATagg                                                                                                                                               >  2:234128/1‑249 (MQ=255)
                                                                                                    tGCACTTCCAGGGTGTTGGTCGAGGGTTCATTGTCCCAGTTATAGATGTCGTTGTAGAGAATTTCCCGATGCACCGGACTGCCTGCTTTGAGCATTAACCGTGACAGCAGAGCATATTCTTTGGGCGTCAGAATCAACTCTTCACCGTCCATCCATACCTGACGGCGACCCATGTTCAGCGTCAGATTGCCAACAATCAGCTCACTTTCGCCCTGATTATTATGGCGTCGTAGCAGGGCTAGGATACGGGc                                                                                                                                            >  1:162880/1‑251 (MQ=255)
                                                                                                    tGCACTTCCAGGGGGTTGGTCGAGGGTTCATTGTCCCAGTTATAGATGTCGTTGTAGAGAATTTCCCGATGCACCGGACTGCCTGCTTTGAGCATTAACCGTGACAGCAGAGCATATTCTTTGGGCGTCAGAATCAACTCTTCACCGTCCATCCATACCTGACGGCGACCCATGTTCAGCGTCAGATTGCCAACAATCAGCTCACTTTCGCCCTGATTATTATGGCGTCGTAGCAGGGCGCGGATACGGGc                                                                                                                                            >  2:354271/1‑251 (MQ=255)
                                                                                                             aGGGTGTTGGTCGAGGGTTCATTGTCCCAGTTATAGATGTCGTTGTAGAGAATTTCCCGATGCACCGGCCTGCCTGCTTTGAGCATTAACCGTGACAGCAGAGCATATTCTTTGGGCGTCAGAATCAACTCTTCACCGTACATCCATACCTGACGGCGCCCCATGTTCAGCGTCAGATTGCCAACAATCAGCTCACttt                                                                                                                                                                                       >  2:218326/1‑199 (MQ=255)
                                                                                                             aGGGTGTTGGTCGAGGGTTCATTGTCCCAGTTATAGATGTCGTTGTAGAGAATTTCCCGATGCACCGGACTGCCTGCTTTGAGCATTAACCGTGACAGCAGAGCATATTCTTTGGGCGTCAGAATCAACTCTTCACCGTCCATCCATACCTGACGGCGACCCATGTTCAGCGTCAGATTGCCAACAATCAGCTCACttt                                                                                                                                                                                       <  1:218326/199‑1 (MQ=255)
                                                                                                                                cATTGTCCCAGTTATAGATGTCGTTGTAGAGAATTTCCCGATGCACCGGACTGCCTGCTTTGAGCATTAACCGTGACAGCAGAGCATATTCTTTGGGCGTCAGAATCAACTCTTCACCGTCCATCCATACCTGACGGCGACCCATGTTCAGCGTCAGATTGCCAACAATCAGCTCACTTTCGCCCTGATTATTATGGCGTCGTAGCAGGGCGCGGATACGGGCATGTa                                                                                                                                       <  1:399945/228‑1 (MQ=255)
                                                                                                                                cATTGTCCCAGTTATAGATGTCGTTGTAGAGAATTTCCCGATGCACCGGACTGCCTGCTTTGAGCATTAACCGTGACAGCAGAGCATATTCTTTGGGCGTCAGAATCAACTCTTCACCGTCCATCCATACCTGACGGCGACCCATGTTCAGCGTCAGATTGCCAACAATCAGCTCACTTTCGCCCTGATTATTATGGCGTCGTAGCAGGGCGCGGATACGGGCATGTa                                                                                                                                       >  2:399945/1‑228 (MQ=255)
                                                                                                                                    gTCCCAGTTATAGATGTCGTTGTAGAGAATTTCCCGATGCACCGGACTGCCTGCTTTGAGCATTAACCGTGACAGCCGAGCATATTCTTTGGGCGTCAGAATCAACTCTTCACCGTCCATCCATACCTGACGGCGACCCATGTTCAGCGTCAGATTGCCAACAATCAGCTCACTTTCGCCCTGATTATTATGGCGTCGTAGCAGGGCGCGGATACGGGCATGTCACTCTTCCAGCGCAAAAGGCTTCAcc                                                                                                             >  2:139285/1‑250 (MQ=255)
                                                                                                                                    gTCCCAGTTATAGATGTCGTTGTAGAGAATTTCCCGATGCACCGGACTGCCTGATTTGAGCATTAACCGTGACAGCAGAGCATATTCTTTGGGCGTCAGAATCAACTCTTCACCGTCCATCCATACCTGACGGCGACCCATGTTCAGCGTCAGATTGCCAACAATCAGCTCACTTTCGCCCTGATTATTATGGCGTCGTAGCAGGGCGCGGATACGGGCATGTAACTCTTCCAGCGCAAAAGGCTTCAcc                                                                                                             >  1:60109/1‑250 (MQ=255)
                                                                                                                                                        ggTAGAGAATTTCCCGCTGCACAGGACTGCCTGCTTTGAGCATTAACCGTTACAGCAGAGCATATTCTTTGGGCGTCAGAATCAACTCTTCACCGTCCTTCCATACCTGACGGCGACCCATGTTCAGCGTCAGATTGCCAACAATCAGCCCACTTTCGCCCTGATTATTCTGGCGTCGTAGCAGGGCGCTGATACGGGCATGTAACTCTTCCAGCGCACAAGGCTTCACCAGCTCGTAGGCGGCACCGaa                                                                                         >  2:69622/2‑249 (MQ=255)
                                                                                                                                                              gaATTTCCCGATTGCACGGCCTGCCTGCTTTGAGCATTACCCGTGTCATCAGAGCATATTCTTTGGGCGTCAGAATCAACTCTTCTCCGTCCATCCATACCTGACGGCGACCCATGTTCAGCGTCAGATTGCCAACAATCAGCTCACTTGCGCCCTGATTATTATGGCGTCGTAGCAGGGCGCGGATACGGGCATGTAACTCTTCCAGCGCAAAAGGCTTCACCAGATAGTCGTCGGCCCCGACATCCAGc                                                                                  <  2:68982/251‑1 (MQ=255)
                                                                                                                                                              aaaTTTCCCGATGCACCGGACTGCCTGCTTTGAGCATTAACCGTGACAGCAGAGCATATTCTTTGGGCGTCAGAATCAACTCTTCACCGTCCATCCATACCTGACGGCGACCCATGTTCAGCGTCAGATTGCCAACAATCAGCTCACTTTCGCCCTGATTATTATGGCGTCGTAGCAGGGCGCGGATACGGGCATGTAACTCTTCCAGCGCAAAAGGCTTCACCAGATAGTCGTCGGCACCGACATCCAGc                                                                                  <  1:46371/250‑1 (MQ=255)
                                                                                                                                                                    cccGATGCACCGGACTGCCTGCTTTGAGCATTAACCGTGACAGCAGAGCATATTCTTTGGGCGTCAGAATCAACTCTTCACCGTCCATCCATACCTGACGGCGACCCATGTTCAGCGTCAGATTGCCAACAATCAGCTCACTTTCGCCCTGATTATTATGGCGTCGTAGCAGGGCGCGGATACGGGCATGTAACTCTTCCAGCGCAAAAGGCTTCACCAGATAGTCGTCGGCACCGACATCCAGCCCGGCg                                                                            >  2:341257/1‑251 (MQ=255)
                                                                                                                                                                              ccGGACTGCCTGCTTTGAGCATTAACCGTGACAGCAGAGCATATTCTTTGGGCGTCAGAATCAACTCTTCACCGTCCATCCATACCTGACGGCGACCCATGTTCAGCGTCAGATTGCCAACAATCAGCTCACTTTCGCCCTGATTATTATGGCGTCGTAGCAGGGCGCGGATACGGGCATGTAACTCTTCCAGCGCAAAAGGCTTCACCAGATAGTCGTCGGCACCGACATCCAGCCCGTCGATTTTGTCg                                                                   <  1:98362/251‑1 (MQ=255)
                                                                                                                                                                                        tgcTTTGAGCATTAACCGTGACAGCAGAGCATATTATTTGGGCGTCAGAATCAACTCTTCACCGTCCATCCATACCTGACGGCGACCCATGTTCAGCGTCAGATTGCCAACAATCAGCTCACTTTCGCCCTGATTATTATGGCGTCGTAGCAGGGCGCGGATACGGGCATGTAACTCTTCCAGCGCAAAAGGCTTCACCAGATAGTCGTCGGCACCGCCATCCAGCCCGGCGATTTTGTCGGTCAGCGTAt                                                         >  1:198870/1‑251 (MQ=255)
                                                                                                                                                                                         gcTTTGAGCATTAACCGTGACAGCAGAGCATATTCTTTGGGCGTCAGAATCAACTCTTCACCGTCCATCCATACCTGACGGCGACCCATGTTCAGCGTCAGATTGCCAACAATCAGCTCACTTTCGCCCTGATTATTATGGCGTCGTAGCAGGGCGCGGATACGGGCATGTAACTCTTCCAGCGCAAAAGGCTTCACCAGATAGTCGTCGGCACCGACATCCAGCCCGGCGATTTTGTCGGTCAGCGTATc                                                        <  1:105594/251‑1 (MQ=255)
                                                                                                                                                                                            ttGAGCATTAACCGTGACAGCAGAGCATATTCTTTGGGCGTCAGAATCAACTCTTCACCGTCCATCCATACCTGACGGCGACCCATGTTCAGCGTCAGATTGCCAACAATCAGCTCACTTTCGCCCTGATTATTATGGCGTCGTAGCAGGGCGCGGATACGGGCATGTAACTCTTCCAGCGCAAAAGGCTTCACCAGATAGTCGTCGGCACCGACATCCAGCCCGGCGATTTTGTCGGTCAGCGTATcgcg                                                     >  2:112012/1‑251 (MQ=255)
                                                                                                                                                                                            ttGAGCATTAACCGTGACAGCAGAGCATATTCTTTGGGCGTCAGAATCAACTCTTCACCGTCCATCCATACCTGACGGCGACCCATGTTCAGCGTCAGATTGCCAACAATCAGCTCACTTTCGCCCTGATTATTATGGCGTCGTAGCAGGGCGCGGATACGGGCATGTAACTCTTCCAGCGCAAAAGGCTTCACCAGATAGTCGTCGGCACCGACATCCAGCCCGGCGATTTTGTCGGTCAGCGTATcg                                                       <  1:110020/249‑1 (MQ=255)
                                                                                                                                                                                                 cATTAACCGTGACAGCAGAGCATATTCTTTGGGCGTCAGAATCAACTCTTCACCGTCCATCCATACCTGACGGCGACCCATGTTCAGCGTCAGATTGCCAACAATCAGCTCACTTTCGCCCTGATTATTATGGCGTCGTAGCAGGGCGCGGATACGGGCATGTAACTCTTCCAGCGCAAAAGGCTTCACCAGATAGTCGTCGGCACCGACATCCCACCCGGCGATTTTGTCGGTCAGCGTATCGCGAGCgg                                                >  2:371884/1‑251 (MQ=255)
                                                                                                                                                                                                    tAACCGTGACAGCAGAGCATATTCTTTGGGCGGCAGAATCAACTCCTCACCGTCCATCCATACCTGACGGCGACCCATGTTCAGCGTCAGATTGCCAACAATCAGCTCACTTTCGCCCTGATTATTATGGCGTCGTAGCAGGGCGCGGATACGGGCATGTAACTCTTCCAGCGCAAAAGTCTTCACCAGATAGTCGTCGGCAc                                                                                             <  1:168410/203‑1 (MQ=255)
                                                                                                                                                                                                    tAACCATGACAGCAGAGCCTATTCTTTGGGCGTCAGAATCAACTCTTCACCGTCCATCCATACCTGACGGCGACCCATGTTCAGCGTCAGATTTCCAACAATCCGCTCACTTTCGCCCTGATTATTATGGCTTCGTAACAGGGCGCGGATACGGGCATGTAACTCTTCCAGCGCAACAGGCTTCACCAGATAGTCGTCGGCAc                                                                                             >  2:168410/1‑203 (MQ=255)
                                                                                                                                                                                                                 agagCATATTCTTTGGGCGTCAGAATCAACTCTTCACCGTCCATCCATACCTGACGGCGACCCATGTTCAGCGTCAGATTGCCAACAATCAGCTCACTTTCGCCCTGATTATTATGGCGTCGTAGCAGGGCGCGGCTACGGGCATGTAACTCTTCCAGCGCAACAGGCTTCACCAGATAGACGTCGGCACCGACATCCAGCCCGGCGATTTTGTCGGTCAGCGTATCGCGAGCGGTTAGGTTCAGTACCgg                                >  2:417590/1‑251 (MQ=255)
                                                                                                                                                                                                                   agCATATTCACTGGGCGTCAGAATCAACTCTTCACCGTCCATCCATACCTGACGGCGACCCATGATCAACGTCAGATTGCCACCCATCAGCTCACTTTCGCCCTGATTATTATGGCGTCGCAGCAGGGCGCGGATACGCGCCTGTCACTCTTCCAGCGCAACAGGCTTCACCAGATAGCCGTCGGCACCGACATCCAGCCCGGCGATTTTGTCGGTCAGCGTCTCGCGAGCGGTGAGGATCAGTACCGGcg                              >  2:389715/1‑250 (MQ=255)
                                                                                                                                                                                                                     cATATTCTTTGGGCGTCAGAATCAACTCTTCCCCGTCCATCCATACCTGACGGCGACCCATGTTCAGCGTCAGATTGCCAACAATCAGCTCACTTTCGCCCTGATTATTATGGCGTCGTAGCAGGGCGCGGATACGGGCATGTAACTCTTCCAGCGCAAAAGGCTTCACCAGATAGTCGTCGGCACCGACATCCAGCCCGGCGATTTTGTCGGTCAGCGTATCGCGAGCGGTGAGGATCAGTACCGGCAgg                            <  1:167137/251‑1 (MQ=255)
                                                                                                                                                                                                                         ttCTTTGGGCATCAGAAGCAACTCTTCCACTTCCAGCCATGCCTTGCGGCGCCCTATGTTCAGCGTCAGATTGCCAACATTCAGCTCCCTTTCGCCCTTATTATTATGGCGTCGTAGCAGGGCGCGTATACGGGCATGTAACTCTTCCAGAGCAAAAGGCTTCACCAGATAGTCGTCGGCACCGACATCCAGCCCGGCGATTTTGTCGGTCAGCGTATCGCGAGCGGTGAGGATCAGTACCGGCAGGGta                         <  2:331349/250‑1 (MQ=255)
                                                                                                                                                                                                                              tGGGCTTCAGAATCAACTCTTCACCGTCCATCTATAACTGACGGCGCCCCATGTTCAGCGTCCGATTGCCAACAATCAGCTCCCTTGCGCCCTGATTATTCTGGCGTCGTAGCAGGGCGCGGATACGGGCATGTAACTCTTCCAGCGCAAAAGGCTTCACCATGTAGTCGTCGGCACCGACATCCGGCACTGCGATTTTGTCGGTCAGCGCATAGAGAGCGGTGAGGCCCGGGACCGGCAGGGGAAAtttt                   >  2:120912/1‑251 (MQ=255)
                                                                                                                                                                                                                                    tCAGAATCAACTCTTCACCGTCCATCCATACCTGACGGCGACCCATGCTCAGCGTCAGATTGCCAACAATCAGCTCACTTTCGCCCTGATTATTATGGCGTCGTAGCAGGGCGCGGATACGGGCATGTAACTCTTCCAGCGCAAAAGGCTTCACCAGATAGTCGTCGGCACCGACATCCAGCCCGGCGATTTTGTCGGTCAGCGTATCGCGAGAGGCGAGGATCAGTACCGGCAGGGTATATTTTTTATGc             >  2:8383/1‑251 (MQ=255)
                                                                                                                                                                                                                                            aaCTCTTCACCGTCCATCCATACCTGACGGCGACCCATGTTCAGCGTCAGATTGCCAACAATCAGCTCACTTTCGCCCTGatta                                                                                                                                                                            >  1:274652/1‑84 (MQ=255)
                                                                                                                                                                                                                                            aaCTCTTCACCGTCCATCCATACCTGACGGCGACCCATGTTCAGCGTCAGATTGCCAACAAGCAGCTCACTTTCGCCCTGatta                                                                                                                                                                            <  2:274652/84‑1 (MQ=255)
                                                                                                                                                                                                                                               tctTCACCGTCCATCCATACCTGACGGCGACCCATGTTCAGCGTCAGATTGCCAACAATCAGCTCACTTTCGCCCTGATTATTATGGCGTCGTAGCAGGGCGCGGATACGGGCATGTAACTCTTCCAGCGCAAAAGGCTTCACCAGATAGTCGTCGGCACCGACATCCAGCCCGGCGATTTTGTCGGTCAGCGTATCGCGAGCGGTGAGGATCAGTACCGGCAGGGTATATTTTTTCTGCCGGATACGGGc  >  1:248821/1‑251 (MQ=255)
                                                                                                                                                                                                                                                        |                                                                                                                                                                                                                                                  
CAGATTCAATTAGTTTTCCTCATTCGCGACCAGCATATAGCCAAAGCCGCGCACGGTGCGGATACGGGCTTTGCCCACTTTGTCGCGCAGATTGTGGATATGCACTTCCAGGGTGTTGGTCGAGGGTTCATTGTCCCAGTTATAGATGTCGTTGTAGAGAATTTCCCGATGCACCGGACTGCCTGCTTTGAGCATTAACCGTGACAGCAGAGCATATTCTTTGGGCGTCAGAATCAACTCTTCACCGCCCATCCATACCTGACGGCGACCCATGTTCAGCGTCAGATTGCCAACAATCAGCTCACTTTCGCCCTGATTATTATGGCGTCGTAGCAGGGCGCGGATACGGGCATGTAACTCTTCCAGCGCAAAAGGCTTCACCAGATAGTCGTCGGCACCGACATCCAGCCCGGCGATTTTGTCGGTCAGCGTATCGCGAGCGGTGAGGATCAGTACCGGCAGGGTATATTTTTTCTGCCGGATACGGGC  >  NC_000913/4333273‑4333761

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