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

CAGATTCAATTAGTTTTCCTCATTCGCGACCAGCATATAGCCAAAGCCGCGCACGGTGCGGATACGGGCTTTGCCCACTTTGTCGCGCAGATTGTGGATATGCACTTCCAGGGTGTTGGTCGAGGGTTCATTGTCCCAGTTATAGATGTCGTTGTAGAGAATTTCCCGAT‑GCACCGGACTGCCTGCTTTGAGCATTAACCGTGACAGCAGAGCATATTCTTTGGGCGTCAGAATCAACTCTTCACCGCCCATCCATACCTGACGGCGACCCATGTTCAGCGTCAGATTGCCAACAAT‑CAGCTCACTTTCGCCCTGATTATTATGGCGTCGTAGCAGGGCGCGGATACGGGCATGTAACTCTTCCAGCGCAAAAGGCTTCACCAGATAGTCGTCGGCACCGACATCCAGCCCGGCGATTTTGTCGGTCAGCGTATCGCGAGCGGTGAGGATCAGTACCGGCAGGGTATATTTTTTCTGCCGGATACGGGC  >  NC_000913/4333273‑4333761
                                                                                                                                                                                                                                                        |                                                                                                                                                                                                                                                  
cAGATTCAATTAGTTTTCCTCATTCGCGACCAGCATATAGCCAAAGCCGCGCACGGTGCGGATACGGGCTTTGCCCACTTTGTCGCGCAGATTGTGGATATGCACTTCCAGGGTGTTGGTCGAGGGTTCATTGTCCCAGTTATAGATGTCGTTGTAGAGAATTTCCCGAT‑GCACCGGACTGCCTGCTTTGAGCATTAACCGTGACAGCAGAGCATATTCTTTGGGCGTCAGAATCAACTCTTCACCGTcca                                                                                                                                                                                                                                                 <  1:222856/251‑1 (MQ=255)
  gATTCAATTAGTTTTCCTCATTCGCGACCAGCATATAGCCAAAGCCGCGCACGGTGCGGATACGGGCTTTGCCCACTTTGTCGCGCAGATTGTGGATATGCACTTCCAGGGTGTTGGTCGAGGGTTCATTGTCCCAGTTATAGATGTCGTTGTAGAGAATTTCCCGAT‑GCACCGGACTGCCTGCTTTGAGCATTAACCGTGACAGCAGAGCATATTCTTTGGGCGTCAGAATCAACTCTTCACCGTccatc                                                                                                                                                                                                                                               >  2:105594/1‑251 (MQ=255)
    ttCAATTAGTTTTCCTCATTCGCGACCAGCATATAGCCAAAGCCGCGCACGGTGCGGATACGGGCTTTGCCCACTTTGTCGCGCAGATTGTGGATATGCACTTCCAGGGTGTTGGTCGAGGGTTCATTGTCCCAGTTATAGATGTCGTTGTAGAGAATTTCCCGAT‑GCACCGGACTGCCTGCTTTGAGCATTAACCGTGACAGCAGAGCATATTCTTTGGGCGTCAGAATCAACTCTTCACCGTccatcc                                                                                                                                                                                                                                              >  1:165942/1‑250 (MQ=255)
       aaTTAGTTTTCCTCATTCGCGACCAGCATATAGCCAAAGCCGCGCACGTTGCGGATACGGGCTTTGCCCACTTTGTCGCGCAGATTGTGGATATGCACTTCCAGGGTGTTGGTCGAGGGTTCATTGTCCCAGTTATAGATGTCGTTGTAGAGAATTTCCCGAT‑GCACCGGACTGCCTGCTTTGAGCATTAACCGTGACAGCCGAGCATATTCTTTGGGCGTCAGAATCAACTCTTCACCGTCCATCCATAc                                                                                                                                                                                                                                          >  2:110020/1‑251 (MQ=255)
       aaTTAGTTTTCCTCATTCGCGACCAGCATATAGCCAAAGCCGCGCACGGTGCGGATACGGGCTTTGCCCACTTTGTCGCGCAGATTGTGGATATGCACTTCCAGGGTGTTGGTCGAGGGTTCATTGTCCCAGTTATAGATGTCGTTGTAGAGAATTTCCCGAT‑GCACCGGACTGCCTGCTTTGAGCATTAACCGTGACAGCAGAGCATATTCTTTGGGCGTCAGAATCAACTCTTCACCGTCCATCCATAc                                                                                                                                                                                                                                          >  1:305027/1‑251 (MQ=255)
          tAGTTTTCCTCATTCGCGACCAGCATATAGCCAAAGCCGCGCACGGTGCGGATACGGGCTTTGCCCACTTTGTCGCGCAGATTGTGGATATGCACTTCCAGGGTGTTGGTCGAGGGTTCATTGTCCCAGTTATAGATGTCGTTGTAGAGAATTTCCCGAT‑GCACCGGACTGCCTGCTTTGAGCATTAACCGTGACAGCAGAGCATATTCTTTGGGCGTCAGAATCAACTCTTCACCGTCCATCCATACCt                                                                                                                                                                                                                                        <  1:206198/250‑1 (MQ=255)
                 ccTCATTCGCGACCAGCATATAGCCAAAGCCGCGCACGGTGCGGATACGGGCTTTGCCCACTTTGTCGCGCAGATTGTGGATATGCACTTCCAGGGTGTTGGTCGAGGGTTCATTGTCCCAGTTATAGATGTCGTTGTAGAGAATTTCCCGAT‑GCACCGGACTGCCTGCTTTGAGCATTAACCGTGACAGCAGAGCATATTCTTTGGGCGTCAGAATCAACTCTTCACCGCCCATCCATACCTGACGGCg                                                                                                                                                                                                                                 <  2:186789/250‑1 (MQ=255)
                            aCCAGCATATAGCCAAAGCCGCGCACGGTGCGGATACGGGCTTTGCCCACTTTGTCGCGCAGATTGTGGATATGCACTTCCAGGGTGTTGGTCGAGGGTTCATTGTCCCAGTTATAGATGTCGTTGTAGAGAATTTCCCGAT‑GCACCGGACTGCCTGCTTTGAGCATTAACCGTGACAGCAGAGCATATTCTTTGGGCGTCAGAATCAACTCTTCACCGTCCATCCATACCTGACGGCGACCCATGTTc                                                                                                                                                                                                                       >  1:67404/1‑249 (MQ=255)
                                               cgcgCACGGTGCGGATACGGGCTTTGCCCACTTTGTCGCGCAGATTGTGGATATGCACTTCCAGGGTGTTGGTCGAGGGTTCATTGTCCCAGTTATAGATGTCGTTGTAGAGAATTTCCCGAT‑GCACCGGACTGCCTGCTTTGAGCATTAACCGTGACAGCAGAGCATATTCTTTGGGCGTCAGAATCAACTCTTCACCGTCCATCCATACCTGACGGCGACCCATGTTCAGCGTCAGATTGCCAACAAT‑c                                                                                                                                                                                                 <  1:147560/251‑1 (MQ=255)
                                                      ggTGCGGATACGGGCTTGGCCCACTTTGTCGCGCAGATTGGGGATATGCACTTCCAGGGTGTTGGTCGAGGGTTCATTGTCCCAGTTATAGATGTCGTTGTAGAGAATTTCCCGAT‑GCACCGGACTGCCTGCTTTGAGCATTAACCGTGACAGCAGAGCATATTCTTTGGGCGTCAGAATCAACTCTTCACCGCCCATCCATCCCTGACGGCGACCCATGTTCAGCGTCAGATTGCCAACAAT‑CAGCTCAc                                                                                                                                                                                          <  2:415263/251‑1 (MQ=255)
                                                        tGCGGATACGGGCTTTGCCCACTTTGTCGCGCAGATTGTGGATATGCACTTCCAGGGTGTTGGTCGAGGGTTCATTGTCCCAGTTATAGATGTCGTTGTAGAGAATTTCCCGAT‑GCACCGGACTGCCTGCTTTGAGCATTAACCGTGACAGCAGAGCATATTCTTTGGGCGTCAGAATCAACTCTTCACCGTCCATCCATACCTGACGGCGACCCATGTTCAGCGTCAGATTGCCAACAAT‑CAGCTCACtt                                                                                                                                                                                        >  1:211237/1‑251 (MQ=255)
                                                                        gCCCACTTTGTCGCGCAGATTGTGGATATGCACTTCCAGGGTGTTGGTCGAGGGTTCATTGTCCCAGTTATAGATGGCGTTGTAGAGAATTTCCCGAT‑GCACCGGACTGCCTGCTTTGAGCATTAACCGTGACAGCAGAGCATATTCTTTGGGCGTCAGAATCAACTCTTCACCGTCCATCCATACCTGACGGCGACCCATGTTCAGCGTCAGATTGCCAACAAT‑CAGCTCACTTTCGCCCTGAttatta                                                                                                                                                                         <  1:221936/250‑1 (MQ=255)
                                                                                tGTCGCGCAGATTGTGGATATGCACTTCCAGGGTGTTGGTCGAGGGTTCATTGTACCAGTTATAGATGACGTTGTAGAGAATTTCCCGCT‑GCACCGGACTGCCTGCTTTGCGCATTAACCGTGACAGCAGAGCATATTCTTTGGGCGTCAGAATCAACTCTTCACCGCCCATCCATACCTGCCGGCGACCCATGTTCAGCGTCAGATTGCCAACAAT‑CAGCTCACTTTCGCCCTGATTATTATGGcgtcgt                                                                                                                                                                >  1:391269/1‑251 (MQ=255)
                                                                                     cgcAGATTGTGGATATGCACTTCCAGGGTGTTGGTCGAGGGTTCATTGTCCCAGTTATAGATGTCGTTGTAGAGAATTCCCCGCT‑GCACCGGACTGCCTGCTTTGAGCATTAACCGTGACAGCAGAGCATATTCTGTGGGCGTCAGAATCAACTCTTCACCGTCCATCCATACCTGACGGCGACCCATGTTCAGCGTCAGATTGCCAACAAT‑CAGCTCACTTTCGCACTGATTATTATGGCGTCGTAGCAg                                                                                                                                                           <  2:145431/251‑1 (MQ=255)
                                                                                      gcAGATTGTGGATATGCACTTCCAGGGTGTTGGTCGAGGGTTCATTGTCCCAGTTATAGATGTCGTTGTAGAGAATTTGCCGAT‑GCACCGGACTGCCTGCTTTGTGCATTAATCGAGACAGCAGAGCATATTCTTTGGGCGTCAGAATCAAATCTTAACAGTCCCTCCATACCTAACGGCGTCCCATGTTCAGCGTAAGATTGACA‑CAATTCAGATCACTTTCGTCCTGTTTATTATAGCGTCGTAGccgg                                                                                                                                                          >  2:294916/1‑248 (MQ=255)
                                                                                        aGATTGTGGATATGCACTTCCAGGGTGTCGGTAGAGGGTTCATTGTCCCAGTTATAGATGTCGTTGTAGAGAATTTCCCGAT‑GCACCGGACTGCCTGCTTTGAGCATTAACCGTGACAGCAGAGCATATTCTTTGGGCGTCAGAATCAACTCTTCACCGCCCATCCATACCTGACGGCGACCCATGTTCAGCGTCAGATTGCCAACAAT‑CAGCTCACTTTCGCCCTGATTATTATGGCGTCGTAGCAGGgc                                                                                                                                                        <  1:208630/251‑1 (MQ=255)
                                                                                             gtgGATATGCACTTCCAGGGTGTTGGTCGAGGGTTCATTGTCCCAGTTATAGATGTCGTTGTAGAGAATTTCCCGAT‑GCACCGGACTGCCTGCTTTGAGCATTAACCGTGACAGCAGAGCATATTCTTTGGGCGTCAGAATCAACTCTTCACCGTCCATCCATACCTGACGGCGACCCATGTTCAGCGTCAGATTGCCAACAAT‑CAGCTCACTTTCGCCCTGATTATTATGGCGTCGTAGCAGGGCGCgg                                                                                                                                                    >  1:409463/1‑250 (MQ=255)
                                                                                                 atatGCACTTCCAGGGTGTTGGTCGAGGGTTCATTGTCCCAGTTGTAGATGTCGTTGTAGAGAATTTCGCGAT‑GCACCGGACTGCCTGCTTTGAGCATTAAACGTGACAGCAGAGCATATTCTTTGGGCGTCAGAATCAACTCTTCACCGTCCATCCATACCTGACGGCGACCCATGTTCAGCGTCAGATTGCCAACAAT‑CAGCTCACTTTCGCCCTGATTATTATGGCGTCGTAGCAGGGCGCGGATagg                                                                                                                                               >  2:234128/1‑249 (MQ=255)
                                                                                                    tGCACTTCCAGGGTGTTGGTCGAGGGTTCATTGTCCCAGTTATAGATGTCGTTGTAGAGAATTTCCCGAT‑GCACCGGACTGCCTGCTTTGAGCATTAACCGTGACAGCAGAGCATATTCTTTGGGCGTCAGAATCAACTCTTCACCGTCCATCCATACCTGACGGCGACCCATGTTCAGCGTCAGATTGCCAACAAT‑CAGCTCACTTTCGCCCTGATTATTATGGCGTCGTAGCAGGGCTAGGATACGGGc                                                                                                                                            >  1:162880/1‑251 (MQ=255)
                                                                                                    tGCACTTCCAGGGGGTTGGTCGAGGGTTCATTGTCCCAGTTATAGATGTCGTTGTAGAGAATTTCCCGAT‑GCACCGGACTGCCTGCTTTGAGCATTAACCGTGACAGCAGAGCATATTCTTTGGGCGTCAGAATCAACTCTTCACCGTCCATCCATACCTGACGGCGACCCATGTTCAGCGTCAGATTGCCAACAAT‑CAGCTCACTTTCGCCCTGATTATTATGGCGTCGTAGCAGGGCGCGGATACGGGc                                                                                                                                            >  2:354271/1‑251 (MQ=255)
                                                                                                             aGGGTGTTGGTCGAGGGTTCATTGTCCCAGTTATAGATGTCGTTGTAGAGAATTTCCCGAT‑GCACCGGCCTGCCTGCTTTGAGCATTAACCGTGACAGCAGAGCATATTCTTTGGGCGTCAGAATCAACTCTTCACCGTACATCCATACCTGACGGCGCCCCATGTTCAGCGTCAGATTGCCAACAAT‑CAGCTCACttt                                                                                                                                                                                       >  2:218326/1‑199 (MQ=255)
                                                                                                             aGGGTGTTGGTCGAGGGTTCATTGTCCCAGTTATAGATGTCGTTGTAGAGAATTTCCCGAT‑GCACCGGACTGCCTGCTTTGAGCATTAACCGTGACAGCAGAGCATATTCTTTGGGCGTCAGAATCAACTCTTCACCGTCCATCCATACCTGACGGCGACCCATGTTCAGCGTCAGATTGCCAACAAT‑CAGCTCACttt                                                                                                                                                                                       <  1:218326/199‑1 (MQ=255)
                                                                                                                                cATTGTCCCAGTTATAGATGTCGTTGTAGAGAATTTCCCGAT‑GCACCGGACTGCCTGCTTTGAGCATTAACCGTGACAGCAGAGCATATTCTTTGGGCGTCAGAATCAACTCTTCACCGTCCATCCATACCTGACGGCGACCCATGTTCAGCGTCAGATTGCCAACAAT‑CAGCTCACTTTCGCCCTGATTATTATGGCGTCGTAGCAGGGCGCGGATACGGGCATGTa                                                                                                                                       <  1:399945/228‑1 (MQ=255)
                                                                                                                                cATTGTCCCAGTTATAGATGTCGTTGTAGAGAATTTCCCGAT‑GCACCGGACTGCCTGCTTTGAGCATTAACCGTGACAGCAGAGCATATTCTTTGGGCGTCAGAATCAACTCTTCACCGTCCATCCATACCTGACGGCGACCCATGTTCAGCGTCAGATTGCCAACAAT‑CAGCTCACTTTCGCCCTGATTATTATGGCGTCGTAGCAGGGCGCGGATACGGGCATGTa                                                                                                                                       >  2:399945/1‑228 (MQ=255)
                                                                                                                                    gTCCCAGTTATAGATGTCGTTGTAGAGAATTTCCCGAT‑GCACCGGACTGCCTGCTTTGAGCATTAACCGTGACAGCCGAGCATATTCTTTGGGCGTCAGAATCAACTCTTCACCGTCCATCCATACCTGACGGCGACCCATGTTCAGCGTCAGATTGCCAACAAT‑CAGCTCACTTTCGCCCTGATTATTATGGCGTCGTAGCAGGGCGCGGATACGGGCATGTCACTCTTCCAGCGCAAAAGGCTTCAcc                                                                                                             >  2:139285/1‑250 (MQ=255)
                                                                                                                                    gTCCCAGTTATAGATGTCGTTGTAGAGAATTTCCCGAT‑GCACCGGACTGCCTGATTTGAGCATTAACCGTGACAGCAGAGCATATTCTTTGGGCGTCAGAATCAACTCTTCACCGTCCATCCATACCTGACGGCGACCCATGTTCAGCGTCAGATTGCCAACAAT‑CAGCTCACTTTCGCCCTGATTATTATGGCGTCGTAGCAGGGCGCGGATACGGGCATGTAACTCTTCCAGCGCAAAAGGCTTCAcc                                                                                                             >  1:60109/1‑250 (MQ=255)
                                                                                                                                                        ggTAGAGAATTTCCCGCT‑GCACAGGACTGCCTGCTTTGAGCATTAACCGTTACAGCAGAGCATATTCTTTGGGCGTCAGAATCAACTCTTCACCGTCCTTCCATACCTGACGGCGACCCATGTTCAGCGTCAGATTGCCAACAAT‑CAGCCCACTTTCGCCCTGATTATTCTGGCGTCGTAGCAGGGCGCTGATACGGGCATGTAACTCTTCCAGCGCACAAGGCTTCACCAGCTCGTAGGCGGCACCGaa                                                                                         >  2:69622/2‑249 (MQ=255)
                                                                                                                                                              gaATTTCCCGATTGCAC‑GGCCTGCCTGCTTTGAGCATTACCCGTGTCATCAGAGCATATTCTTTGGGCGTCAGAATCAACTCTTCTCCGTCCATCCATACCTGACGGCGACCCATGTTCAGCGTCAGATTGCCAACAAT‑CAGCTCACTTGCGCCCTGATTATTATGGCGTCGTAGCAGGGCGCGGATACGGGCATGTAACTCTTCCAGCGCAAAAGGCTTCACCAGATAGTCGTCGGCCCCGACATCCAGc                                                                                  <  2:68982/251‑1 (MQ=255)
                                                                                                                                                              aaaTTTCCCGAT‑GCACCGGACTGCCTGCTTTGAGCATTAACCGTGACAGCAGAGCATATTCTTTGGGCGTCAGAATCAACTCTTCACCGTCCATCCATACCTGACGGCGACCCATGTTCAGCGTCAGATTGCCAACAAT‑CAGCTCACTTTCGCCCTGATTATTATGGCGTCGTAGCAGGGCGCGGATACGGGCATGTAACTCTTCCAGCGCAAAAGGCTTCACCAGATAGTCGTCGGCACCGACATCCAGc                                                                                  <  1:46371/250‑1 (MQ=255)
                                                                                                                                                                    cccGAT‑GCACCGGACTGCCTGCTTTGAGCATTAACCGTGACAGCAGAGCATATTCTTTGGGCGTCAGAATCAACTCTTCACCGTCCATCCATACCTGACGGCGACCCATGTTCAGCGTCAGATTGCCAACAAT‑CAGCTCACTTTCGCCCTGATTATTATGGCGTCGTAGCAGGGCGCGGATACGGGCATGTAACTCTTCCAGCGCAAAAGGCTTCACCAGATAGTCGTCGGCACCGACATCCAGCCCGGCg                                                                            >  2:341257/1‑251 (MQ=255)
                                                                                                                                                                              ccGGACTGCCTGCTTTGAGCATTAACCGTGACAGCAGAGCATATTCTTTGGGCGTCAGAATCAACTCTTCACCGTCCATCCATACCTGACGGCGACCCATGTTCAGCGTCAGATTGCCAACAAT‑CAGCTCACTTTCGCCCTGATTATTATGGCGTCGTAGCAGGGCGCGGATACGGGCATGTAACTCTTCCAGCGCAAAAGGCTTCACCAGATAGTCGTCGGCACCGACATCCAGCCCGTCGATTTTGTCg                                                                   <  1:98362/251‑1 (MQ=255)
                                                                                                                                                                                        tgcTTTGAGCATTAACCGTGACAGCAGAGCATATTATTTGGGCGTCAGAATCAACTCTTCACCGTCCATCCATACCTGACGGCGACCCATGTTCAGCGTCAGATTGCCAACAAT‑CAGCTCACTTTCGCCCTGATTATTATGGCGTCGTAGCAGGGCGCGGATACGGGCATGTAACTCTTCCAGCGCAAAAGGCTTCACCAGATAGTCGTCGGCACCGCCATCCAGCCCGGCGATTTTGTCGGTCAGCGTAt                                                         >  1:198870/1‑251 (MQ=255)
                                                                                                                                                                                         gcTTTGAGCATTAACCGTGACAGCAGAGCATATTCTTTGGGCGTCAGAATCAACTCTTCACCGTCCATCCATACCTGACGGCGACCCATGTTCAGCGTCAGATTGCCAACAAT‑CAGCTCACTTTCGCCCTGATTATTATGGCGTCGTAGCAGGGCGCGGATACGGGCATGTAACTCTTCCAGCGCAAAAGGCTTCACCAGATAGTCGTCGGCACCGACATCCAGCCCGGCGATTTTGTCGGTCAGCGTATc                                                        <  1:105594/251‑1 (MQ=255)
                                                                                                                                                                                            ttGAGCATTAACCGTGACAGCAGAGCATATTCTTTGGGCGTCAGAATCAACTCTTCACCGTCCATCCATACCTGACGGCGACCCATGTTCAGCGTCAGATTGCCAACAAT‑CAGCTCACTTTCGCCCTGATTATTATGGCGTCGTAGCAGGGCGCGGATACGGGCATGTAACTCTTCCAGCGCAAAAGGCTTCACCAGATAGTCGTCGGCACCGACATCCAGCCCGGCGATTTTGTCGGTCAGCGTATcgcg                                                     >  2:112012/1‑251 (MQ=255)
                                                                                                                                                                                            ttGAGCATTAACCGTGACAGCAGAGCATATTCTTTGGGCGTCAGAATCAACTCTTCACCGTCCATCCATACCTGACGGCGACCCATGTTCAGCGTCAGATTGCCAACAAT‑CAGCTCACTTTCGCCCTGATTATTATGGCGTCGTAGCAGGGCGCGGATACGGGCATGTAACTCTTCCAGCGCAAAAGGCTTCACCAGATAGTCGTCGGCACCGACATCCAGCCCGGCGATTTTGTCGGTCAGCGTATcg                                                       <  1:110020/249‑1 (MQ=255)
                                                                                                                                                                                                 cATTAACCGTGACAGCAGAGCATATTCTTTGGGCGTCAGAATCAACTCTTCACCGTCCATCCATACCTGACGGCGACCCATGTTCAGCGTCAGATTGCCAACAAT‑CAGCTCACTTTCGCCCTGATTATTATGGCGTCGTAGCAGGGCGCGGATACGGGCATGTAACTCTTCCAGCGCAAAAGGCTTCACCAGATAGTCGTCGGCACCGACATCCCACCCGGCGATTTTGTCGGTCAGCGTATCGCGAGCgg                                                >  2:371884/1‑251 (MQ=255)
                                                                                                                                                                                                    tAACCGTGACAGCAGAGCATATTCTTTGGGCGGCAGAATCAACTCCTCACCGTCCATCCATACCTGACGGCGACCCATGTTCAGCGTCAGATTGCCAACAAT‑CAGCTCACTTTCGCCCTGATTATTATGGCGTCGTAGCAGGGCGCGGATACGGGCATGTAACTCTTCCAGCGCAAAAGTCTTCACCAGATAGTCGTCGGCAc                                                                                             <  1:168410/203‑1 (MQ=255)
                                                                                                                                                                                                    tAACCATGACAGCAGAGCCTATTCTTTGGGCGTCAGAATCAACTCTTCACCGTCCATCCATACCTGACGGCGACCCATGTTCAGCGTCAGATTTCCAACAAT‑CCGCTCACTTTCGCCCTGATTATTATGGCTTCGTAACAGGGCGCGGATACGGGCATGTAACTCTTCCAGCGCAACAGGCTTCACCAGATAGTCGTCGGCAc                                                                                             >  2:168410/1‑203 (MQ=255)
                                                                                                                                                                                                                 agagCATATTCTTTGGGCGTCAGAATCAACTCTTCACCGTCCATCCATACCTGACGGCGACCCATGTTCAGCGTCAGATTGCCAACAAT‑CAGCTCACTTTCGCCCTGATTATTATGGCGTCGTAGCAGGGCGCGGCTACGGGCATGTAACTCTTCCAGCGCAACAGGCTTCACCAGATAGACGTCGGCACCGACATCCAGCCCGGCGATTTTGTCGGTCAGCGTATCGCGAGCGGTTAGGTTCAGTACCgg                                >  2:417590/1‑251 (MQ=255)
                                                                                                                                                                                                                   agCATATTCACTGGGCGTCAGAATCAACTCTTCACCGTCCATCCATACCTGACGGCGACCCATGATCAACGTCAGATTGCCACCCAT‑CAGCTCACTTTCGCCCTGATTATTATGGCGTCGCAGCAGGGCGCGGATACGCGCCTGTCACTCTTCCAGCGCAACAGGCTTCACCAGATAGCCGTCGGCACCGACATCCAGCCCGGCGATTTTGTCGGTCAGCGTCTCGCGAGCGGTGAGGATCAGTACCGGcg                              >  2:389715/1‑250 (MQ=255)
                                                                                                                                                                                                                     cATATTCTTTGGGCGTCAGAATCAACTCTTCCCCGTCCATCCATACCTGACGGCGACCCATGTTCAGCGTCAGATTGCCAACAAT‑CAGCTCACTTTCGCCCTGATTATTATGGCGTCGTAGCAGGGCGCGGATACGGGCATGTAACTCTTCCAGCGCAAAAGGCTTCACCAGATAGTCGTCGGCACCGACATCCAGCCCGGCGATTTTGTCGGTCAGCGTATCGCGAGCGGTGAGGATCAGTACCGGCAgg                            <  1:167137/251‑1 (MQ=255)
                                                                                                                                                                                                                         ttCTTTGGGCATCAGAAGCAACTCTTCCACTTCCAGCCATGCCTTGCGGCGCCCTATGTTCAGCGTCAGATTGCCAACATT‑CAGCTCCCTTTCGCCCTTATTATTATGGCGTCGTAGCAGGGCGCGTATACGGGCATGTAACTCTTCCAGAGCAAAAGGCTTCACCAGATAGTCGTCGGCACCGACATCCAGCCCGGCGATTTTGTCGGTCAGCGTATCGCGAGCGGTGAGGATCAGTACCGGCAGGGta                         <  2:331349/250‑1 (MQ=255)
                                                                                                                                                                                                                              tGGGCTTCAGAATCAACTCTTCACCGTCCATCTATAACTGACGGCGCCCCATGTTCAGCGTCCGATTGCCAACAAT‑CAGCTCCCTTGCGCCCTGATTATTCTGGCGTCGTAGCAGGGCGCGGATACGGGCATGTAACTCTTCCAGCGCAAAAGGCTTCACCATGTAGTCGTCGGCACCGACATCCGGCACTGCGATTTTGTCGGTCAGCGCATAGAGAGCGGTGAGGCCCGGGACCGGCAGGGGAAAtttt                   >  2:120912/1‑251 (MQ=255)
                                                                                                                                                                                                                                    tCAGAATCAACTCTTCACCGTCCATCCATACCTGACGGCGACCCATGCTCAGCGTCAGATTGCCAACAAT‑CAGCTCACTTTCGCCCTGATTATTATGGCGTCGTAGCAGGGCGCGGATACGGGCATGTAACTCTTCCAGCGCAAAAGGCTTCACCAGATAGTCGTCGGCACCGACATCCAGCCCGGCGATTTTGTCGGTCAGCGTATCGCGAGAGGCGAGGATCAGTACCGGCAGGGTATATTTTTTATGc             >  2:8383/1‑251 (MQ=255)
                                                                                                                                                                                                                                            aaCTCTTCACCGTCCATCCATACCTGACGGCGACCCATGTTCAGCGTCAGATTGCCAACAAT‑CAGCTCACTTTCGCCCTGatta                                                                                                                                                                            >  1:274652/1‑84 (MQ=255)
                                                                                                                                                                                                                                            aaCTCTTCACCGTCCATCCATACCTGACGGCGACCCATGTTCAGCGTCAGATTGCCAACAAG‑CAGCTCACTTTCGCCCTGatta                                                                                                                                                                            <  2:274652/84‑1 (MQ=255)
                                                                                                                                                                                                                                               tctTCACCGTCCATCCATACCTGACGGCGACCCATGTTCAGCGTCAGATTGCCAACAAT‑CAGCTCACTTTCGCCCTGATTATTATGGCGTCGTAGCAGGGCGCGGATACGGGCATGTAACTCTTCCAGCGCAAAAGGCTTCACCAGATAGTCGTCGGCACCGACATCCAGCCCGGCGATTTTGTCGGTCAGCGTATCGCGAGCGGTGAGGATCAGTACCGGCAGGGTATATTTTTTCTGCCGGATACGGGc  >  1:248821/1‑251 (MQ=255)
                                                                                                                                                                                                                                                        |                                                                                                                                                                                                                                                  
CAGATTCAATTAGTTTTCCTCATTCGCGACCAGCATATAGCCAAAGCCGCGCACGGTGCGGATACGGGCTTTGCCCACTTTGTCGCGCAGATTGTGGATATGCACTTCCAGGGTGTTGGTCGAGGGTTCATTGTCCCAGTTATAGATGTCGTTGTAGAGAATTTCCCGAT‑GCACCGGACTGCCTGCTTTGAGCATTAACCGTGACAGCAGAGCATATTCTTTGGGCGTCAGAATCAACTCTTCACCGCCCATCCATACCTGACGGCGACCCATGTTCAGCGTCAGATTGCCAACAAT‑CAGCTCACTTTCGCCCTGATTATTATGGCGTCGTAGCAGGGCGCGGATACGGGCATGTAACTCTTCCAGCGCAAAAGGCTTCACCAGATAGTCGTCGGCACCGACATCCAGCCCGGCGATTTTGTCGGTCAGCGTATCGCGAGCGGTGAGGATCAGTACCGGCAGGGTATATTTTTTCTGCCGGATACGGGC  >  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: ‑