Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 122,761 C→T H224Y (CAT→TAT)  pdhR → DNA‑binding transcriptional dual regulator PdhR

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913122,7610CT97.7% 136.7 / NA 43H224Y (CAT→TAT) pdhRDNA‑binding transcriptional dual regulator PdhR
Reads supporting (aligned to +/- strand):  ref base C (0/0);  major base T (14/28);  minor base A (1/0);  total (15/28)
Fisher's exact test for biased strand distribution p-value = 3.49e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.82e-01

GATCTGGACGCGGAATCAAACGCCGTACTCCAGTATCAGATTGCCGTCACCGAAGCGGCCCACAATGTGGTTCTGCTTCATCTGCTAAGGTGTATGGAGCCGATGTTGGCCCAGAATGTCCGCCAGAACTTCGAATTGCTCTATTCGCGTCGCGAGATGCTGCCGCTGGTGAGTAGTCACCGCACCCGCATATTTGAAGCGATTATGGCCGGTAAGCCGGAAGAAGCGCGCGAAGCATCGCATCGCCATCTGGCCTTTATCGAAGAAATTTTGCTCGACAGAAGTCGTGAAGAGAGCCGCCGTGAGCGTTCTCTGCGTCGTCTGGAGCAACGAAAGAATTAGTGATTTTTCTGGTAAAAATTATCCAGAAGATGTTGTAAATCAAGCGCATATAAAAGCGCGGCAACTAAACGTAGAACCTGTCTTATTGAGCTTTCCGGCGAGAGTTCAATGGGACAGGTTCCAGAAAACTCAACGTTATTAGATAGATAA  >  NC_000913/122515‑123006
                                                                                                                                                                                                                                                      |                                                                                                                                                                                                                                                     
