Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,519,818 C→T E296K (GAA→AAA)  hofQ ← DNA utilization protein HofQ

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,519,8180CT98.3% 199.2 / NA 60E296K (GAA→AAA) hofQDNA utilization protein HofQ
Reads supporting (aligned to +/- strand):  ref base C (0/0);  major base T (35/24);  minor base G (0/1);  total (35/25)
Fisher's exact test for biased strand distribution p-value = 4.17e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CAACGGTACGCTATCCTGACCCGATTTATTTTTACGGGTAAAAATGCCGCCCAGCGCCAACGTTTCTCCGCTTTTGACCTCGACCTGCGTTTCGATCTCCTGCTTATCAATCGCCAGCACTTCGCCATCGGCCTGCTGTAGCACCTGCCCCGGAACGTTCTGGCTGATGTGTAATTTCAGCCGGATGCGACCTTTTTGTAACACCGTGGGCGTGACCTCCATCCCCAGGACGGCCTCTTTAAATTCCACCGACGTCGCGCCACTTTCCCCGCTGGAAACCTGATATGGAATTTCGCTCCCCTGTTTAATGCTGGCAGGCTGAAGATGTGAGGCCAGCAGACGCGGACTGGCGATAATATCCAGCTGCTGTTTTTGTTCGAGCGCGGAAAGCTCAAGATCCAGCAAGCGTCCGTTGATGCGCCCAATGTTAAAACCGACATGCGTTGTCGCCGTCGCTACGGAGAGGTCGCTACCAAGCGTG  >  NC_000913/3519573‑3520053
                                                                                                                                                                                                                                                     |                                                                                                                                                                                                                                           
cAACGGTACGCTATCCTGACCCGATTTATTTTTACGGGTAAAAATGCCGACCAGCGCCAACGTTTCTCCGCTTTTGACCTCGACCTGCGTTTCGATCTCCTGCTTATCAATCGCCAGCACTTCGCCATCGGCCTGCTGTAGCACCTGCCCCGGAACGTTCTGGCTGATGTGTAATTTCAGCCGGATGCGACCTTTTTGTAACACCGTGGGCGTGACCTCCATCCCCAGGACGGCCTCTTTAAATTTCACCg                                                                                                                                                                                                                                        >  1:372055/1‑251 (MQ=255)
   cGGTACGCTATCCTGACCCGATTTATTTTTACGGGTAAAAATGCCGCCCAGCGCCAACGTTTCTCCGCCTTTGACCTCGACCTGCGTTTCGATCTCCTGCTTATCCATCGCCAGCACTTCGCCATCGGCCTGCTGTAGCACCTGCCCCGGAACGTTCTGGCTGATGTGTAATTTCAGCCGGATGCGACCTTTTTGTAACACCGTGGGCGTGACCTCCATCCCCAGGACGGCCTCTTTACATTTCACCGACg                                                                                                                                                                                                                                     >  1:32732/1‑251 (MQ=255)
            aTCCTGACCCGATTTATTTTTACGGGTAAAAATGCCGCCCAGCGCCAACGTTTCTCCGCTTTTGACCTCGACCTGCGTTTCGATCTCCTGCTTATCAATCGCCAGCACTTCGCCATCGGCCTGCTGTAGCACCTGCCCCGGAACGTTCTGGCTGATGTGTAATTTCAGCCGGATGCGACCTTTTTGTAACACCGTGGGCGTGACCTCCATCCCCAGGACGGCCTCTTTAAATTTCACCGACGTCGCGCCAc                                                                                                                                                                                                                            <  1:37875/251‑1 (MQ=255)
              ccTGACCCGATTTATTTTTACGGGTAAAAATGCCGCCCAGCGCCAACGTTTCTCCGCTTTTGACCTCGACCTGCGTTTCGATCTCCTGCTTATCAATCGCCAGCACTTCGCCATCGGCCTGCTGTAGCACCTGCCCCGGAACGTTCTGGCTGATGTGTAATTTCAGCCGGATGCGACCTTTTTGTAACACCGTGGGCGTGACCTCCATCCCCAGGACGGCCTCTTTAAATTTCACCGACGTCGCGCCACtt                                                                                                                                                                                                                          <  1:132078/251‑1 (MQ=255)
               cTGACCCGATTTATTTTTACGGGTAAAAATGCCGCCCAGCGCCAACGTTTCTCCGCTTTTGACCTCGACCTGCGTTTCGATCTCCTGCTTATCAATCGCCAGCACTTCGCCATCGGCCTGCTGTAGCACCTGCCCCGGAACGTTCTGGCTGATGTGTAATTTCAGCCGGATGCGACCTTTTTGTAACACCGTGGGCGTGACCTCCATCCCCAGGACGGCCTCTTTAAATTTCACCGACGTCGCGCCACttt                                                                                                                                                                                                                         >  1:285477/1‑251 (MQ=255)
                  aCCCGATTTATTTTTACGGGTAAAAATGCCGCCCAGCGCCACCGTTTCTCCGCTTTTGACCTCGAACTGCGTTTCGATCTCCTGCTTATCAATCGCCAGCACTTCGCCATCGGCCTGCTGTAGCACCTGCCCCGGCACGTTCTCGCTGTTGTGTACTTTCAGCCGGATGCGACCTTTTTGTAACACCGTGGGCGTGACCTCCATCCCCAGGCCGGCCTCTTTACATTTCACCGCCGTCGCGCCCCTTTccc                                                                                                                                                                                                                      >  2:22724/1‑251 (MQ=255)
                  aCCCGATTTATTTTTACGGGTAAAAATGCCGCCCAGCGCCAACGTTTCTCCGCTTTTGACCTCGACCTGCGTTTCGATCTCCTGCTTATCAATCGCCAGCACTTCGCCACCGGCCTGCTGTAGCACCTGCCCCGGAACGTTCTGGCTGATGTGTAATTTCAGCCGGATGCGACCTTTTTGTAACACCGTGGGCGTGACCTCCATCCCCAGGACGGCCTCTTTAAATTTCACCGACGTCGCGCCACTTTccc                                                                                                                                                                                                                      >  1:245070/1‑251 (MQ=255)
                                 aCGGGTAAAAATGCCGCCCAGCGCCAACGTTTCTCCGCTTTTGACCTCGACCTGCGTTTCGATCTCCTGCTTATCAATCGCCAGCACTTCGCCATCGGCCTGCTGTAGCACCTGCCCCGGAACGTTCTGGCTGATGTGTAATTTCAGCCGGATGCGACCTTTTTGTAACACCGTGGGCGTGACCTCCATCCCCAGGACGGCCTCTTTAAATTTCACCGACGTCGCGCCACTTTCCCCGCTGGAAACCTGat                                                                                                                                                                                                       <  1:378646/251‑1 (MQ=255)
                                  cGGGTAAAAATGCCGCCCAGCGCCAACGTTTCTCCGCTTTTGACCTCGACCTGCGTTTCGATCTCCTGCTTATCAATCGCCAGCACTTCGCCATCGGCCTGCTGTAGCACCTGCCCCGGAACGTTCTGGCTGATGTGTAATTTCAGCCGGATGCGACCTTTTTGTAACACCGTGGGCGTGACCTCCATCCCCAGGACGGCCTCTTTAAATTTCACCGACGTCGCGCCACTTTCCCCGCTGGAAACCTGata                                                                                                                                                                                                      <  1:395625/251‑1 (MQ=255)
                                           aTGCCGCCCAGCGCCAACGTTTCTCCGCTTTTGACCTCGACCTGCGTTTCGATCTCCTGCTTATCAATCGCCAGCACTTCGCCATCGGCCTGCTGTAGCACCTGCCCCGGAACGTTCTGGCTGATGTGTAATTTCAGCCGGATGCGACCTTTTTGTAACACCGTGGGCGTGACCTCCATCCCCAGGACGGCCTCTTTAAATTTCACCGACGTCGCGCCACTTTCCCCGCTGGAAACCTGATATGGAATTTc                                                                                                                                                                                             >  1:70680/1‑251 (MQ=255)
                                                  ccAGCGCCAACGTTTCTCCGCTTTTGACCTCGACCTTCGTTTCGATCTCCTGCTTATCAATCGCCAGCACTTCGCCATCGGCCTGCTGTAGCACCTGCCCCGGAACGCTCTGGCTAATGTGTAATTTCAGCCGGATGCGACCTTTTTATAACACCGTGGGCGTGACGTCCATCCCCTGGACGGCCTCTTTAACTTTCACCGACGTCGCGCCACTTTCCCCGCTGGAAACCTGATATGGACTTTCGCTcccc                                                                                                                                                                                      >  1:221501/1‑251 (MQ=255)
                                                  ccAGCGCCAACGTTTCTCCGCTTTTGACCTCGACCTGCGTTTCGATCTCCTGCTTATCAATCGCCAGCACTTCGCCATCGGCCTGCTGTAGCACCTGCCCCGGAACGTTCTGGCTGATGTGTAATTTCAGCCGGATGCGACCTTTTTGTAACACCGTGGGCGTGACCTCCATCCCCAGGACGGCCTCTTTAAATTTCACCGACGTCGCGCCACTTTCCCCGCTGGAAACCTGATATGGAATTTCGCTcccc                                                                                                                                                                                      >  1:191657/1‑251 (MQ=255)
                                                  ccAGCGCCAACGTTTCTCCGCTTTTGACCTCGACCTGCGTTTCGATCTCCTGCTTATCAATCGCCAGCACTTCGCCATCGGCCTGCTGTAGCACCTGCCCCGGAACGTGCTGGCTGATGTGTAATTTCAGCCGGATGCGACCTTGTTGTAACACCGTGGGCGTGACCTCCATCCCCAGGACTGCCTCTTTAAATTTCACCGACGTCGCGCCACTTTCCCCGCTGGAAACCTGATATGGAATTTCGCTcccc                                                                                                                                                                                      <  1:273006/251‑1 (MQ=255)
                                                       gcCAACGTTTCTCCGCTTTTGACCTCGACCTGCGTTTAGATCTCCTGCTTATCAATCGCCAGCACTTCGCCATCGGCCTGCTGTAGCACCTGCCCCGGAACGTTCTGGCTGATGTGTAATTTCAGCCGGATGCGACCTTTTTGTAACACCGTGGGCGTGACCTCCATCCCCAGGACGGCCTCTTTAAATTTCACCGACGTCGCGCCACTTTCCCCGCTGGAAACCTGATATGGACTTTCGCTCCCCTGttt                                                                                                                                                                                 >  2:364266/1‑251 (MQ=255)
                                                          aaCGTTTCTCCGCTTTTGACCTCGACCTGCGTTTCGATCTCCTGCTTATCAATCGCCAGCACTTCGCCATCGGCCTGCTGTAGCACCTGCCCCGGAACGTTCTGGCTGATGTGTAATTTCAGCCGGATGCGACCTTTTTGTAACACCGTGGGCGTGACCTCCATCCCCAGGACGGCCTCTTTAAATTTCACCGACGTCGCGCCACTTTCCCCGCTGGCAACCTGATATGGAATTTCGCTCCCCTGTTTAAt                                                                                                                                                                              >  1:226061/1‑251 (MQ=255)
                                                           aCGTTTCTCCGCTTTTGACCTCGACCTGCGTTTCGATCTCCTGCTTATCAATCGCCAGCACTTCGCCATCGGCCTGCTGTAGCACCTGCCCCGGAACGTTCTGGCTGATGTGTAATTTCAGCCGGATGCGACCTTTTTGTAACACCGTGGGCGTGACCTCCATCCCCAGGACGGCCTCTTTAAATTTCACCGACGTCGCGCCACTTTCCCCGCTGGAAACCTGATATGGACTTTCGCTCCCCTGTTTAATg                                                                                                                                                                             >  2:143397/1‑251 (MQ=255)
                                                                        tttGACCTCGACCTGCGTTTCGATCTCCTGCTTATCAATCGCCAGCACTTCGCCATCGGCCTGCTGTAGAACCTGCCCCGGAACGTTCTGGCTGATGTGTAATTTCAGCCGGATGCGACCTTTTTGTAACACCGTGGGCGTGACCTCCATCCCCAGGACGGCCTCTTTAAATTTCACCGACGTCGCGCCACTTTCCCCGCTGGAAACCTGATATGGAATTTCGCTCCCCTGTTTAATGCTGGCAGGCTGa                                                                                                                                                                 >  1:374143/1‑250 (MQ=255)
                                                                         ttGACCTGGACCTGCGTTTCGATCTCCTGCTTATCAATCGCCAGCACTTCGCCATCGGCCTGCTGTAGCACCTGCCCCTGAACGTTCTGGCTGATGTGTAATTCCAGCCGGATGCGACCTTTTTGTAACACCGTTGGCGTGACCTCCATCCCCAGGACGGCCACTTTAAATTTCACCGCCTTCGCGCCACTTTCCCCGCTGGAAACCTGAGATGGAATTTCGCTCCCCTGTTTAATGCTTGCAGGCTGaa                                                                                                                                                                <  2:160329/250‑1 (MQ=255)
                                                                                cGACCTGCGTTTCGATCTCCTGCTTATCAATCGCCAGCACTTCGCCATCGGCCTGCTGTAGCACCTGCCCCGGAACGTTCTGGCTGATGTGTAATTTCAGCCGGATGCGACCTTTTTGTAACACCGTGGGCGTGACCTCCATCCCCAGGACGGCCTCTTTAAATTTCACCGACGTCGCGCCACTTTCCCCGCTGGAAACCTGATATGGAATTTCGCTCCCCTGTTTAATGCTGGCAGGCTGAAGATGTGAg                                                                                                                                                        <  1:321625/251‑1 (MQ=255)
                                                                                 gACCTGCGTTTCGATCTCCTGCTTATCAATCGCCAGCACTTCGCCATCGGCCTGCTGTAGCACCTGCCCCGGAACGTTCTGGCTGATGTGTAATTTCAGCCGGATGCGACCTTTTTGTAACACCGTGGGCGTGACCTCCATCCCCAGGACGGCCTCTTTAAATTTCACCGACGTCGCGCCACTTTCCCCGCTGGAAACCTGATATGGAATTTCGCTCCCCTGTTTAATGCTGGCAGGCTGAAGATGTGAgg                                                                                                                                                       >  1:271881/1‑251 (MQ=255)
                                                                                  aCCTGCGTTTCGATCTCCTGCTTATCAATCGCCAGCACTTCGCCATCGGCCTGCTGTAGCACCTGCCCCGGACCGTTCTGGCTGATGTGTAATTTCAGCCGGATGCGACCTTTTTGTAACACCGTGGGCGTGACCTCCATCCCCAGGACGGCCTCTTTAAATTTCACCGACGTCGCGCCACTTTCCCCGCTGGAAACCTGATATGGAATTTCGCTCCCCTGTTTAATGCTGGCCGGCTGAAGATGTGAGGc                                                                                                                                                      >  1:386598/1‑251 (MQ=255)
                                                                                   ccTGCGTTTCGATCTCCTGCTTATCAATCGCCAGCACTTCGCCATCGGCCTGCTGTAGCACCTGCCCCGGAACGTTCTGGCTGATGTGTAATTTCAGCCGGATGCGACCTTTTTGTAACACCGTGGGCGTGACCTCCATCCCCAGGACGGCCTCTTTAAATTTCACCGACGTCGCGCCACTTTcccc                                                                                                                                                                                                                     >  1:214332/1‑187 (MQ=255)
                                                                                   ccTGCGTTTCGATCTCCTGCTTATCAATCGCCAGCACTTCGCCATCGGCCTGCTGTAGCACCTGCCCCGGAACGTTCTGGCTGATGTGTAATTTCAGCCGGATGCGACCTTTTTGTAACACCGTGGGCGTGACCTCCATCCCCAGGACGGCCTCTTTAAATTTCACCGACGTCGCGCCACTTTcccc                                                                                                                                                                                                                     <  2:214332/187‑1 (MQ=255)
                                                                                    cggCGTTTCGATCTACTGCTTATCACTCGCCAGCACTTCGCCATCGGCCTGCTGTAGCACCTGCCCCGGAACGTTCTGGCTGATGTGTAATTTCAGCCGGATGCGACCTTTTTGTAACACCGTGGGCGTGACCTCCATCCCCAGGACGGCCTCTTTAAATTTCACCGACGTCGCGCCACTTTCCCCGCTGGAAACCTGATATGGAATTTCGCTCCCCTGTTTAATGCTGGCAGGCTGAAGATGTGAGGCca                                                                                                                                                    <  2:2640/249‑1 (MQ=255)
                                                                                    cTGCGTTTCGATCTCCTGCTTATCAATCGCCAGCACTTCGCCATCGGCCTGCTGTAGCACCTGCCCCGGAACGTTCTGGCTGATGTGTAATTTCAGCCGGATGCGACCTTTTTGTAACACCGTGGGCGTGACCTCCATCCCCAGGACGGCCTCTTTAAATTTCACCGACGTCGCGCCACTTTCCCCGCTGGAAACCTGATATGGAATTTCGCTCCCCTGTTTAATGCTGGCAGGCTGAAGATGTGAGGcc                                                                                                                                                     >  1:297005/1‑250 (MQ=255)
                                                                                    cTGCGTTTCGATCTCCTGCTTATCAATCGCCAGCACTTCGCCATCGGCCTGCTGTAGCACCTGCCCCGGAACGTTCTGGCTGATGTGTAATTTCAGCCGGATGCGACCTTTTTGTAACACCGTGGGCGTGACCTCCATCCCCAGGACGGCCTCTTTAAATTTCACCGACGTCGCGCCACTTTCCCCGCTGGAAACCTGATATGGAATTTCGCTCCCCTGTTTAATGCTGGCAGGCTGAAGATGTGAGGcc                                                                                                                                                     >  1:95193/1‑250 (MQ=255)
                                                                                                  ccTGCTTATCAATCGCCAGCACTTCTCCATCGGCCTGCTGTAGCACCTGCCCCGGAACGTTCTGGCTGATGTGTAATTTCAGCCGGATGCGACCTTTTTGTAACACCGTGGGCGTGACCTCCATCCCCAGGACGGCCTCTTTAAATTTCACCGACGGCGCGCCACTGTCGCCTCTGGAAACCTGAGATGGAAGTTCGCTCCCCTGGTTAATGCTGGCAGGCTGAAGATGGGAGGCCAGCATAcgcg                                                                                                                                           <  2:404540/246‑1 (MQ=255)
                                                                                                   cggCTTAGCAAGCGCCAGCACTGCGCCATCGGCCGGTTGTAGTACCTGCCCCGGAACGTTCGGGCTGAGGGTTAATTTCAGCCGGATGCGACCTTTTTGTAACACCGTGGGCGTGACCACCATCCCCAGGACGGCCTCTTTAAATGTCACCGACGTCGCGCCATTTTCCCTGCTGGAAATCTGATATGGAAATTCGCTCCCCGGTTTAAGGCTGGCAGGCTGAAGATGTGAGGCCAGCAGACGCggccggg                                                                                                                                     <  2:401833/249‑6 (MQ=255)
                                                                                                    tGCTTCTCAACCGCCAGCACTTCGCCATCGGCCTGCTGTAGCACCTGCCCCGGACCGTTCTGGCTGATGTGTAATTTCAGCCGGATGCGACCTTTTTGTAACACCGTGGGCGTGACCTCCATCCCCAGGACGGCCTCTTTAAATTTCACCGACGTCGCGCCACTTTCCCCGCTGGAAACCTGATATGGAATTTCGCTCCCCTGTTTAATGCTGGCAGGCTGAAGATGTGAGGCCAGCAGACGCGGACTGGc                                                                                                                                    >  2:397876/1‑251 (MQ=255)
                                                                                                    tGCTTATCAATCGCCAGCACTTCGCCATCGGCCTGCTGTAGCACCTGCCCCGGAACGTTCTGGCTGATGTGTAATTTCAGCCGGATGCGCCCTTTTTGTAACACCGTGGGCGTGACCTCCATCCACAGCACGGCCTCTTTAAATTTCACCGCCGTCGCGACACTTTCCCCGCTGGACACATGATATGGAATTTCGCTCCCCTGTTTACTGCTGGCAGGCTGCAGATGAGAGGCCAGCAGACGCGGCCTGGc                                                                                                                                    >  1:417997/1‑251 (MQ=255)
                                                                                                     gCTTATCAATCGCCAGCACTGCGCCATCGGCCTGCTGTAGCACCTGCCCCGGAACGTTCTGGCTGATGTGTAATGTCAGCCGGATGCGACCTTTTTGTAACACCGTGGGCGTGACCTCCAGCCCCAGGACGGCCTCTTTAAATTTCACCGACGTCGCGCCACTTTCCCCGCTGGAAACCTGAGATGGAATTTCGCTCCCCTGGTTAATGCTGGCAGGCTGAAGATGTGAGGCCAGCAGACGCGGACTGGc                                                                                                                                    <  1:380912/250‑1 (MQ=255)
                                                                                                        tagcAATCTCCAGCGCTTCGCCATCGGCCTGCTGTAGCACCTGCCCCGGAACGTTCTGGCTGATGTGTAATTTCAGCCGTATGCGACCTTTTTGTAACACCGTGGGCGTGACCTCCATCCCCAGGACGGCCTCTTTAAATTTCACCGACGTCGCGCCACTTTCCCCGCGGGAAACCTGATATGGAATTGCGCTCCCCTGTTTAATGCTTGCAGGCTGAAGATGTGAGGCCAGCAGACGCGTACTGGCGa                                                                                                                                  <  2:197382/246‑1 (MQ=255)
                                                                                                         aTCAATCGCCAGCACTTCGCCAGCGGCCTGCTGTAGCACCTGCCCCGGAACGTTCTGGCTGAGGTGTAATTTCAGCCGGATGCGACCTTTTTGTAACACCGTGGGCGTGACCTCCATCCCCAGGACGGCCTCTTTAAATTGCACCGACGTCGCGCCACTTTCCCCGCTGGAAACCTGATATGGAATTTCGCTCCCCTGTTTAATGCTGGCAGGCTGAAGATGTGAGTCCAGCAGACGCGGACTGGCGataa                                                                                                                               <  2:285477/251‑1 (MQ=255)
                                                                                                            aaTCGCCAGCACTTCGCCATCGGCCTGCTGTAGCACCTGCCCCGGAACGTTCTGGCTGATGTGTAATTTCAGCCGGATGCGACCTTTTTGTAACACCGTGGGCGTGACCTCCATCCCCAGGACGGCCTCTTTAAATTTCACCGACGTCGCGCCACTTTCCCCGCTGGAAACCTGATATGGAATTTCGCTCCCCTGTTTAATGCTGGCAGGCTGAAGATGTGAGGCCAGCAGACGCGGACTGGCGATAatat                                                                                                                            >  1:326289/1‑251 (MQ=255)
                                                                                                               cGCCAGCACTTCGCCATCGGCCTGCTGTAGCACCTGCCCCGGAACGTTCTGGCTGATGTGTAATTTCAGCCGGATGCGACCTTTTTGTAACACCGTGGGCGTGACCTCCATCCCCAGGACGGCCTCTTTAAATTTCACCGACGTCGCGCCACTTTCCCCGCTGGAAACCTGATATGGAATTTCGCTCCCCTGTTTAATGCTGGCAGGCTGAAGATGTGAGGCCAGCAGACGCGGACTGGCGATAATATcc                                                                                                                          >  2:208519/1‑250 (MQ=255)
                                                                                                                  cAGCACTTCGCCATCGGCCTGCTGTAGCACCTGCACCGGAACGTTCCGGCTGATGTGTAATTTCAGCCGGATGCGACCTTTTTGTAACACCGTGGGCGTGACCTCCATCCCCAGGACGGCCTCTTTAAATTTCACCGACGTCGCGCCACTTTCCCCGCTGGAAACCTGATATGGAATTTCGCTCCCCTGTTTAATGCTGGCAGGCTGAAGATGTGAGGCCAGCAGACGCGGACTGGCGATAATATACAgct                                                                                                                      >  1:132899/1‑251 (MQ=255)
                                                                                                                         gcgccaatggCCTGCTTTAGCACCTGCCCCGGAACTTTCTGGCTGATGTGTAATTTCAGCCGGATGCGACCTTTTTGTCACACCGTGGGCGTTACCTCCATCCCCAGGACGGCCTCTTTAAATTTCACCGACGTCGCGCCACTTTCCCCGTTGGAAACCTGATATGGCATTTCGCTCCCCTGTTTAATGCTGGCAGGCTGAAGTTGTGAGGACAGCAGACGCGGACTGGCGATAATCTCCAGCTGCTGttg                                                                                                               <  2:226061/243‑2 (MQ=255)
                                                                                                                          cGCCATCGGCCTGCTGTAGCACCTGCCCCGGAACGTTCTGGCTGATGTGTAATTTCAGCCGGATGCGACCTTTTTGTAACACCGTGGGCGTGACCTCCATCCCCAGGACGGCCTCTTTAAATTTCACCGACGTCGCGCCACTTTCCCCGCTGGAAACCTGATATGGAATTTCGCTCCCCTGTTTAATGCTGGCAGGCTGAAGATGTGAGGCCAGCAGACGCGGACTGGCGATAATATCCAGCTGCTGttt                                                                                                               <  1:359774/250‑1 (MQ=255)
                                                                                                                           gCCATCGGCCTGCTGTAGCACCTGCCCCGGAACGTTCTGGCTGATGTGTAATTTCAGCCGGATGCGACCTTTTTGTAACACCGTGGGCGTGACCTCCATCCCCAGGACGGCCTCTTTAAATTTCACCGACGTCGCGCCACTTTCCCCGCTGGAAACCTGATATGGAATTTCGCTCCCCTGTTTAATGCTGGCAGGCTGAAGATGTGAGGCCAGCAGACGCGGACTGGCGATAATATCCAGCTGCTGttttt                                                                                                             <  1:339672/251‑1 (MQ=255)
                                                                                                                               ccGGCCTGCTGTAGCACCTGCCCCGGAACGTTCTGGCTGATGTGTAATTTCAGCCGGATGCGACCTTTTTGTAACACCGTGGGCGTGACCTCCATCCCCAGGACGGCCTCTTTAAATTTCACCGACGTCGCGCCACTTTCCCCGCTGGAAACCTGATATGGAATTTCGCTCCCCTGTTTAATGCTGGCAGGCTGAAGATGTGAGGCCAGCAGACGCGGACTGGCGATAATATCCAGCTGCTGTTTTTGTTc                                                                                                         <  2:154263/250‑1 (MQ=255)
                                                                                                                                 ggCCTGCTGTAGCACCTGCCCCGGAACGTTCTGGCTGATGTGTAATTTCAGCCGGATGCGACCTTTTTGTAACACCGTGGGCGTGACCTCCATCCCCAGGACGGCCTCTTTAAATTTCACCGACGTCGCGCCACTTTCCCCGCTGGCAACCTGATATGGAATTTCGCTCCCCTGTTTAATGCTGGCAGGCTGAAGATGTGAGGCCAGCAGACGCGGACTGGCGATAATATCCAGCTGCTGTTTTTGTTCg                                                                                                        >  2:3692/1‑250 (MQ=255)
                                                                                                                                     tgctgTAGCACCTGCCCCGGAACGTTCTGGCTGATGTGTAATTTCAGCCGGATGCGACCTTTTTGTAACACCGTGGGCGTGACCTCCATCCCCAGGACGGCCTCTTTAAATTTCACCGACGTCGCGTCACTTTCCCCGCTGGAAACCTGATATGGAATTTCGCTCCCCTGTTTAATGCTGGCAGGCTGAAGATGTGAGGCCAGCAGACGCGGACTGGCGATAATATCCAGCTGCTGTTTTTGTTCGAgcgc                                                                                                   >  2:360705/1‑251 (MQ=255)
                                                                                                                                                  gCCCCGGAACGTTCTGGCTGATGTGTAATTTCAGCCGGATGCGACCTTTTTGTAACACCGTGGGCGTGACCTCCATCCCCAGGACGGCCTCTTTAAATTTCACCGACGTCGCGCCACTTTCCCCGCTGGAAACCTGATATGGAATTTCGCTCCCCTGTTTAATGCTGGCAGGCTGAAGATGTGAGGCCAGCAGACGCGGACTGGCGATAATATCCAGCTGCTGTTTTTGTTCGAGCGCGGAAAGCTCAAg                                                                                       <  2:95193/250‑1 (MQ=255)
                                                                                                                                                       ggAACGTTCTGGCTGATGTGTACTTTCAGCCGGATGCGACCTTTTTGTAACACCGTGGGCGTGACCTCCATCCCCAGGACGGCCTCTTTAAATTTCACCGACGTCGCGCCACTTTCCCCGCTGGAAACCTGATATGGAATTTCGCTCCCCTGTTTAATGCTGGCAGGCTGAAGATGTGAGGCCAGCAGACGCGGACTGGCGATAATATCCAGCTGCTGTTTTTGTTCGAgc                                                                                                     >  2:374003/1‑231 (MQ=255)
                                                                                                                                                       ggAACGTTCTGGCTGATGTGTAATTTCAGCCGGATGCGACCTTTTTGTAACACCGTGGGCGTGACCTCCATCCCCAGGACGGCCTCTTTAAATTTCACCGACGTCGCGCCACTTTCCCCGCTGGAAACCTGATATGGAATTTCGCTCCCCTGTTTAATGCTGGCAGGCTGAAGATGTGAGGCCAGCAGACGCGGACTGGCGATAATATCCAGCTGCTGTTTTTGTTCGAgc                                                                                                     <  1:374003/231‑1 (MQ=255)
                                                                                                                                                            ggtCTGGCTGATGTGTAATTTCAGCCGGATGCGACCTTTTTGTAACACCGTGGGCGTGACCTCCATCCCCAGGACGGCCTCTTTAAATTTCACCGACGTCGCGCCACTTTCCCCGCTGGAAACCTGATATGGAATTTCGCTCCCCTGTTTAATGCTGGCAGGCTGAAGATGTGAGGCCAGCAGACGCGGACTGGCGATAATATCCAGCTGCTGTTTTTGTTCGAGCGCGGAAAGCTCAAGATCCAGCAAGc                                                                            <  2:316079/249‑1 (MQ=255)
                                                                                                                                                                        gtgtAATTTCAGCCGGATGCGACCTTTTTGTAACACCGTGGGCGTGACCTCCATCCCCAGGACGGCCTCTTTAAATTTCACCGACGTCGCGCCACTTTCCCCGCTGGAAACCTGATATGGAATTTCGCTCCCCTGTTTAATGCTGGCAGGCTGAAGATGTGAGGCCAGCAGACGCGGACTGGCGATAATATCCAGCTGCTGTTTTTGTTCGAGCGCGGAAAGCTCAAGATCCAGCAAGCGTCCGTTGATgc                                                                >  2:346185/1‑251 (MQ=255)
                                                                                                                                                                           tAATTTCAGCCGGATGCGACCTTTTTGTAACACCGTGGGCGTGACCTCCATCCCCAGGACGGCCTCTTTAAATTTCACCGACGTCGCGCCACTTTCCCCGCTGGAAACCTGATATGGAATTTCGCTCCCCTGTTTAATGCTGGCAGGCTGAAGATGTGAGGCCAGCAGACGCGGACTGGCGATAATATCCAGCTGCTGTTTTTGTTCGAGCGCGGAAAGCTCAAGATCCAGCAAGCGTCCGTTGATGCGcc                                                             >  2:120174/1‑251 (MQ=255)
                                                                                                                                                                                           cGACCTTTTTTTAACACCGTGGGCGTGACCTCCATCCCCAGGACGGCCTCTTTAAATTTCACCTACGTCGCGCCACTTTCCCCGCTGGAATCCTGATATGGAATTTCGCTCCCCTGTTTAATGCTGGCAGGCTGAAGATGTGAGGCCAGCAGACGCGGACTGGCGATAATATCCAGCTGCTGTTTTGGTTCGAGCGCGGAAAGCTCAAGATCCATCAAGCGTCCGTTGATGCGCCCAACTTTAAAACCGAc                                             <  2:306032/251‑1 (MQ=255)
                                                                                                                                                                                              ccTTTTTGTAACACCGTGGGCGTGACCTCCATCCCCAGGACGGCCTCTTTAAATTTCACCGACGTCGCGCCACTTTCCCCGCTGGAAACCTGATATGGAATTTCGCTCCCCTGTTTAATGCTGGCAGGCTGAAGATGTGAGGCCAGCAGACGCGGACTGGCGATAATATCCAGCTGCTGTTTTTGTTCGAGCGCGGAAAGCTCAAGATCCAGCAAGCGTCCGTTGATGCGCCCAATGTTAAAACCGACATg                                          >  1:389176/1‑251 (MQ=255)
                                                                                                                                                                                               cTTTTTGTAACACCGTGGGCGTGACCTCCATCCCCAGGACGGCCTCTTTCAATTTCACCGCCGTCGCGCCACTTTCCCCGCTGGCAACCTGATATGGACTTTCGCTCCCCTGTTTAATGCTGGCAGGCTGAAGATGTGAAGCCAGCAGACGCGGACTGGCGATAATATCCAGCTACTGTTTTTGTTCGAGCGCGGAAAGCTCAAGATCCAGCACGCGTCCGTTGATGCGCCCAATGTTAAAACCGCCATGc                                         >  1:401256/1‑251 (MQ=255)
                                                                                                                                                                                                        acacCTTGGGCGTGACCTCCATCCCCAGGACGGCCGCTTTAAATTTCTCCGACGTCGCGCCACTTTCCCCGCTGGAAACCGTATATGGAATTTCGCTCCCCTGTTTAATGCTGGCAGGCTGAAGATGTGAGGCCAGCAGACGCGGACTGGCGATAATATCCAGCTGCTGTTTTTGTTCGAGCGCTGAAGTGTCAAGATCCAGCAAGCGTCCGTTGATGCGCGCAATGTTAAAACCGACATGCGTTGTCGcc                                <  2:197201/251‑1 (MQ=255)
                                                                                                                                                                                                            cGTGGGCGTGACCTCCATCCCCAGGACGGCCTCTTTAAATTTCACCGACGTCGCGCCACTTTCCCCGCTGGAAACCTGATATGGAATTTCGCTCACCTGTTTAATGCTGGCAGGCTGACGATGTGAGGCCAGCAGACGCGGACTGGCGATGATATCCAGCTGCTGTTTTTGTTCGAGCGCGGAAAGCTCAAGATCCAGCAAGCGTCCGTTGATGCGCCCAATGTTAAAACCGACATGCGTTGTCGCCGTCg                            >  1:254129/1‑251 (MQ=255)
                                                                                                                                                                                                             gTGGGCGTGACCTCCATCCCCAGGACGGCCTCTTTAAATTTCACCGACGTCGCGCCACTTTCCCCGCTGGAAACCTGATATGGAATTTCGCTCCCCTGTTTAATGCTGGCAGGCTGAAGATGTGAGGCCAGCAGACGCGGACTGGCGATAATATCCAGCTGCTGTTTTTGTTCGAGCGCGGAAAGCTCAAGATCCAGCAAGCGTCCGTTGATGCGCCCAATGTTAAAACCGACATGCGTTGTCGCCGTCGc                           >  1:146702/1‑251 (MQ=255)
                                                                                                                                                                                                                  cGTGACCTCCCTCCCCAGGACGGCCTCTTTAAATTTCACCGCCGTCGCGCCACTTTCCCCGCTGGAAACCTGATATGGAATTTCGCTCCCCTGTTTAATGCTGGCAGGCTGAAGATGTGAGGCCAGCAGACGCGGCCTGGCGATACTATCCAGCTGCTGTTTTTGTTCGAGCGCGGAAAGCTcca                                                                                        >  2:164985/1‑183 (MQ=255)
                                                                                                                                                                                                                  cGTGACCTCCATCCCCAGGACGGCCTCTTTAAATTTCACCGACGTCGCGCCACTTTCCCCGCTGGAAACCTGATATGGAATTTCGCTCCCCTGTTTAATGCTGGCAGGCTGAAGATGTGAGGCCAGCAGACGCGGACTGGCGATAATATCCAGCTGCTGTTTTTGTTCGAGCGCGGAAAGCTCaa                                                                                        <  1:164985/185‑1 (MQ=255)
                                                                                                                                                                                                                   gTGACCTCCATCCCCAGAACGGCCTCTTTAAATTTCACCGACGTCGCGCCACTTTCCCCGCTGGAAACCTGATATGGAATTTCGCTCCCCTGTTTAATGCTGGCAGGCTGAAGATGTGAGGCCAGCAGACGCGGACTGGCGATAATATCCAGCTGCTGTTTTTGTTCGAGCGCGGAAAGCTCAAGATCCAGCAAGCGTCCGTTTATGCGCCCAATGTTAAAACCGACATGCGTTGTCGCCGTCGCTACgg                      >  2:229596/1‑250 (MQ=255)
                                                                                                                                                                                                                     gACCTCCATCACAAGTACGAACTCTTTAAATTTCACAGACGTCACACCCCTTTACCCGCTGAAACCCTGATATGGAATTTCGCTCCCCTGTTTAATGCTGGAAGGCTGAAGATTTGAGGCCAGCACACGCGCACTGGCTATAATATCCAGCTGCTGTTTTTGTTCGAGCGCGGAACGCTCAAGCTCCAGCACACGTCCTTTGATCTGCCAAATGTTAACACCGACATTCGTTGTCGCCGTCGCTACgg                      >  2:32814/1‑248 (MQ=255)
                                                                                                                                                                                                                                     aCGGCCTCTTTAAATTTCACCGACGTCGCGTCACTTTCCCCGCTGGAAACCTGATATGGAATTTCGCTCCCCTGTTTAATGCTGGCAGGCTGAAGATGTGAGGCCAGCAGACGCGGACTGGCGATAATATCCAGCTGCTGTTTTTGTTCGAGCGCGGAAAGCTCAAGATCCAGCAAGCGTCCGTTGATGCGCCCAATGTTAAAACCGACATGCGTTGTCGCCGTCGCTACGGAGAGGTCGCTACCAAGCgt   <  1:360705/251‑1 (MQ=255)
                                                                                                                                                                                                                                      cGGCCTCTTTAAATTTCACCGACGTCGCGCCACTTTCCCCGCTGGAAACCTGATATGGAATTTCGCTCCCCTGTTTAATGCTGGCAGGCTGAAGATGTGAGGCCAGCAGACGCGGACTGGCGATAATATCCAGCTGCTGTTTTTGTTCGAGCGCGGAAAGCTCAAGATCCAGCAAGCGTCCGTTGATGCGCCCAATGTTAAAACCGACATGCGTTGTCGCCGTCGCTACGGAGAGGTCGCTACCAAGCgtg  >  1:418806/1‑251 (MQ=255)
                                                                                                                                                                                                                                                     |                                                                                                                                                                                                                                           
CAACGGTACGCTATCCTGACCCGATTTATTTTTACGGGTAAAAATGCCGCCCAGCGCCAACGTTTCTCCGCTTTTGACCTCGACCTGCGTTTCGATCTCCTGCTTATCAATCGCCAGCACTTCGCCATCGGCCTGCTGTAGCACCTGCCCCGGAACGTTCTGGCTGATGTGTAATTTCAGCCGGATGCGACCTTTTTGTAACACCGTGGGCGTGACCTCCATCCCCAGGACGGCCTCTTTAAATTCCACCGACGTCGCGCCACTTTCCCCGCTGGAAACCTGATATGGAATTTCGCTCCCCTGTTTAATGCTGGCAGGCTGAAGATGTGAGGCCAGCAGACGCGGACTGGCGATAATATCCAGCTGCTGTTTTTGTTCGAGCGCGGAAAGCTCAAGATCCAGCAAGCGTCCGTTGATGCGCCCAATGTTAAAACCGACATGCGTTGTCGCCGTCGCTACGGAGAGGTCGCTACCAAGCGTG  >  NC_000913/3519573‑3520053

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: