Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 1,667,887 Δ1 bp coding (678/1221 nt) mlc ← DNA‑binding transcriptional repressor Mlc

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,667,8870C.97.2% 313.7 / ‑1.6 72coding (678/1221 nt)mlcDNA‑binding transcriptional repressor Mlc
Reads supporting (aligned to +/- strand):  ref base C (2/0);  new base . (26/44);  total (28/44)
Fisher's exact test for biased strand distribution p-value = 1.48e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.87e-01

CGCCGCGCAATGCCGCCTGACACAATGAGTCCACGGTTAACGGTTGTCCATGTAACATCGAGCTCATGGATTGATTAAGACGCAGCTGTGCCAGCTCAAGAATACTGTCCACGCTGGCGATGGTTTCGAGGCAGCCGTGATTCCCGCAATAACAGCGTTTCCCATACGGGTCGACCTGTGTGTGGCCTATTTCCACGAGACTACTGCTGCCTGCGTGTAGCAGATGACCATCGGTAATGACGCCCGCCCCCACGTTGTGATCGATAACCACCTGAATCACATCGCGCGCCCCGCGTGAGGCACCAAACAAGGCCTCTGCCATCGTCCATGCGCTGATATCATGCTGAATATAAACCGGAACGCCGGTATGCTGCTCCAGCGCCTCGCCGAGCGGCATCTCTTTTACATCCTCGTAGAACGGCATGCGATGTACAATACCATTTTCCGTATCAATAATTCCCGGCAAGGTTATGGCAATCGAAGTTAGACGC  >  NC_000913/1667643‑1668133
                                                                                                                                                                                                                                                    |                                                                                                                                                                                                                                                      
cgccgcGCAATGCCGCCTGACACAATGAGTCCACGGTTAACGGTTGTCCATGTAACATCGAGCTCATGGATTGATTAAGACGCAGCTGTGCCAGCTCAAGAATACTGTCCACGCTGGCGATGGTTTCGAGGCAGCCGTGATTCCCGCAATAACAGCGTTTCCCATACGGGTCGACCTGTGTGTGGCCTATTTCCACGAGACTACTGCTGCCTGCGTGTAGCAGATGACCATCGGTAATGACGCCGccccc                                                                                                                                                                                                                                                  <  1:334176/250‑1 (MQ=255)
      gcAATGCCGCCTGACACAATGAGTCCACGGTTAACGGTTGTCCATGTAACATCGAGCTCATGGATTGATTAAGACGCAGCTGTGCCAGCTCAAGAATACTGTCCACGCTGGCGATGGTTTCGAGGCAGCCGTGATTCCCGCAATAACAGCGTTTCCCATACGGGTCGACCTGTGTGTGGCCTATTTCCACGAGACTACTGCTGCCTGCGTGTAGCAGATGACCATCGGTAATGACGCCGCCCCCACGTtgt                                                                                                                                                                                                                                           <  2:421596/251‑1 (MQ=255)
            ccgccTGACACAATGAGTCCACGGTTAACGGTTGTCCATGTAACATCGAGCTCATGGATTGATTAAGACGCAGCTGTGCCAGCTCAAGAATACTGTCCACGCTGGCGATGGTTTCGAGGCAGCCGTGATTCCCGCAATAACAGCGTTTCCCATACGGGTCGACCTGTGTGTGGCCTATTTCCACGAGACTACTGCTGCCTGCGTGTAGCAGATGACCATCGGTAATGACGCCGCCCCCACGTTGTGATCg                                                                                                                                                                                                                                      >  1:327649/1‑250 (MQ=255)
              gccTGACACAATTAGTCCACGGTTAACGGTTGTCCAGGTAACATCGAGCTCATGGATTGATTAAGACGCAGCTGTGCCAGCTCAAGAAGACTGTCCACGCTGGCGATGGTTTCGAGGCAGCCGTGATTCCCGCAATAACAGCGTTTCCCATACGGGTCGATCGGTGTGTGGCCTAGTTCCACGAGACTACTGCTGCCTGCGTGTATCAGATGACCATCGGTAATGACGCCGCTCCCACGTTGTGATCGATa                                                                                                                                                                                                                                   <  2:285634/251‑1 (MQ=255)
              gccTGACACAATGAGTCCACGGTTAACGGTTGTCCATGTAACATCGAGCTCATGGATTGATTAAGACGCAGCTGTGCCAGCTCAAGAATACTGTCCACGCTGGCGATGGTTTCGAGGCAGCCGTGATTCCCGCAATAACAGCGTTTCCCATACGGGTCGACCTGTGTGTGGCCTATTTCCACGAGACTACTGCTGCCTGCGTGTAGCAGATGACCATCGGTAATGACGCCGCCCCCACGTTGTGATCGATa                                                                                                                                                                                                                                   <  1:124248/251‑1 (MQ=255)
               ccTGACACAATGAGTCCACGGTTAACGGTTGTCCATGTAACATCGAGCTCATGGATTGATTAAGACGCAGCTGTGCCAGCTCAAGAATACTGTCCACGCTGGCGATGGTTTCGAGGCAGCCGTGATTCCCGCAATAACAGCGTTTCCCATACGGGTCGACCTGTGTGTGGCCTATTTCCACGAGACTACTGCTGCCTGCGTGTAGCAGATGACCATCGGTAATGACGCCGCCCCCACGTTGTGATCGATaa                                                                                                                                                                                                                                  <  1:459683/251‑1 (MQ=255)
               ccTGACACAATGAGTCCACGGTTAACGGTTGTCCATGTAACATCGAGCTCATGGATTGATTAAGACGCAGCTGTGCCAGCTCAAGAATACTGTCCACGCTGGCGATGGTTTCGAGGCAGCCGTGATTCCCGCAATAACAGCGTTTCCCATACGGGTCGACCTGTGTGTGGCCTATTTCCACGAGACTACTGCTGCCTGCGTGTAGCAGATGACCATCGGTAATGACGCCGCCCCCACGTTGTGATCGATaa                                                                                                                                                                                                                                  <  2:519579/251‑1 (MQ=255)
                        aggAGTCCACGGTTAACGGTTGTCCATGTAACGTCGAGCTCATGGATTGATTAAGACGCAGCTGTGCCAGCTCAAGAATACTGTCCACGCTGGCGATGGTTTCGAGGCAGCCGTGATTCCCGCAATAACAGCGTTTCCCATACGGGTCGACCTGTGTGTGGCCTATTTCCACGAGACTACTGCTGCCTGTGTGTAGCAAATGACCATCGGTAAGGACGCCGCCCCCTCGTTGTGATCGATAACCCCCTGaa                                                                                                                                                                                                                         <  2:477974/249‑1 (MQ=255)
                        aTGAGTCCACGGTTAACGGTTGTCCATGTAACATCGAGCTCATGGATTGATTAAGACGCAGCTGTGCCAGCTCAAGAATACTGTCCACGCTGGCGATGGTTTCGAGGCAGCCGTGATTCCCGCAATAACAGCGTTTCCCATACGGGTCGACCTGTGTGTGGCCTATTTCCACGAGACTACTGCTGCCTGCGTGTAGCAGATGACCATCGGTAATGACGCCGCCCCCACGTTGTGATCGATAACCACCTGaa                                                                                                                                                                                                                         <  1:594435/251‑1 (MQ=255)
                               cACGGTTAACGGTTGTCCATGTAACATCGAGCTCATGGATTGATTAAGACGCAGCTGTGCCAGCTCAAGAATACTGTCCACGCTGGCGATGGTTTCGAGGCAGCCGTGATTCCCGCAATAACAGCGTTTCCCATACGGGTCGACCTGTGTGTGGCCTATTTCCACGAGACTACTGCTGCCTGCGTGTAGCAGATGACCATCGGTAATGACGCCGCCCCCACGTTGTGATCGATAACCACCTGAATCACATc                                                                                                                                                                                                                  >  1:223293/1‑251 (MQ=255)
                                aCGGTTAACGGTTGTCCATGTAACATCGAGCTCATGGATTGATTAAGACGCAGCTGTGCCAGCTCAAGAATACTGTCCACGCTGGCGATGGTTTCGAGGCAGCCGTGATTCCCGCAATAACAGCGTTTCCCATACGGGTCGACCTGTGTGTGGCCTATTTCCACGAGACTACTGCTGCCTGCGTGTAGCAGATGACCATCGGTAATGACGCCGCCCCCACGTTGTGATCGATAACCACCTGAATCACATcg                                                                                                                                                                                                                 >  1:30256/1‑251 (MQ=255)
                                 cGGTTAACGGTTGTCCATGTAACATCGAGCTCATGGATTGATTAAGACGCAGCTGTGCCAGCTCAAGAATACTGTCCACGCTGGCGATGGTTTCGAGGCAGCCGTGATTCCCGCAATAACAGCGTTTCCCATACGGGTCGACCTGTGTGTGGCCTATTTCCACGAGACTACTGCTGCCTGCGTGTAGCAGATGACCATCGGTAATGACGCCGCCCCCACGTTGTGATCGATAACCACCTGAATCACATcgc                                                                                                                                                                                                                >  2:68470/1‑251 (MQ=255)
                                 cGGTTAACGGTTGTCCATGTAACATCGAGCTCATGGATTGATTAAGACGCAGCTGTGCCAGCTCAAGAATACTGTCCACGCTGGCGATGGTTTCGAGGCAGCCGTGATTCCCGCAATAACAGCGTTTCCCATACGGGTCGACCTGTGTGTGGCCTATTTCCACGAGACTACTGCTGCCTGCGTGTAGCAGATGACCATCGGTAATGACGCCGCCCCCACGTTGTGATCGATAACCACCTGAATCACATcgc                                                                                                                                                                                                                >  1:313189/1‑251 (MQ=255)
                                      aaCGGTTGTCCATGTAACATCGAGCTCATGGATTGATTAAGACGCAGCTGTGCCAGCGCAAGAATACTGTCCACGCTGGCGATGGTTTCGAGGCAGCCGTGATTCCCGCAATAACAGCGTTTCCCATACGGGTCGACCTGTGTGTGGCCTATTTCCACGAGACTACTGCTGCCTGCGTGTAGCAGATGACCATCGGTAATGACGCCGCCCCCACGTTGTGATCGATAACCACCTGAATCACATCGCGCGcc                                                                                                                                                                                                           <  1:445286/251‑1 (MQ=255)
                                            tGTCCATGTAACATCGAGCTCATGGATTGATTAAGACGCAGCTGTGCCAGCTCAAGAATACTGTCCACGCTGGCGATGGTTTCGAGGCAGCCGTGATTCCCGCAATAACAGCGTTTCCCATACGGGTCGACCTGTGTGTGGCCTATTTCCACGAGACTACTGCTGCCTGCGTGTAGCAGATGACCATCGGTAATGACGCCGCCCCCACGTTGTGATCGATAACCACCTGAATCACATCGCGCGCCCCGCGt                                                                                                                                                                                                     >  2:506580/1‑251 (MQ=255)
                                            tGTCCATGTAACATCGAGCTCATGGATTGATTAAGACGCAGCTGTGCCAGCTCAAGAATACTGTCCACGCTGGCGATGGTTTCGAGGCAGCCGTGATTCCCGCAATAACAGCGTTTCCCATACGGGTCGACCTGTGTGTGGCCTATTTCCACGAGACTACTGCTGCCTGCGTGTAGCAGATGACCATCGGGAATGACGCCGCCCCCACGTTGTGATCGATAACCACCTGAATCACATCGCGCGCCCcgcg                                                                                                                                                                                                      <  2:327649/250‑1 (MQ=255)
                                            ggTCCATGTAACATCGAGCTCATGGATTGATTAAGACGCAGCTGTGCCAGCTCAAGAATACTGTCCACGCTGGCGATGGTTTCGAGGCAGCCGTGATTCCCGCAATAACAGCGTTTCCCATACGGGTCGACCTGTGTGTGGCCTATTCCCACGAGACTACTGCTGCCTGCGTGTAGCAGATGACCATCGGTAATGACGCCGCCCCCACGTTGTGATCGATAACCACCTGAATCACATCGCGCGCCCCGCGt                                                                                                                                                                                                     <  2:284237/250‑1 (MQ=255)
                                             gTCCATGTAACATCGAGCTCATGGATTGATTAAGACGCAGCTGTGCCAGCTCAAGAATACTGTCCACGCTGGCGATGGTTTCGAGGCAGCCGTGATTCCCGCAATAACAGCGTTTCCCATACGGGTCGACCTGTGTGTGGCCTATTTCCACGAGACTACTGCTGCCTGCGTGTAGCAGATGACCATCGGTAATGACGCCGCCCCCACGTTGTGATCGATAACCACCTGAATCACATCGCGCGCCCCGCGTg                                                                                                                                                                                                    <  2:313189/251‑1 (MQ=255)
                                             gTCCATGTAACATCGAGCTCATGGATTGATTAAGACGCAGCTGTGCCAGCTCAAGAATACTGTCCACGCTGGCGATGGTTTCGAGGCAGCCGTGATTCCCGCAATAACAGCGTTTCCCATACGGGTCGACCTGTGTGTGGCCTATTTCCACGAGACTACTGCTGCCTGCGTGTAGCAGATGACCATCGGTAATGACGCCGCCCCCACGTTGTGATCGATAACCACCTGAATCACATCGCGCGCCCCGCGTg                                                                                                                                                                                                    <  2:390504/251‑1 (MQ=255)
                                                cATGTAACATCGAGCTCATGGATTGATTAAGACGCAGCTGTGCCAGCTCAAGAATACTGTCCACGCTGGCGATGGATTCGAGGCAGCCGTGATTCCCGCAATAACAGCGTTTCCCATACGGGTTGACCTGTGTGTGGCCTATTTCCACGAGACTACTGCTGCCTGCGTGTAGCAGATGACCATCGGTAATGACGCCGCCCCCACGTTGTGATCGATAACCACCTGAATCACATCGCGCGCCCCGCGTGAgg                                                                                                                                                                                                 >  1:264400/1‑251 (MQ=255)
                                                   gTAACATCGAGCTCATGGATTGATTAAGACGCAGCTGTGCCAGCTCAAGAATACTGTCCACGCTGGCGATGGTTTCGAGGCAGCCGTGATTCCCGCAATAACAGCGTTTCCCATACGGGTCGACCTGTGTGTGGCCTATTTCCACGAGACTACTGCTGCCTGCGTGTAGCAGATGACCATCGGTAATGACGCCGCCCCCACGTTGTGATCGATAACCACCTGAATCACATCGCGCGCCCCGCGTGAGGCAc                                                                                                                                                                                              <  1:601424/251‑1 (MQ=255)
                                                     aaCATCGAGCTCATGGATTGATTAAGACGCAGCTGTGCCAGCTCAAGAATACTGTCCACGCTGGCGATGGTTTCGAGGCAGCCGTGATTCCCGCAATAACAGCGTTTCCCATACGGGTCGACCTGTGTGTGGCCTATTTCCACGAGACTACTGCTGCCTGCGTGTAGCAGATGACCATCGGTAATGACGCCGCCCCCACGTTGTGATCGATAACCACCTGAATCACATCGCGCGCCCCGCGTGAGGCAcc                                                                                                                                                                                             >  2:398449/1‑250 (MQ=255)
                                                      aCATCGAGCTCATGGATTGATTAAGACGCAGCTGTGCCAGCTCAAGAATACTGTCCACGCTGGCGATGGTTTCGAGGCAGCCGTGATTCCCGCAATAACAGCGTTTCCCATACGGGTCGACCTGTGTGTGGCCTATTTCCACGAGACTACTGCTGCCTGCGTGTAGCAGATGACCATCGGTAATGACGCCGCCCCCACGTTGTGATCGATAACCACCTGAATCACATCGCGCGCCCCGCGTGAGGCACCa                                                                                                                                                                                            >  2:112033/1‑250 (MQ=255)
                                                          cGAGCTCATGGATTGATTAAGACGCAGCTGTGCCAGCTCAAGAATACTGTCCACGCTGGCGATGGTTTCGAGGCAGCCGTGATTCCCGCAATAACAGCGTTTCCCATACGGGTCGACCTGTGTGTGGCCTATTTCCACGAGACTACTGCTGCCTGCGTGTAGCAGATGACCATCGGTAATGACGCCGCCCCCACGTTGTGATCGATAACCACCTGAATCACATCGCGCGCCCCGCGTGAGGCACCAAACa                                                                                                                                                                                        >  1:209499/1‑250 (MQ=255)
                                                                 aTGGATTGATTAAGACGCAGCTGTGCCAGCTCAAGAATACTGTCCACGCTGGCGATGGTTTCGAGGCAGCCGTGATTCCCGCAATAACAGCGTTTCCCATACGGGTCGACCTGTGTGTGGCCTATTTCCACGAGACTACTGCTGCCTGCGTGTAGCAGATGACCATCGGTAATGACGCCGCCCCCACGTTGTGATCGATAACCACCTGAATCACATCGCGCGCCCCGCGTGAGGCACCAAACAAGGCctct                                                                                                                                                                                <  2:465780/251‑1 (MQ=255)
                                                                 aTGGATTGATTAAGACGCAGCTGTGCCAGCTCAAGAATACTGTCCACGCTGGCGATGGTTTCGAGGCAGCCGTGATTCCCGCAATAACAGCGTTTCCCATACGGGTCGACCTGTGTGTGGCCTATTTCCACGAGACTACTGCTGCCTGCGTGTAGCAGATGACCATCGGTAATGACGCCGCCCCCACGTTGTGATCGATAACCACCTGAATCACATCGCGCGCACCGCGTGAGGCACCAAACAAGGCctct                                                                                                                                                                                >  2:277381/1‑251 (MQ=255)
                                                                       tgattAAGACGCAGCTGTGCCAGCTCAAGAATACTGTCCACGCTGGCGATGGTTTCGAGGCAGCCGTGATTCCCGCAATAACAGCGTTTCCCATACGGGTCGACCTGTGTGTGGCCTATTTCCACGAGACTACTGCTGCCTGCGTGTAGCAGATGACCATCGGTAATGACGCCGCCCCCACGTTGTGATCGATAACCACCTGAATCACATCGCGCGCCCCGCGTGAGGCACCAAACAAGGCCTCTGCCATc                                                                                                                                                                          >  1:559055/1‑251 (MQ=255)
                                                                         attAAGACGCAGCTGTGCCAGCTCAAGAATACTGTCCACGCTGGCGATGGTTTCGAGGCAGCCGTGATTCCCGCAATAACAGCGTTTCCCATACGGGTCGACCTGTGTGTGGCCTATTTCCACGAGACTACTGCTGCCTGCGTGTAGCAGATGACCATCGGTAATGACGCCGCCCCCACGTTGTGATCGATAACCACCTGAATCACATCGCGCGCCCCGCGTGAGGCACCAAACAAGGCCTCTGCCATCGt                                                                                                                                                                        <  2:215839/251‑1 (MQ=255)
                                                                             aGACGCAGCTGTGCCAGCTCAAGAATACTGTCCACGCTGGCGATGGTTTCGAGGCAGCCGTGATTCCCGCAATAACAGCGTTTCCCATACGGGTCGACCTGTGTGTGGCCTATTTCCACGAGACTACTGCTGCCTGCGTGTAGCAGATGACCATCGGTAATGACGCCGCCCCCACGTTGTGATCGATAACCACCTGAATCACATCGCGCGCCCCGCGTGAGGCACCAAACAAGGCCTCTGCCATCGTCCAt                                                                                                                                                                    <  2:494735/251‑1 (MQ=255)
                                                                             aGACGCAGCTGTGCCAGCTCAAGAATACTGTCCACGCTGGCGATGGTTTCGAGGCAGCCGTGATTCCCGCAATAACAGCGTTTCCCATACGGGTCGACCTGTGTGTGGCCTATTTCCACGAGACTACTGCTGCCTGCGTGTAGCAGATGACCATCGGTAATGACGCCGCCCCCACGTTGTGATCGATAACCACCTGAATCACATCGCGCGCCCCGCGTGAGGCACCAAACAAGGCCTCTGCCATCGTCCAt                                                                                                                                                                    >  2:54173/1‑251 (MQ=255)
                                                                              gACGCAGCTGTGCCAGCTCAAGAATACTGTCCACGCTGGCGATGGTTTCGAGGCAGCCGTGATTCCCGCAATAACAGCGTTTCCCATACGGGTCGACCTGTGTGTGGCCTATTTCCACGAGACTACTGCTGCCTGCGTGTAGCAGATGACCATCGGTAATGACGCCGCCCCCACGTTGTGATCGATAACCACCTGAATCACATCGCGCGCCCCGCGTGAGGCACCAAACAAGGCCTCTGCCATCGTCCATg                                                                                                                                                                   >  2:69610/1‑251 (MQ=255)
                                                                                 gCAGCTGTGCCAGCTCAAGAATACTGTCCACGCTGGCGATGGTTTCGAGGCAGCCGTGCTTCCCGCAATAACAGCGTTTCCCATACGGGTCGACCTGTGTGTGGCCTATTTCCACGAGACTACTGCTGCCTGCGTGTAGCAGATGACCATCGGTAATGACGCCGCCCCCACGTTGTGATCGATAACCACCTGAATCACATCGCGCGCCACGCGTGAGGCACCAAACAAGGCCTCTGCCATCGTCCATgcgc                                                                                                                                                                >  2:457622/1‑251 (MQ=255)
                                                                                 gCAGCTGTGCCAGCGCAAGAATACTGTCCACGCTGGCGATGGTTTCGAGGCAGCCGTGATTCCCGCAATAACAGCGTTTCCCATACGGGTCGACCTGTGTGTGGCCTATTTCCACGAGACTACTGCTGCCTGCGTGTAGCAGATGACCATCGGTAATGACGCCGCCCCCACGTTGTGATCGATAACCACCTGAATCACATCGCGCGCCCCGCGGGAGGCACCAAACAAGGCCTCTGCCATCGTCCATgcgc                                                                                                                                                                <  2:455262/251‑1 (MQ=255)
                                                                                            aGCTCAAGAATACTGTCCACGCTGGCGATGGTTTCGAGGCAGCCGTGATTCCCGCAATAACAGCGTTTCCCATACGGGTCGACCTGTGTGTGGCCTATTTCCACGAGACTACTGCTGCCTGCGTGTAGCAGATGACCATCGGTAATGACGCCGCCCCCACGTTGTGATCGATAACCACCTGAATCACATCGCGCGCCCCGCGTGAGGCACCAAACAAGGCCTCTGCCATCGTCCATGCGCTGATATCATGc                                                                                                                                                     >  1:485980/1‑251 (MQ=255)
                                                                                                cAAGAATACTGTCCACGCTGGCGATGGTTTCGAGGCAGCCGTGATTCCCGCAATAACAGCGTTTCCCATACGGGTCGACCTGTGTGTGGCCTATTTCCACGAGACTACTGCTGCCTGCGTGTAGCAGATGACCATCGGTAATGACGCCGCCCCCACGTTGTGATCGATAACCACCTGAATCACATCGCGCGCCCCGCGTGAGGCACCAAACAAGGCCTCTGCCATCGTCCATGCGCTGATATCATGCTGa                                                                                                                                                  >  2:485625/1‑250 (MQ=255)
                                                                                                cAAGAATACTGTCCACGCTGGCGATGGTTTCGAGGCAGCCGTGATTCCCGCAATAACAGCGTTTCCCATACGGGTCGACCTGTGTGTGGCCTATTTCCACGAGACTACTGCTGCCTGCGTGTAGCAGATGACCATCGGTAATGACGCCGCCCCCACGTTGTGATCGATAACCACCTGAATCACATCGCGCGCCCCGCGTGAGGCACCAAACAAGGCCTCTGCCATCGTCCATGCGCTGATATCATGCTGa                                                                                                                                                  <  1:23461/250‑1 (MQ=255)
                                                                                                aaaGAATACTGTCCACGCTGGCGATGGTTTCGAGGCAGCCGTGATTCCCGCAATAACAGCGTTTCCCATACGGGTCGACCTGTGTGTGGCCTATTTCCACGAGACTACTGCTGCCTGCGTGTAGCAGATGACCATCGGTAATGACGCCGCCCCCACGTTGTGATCGATAACCACCTGAATCACATCGCGCGCCCCGCGTGAGGCACCAAACAAGGCCTCTGCCATCGTCCATGCGCTGATATCATGCTGaa                                                                                                                                                 <  2:478847/250‑1 (MQ=255)
                                                                                                       aCTGTCCACGCTGGCGATGGTTTCGAGGCAGCCGTGATTCCCGCAATAACAGCGTTTCCCATACGGGTCGACCTGTGTGTGGCCTATTTCCACGAGACTACTGCTGCCTGCGTGTAGCAGATGACCATCGGTAATGACGCCGCCCCCACGTTGTGATCGATAACCACCTGAATCACATCGCGCGCCCCGCGTGAGGCACCAAACAAGGCCTCTGCCATCGTCCATGCGCTGATATCATGCTGAATATAAAc                                                                                                                                          <  1:51266/251‑1 (MQ=255)
                                                                                                            ccccGCTGGCGATGGTTTCGGGGCAGCCGTGATTCCCGCAATAACAGCGTTTCCCATACGGGTCGACCTGTGTGGGGCCTATTTCCACGGGACTGCTTCTGCCTGCGTGTAGCAGATGACCATCGGTAATGACGCCGCCCCCACGTTGTGATCGATAACCACCTGAATCACATCGCGCGCCCCGCGTGAGGCCCCAAACAAGGCCTCTGCCATCGTCCATGCGCTGATATCATGCTGAATATAAACCGGa                                                                                                                                      <  2:348106/247‑1 (MQ=255)
                                                                                                            ccACGCTGGCGATGGTTTCGAGGCAGCCGTGATTCCCGCAATAACAGCGTTTCCCATACGGGTCGACCTGTGTTTGGCCTATTTCCACGAGACTACTGCTGCCTGCGTGTAGCAGATGACCATCGGTAATGACGCCGCCCCCTCTTTGTGATCGATAACCACCTGAATCACATCGCGCGCCCCTCGTGAGGCACCAATCAAGTCCTCTGCCATCGTCCATGTGCTGATATCATGCTGAATATAAATCggta                                                                                                                                     <  2:98996/251‑3 (MQ=255)
                                                                                                                   ggCGATGGTTTCGAGGCAGCCGTGATTCCCGCAATAACAGCGTTTCCCATACGGGTCGACCTGTGTGTGGCCTATTTCCACGAGACTACTGCTGCCTGCGTGTAGCAGATGACCATCGGTAATGACGCCGCCCCCACGTTGTGATCGATAACCACCTGAATCACATCGCGCGCCCCGCGTGAGGCACCAAACAAGGCCTCTGCCATCGTCCATGCGCTGATATCATGCTGAATATAAACCGGAACGCCgg                                                                                                                               <  1:389501/250‑1 (MQ=255)
                                                                                                                      gATGGTTTCGAGGCAGCCGTGATTCCCGCAATAACAGCGTTTCCCATACGGGTCGACCTGTGTGTGGCCTATTTCCACGAGACTACTGCTGCCTGCGTGTAGCAGATGACCATCGGTAATGACGCCGCCCCCACGTTGTGATCGATAACCACCTGAATCACATCGCGCGCCCCGCGTGAGGCACCAAACAAGGCCTCTGCCATCGTCCATGCGCTGATATCATGCTGAATATAAACCGGAACGCCGGTAtg                                                                                                                           >  1:414240/1‑251 (MQ=255)
                                                                                                                       aTGGTTTCGAGGCAGCCGTGATTCCCGCAATAACAGCGTTTCCCATACGGGTCGACCTGTGTGTGGCCTATTTCCACGAGACTACTGCTGCCTGCGTGTAGCAGATGACCATCGGTAATGACGCCGCCCCCACGTTGTGATCGATAACCACCTGAATCACATCGCGCGCCCCGCGTGAGGCACCAAACAAGGCCTCTGCCATCGTCCATGCGCTGATATCATGCTGAATATAAACCGGAACGCCGGTAtgc                                                                                                                          >  2:307130/1‑251 (MQ=255)
                                                                                                                        tGGTTTCGAGGCAGCCGTGATTCCCGCAATAACAGCGTTTCCCATACGGGTCGACCTGTGTGTGGCCTATTTCCACGAGACTACTGCTGCCTGCGTGTAGCAGATGACCATCGGTAATGACGCCGCCCCCACGTTGTGATCGATAACCACCTGAATCACATCGCGCGCCCCGCGTGAGGCACCAAACAAGGCCTCTGCCATCGTCCATGCGCTGATATCATGCTGAATATAAACCGGAACGCCGGTAtgct                                                                                                                         >  2:635115/1‑251 (MQ=255)
                                                                                                                                 ggCAGCCGTGATTCCCGCAATAACAGCGTTTCCCATACGGGTCGCCCTGTGTGTGGCCTATTTCCACGAGACTACTGCTGCCTGCGTGTAGCAGATGACCATCGGTAATGACGCCGCCCCCACGTTGTGATCGATAACCACCTGAATCACATCGCGCGCCCCGCGTGAGGCACCAAACAAGGCCTCTGCCATCGTCCATGCGCTGATATCATGCTGAATATAAACCGGAACGCCGGTATGCTGCTCCAgcg                                                                                                                <  1:635115/251‑1 (MQ=255)
                                                                                                                                  gCAGCCGTGATTCCCGCAATAACAGCGTTTCCCATACGGGTCGACCTGTGTGTGGCCTATTTCCACGAGACTACTGCTGCCTGCGTGTAGCAGATGACCATCGGTAATGACGCCGCCCCCACGTTGTGATCGATAACCACCTGAATCACATCGCGCGCCCCGCGTGAGGCACCAAACAAGGCCTCTGCCATCGTCCATGCGCTGATATCATGCTGAATATAAACCGGAACGCCGGTATGCTGCTCCAgcgc                                                                                                               <  2:414240/251‑1 (MQ=255)
                                                                                                                                     gCCGTTATTCCCGCAATAACAGCGTTTCCCATTCGGGTCGACCTGTGTGTGGCCTATTTCCACGAGACTGCTGCTGCCTGCGTGTAGCAGATGACCATCGGTAATGACGCCGCCCCCACGTTGTGATCGATAACCACCCGAATCACATCGCGCGCCCCGCGTGAGGCACCAAACAAGGCCTCTGCCATCGTCCATGCGCTGATATCATGCTGAATATAAACCGGAACGCCGGTATGCTGCTCCAGCGCCTc                                                                                                            <  2:63385/251‑1 (MQ=255)
                                                                                                                                        gTGATTCCCGCAATAACAGCGTTTCCCATACGGGTCGACCTGTGTGTGGCCTATTTCCACGAGACTACTGCTGCCTGCGTGTAGCAGATGACCATCGGTAATGACGCCGCCCCCACGTTGTGATCGATAACCACCTGAATCACATCGCGCGCCCCGCGTGAGGCACCAAACAAGGCCTCTGCCATCGTCCATGCGCTGATATCATGCTGAATATAAACCGGAACGCCGGTATGCTGCTCCAGCGCCTCGcc                                                                                                         <  1:3059/251‑1 (MQ=255)
                                                                                                                                           aTTCCCGCAATAACAGCGTTTCCCATACGGGTCGACCTGTGTGTGGCCTATTTCCACGAGACTACTGCTGCCTGCGTGTAGCAGATGACCATCGGTAATGACGCCGCCCCCACGTTGTGATCGATAACCACCTGAATCACATCGCGCGCCCCGCGTGAGGCACCAAACAAGGCCTCTGCCATCGTCCATGCGCTGATATCATGCTGAATATAAACCGGAACGCCGGTATGCTGCTCGAGCGCCTCGCCGAg                                                                                                      <  1:202543/251‑1 (MQ=255)
                                                                                                                                                  cAATAACAGCGTTTCCCATACGGGTCGACCTGTGTGTGGCCTATTTCCACGAGACTACTGCTGCCTGCGTGTAGCAGATGACCATCGGTAATGACGCCGCCCCCACGTTGTGATCGATAACCACCTGAATCACATCGCGCGCCCCGCGTGAGGCACCAAACAAGGCCTCTGCCATCGTCCATGCGCTGATATCATGCTGAATATAAACCGGAACGCCGGTATGCTGCTCCAGCGCCTCGCCGAGCGGCAtc                                                                                               >  1:334621/1‑251 (MQ=255)
                                                                                                                                                   aaTAACAGCGTTTCCCATACGGGTCGACCTGTGTGTGGCCTATTTCCACGAGACTACTGCTGCCTGCGTGTAGCAGATGACCATCGGTAATGACGCCGCCCCCACGTTGTGATCGATAACCACCTGAATCACATCGCGCGCCCCGCGTGAGGCACCAAACAAGGCCTCTGCCATCGTCCATGCGCTGATATCATGCTGAATATAAACCGGAACGCCGGTATGCTGCTCCAGCGCCTCGCCGAGCGGCAtct                                                                                              >  1:426280/1‑251 (MQ=255)
                                                                                                                                                      aaCAGCGTTTCCCATACGGGTCGACCTGTGTGTGGCCTATTTCCACGAGACTACTGCTGCCTGCGTGTAGCAGATGACCATCGGTAATGACGCCGCCCCCACGTTGTGATCGATAACCACCTGAATCACATCGCGCGCCCCGCGTGAGGCACCAAACAAGGCCTCTGCCATCGTCCATGCGCTGATATCATGCTGAATATAAACCGGAACGCCGGTATGCTGCTCCAGCGCCTCGCCGAGCGGCATctct                                                                                            <  2:541703/250‑1 (MQ=255)
                                                                                                                                                      aaCAGCGTTTCCCATACGGGTCGACCTGTGTGTGGCCTATTTCCACGAGACTACTGCTGCCTGCGTGTAGCAGATGACCATCGGTAATGACGCCGCCCCCACGTTGTGATCGATAACCACCTGAATCACATCGCGCGCCCCGCGTGAGGCACCAAACAAGGCCTCTGCCATCGTCCATGCGCTGATATCATGCTGAATATAAACCGGAACGCCGGTATGCTGCTCCAGCGCCTCGCCGAGCGGCATctct                                                                                            <  1:307130/250‑1 (MQ=255)
                                                                                                                                                        cAGCGTTTCCCATACGGGTCGACCTGTGTGTGGCCTATTTCCACGAGACTCCTGCTGCCTGCGTGTAGCAGATGACCATCGGTAATGACGCCGCCCCCACGTTGTGATCGATAACCACCTGAATCACATCGCGCGCCCCGCGTGAGGCACCAAACAAGGCCTCTGCCATCGTCCATGCGCTGATATCATGCTGAATATAAACCGGAACGCCGGTATGCTGCTCCAGCGCCTCGCCGAGCGGCATCTCtttt                                                                                         <  2:438784/251‑1 (MQ=255)
                                                                                                                                                           cGTTTCCCATACGGGTCGACCTGTGTGTGGCCTATTTCCACGAGACTACTGCTGCCTGCGTGTAGCAGATGACCATCGGTAATGACGCCGCCCCCACGTTGTGATCGATAACCACCTGAATCACATCGCGCGCCCCGCGTGAGGCACCAAACAAGGCCTCTGCCATCGTCCATGCGCTGATATCATGCTGAATATAAACCGGAACGCCGGTATGCTGCTCCAGCGCCTCGCCGAGCGGCATCTCTTTTAc                                                                                       >  2:307638/1‑250 (MQ=255)
                                                                                                                                                               tCCCATACGGGTCGACCTGTGTGTGGCCTATTTCCACGAGACTACTGCTGCCTGCGTGTAGCAGATGACCATCGGTAATGACGCCGCCCCCACGTTGTGATCGATAACCACCTGAATCACATCGCGCGCCCCGCGTGAGGCACCAAACAAGGCCTCTGCCATCGTCCATGCGCTGATATCATGCTGAATATAAACCGGAACGCCGGTATGCTGCTCCAGCGCCTCGCCGAGCGGCATCTCTTTTACATcc                                                                                   <  1:594776/250‑1 (MQ=255)
                                                                                                                                                                 ccATACGGGTCGACCTGTGTGTGGCCTATTTCCACGAGACTACTGCTGCCTGCGTGTAGCAGATGACCATCGGTAATGACGCCGCCCCCACGTTGTGATCGATAACCACCTGAATCACATCGCGCGCCCCGCGTGAGGCACCAAACAAGGCCTCTGCCATCGTCCATGCGCTGATATCATGCTGAATATAAACCGGAACGCCGGTATGCTGCTCCAGCGCCTCGCCGAGCGGCATCTCTTTTACATCCTCg                                                                                <  1:277381/251‑1 (MQ=255)
                                                                                                                                                                      cGGGTCGACCTGTGTGTGGCCTATTTCCACGAGACTACTGCTGCCTGCGTGTGGCAGGTGACCATCGGTAATGACACCGCCCCCACGTTGTGATCGATAACCACCTGAATCACATCGCGCGCCCCGCGTGAAGCACCAAACAAGACCTCTGCCATCGTCCATGCGCTGATATCATGCTGAATATAAACCGGAACGCCGGTATGCTGCTCCAGCGCCTCGCCGAGCGGCATCTCTTTTACATCCTCGTAGa                                                                            >  1:452126/1‑250 (MQ=255)
                                                                                                                                                                        gaacGACCTGAGTGTGGCCTATTTCCACGAGACTACTGATTCCTGCGTGTAGCAGATGCCCATCGGTAATAACGCCCCCCCCACGTTGTGATCGATAACCAACTGAATCACATCGCGAGCCCCGCGTGACGCACCAAACAAGGCCTCTGCCATCGTCCATGCGCTGATATCATGCTGAATATAAACCCGAACGCCGGTATGCTGCTCCAGCGCCTCGCCGAGCGGCCTCTCTTTTACATCCTAGTAGAACg                                                                         >  2:450860/4‑251 (MQ=255)
                                                                                                                                                                         gTCGACCTGTGTGTGGCCTATTTCCACGAGACTACTGCTGCCTGCGTGTAGCAGATGACCATCGGTAATGACGCCGCCCCCACGTTGTGATCGATAACCACCTGAATCACATCGCGCGCCCCGCGTGAGGCACCAAACAAGGCCTCTGCCATCGTCCATGCGCTGATATCATGCTGAATATAAACCGGAACGCCGGTATGCTGCTCCAGCGCCTCGCCGAGCGGCATCTCTTTTACATCCTCGTAGAACgg                                                                        >  1:78218/1‑251 (MQ=255)
                                                                                                                                                                           cGACCTGTGTGTGGCCTATTTCCACGAGACTACTGCTGCCTGCGTGTAGCAGATGACCATCGGTAATGACGCCGCCCCCACGTTGTGATCGATAACCACCTGAATCACATCGCGCGCCCCGCGTGAGGCACCAAACAAGGCCTCTGCCATCGTCCATGCGCTGATATCATGCTGAATATAAACCGGAACGCCGGTATGCTGCTCCAGCGCCTCGCCGAGCGGCATCTCTTTTACATCCTCGTAGAACGGc                                                                       >  1:102995/1‑250 (MQ=255)
                                                                                                                                                                                 gtgtgtgGCCTATTTCCACGAGACTACTGCTGCCTGCGTGTAGCAGATGACCATCGGTAATGACGCCGCCCCCACGTTGTGATCGATAACCACCTGAATCACATCGCGCGCCCCGCGTGAGGCACCAAACAAGGCCTCTGCCATCGTCCATGCGCTGATATCATGCTGAATATAAACCGGAACGCCGGTATGCTGCTCCAGCGCCTCGCCGAGCGGCATCTCTTTTACATCCTCGTAGAACGGCATGCg                                                                  <  1:68470/249‑1 (MQ=255)
                                                                                                                                                                                       tgCCTATTTCCACGAGACTACGGCTGCCTGCGTGTAGCAGATGACCATCGGTAATGACGCCGCCCCCACGTTGTGATCGATAACCACCTGAATCACATCGCGCGCCCCGCGTGAGGCACCAAACAAGGCCTCTGCCATCGTCCATGCGCTGATATCATGCTGAATATAAACCGGAACGCCGGTATGCTGCTCCAGCGCCTCGCCGAGCGGCATCTCTTTTACATCCTCGTAGAACGGCATGCGATGTACa                                                           <  2:544103/249‑1 (MQ=255)
                                                                                                                                                                                            aTTTCCACGAGACTACTGCTGCCTGCGTGTAGCAGATGACCATCGGTAATGACGCCGCCCCCACGTTGTGATCGATAACCACCTGAATCACATCGCGCGCCCCGCGTGAGGCACCAAACAAGGCCTCTGCCATCGTCCATGCGCTGATATCATGCTGAATATAAACCGGAACGCCGGTATGCTGCTCCAGCGCCTCGCCGAGCGGCATCTCTTTTACATCCTCGTAGAACGGCATGCGATGTACAATAcc                                                      <  1:506580/250‑1 (MQ=255)
                                                                                                                                                                                                ccACGAGACTACTGCTGCCTGCGTGTAGCAGATGACCATCGGTAATGACGCCGCCCCCACGTTGTGATCGATAACCACCTGAATCACATCGCGCGCCCCGCGTGAGGCACCAAACAAGGCCTCTGCCATCGTCCATGCGCTGATATCATGCTGAATATAAACCGGAACGCCGGTATGCTGCTCCAGCGCCTCGCCGAGCGGCATCTCTTTTACATCCTCGTAGAACGGCATGCGATGTACAATACCAttt                                                  <  1:598521/250‑1 (MQ=255)
                                                                                                                                                                                                ccACGAGACTACTGCTGCCTGCGTGTAGCAGATGACCATCGGTAATGACGCCGCCCCCACGTTGTGATCGATAACCACCTGAATCACATCGCGCGCCCCGCGTGAGGCACCAAACAAGGCCTCTGCCAGCGTCCATGCGCTGATATCATGCTGAATATAAACCGGAACGCCGGTATGCTGCTCCAGCGCCTCGCCGAGCGGCATCTCTTTTACATCCTCGTAGAACGGCATGCGATGTACAATACCAttt                                                  <  2:102995/250‑1 (MQ=255)
                                                                                                                                                                                                     agaCTACTGCTGCCTGCGTGTAGCAGATAACCATCGGTAATGACGCCGCCCCCACGTTGTGATCGATAACCACCTGAATCACATCGCGCGCCCCGCGTGAGGCACCAAACAAGGCCTCTGCCATCTTCCATGCGCTGATATCATTCTGAATATAATCCGGATCGCCGGTATGCGGCTCCAGCGCCTCGCCGAGCGTCATCTCTTTTACAGCCTCGTAGAACGGCATGCGATGTACAATACCATTTTCCGTa                                            <  2:264400/251‑1 (MQ=255)
                                                                                                                                                                                                          actGCTGCCTGCGTGTAGCAGATGACCATCGGTAATGACGCCGCCCCCACGTTGTGATCGATAACCACCTGAATCACATCGCGCGCCCCGCGTGAGGCACCAAACAAGGCCTCTGCCATCGTCCATGCGCTGATATCATGCTGAATATAAACCGGAACGCCGGTATGCTGCTCCAGCGCCTCGCCGAGCGGCATCTCTTTTACATCCTCGTAGAACGGCATGCGATGTACAATACCATTTTCCGTATCa                                         <  1:620840/249‑1 (MQ=255)
                                                                                                                                                                                                                   tgcGTGTAGCAGATGACCATCGGTAATGACGCCGCCCCCACGTTGTGATCGATAACCACCTGAATCACATCGCGCGCCCCGCGTGAGGCACCAAACAAGGCCTCTGCCATCGTCCATGCGCTGATATCATGCTGAATATAAACCGGAACGCCGGTATGCTGCTCCAGCGCCTCGCCGAGCGGCATCTCTTTTACATCCTCGTAGAACGGCATGCGATGTACAATACCATTTTCCGTATCAATAATTCCCgg                              >  1:46717/1‑251 (MQ=255)
                                                                                                                                                                                                                    gcGTGTAGCAGAAGACTATAGGTAATGACGCCGCCCTCACGTAGTGATCGATAACCACCTGAATCACATCGCGCGCTCCGCGTGAGGCACCAAACAAGGCCTCTGCCATCGTCCATGCGCTGATATCATGCTGAATATAAACCGGAACGCCGGTATGCTGCTCCAGCGCCTCGCCGAGCGGCATCTCTTTTACATCCTCGTAGAACGGCATGCGATGTACAATACCATTTTCCGTATCAATAATTCCCGGc                             <  2:452126/251‑1 (MQ=255)
                                                                                                                                                                                                                              gATGACCATCGGTAATGACGCCGCCCCCACGTTGTGATCGATAACCACCTGAATCACATCGCGCGCCCCGCGTGAGGCACCAAACAAGGCCTCTGCCATCGTCCATGCGCTGATATCATGCTGAATATAAACCGGAACGCCGGTATGCTGCTCCAGCGCCTCGCCGAGCGGCATCTCTTTTACATCCTCGTAGAACGGCATGCGATGTACAATACCATTTTCCGTATCAATAATTCCCGGCAAGGTTATgg                   <  2:559055/251‑1 (MQ=255)
                                                                                                                                                                                                                                   ccATCGGTAATGACGCCGCCCCCACGTTGTGATCGATAACCACCTGAATCACATCGCGCGCCCCGCGTGAGGCACCAAACAAGGCCTCTGCCATCGTCCATGCGCTGATATCATGCTGAATATAAACCGGAACGCCGGTATGCTGCTCCAGCGCCTCGCCGAGCGGCATCTCTTTTACATCCTCGTAGAACGGCATGCGATGTACAATACCATTTTCCGTATCAATAATTCCCGGCAAGGTTATGGCAATc              <  1:54173/251‑1 (MQ=255)
                                                                                                                                                                                                                                   ccATCGGTAATGACGCCGCCCCCACGTTGTGATCGATAACCACCTGAATCACATCGCGCGCCCCGCGTGAGGCACCAAACAAGGCCTCTGCCATCGTCCATGCGCTGATATCATGCTGAATATAAACCGGAACGCCGGTATGCTGCTCCAGCGCCTCGCCGAGCGGCATCTCTTTTACATCCTCGTAGAACGGCATGCGATGTACAATACCATTTTCCGTATCAATAATTCCCGGCAAGGTTAGGGCAATc              <  1:514281/251‑1 (MQ=255)
                                                                                                                                                                                                                                                acgccgcccCCACGTTGTGATCGATAACCACCTGAATCACATCGCGCGCCCCGCGTGAGGCACCAAACAAGGCCTCTGCCATCGTCCATGCGCTGATATCATGCTGAATATAAACCGGAACGCCGGTATGCTGCTCCAGCGCCTCGCCGAGCGGCATCTCTTTTACATCCTCGTAGAACGGCATGCGATGTACAATACCATTTTCCGTATCAATAATTCCCGGCAAGGTTATGGCAATCGAAGTTAGACGc  >  1:466436/4‑251 (MQ=255)
                                                                                                                                                                                                                                                    |                                                                                                                                                                                                                                                      
CGCCGCGCAATGCCGCCTGACACAATGAGTCCACGGTTAACGGTTGTCCATGTAACATCGAGCTCATGGATTGATTAAGACGCAGCTGTGCCAGCTCAAGAATACTGTCCACGCTGGCGATGGTTTCGAGGCAGCCGTGATTCCCGCAATAACAGCGTTTCCCATACGGGTCGACCTGTGTGTGGCCTATTTCCACGAGACTACTGCTGCCTGCGTGTAGCAGATGACCATCGGTAATGACGCCCGCCCCCACGTTGTGATCGATAACCACCTGAATCACATCGCGCGCCCCGCGTGAGGCACCAAACAAGGCCTCTGCCATCGTCCATGCGCTGATATCATGCTGAATATAAACCGGAACGCCGGTATGCTGCTCCAGCGCCTCGCCGAGCGGCATCTCTTTTACATCCTCGTAGAACGGCATGCGATGTACAATACCATTTTCCGTATCAATAATTCCCGGCAAGGTTATGGCAATCGAAGTTAGACGC  >  NC_000913/1667643‑1668133

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