ggtCTGGACGCGGAATCAAACGCCGTACTCCAGTATCAGATTGCCGTCACCGAAGCGGCCCACAATGTGGTTCTGCTTCATCTGCTAAGGTGTATGGAGCCGATGTTGGCCCAGAATGTCCGCCAGAACTTCGAATTGCTCTATTCGCGTCGCGAGATGCTGCCGCTGGTGAGTAGTCACCGCACCCGCATATTTGAAGCGATTATGGCCGGTAAGCCGGAAGAAGCGCGCGAAGCATCGCATCGCTATCt                                                                                                                                                                                                                                                   <  2:252700/249‑1 (MQ=255)
     ggACGCGGAATCAAACGCCGTACTCCAGTATCAGATTGCCGTCACCGAAGCGGCCCACAATGTGGTTCTGCTTCATCTGCTAAGGTGTATGGAGCCGATGTTGGCCCAGAATGTCCGCCAGAACTTCGAATTGCTCTATTCGCGTCGCGAGATGCTGCCGCTGGTGAGTAGTCACCGCACCCGCATATTTGAAGCGATTATGGCCGGTAAGCCGGAAGAAGCGCGCGAAGCATCGCATCGCTATc                                                                                                                                                                                                                                                    <  2:87518/245‑1 (MQ=255)
     ggACGCGGAATCAAACGCCGTACTCCAGTATCAGATTGCCGTCACCGAAGCGGCCCACAATGTGGTTCTGCTTCATCTGCTAAGGTGTATGGAGCCGATGTTGGCCCAGAATGTCCGCCAGAACTTCGAATTGCTCTATTCGCGTCGCGAGATGCTGCCGCTGGTGAGTAGTCACCGCACCCGCATATTTGAAGCGATTATGGCCGGTAAGCCGGAAGAAGCGCGCGAAGCATCGCATCGCTATCt                                                                                                                                                                                                                                                   <  1:178737/246‑1 (MQ=255)
     ggACGCGGAATCAAACGCCGTACTCCAGTATCAGATTGCCGTCACCGAAGCGGCCCACAATGTGGTTCTGCTTCATCTGCTAAGGTGTATGGAGCCGATGTTGGCCCAGAATGTCCGCCAGAACTTCGAATTGCTCTATTCGCGTCGCGAGATGCTGCCGCTGGTGAGTAGTCACCGCACCCGCATATTTGAAGCGATTATGGCCGGTAAGCCGGAAGAAGCGCGCGAAGCATCGCATCGCTATCt                                                                                                                                                                                                                                                   <  2:407079/246‑1 (MQ=255)
                 aaaCGCCGTACTCCAGTATCAGATTGCCGTCACCGAAGCGGCCCACAATGTGGTTCTGCTTCATCTGCTAAGGTGTATGGAGCCGATGTTGGCCCAGAATGTCCGCCAGAACTTCGAATTGCTCTATTCGCGTCGCGAGATGCTGCCGCTGGTGAGTAGTCACCGCACCCGCATATTTGAAGCGATTATGGCCGGTAAGCCGGAAGAAGCGCGCGAAGCATCGCATCGCTATCTGGCCTTTATCGAAGaaa                                                                                                                                                                                                                                  <  2:244355/251‑1 (MQ=255)
                            tccagtatccccttGCCGTCAACGAACCGGCCCACAATGTGGTTCTGCTTAATCTGCTCAGGTGTATGGAGCCGATGTTGGCCCAGAATGTCCGCCAGCACTTCGAATTGCTCTATTCGCGTCGCGAGATGCTGCCGCTGGTGAGTAGTCACCGCACCCGCATATTTGAAGCGATTATGGCCGGTAAGCCGGAAGAAGCGCGCGAAGCATCGCATCGCTATCTGGCCTTTATCGAAGAAATTTTGCTCGAc                                                                                                                                                                                                                       >  2:68373/13‑251 (MQ=255)
                               agtttCAGATTGCCTTCACCGAAGCGGCCCACAATGGGGTTCTGCTTCATCTTCTAAGGTGTATGGAGCCTATGTTGGCCCAGAATGTCCGCCAGAACTGCGAATTGCTCTATTCGCGTCGCGAGATGCTGCCGCTGGTGAGTAGTCACCGCACCCGCATATTTGAAGCGATTATGGCCGGTAAGCCGGAAGAAGCGCGCGAAGCATCGCATCGCTATCTGGCCTTTATCGAAGAAATTTTGCTCGACAGa                                                                                                                                                                                                                    <  2:182341/247‑1 (MQ=255)
                                      gATTGCCGTCACCGAAGCGGCCCACAATGTGGTTCTGCTTCATCTGCTAAGGTGTATGGAGCCGATGTTGGCCCAGAATGTCCGCCAGAACTTCGAATTGCTCTATTCGCGTCGCGAGATGCTGCCGCTGGTGAGTAGTCACCGCACCCGCATATTTGAAGCGATTATGGCCGGTAAGCCGGAAGAAGCGCGCGAAGCATCGCATCGCTATCTGGCCTTTATCGAAGAAATTTTGCTCGACAGAAGTCGTg                                                                                                                                                                                                             >  1:199220/1‑251 (MQ=255)
                                              tCACCGAAGCGGCCCACAATGTGGTTCTGCTTCATCTGCTAAGGTGTATGGAGCCGATGTTGGCCCAGAATGTCCGCCAGAACTTCGAATTGCTCTATTCGCGTCGCGAGATGCTGCCGCTGGTGAGTAGTCACCGCACCCGCATATTTGAAGCGATTATGGCCGGTAAGCCGGAAGAAGCGCGCGAAGCATCGCATCGCTATCTGGCCTTTATCGAAGAAATTTTGCTCGACAGAAGTCGTGAAGAGAgc                                                                                                                                                                                                     >  1:29648/1‑251 (MQ=255)
                                               cACCGAAGCGGCCCACAATGTGGTTCTGCTTCATCTGCTAAGGTGTATGGAGCCGATGTTGGCCCAGAATGTCCGCCAGAACTTCGAATTGCTCTATTCGCGTCGCGAGATGCTGCCGCTGGTGAGTAGTCACCGCACCCGCATATTTGAAGCGATTATGGCCGGTAAGCCGGAAGAAGCGCGCGAAGCATCGCATCGCTATCTGGCCTTTATCGAAGAAATTTTGCTCGACAGAAGTCGTGAagagatcc                                                                                                                                                                                                    <  2:336639/251‑4 (MQ=255)
                                                          cccACAATGTGGTTCTGCTTCATCTGCTAAGGTGTATGGAGCCGATGTTGGCCCAGAATGTCCGCCAGAACTTCGAATTGCTCTATTCGCGTCGCGAGATGCTGCCGCTGGTGAGTAGTCACCGCACCCGCATATTTGAAGCGATTATGGCCGGTAAGCCGGAAGAAGCGCGCGAAGCATCGCATCGCTATCTGGCCTTTATCGAAGAAATTTTGCTCGACAGAAGTCGTGAAGAGAGCCGCCGTGAGCGt                                                                                                                                                                                         <  2:211378/251‑1 (MQ=255)
                                                            caaaaTGTGGTTCTGCTTCATCTGCTAAGTTGTATGGGGCCGATGTTGGCCCAGAAGGTCCGCCAGAACTTCGAATTGCTCTATTCGCGTCGCGAGATGCTGCCGCGGGGGAGTAGTCACCGCACCCGCATATTTGAAGCGATTATGGCCGGTAGGCCGGAAGAAGCGCGCGAAGCATCGCATCGCTATCTGGCCTTTATCGAAGAAATTTTGCTCGACAGAAGTCGTGAAGAGAGCCGCCGTGAGCGTtc                                                                                                                                                                                       <  2:464681/248‑1 (MQ=255)
                                                                     gTTCTGCTTCATCTGCTAAGGTGTATGGAGCCGATGTTGGCCCAGAATGTCCGCCAGAACTTCGAATTGCTCTATTCGCGTCGCGAGATGCTGCCGCTGGTGAGTAGTCACCGCACCCGCATATTTGAAGCGATTATGGCCGGTAAGCCGGAAGAAGCGCGCGAAGCATCGCATCGCTATCTGGCCTTTATCGAAGAAATTTTGCTCGACAGAAGTCGTGAAGAGAGCCGCCGTGAGCGTTCTCTGcgtcg                                                                                                                                                                              >  2:131218/1‑251 (MQ=255)
                                                                          gctgcATCTGCTAAGGTGTATGGCCCCTAGGTGGGCCCTGAATGTCCGCCAGACCTTCGAATTGCCCTATTCGCGTCGCGAGATGCTGTCGCTGGTGAGTAGTCACCGCACCCGCATATTTGAAGCGATTATGCCCGGTAAGCCGGAAGAAGCGCGCGAAGCATCGCATCGTTATCTGGCCTTTATCGAAAAAATTTTGCTCGACAGAAGTCTTGAAGAGAGCCGCCGTGAGCGTTCTCTGCGTCGTCTg                                                                                                                                                                          <  2:281661/246‑1 (MQ=255)
                                                                          gCTTCATCTGCTAAGTTGTATGGGGCCGATGTTGGCCCAGAATGTCCGCCAGAACTTCGAATTGCTCTATGCGCGGCGCGAGATGCTGCCGCTGGTGAGTAGTCCCCGCACCCGCATATTTGAAGCGATTATGGCCGGTAAGCCGGAAGAAGCCCGCGAAGCATCGCATCGCTAGCTGGCCTTTATCGAAGAAATTTTGCTCGACAGAAGTCGTGAAGAGAGCCGCCGTGAGCGTTCTCTGCGTCGTCTg                                                                                                                                                                          <  2:125370/250‑1 (MQ=255)
                                                                          gCTTCATCTGCTAAGGTGTATGGAGCCGATGTTGGCCCAGAATGTCCGCCAGAACTTCGAATTGCTCTATTCGCGTCGCGAGAGGCTGCCGCTGGTGAGTAGTCACCGCACCCGCATATTTGAAGCGATTATGGCCGGTAAGCCGGAAGAAGCGCGCGAAGCATCGCATCGCTATCTGGCCTTTATCGAAGAAATTTTGCTCGACAGAAGTCGTGAAGAGAGCCGCCGTGAGCGTTCTCTGCgtcgtc                                                                                                                                                                            <  2:464497/248‑1 (MQ=255)
                                                                          gCTTCATCTGCTAAGGTGGAGGGAGCCGGTTTTGGCCCAGAATGTCCGCCAGAACTGCGAATTGCTCTATTCGCGTAGCGAGATGCTGCCGCTGGTGAGTAGTCACCGCACCCGCATATTTGAAGCGATTATGGCCGGTAAGCCGGAAGAAGCGCGCGAAGCATCGCATCGCTATCTGGCCTTTATCGAAGAAATTTTGCTCGACAGAAGTCGTGAAGAGAGCCGCCGTGAGCGTTCTCTGCGTCGTCTg                                                                                                                                                                          <  2:281688/250‑1 (MQ=255)
                                                                              cATCTGCTAAGGTGTATGGAGCCGATGTTGGCCCAGAATGTCCGCCAGAACTTCGAATTGCTCTATTCGCGTCGCGAGATGCTGCCGCTGGTGAGTAGTCACCGCACCCGCATATTTGAAGCGATTATGGCCGGTAAGCCGGAAGAAGCGCGCGAAGCATCGCATCGCTATCTGGCCTTTATCGAAGAAATTTTGCTCGACAGAAGTCGTGAAGAGAGCCGCCGTGAGCGTTCTCTGCGTCGTCTGGAGc                                                                                                                                                                      >  2:50185/1‑250 (MQ=255)
                                                                                 cTGCTAAGGAGTATGGAGCCGATGTTGGCCCAGACTGTTTGCCAGAAATTAGAAATGCTCTATTCGCGTCGCGAGCTGCTGCCGCTGGTGAGTAGTCACAGCACGCGCATATTTGAAGCGATTATGGCCGGCAAGCCGGAAGAAGCGCGCGAAGCATCGCATAGCTATCTGGCCTTTATCGAAGAAATTTTGCTCGACAGAAGTCGTGCAGAGAGACGCCGTGAGCGTTCTCTGCGTCGTCTGGAGCAACg                                                                                                                                                                  >  1:306113/1‑251 (MQ=255)
                                                                                                                  aaTGTCCGCCAGAACTTCGAATTGCTCTATTCGCGTCGCGAGATGCTGCCGCTGGAGAGTAGTCACCGCACCCGCATATTTGAAGCGATTATGGCCGGTAAGCCGGAAGAAGCGCGCGAAGCATCGCATCGCTATCTGGCCTTTATCGAAGAAATTTTGCTCGACAGAAGTCGTGAAGAGAGCCGCCGTGAGCGTTCTCTGCGTCGTCTGGAGCAACGAAAGAATTAGTGATTTTTCTGGTAAAAATTATc                                                                                                                                 >  2:465872/1‑251 (MQ=255)
                                                                                                                          ccAGAACTTCGAATTGCTCTATTCGCGTCGCGAGATGCTGCCGCTGGTGAGTAGTCACCGCACCCGCATATTTGAAGCGATTATGGCCGGTAAGCCGGAAGAAGCGCGCGAAGCATCGCATCGCTATCTGGCCTTTATCGAAGAAATTTTGCTCGACAGAAGTCGTGAAGAGAGCCGCCGTGAGCGTTCTCTGCGTCGTCTGGAGCAACGAAAGAATTAGTGATTTTTCTGGTAAAAATTATCCAGAAGAt                                                                                                                         <  1:381617/251‑1 (MQ=255)
                                                                                                                                    gAATTGCTCTATTCGCGTCGCGAGATGCTGCCGCTGGTGAGTAGTCACCGCACCCGCATATTTGAAGCGATTATGGCCGGTAAGCCGGAAGAAGCGCGCGAAGCATCGCATCGCTATCTGGCCTTTATCGAAGAAATTTTGCTCGACAGAAGTCGTGAAGAGAgccgc                                                                                                                                                                                                  <  1:174057/168‑1 (MQ=255)
                                                                                                                                    gAATTGCTCTATTCGCGTCGCGAGATGCTGCCGCTGGTGAGTAGTCACCGCACCCGCATATTTGAAGCGATTATGGCCGGTAAGCCGGAAGAAGCGCGCGAAGCATCGCATCGCTATCTGGCCTTTATCGAAGAAATTTTGCTCGACAGAAGTCGTGAAGAGAgccgc                                                                                                                                                                                                  >  2:174057/1‑168 (MQ=255)
                                                                                                                                         gCTCTATTCGCGTCGCGAGATGCTGCCGCTGGTGAGTAGTCACCGCACCCGCATATTTGAAGCGATTATGGCCGGTAAGCCGGAAGAAGCGCGCGAAGCATCGCATCGCTATCTGGCCTTTATCGAAGAAATTTTGCTCGACAGAAGTCGTGAAGAGAGCCGCCGTGAGCGTTCTCTGCGTCGTCTGGAGCAACGAAAGGATTAGTGATTTTTCTGGTACAAATTATCCAGAAGATGTTGTAACTCAAgcg                                                                                                          >  2:370938/1‑251 (MQ=255)
                                                                                                                                           tctATTCGCGTCGCGAGATGCTGCCGCTGGTGAGTAGTCACCGCACCCGCATATTTGAAGAGATTATGGCCGGTAAGCCGGAAGAAGCGCGCGCAGCATCGCATCGCTATCTGGCCTTTATCGAAGAAATTTTGCTCGACCGACGTCGTGAAGAGAGCCGCCGTGAGCGTTCTCTGCGTCGTCTGGAGCAACGAAAGAATTAGTGATTTTTCTGGTAAAAATTATCCAGAAGATGTTGTAAATCAAgcgc                                                                                                         >  1:267561/1‑250 (MQ=255)
                                                                                                                                                 cgcgtcgcgAGATGCTGCCGCTGGTGAGTAGTCACCGCACCCGCATATTTGAAGCGATTATGGCCGGTAAGCCGGAAGAAGCGCGCGAAGCATCGCATCGCTATCTGGCCTTTATCGAAGAAATTTTGCTCGACAGAAGTCGTGAAGAGAGCCGCCGTGAGCGTTCTCTGCGTCGTCTGGAGCAACGAAAGAATTAGTGATTTTTCTGGTAAAAATTATCCAGAAGATGTTGTAAATCAAGCGCATATaa                                                                                                   >  2:114164/1‑250 (MQ=255)
                                                                                                                                                  gagtcgcgAGATGCTGCCGCGGGTGAGTAGTCCCCGCACCCGCATATTTGAAGCGATTATGGCCGGTAAGCCGGAAGAAGCGCGCGAAGCATCGCATCGCTATCTGGCCTTTATCGAAGAAATTTTACTCGACAGAAGTCGTGAAGAGAGCCGCCGTGAGCGTTCTCTGCGTCGTCTGGAGCAACGAAAGAATTAGTGATTTTTCTGGTAAAAATTATCCAGAAGATGTTGTAAATCAAGCGCATATaaaa                                                                                                 <  2:328048/249‑1 (MQ=255)
                                                                                                                                                                 gcCGCTGGTGAGTAGGCACCGCACCCGCATATTTGAAGCGATTATGGCCGGTAAGCCGGAAGAAGCGCGCGAAGCATCGCATCGCTATCTGGCCTTTATCGAAGAAATTTTGCTCGACAGAAGTCGTGAAGAGAGCCGCCGTGAGCGTTCTCTGCGTCGTCTGGAGCAACGAAAGAATTAGTGATTTTTCTGGTAAAAATTATCCAGAAGATGTTGTAAATCAAGCGCATATAAAAGCGCGGCAACTaa                                                                                    <  2:142825/249‑1 (MQ=255)
                                                                                                                                                                      gggTGAGTAGTCACCGCACCCGCATATTTGAAGCGATTATGGCCGGTAAGCCGGAAGAAGCGCGCGAAGCATCGCATCGCTATCTGGCCTTTATCGAAGAAATTTTGCTCGACAGAAGTCGTGAAGAGAGCCGCCGTGAGCGTTCTCTGCGTCGTCTGGAGCAACGAAAGAATTAGTGATTTTTCTGGTAAAAATTATCCAGAAGATGTTGTAAATCAAGCGCATATAAAAGCGCGGCAACTAAACGTAGa                                                                             <  2:196783/250‑1 (MQ=255)
                                                                                                                                                                       ggTGAGTAGTCACCGCACCCGCATATTTGAAGCGATTATGGCCGGTAAGCCGGAAGAAGCGCGCGAAGCATCGCATCGCTATCTGGCCTTTATCGAAGAAATTTTGCTCGACAGAAGTCGTGAAGAGAGCCGCCGTGAGCGTTCTCTGCGTCGTCTGGAGCAACGAAAGAATTAGTGATTTTTCTGGTAAAAATTATCCAGAAGATGTTGTAAATCAAGCGCATATAAAAGCGCGGCAACTAAACGTAGaa                                                                            <  2:137509/251‑1 (MQ=255)
                                                                                                                                                                               gtCACCGCACCCGCATATTTGAAGCGATTATGGCCGGTAAGCCGGAAGAAGCGCGCGAAGCAGCGCATCGCTATCTGGCCTTTATCGAAGAAATTTTGCTCGACAGAAGTCGTGAAGAGAGCCGCCGTGAGCGTTCTCTGCGTCGTCTGGAGCAACGAAAGAATTAGTGATTTTTCTGGTAAAAATTATCCAGAAGATGTTGTAAATCAAGCGCATATAAAAGCGCGGCAACTAAACGTAGAACCTGTCtt                                                                    <  1:18831/251‑1 (MQ=255)
                                                                                                                                                                                             atatTTGAAGCGATTATGGACGGTAAGCCGGAAGCAGCGCGCGAAGCATCGCATCGCTATCTGGCCTTTATCGAAGAAATTTTGCTCGACAGAAGTCGTGAAGAGAGCCGCCGTGAGCGTTCTCTGCGTCGTCTGGAGCAACGAAAGAATTAGTGATTTTTCTGGTAAAAATTATCCAGAAGATGTTGTAAATCAAGCGCATATAAAAGCGCGGCAACTAAACGTAGAACCTGTCTTATTGAGCTTTCCgg                                                      >  2:302611/1‑251 (MQ=255)
                                                                                                                                                                                                         aTTATGGCCGGTAAGCCGGAAGAAGCGCGCGAAGCATCGCATCGCTATCTGGCCTTTATCGAAGAAATTTTGCTCGACAGAAGTCGTGAAGAGAGCCGCCGTGAGCGTTCTCTGCGTCGTCTGGAGCAACGAAAGAATTAGTGATTTTTCTGGTAAAAATTATCCAGAAGATGTTGTAAATCAAGCGCATATAAAAGCGCGGCAACTAAACGTAGAACCTGTCTTATTGAGCTTTCCGGCGAGAGTTCAAt                                          >  1:94814/1‑251 (MQ=255)
                                                                                                                                                                                                           tATGGCCGGTAAGCCGGAAGAAGCGCGCGAAGCATCGCATCGCTATCTGGCCTTTATCGAAGAAATTTTGCTCTACAGAAGTCGTGAAGAGAGCCGCCGTGAGCTTTCTCTGCGTCGTCTGGAGCAACGAAAGAATTAGTGATTTTTCTGGTAAAAATTATCCAGAAGATGTTGTAAAGCAAGCGCATATAAAAGCGCGGCAACTAAACGTAGAACCTTTCTTATTGAGCTTTCCGGCGAGAGTTCAATg                                         <  2:412644/250‑1 (MQ=255)
                                                                                                                                                                                                               gCCGGTAAGCCGGAAGAAGCGCGCGAAGCATCGCATCGCTATCTGGCCTTTATCGAAGAAATTTTGCTCGACAGAAGTCGTGAAGAGAGCCGCCGTGAGCGTTCTCTGCGTCGTCTGGAGCAACGAAAGAATTAGTGATTTTTCTGGTAAAAATTATCCAGAAGATGTTGTAAATCAAGCGCATATAAAAGCGCGGCAACTAAACGTAGAACCTGTCTTATTGAGCTTTCCGGCGAGAGTTCAATGGGACa                                    <  1:50185/251‑1 (MQ=255)
                                                                                                                                                                                                                    aaaGCCGGAAGAAGCGCGCGAAGCATCGCATCGCTATCTGGCATTTATCGAAGAAATTTTGCTCGACAGAATTCGTGAAGAGAGCCGCCGTGAGCGTTCTCTGCGTCGTCTGGAGCAACGAAAGAATTAGTGATTTTTCTGGTAAAAATTATCCAGAAGATGTTGTAAATCAAGCGCATATAAAAGCGCGGCAACTAAACGTAGAACCTGTCTTATTGAGCTTTCCGGCGAGAGTTCAATGGGACAGGTTc                               <  1:313248/250‑1 (MQ=255)
                                                                                                                                                                                                                     aaGCCGGAAGAAGCGCGCGAAGCATCGCATCGCAATCTGGCCTTTATCGCAGAAATTTTGCTCGACAGAAGTCGTGAAGAGACCCGCCGTGAGCGTTCTCTGCGTCGTATGGAGCAACGAAAGAATTAGAGATTTTTCTGGTAAAAATTATCCAGAAGATGTTGTAAATCAAGCGCATATAAAAGCGCGGCAACTAAACGTAGAACCTGTCTTATTGAGCTTTCCGGCGAGAGTTCAATGAGACAGGTTcc                              >  1:218865/1‑251 (MQ=255)
                                                                                                                                                                                                                         cGGAAGAAGCTCGCGAAGCATCGCATCGCTATCTGGCCTTTATCGAATAAATTTTGCTCGACAGAAGTCGTGAAGAGAGCCGCCGTGAGCGTTCTCTGCGTCGTCTGGAGCAACGAAAGAATTAGTGATTTTTCTGGTAAAAATTATCCAGAAGATGTTGTAAATCAAGCGCATATAAAAGCGCGGCAACTAAACGTAGAACCTGTCTTATTGAGCTTGCCGGCGAGAGTTCAATGGGACAGGTTCCAGaa                          <  1:35456/251‑1 (MQ=255)
                                                                                                                                                                                                                         cGGAAGAAGCGCGCGAAGCATCGCATCGCTATCTGGCCTTTATCGAAGAAATTTTGCTCGACAGAAGTCGTGAAGAGAGCCGCCGTGAGCGTTCTCTGCGTCGTCTGGAGCAACGAAAGAATTAGTGATTTTTCTGGTAAAAATTATCCAGAAGATGTTGTAAATCAAGCGCATATAAAAGCGCGGCAACTAAACGTAGAACCTGTCTTATTGAGCTTTCCGGCGAGAGTTCAATGGGACAGGTTCCAGa                           >  2:363769/1‑250 (MQ=255)
                                                                                                                                                                                                                            aagaagCGCGCGAAGCATCGCATCGCTATCTGGCCTTTATCGAAGAAATTTTGCTCGACAGAAGTCGTGAAGAGAGCCGCCGTGAGCGTTCTCTGCGTCGTCTGGAGCAACGAAAGAATTAGTGATTTTTCTGGTAAAAATTATCCAGAAGATGTTGTAAATCAAGCGCATATAAAAGCGCGGCAACTAAACGTAGAACCTGTCTTATTGAGCTTTCCGGCGAGAGTTCAATGGGACAGGTTCCAGAAAAc                       <  2:11739/251‑1 (MQ=255)
                                                                                                                                                                                                                               acgcgcgcGAAGCATCGCATCGTTATCTGGCCTTTATCGAAGAAATTTTGCTCGACAGAAGTCGTGAAGAGAGCCGCCGTGAGCGTTCTCTGCGTCGTCTGGAGCAACGAAAGAATTAGTGATTTTTCTGGTAAAAATTATCCAGAAGATGTTGTAAATCAAGCGCATATAAAAGCGCGGCAACTAAACGTAGAACCTGTCTTATTGAGCTTTCCGGCGAGAGTTCAATGGGACAGGTTCCAGAAAACTCa                    <  1:292285/249‑1 (MQ=255)
                                                                                                                                                                                                                                     gcgAAGCATCGCATCGCTATCTGGCCTTTATCGAAGAAATTTTGCTCGACAGAAGTCGTGAAGAGAGCCGCCGTGAGCGTTCTCTGCGTCGTCTGGAGCAACGAAAGAATTAGTGACTTTTCTGGTAAAAATTATCCAGAAGATGTTGTAAATCAAGCGCATATAAAAGCGCGGCAACTAAACGTAGAACCTGTCTTATTGAGCTTTCCGGCGAGAGTTCAATGGGACAGGTTCCAGAAAACTCAACGtta              <  1:50835/251‑1 (MQ=255)
                                                                                                                                                                                                                                                 atcgcTATCTGGCCTTTATCGAAGAAATTTTGCTCGACAGAAGTCGTGAAGAGAGCCGCCGTGAGCGTTCTCTGCGTCGTCTGGAGCAACGAAAGAATTAGTGATTTTTCTGGTAAAAATTATCCAGAAGATGTTGTAAATCAAGCGCATATAAAAGCGCGGCAACTAAACGTAGAACCTGTCTTATTGAGCTTTCCGGCGAGAGTTCAATGGGACAGGTTCCAGAAAACTCAACGTTATTAGATAGATaa  <  1:461215/251‑1 (MQ=255)
                                                                                                                                                                                                                                                      |                                                                                                                                                                                                                                                     
GATCTGGACGCGGAATCAAACGCCGTACTCCAGTATCAGATTGCCGTCACCGAAGCGGCCCACAATGTGGTTCTGCTTCATCTGCTAAGGTGTATGGAGCCGATGTTGGCCCAGAATGTCCGCCAGAACTTCGAATTGCTCTATTCGCGTCGCGAGATGCTGCCGCTGGTGAGTAGTCACCGCACCCGCATATTTGAAGCGATTATGGCCGGTAAGCCGGAAGAAGCGCGCGAAGCATCGCATCGCCATCTGGCCTTTATCGAAGAAATTTTGCTCGACAGAAGTCGTGAAGAGAGCCGCCGTGAGCGTTCTCTGCGTCGTCTGGAGCAACGAAAGAATTAGTGATTTTTCTGGTAAAAATTATCCAGAAGATGTTGTAAATCAAGCGCATATAAAAGCGCGGCAACTAAACGTAGAACCTGTCTTATTGAGCTTTCCGGCGAGAGTTCAATGGGACAGGTTCCAGAAAACTCAACGTTATTAGATAGATAA  >  NC_000913/122515‑123006

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: