Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,077,161 A→C 7.5% H245Q (CAT→CAG)  PP_2687 ← hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029473,077,1610AC 7.5% 253.8 / 10.5 93H245Q (CAT→CAG) PP_2687hypothetical protein
Reads supporting (aligned to +/- strand):  ref base A (45/41);  new base C (2/5);  total (47/46)
Fisher's exact test for biased strand distribution p-value = 2.67e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.77e-01

CTTAGCCAGGGATGCAACCCTGCATCGAGGATTTCGTCGATTGCGGTGTAGCGCAGGCGCGCGTTCAGTTCGGCGGCCATGCGCTGGGCGGCACGCCCGGTACTGCCCGGCAGGCCTGCCAGGATCAGGTCCAGTTCCTCGATGCAGGCATGCAGCGAGCGCGG‑CACGTCCACCCGCAACAGCAGCAACTCGGACACCGGCCGGGCGCTGGGTGCGGCGCGATAGAGCTCGTTGAAGGCTTCGTACGAGGTCAGCGCGCGCAGCAGGGCGCTCCATTGGTAATAACCGCGGGCGGAGTC  >  NC_002947/3077012‑3077310
                                                                                                                                                     |                                                                                                                                                      
cTTAGCCAGGGATGCAACCCTGCATCGAGGGAGTGGGCGTGTGCGGGGTGGGGCATGCGGGCGTTCAGGTGGGAGCACCTGCTAGGTCAGTCTCGTCCTGTACTGCCCGGCAGGCCTGCCAGGATCAGGTCCAGTTCCTCGATGCAGGCa                                                                                                                                                        >  1:1751237/1‑150 (MQ=255)
 ttAGCCAGGGATGCAACCCTGCATCGAGGATTTCGTCGATTGCGGTGTAGCGCAGGCGCGCGTTCAGTTCGGCGGCCATGCGCTGGGCGGCACGCCCGGTACTGCCCGGCAGGCCTGCCAGGATCAGGTCCTGTTCCTCGATGCAGGCAt                                                                                                                                                       <  1:81694/150‑1 (MQ=255)
 ttAGCCAGGGATGCAACCCTGCATCGAGGATTTCGTCGATTGCGGTGTAGCGCAGGCGCGCGTTCAGTTCGGCGGCCATGCGCTGGGCGGCACGCCCGGTACTGCCCGGCAGGCCTGCCAGGATCAGGTCCAGTTCCTCGATGCAGGCAt                                                                                                                                                       <  2:2365115/150‑1 (MQ=255)
     ccAGGGATGCAACCCTGCATCGAGGATTTCGTCGATTGCGGTGTAGCGCAGGCGCGCGTTCAGTTCGGCGGCCATGCGCTGGGCGGCACGCCCGGTACTGCCCGGCAGGCCTGCCAGGATCAGGTCCAGTTCCTCGATGCAGGCATGCag                                                                                                                                                   >  2:1186889/1‑150 (MQ=255)
     ccAGGGATGCAACCCTGCATCGAGGATTTCGTCGATTGCGGTGTAGCGCAGGCGCGCGTTCAGTTCGGCGGCCATGCGCTGGGCGGCACGCCCGGTACTGCCCGGCAGGCCTGCCAGGATCAGGTCCAGTTCCTCGATGCAGGCATGCag                                                                                                                                                   >  1:202708/1‑150 (MQ=255)
      cAGGGATGCAACCCTGCATCGAGGATTTCGTCGATTGCGGTGTAGCGCAGGCGCGCGTTCAGTTCGGCGGCCATGCGCTGGGCGGCACGCCCGGTACTGCCCGGCAGGCCTGCCAGGATCAGGTCCAGTTCCTCGATGCAGGCATGCagc                                                                                                                                                  >  1:1673198/1‑150 (MQ=255)
        gggATGCAACCCTGCATCGAGGATTTCGTCGATTGCGGTGTAGCGCAGGCGCGCGTTCAGTTCGGCGGCCGTGCGCTGGGCGGTACGCCCGGTACTGCCCGGCAGGCCTGCCAGGATCAGGTCCAGTTCCTCGATGCAGGCATGCagcga                                                                                                                                                <  1:1989655/150‑1 (MQ=255)
         ggATGCAACCCTGCATCGAGGATTTCGTCGATTGCGGTGTAGCGCAGGCGCGCGTTCAGTTCGGCGGCCATGCGCTGGGCGGCACGCCCGGTACTGCCCGGCAGGCCTGCCAGGATCAGGTCCAGTTCCTCGCTGCATGCCTGCatcgag                                                                                                                                               <  1:2297467/150‑1 (MQ=255)
         ggATGCAACCCTGCATCGAGGATTTCGTCGATTGCGGTGTAGCGCAGGCGCGCGTTCAGTTCGGCGGCCATGCGCTGGGCGGCACGCCCGCTATCTTTTGGCAGGCCTGCCAGGATCAGGTCCAGTTCCTCGCTGCATGCCTGCatcgag                                                                                                                                               <  2:2433251/150‑1 (MQ=255)
             gCAACCCTGCATCGAGGATTTCGTCGATTGCGGTGTAGCGCAGGCGCGCGTTCAGTTCGGCGGCCATGCGCTGGGCGGCACGCCCGGTACTGCCCGGCAGGCCTGCCAGGATCAGGTCCAGTTCCTCGATGCAGGCATGCAGCGAGcgcg                                                                                                                                           >  1:1079430/1‑150 (MQ=255)
             gCAACCCTGCATCGAGGATTTCGTCGATTGCGGTGTAGCGCAGGCGCGCGTTCAGTTCGGCGGCCATGCGCTGGGCGGCACGCCCGGTACTGCCCGGCAGGCCTGCCAGGATCAGGTCCAGTTCCTCGATGCAGGCATGCAGCGAGcgcg                                                                                                                                           <  1:1699932/150‑1 (MQ=255)
               aaCCCTGCATCGAGGATTTCGTCGATTGCGGTGTAGCGCAGGCGCGCTTTCAGTTCGGCGGCCATGCGCTGGGCGGCACGCCCGGTACTGCCCGGCAGGCCTGCCAGGATCAGGTCCAGTTCCTCGCTGCATGCCTGCatcgaggaactg                                                                                                                                         <  2:100667/150‑7 (MQ=255)
                   cTGCATCGAGGATTTCGTCGATTGCGGTGTAGCGCAGGCGCGCGTTCAGTTCGGCGGCCATGCGCTGGGCGGCACGCCCGGTACTGCCCGGCAGGCCTGCCAGGATCAGGTCCAGTTCCTCGCTGCATGCCTGCatcga                                                                                                                                                <  2:1669799/139‑1 (MQ=255)
                   cTGCATCGAGGATTTCGTCGATTGCGGTGTAGCGCAGGCGCGCGTTCAGTTCGGCGGCCATGCGCTGGGCGGCACGCCCGGTACTGCCCGGCAGGCCTGCCAGGATCAGGTCCAGTTCCTCGCTGCATGCCTGCatcga                                                                                                                                                >  1:1669799/1‑139 (MQ=255)
                      cagcGAGGATTTCGTCGATTGCGGTGTAGCGCAGGCGGGCGTTCAGTTCGGCGGCCATGCGCTGGGCGGCACGCCCGGTACTGCCCGGCAGGACTGCCAGGATCAGGTCCAGTTCCTCGATGCAGGCATGCAGCGAGCGCGG‑CACGTCCa                                                                                                                                 <  1:764705/147‑1 (MQ=255)
                       aTCGAGGATTTCGTCGATTGCGGTGTAGCGCAGGCGCGCGTTCAGTTCGGCGGCCATGCGCTGGGCGGCACGCCCGGTACTGCCCGGCAGGCCTGCCAGGATCAGGTCCAGTTCCTCGATGCAGGCATGCAGCGAGCGCGG‑CACGTCCAc                                                                                                                                <  2:1673198/150‑1 (MQ=255)
                            ggATTTCGTCGATTGCGGTGTAGCGCAGGCGCGCGTTCAGTTCGGCGGCCATGCGCTGGGCGGCACGCCCGGTACTGCCCGGCAGGCCTGCCAGGATCAGGTCCAGTTCCTCGATGCAGGCATGCAGCGAGCGCGG‑CACGTCCACCCGCa                                                                                                                           >  2:1556581/1‑150 (MQ=255)
                                   gtcgATTGCGGTGTAGCGCAGGCGCGCGTTCAGTTCGGCGGCCATGCGCTGGGCGGCACGCCCGGTACTGCCCGGCAGGCCTGCCAGGATCAGGTCCAGTTCCTCGATGCAGGCATGCAGCGAGCGCGG‑CACGTCCACCCGCAAcagcag                                                                                                                    <  1:2005083/150‑1 (MQ=255)
                                   gtcgATGGCGGTGTAGCGCAGGCGCGCGTTCAGTTCGGCGGCCATGCGCTGGGCGGCACGCCCGGTACTGCCCGGCAGGCCCGCCAGGATCAGGTCCAGTTCCTCGATGCCGGCATGCAGCGAGCGCGG‑CACGTCCACCCGCAAcagcag                                                                                                                    <  2:122710/150‑1 (MQ=255)
                                    tcgATTGCGGTGTAGCGCAGGCGCGCGTTCAGTTCGGCGGCCATGCGCTGGGCGGCACGCCCGGTACTGCCCGGCAGGCCTGCCAGGATCAGGTCCAGTTCCTCGATGCAGGCATGCAGCGAGCGCGG‑CACGTCCACCCGCAACagcagc                                                                                                                   >  2:252897/1‑150 (MQ=255)
                                    tcgATTGCGGTGTAGCGCAGGCGCGCGTTCAGTTCGGCGGCCATGCGCTGGGCGGCACGCCCGGTACTGCCCGGCAGGCCTGCCAGGATCAGGTCCAGTTCCTCGATGCAGGCATGCAGCGAGCGCGG‑CACGTCCACCCGCAACagcagc                                                                                                                   >  1:566191/1‑150 (MQ=255)
                                         tGCGGTGTAGCGCAGGCGCGCGTTCAGTTCGGCGGCCATGCGCTGGGCGGCACGCCCGTTACTGCCCGGCAGGCCTGCCAGGATCAGGTCCAGTTCCTCGATGCAGGCATGCAGCGAGCGCGG‑CACGTCCACCCGCAACAGCAGCAACTc                                                                                                              <  2:566191/150‑1 (MQ=255)
                                                  gcgcAGGCGCGCGTTCAGTTCGGCGGCCATGCGCTGGGCGGCGCGCCCGGTACTGCCCGGCAGGCCTGCCAGGATCAGGTCCAGTTCCTCGATGCAGGCATGCAGCGAGCGCGG‑CACGTCCACCCGCAACAGCAGCAACTCGGACAccgg                                                                                                     >  1:555591/1‑150 (MQ=255)
                                                      aGGCGCGCGTTCAGTTCGGCGGCCATGCGCTGGGCGGCACGCCCGGTACTGCCCGGCAGGCCTGCCAGGATCAGGTCCAGTTCCTCGATGCAGGCATGCAGCGAGCGCGG‑CACGTCCACCCGCAACAGCAGCAACTCGGACAccggccgg                                                                                                 >  2:1951142/1‑150 (MQ=255)
                                                          gcgcgTTCAGTTCGGCGGCCATGCGCTGGGCGGCACGCCCGGTACTGCCCGGCAGGCCTGCCAGGATCAGGTCCAGTTCCTCGATGCAGGCATg                                                                                                                                                      >  1:744628/1‑94 (MQ=255)
                                                          gcgcgTTCAGTTCGGCGGCCATGCGCTGGGCGGCACGCCCGGTACTGCCCGGCAGGCCTGCCAGGATCAGGTCCAGTTCCTCGATGCAGGCATg                                                                                                                                                      <  2:744628/94‑1 (MQ=255)
                                                            gcgTTCAGTTCGGCGGCCATGCGCTGGGCGGCACGCCCGGTACTGCCCGGCAGGCCTGCCAGGATCAGGTCCAGTTCCTCGATGCAGGCATGCAGCGAGCGCGG‑CACGTCCACCCGCAACAGCAGCAACTCGGACACCGGCCGGGCGCTg                                                                                           <  2:264701/150‑1 (MQ=255)
                                                             cgTTCAGTTCGGCGGCCATGCGCTGGGCGGCACGCCCGGTACTGCCCGGCAGGCCTGCCAGGATCAGGTCCAGTTCCTCGATGCAGGCATGCAGCGAGCGCGG‑CACGTCCACCCGCAACAGCAGCAACTCGGACACCGGCCGGGCGCTgg                                                                                          <  1:396175/150‑1 (MQ=255)
                                                             cgTTCAGTTCGGCGGCCATGCGCTGGGCGGCACGCCCGGTACTGCCCGGCAGGCCTGCCAGGATCAGGTCCAGTTCCTCGATGCAGGCATGCAGCGAGCGCGG‑CACGGCCACCCGCAACAGCAGCAACTCGGACACCGGCCGGGCGCTgg                                                                                          <  2:710451/150‑1 (MQ=255)
                                                              gTTCAGTTCGGCGGCCATGCGCTGGGCGGCACGCCCGGTACTGCCCGGCGGGCCTGCCAGGATCAGGTCCAGTTCCTCGATGCAGGCATGCAGCGAGCGCGG‑CACGTCCACCCGCAACAGCAGCAACTCGGACACCGGCCGGGCGCTggg                                                                                         >  1:1879154/1‑150 (MQ=255)
                                                              gTTCAGTTCGGCGGCCATGCGCTGGGCGGCACGCCCGGTACTGCCCGGCAGGCCTGCCAGGATCAGGTCCAGTTCCTCGATGCAGGCATGCAGCGAGCGCGG‑CACGTCCACCCGCAACAGCAGCAACTCGGACACCGGCCGGGCGCTggg                                                                                         >  2:1200447/1‑150 (MQ=255)
                                                               ttCAGTTCGGCGGCCATGCGCTGGGCGGCACGCCCGGTACTGCCCGGCAGGCCTGCCAGGATCAGGTCCAGTTCCTCGATGCAGGCATGCAGCGAGCGCGG‑CACGTCCACCCGCAACAGCAGCAACTCGGACACCGGCCGGGCGCTGGGt                                                                                        >  2:1652528/1‑150 (MQ=255)
                                                               ttCAGTTCGGCGGCCATGCGCTGGGCGGCACGCCCGGTACTGCCCGGCAGGCCTGCCAGGATCAGGTCCAGTTCCTCGATGCAGGCATGCAGCGAGCGCGG‑CACGTCCACCCGCAACAGCAGCAACTCGGACACCGGCCGGGCGCTGGGt                                                                                        >  2:2447998/1‑150 (MQ=255)
                                                                      cggcggCCATGCGCTGGGCGGCACGCCCGGTACTGCCCCGCAGGCCTGCCAGGATCAGGTCCAGTTCCTCGATGCAGGCATGCAGCGAGCGCGG‑CACGTCCACCCGCAACAGCAGCAACTCGGACACCGGCCGGGCGCTGGGTGCGgcgc                                                                                 <  1:623780/150‑1 (MQ=255)
                                                                        gcggcCATGCGCTGGGCGGCACGCCCGGTACTGCCCGGCAGGCCTGCCAGGATCAGGTCCAGTTCCTCGATGCAGGCATGCAGCGAGCGCGG‑CACGTCCACCCGCAACAGCAGCAACTCGGACACCGGCCGGGCGCTGGGTGCGGCGCGa                                                                               <  1:2447998/150‑1 (MQ=255)
                                                                          ggcCATGCGCTGGGCGGCACGCCCGGTACTGCCCGGCAGGCCTGCCAGGATCAGGTCCAGTTCCTCGATGCAGGCATGCAGCGAGCGCGG‑CACGTCCACCCGCAACAGCAGCAACTCGGACACCGGCCGGGCGCTGGGTGCGGCGCGATa                                                                             <  2:955664/150‑1 (MQ=255)
                                                                          ggcCATGCGCTGGGCGGCACGCCCGGTACTGCCCGGCAGGCCTGCCAGGATCAGGTCCAGTTCCTCGATGCAGGCATGCAGCGAGCGCGG‑CACGTCCACCCGCAACAGCAGCAACTCGGACACCGGCCGGGCGCTGGGTGCGGCGCGATa                                                                             <  1:252897/150‑1 (MQ=255)
                                                                                 cgcTGGGCGGCACGCCCGGTACTGCCCGGCAGGCCTGCCAGGATCAGGTCCAGTTCCTCGCTGCATGCCTGCatcgagga                                                                                                                                             <  1:684706/80‑3 (MQ=37)
                                                                                 cgcTGGGCGGCACGCCCGGTACTGCCCGGCAGGCCTGCCAGGATCAGGTCCAGTTCCTCGCTGCATGCCTGCatcgagga                                                                                                                                             >  2:684706/1‑78 (MQ=37)
                                                                                 cgcTGGGCGGCACGCCCGGTACTGCCCGGCAGGCCTGCCAGGATCAGGTCCAGTTCCTCGATGCAGGCATGCAGCGAGCGCGG‑CACGTCCACCCGCAACAGCAGCAACTCGGACACCGGCCGGGCGCTGGGTGCGGCGCGATAGAGCTCg                                                                      <  2:1089615/150‑1 (MQ=255)
                                                                                   cTGGGCGGCACGCCCGGTACTGCCCGGCAGGCCTGCCAGGATCAGGTCCAGTTCCTCGATGCAGGCATGCAGCGAGCGCGG‑CACGTCCACCCGCAACAGCAGCAACTCGGACACCGGCCGGGCGCTGGGTGCGGCGCGATAGAGCTCGtt                                                                    <  2:697433/150‑1 (MQ=255)
                                                                                   cTGGGCGGCACGCCCGGTACTGCCCGGCAGGCCTGCCAGGATCAGGTCCAGTTCCTCGATGCAGGCATGCAGCGAGCGCGG‑CACGTCCACCCGCAACAGCAGCAACTCGGACACCGGCCGGGCGCTGGGTGCGGCGCGATAGAGCTCGtt                                                                    <  2:1068069/150‑1 (MQ=255)
                                                                                    tGGGCGGCACGCCCGGTACTGCCCGGCAGGCCTGCCAGGATCAGGTCCAGTTCCTCGATGCAGGCATGCAGCGAGCGCGG‑CACGTCCACCCGCAACAGCAGCAACTCGGACACCGGCCGGGCGCTGGGTgcg                                                                                     >  2:288263/1‑132 (MQ=255)
                                                                                    tGGGCGGCACGCCCGGTACTGCCCGGCAGGCCTGCCAGGATCAGGTCCAGTTCATCGATGCAGGCATGCAGCGAGCGCGG‑CACGTCCACCCGCAACAGCAGCAACTCGGACACCGGCCGGGCGCTGGGTgcg                                                                                     <  1:288263/132‑1 (MQ=255)
                                                                                      ggcggcACGCCCGGTACTGCCCGGCAGGCCTGCCAGGATCAGGTCCAGTTCCTCGATGCAGGCATGCAGCGAGCGCGG‑CACGTCCACCCGCAACAGCAGCAACTCGGACACCGGCCGGGCGCTGGGTGCGGCGCGATAGAGCTCGTTGaa                                                                 <  2:1934768/150‑1 (MQ=255)
                                                                                          gcACGCCCGGTACTGCCCGGCAGGCCTGCCAGGATCAGGTCCAGTTCCTCGATGCAGGCATGCAGCGAGCGCGG‑CACGTCCACCCGCAACAGCAGCAACTCGGACACCGGCCGGGCGCTGGGTGCGGCGCGATAGAGCTCGTTGAAGGCt                                                             <  2:884071/150‑1 (MQ=255)
                                                                                            aCGCCCGGTACTGCCCGGCAGGCCTGCCAGGATCAGGTCCAGTTCCTCGCTCCAGTCATGCAGCGAGCGCGG‑CACGTCCACCCGCAACAGCAGCAACTCGGACACCGGCCGGGCGCTGGGTGCGGCGCGATAGAGCTCGTTGAAGGCTTc                                                           >  1:2132501/1‑150 (MQ=255)
                                                                                             cGCCCGGTACTGCCCGGCAGGCCTGCCAGGATCAGGTACAGTTCCTCGATGCAAGCATGCAGCGAGCGCGG‑CACGTCCACCCGCAACAGCAGCAACTCGGACACCGGCCGGGCGCTGGGTGCGGCGCGATAGAGCTCGTTGAAGGCTTCg                                                          >  2:1044467/1‑150 (MQ=255)
                                                                                              gCCCGGTACTGCCCGGCAGGCCTGCCAGGATCAGGTCCAGTTCCTCGATGCAGGCATGCAGCGAGCGCGG‑CACGTCCACCCGCAACAGCAGCAACTCGGACACCGGCCGGGCGCTGGGTGCGGCGCGATAGAGCTCGTTGAAGGCTTCGt                                                         >  2:280729/1‑150 (MQ=255)
                                                                                              gCCCGGTACTGCCCGGCAGGCCTGCCAGGATCAGGTCCAGTTCCTCGATGCAGGCATGCAGCGAGCGCGG‑CACGTCCACCCGCAACAGCAGCAACTCGGACACCGGCCGGGCGCTGGGTGCGGCGCGATAGAGCTCGTTGAAGGCTTCGt                                                         <  1:1556581/150‑1 (MQ=255)
                                                                                                  ggTACTGCCCGGCAGGCCTGCCAGGATCAGGTCCAGTTCCTCGATGCAGGCATGCAGCGAGCGCGG‑CACGTCCACCCGCAACAGCAGCAACTCGGACACCGGCCGGGCGCTGGGTGCGGCGCGATAGAGCTCGTTGAAGGCTTCGTACGa                                                     <  2:426236/150‑1 (MQ=255)
                                                                                                  ggTACTGCCCGGCAGGCCTGCCAGGATCAGGTCCAGTTCCTCGATGCAGGCATGCAGCGAGCGCGG‑CACGTCCACCCGCAACAGCAGCAACTCGGACACCGGCCGGGCGCTGGGTGCGGCGCGATAGAGCTCGTTGAAGGCTTCGTACGa                                                     <  1:2234267/150‑1 (MQ=255)
                                                                                                      cTGCCCGGCAGGCCTGCCAGGATCAGGTCCAGTTCCTCGATGCAGGCATGCAGCGAGCGCGG‑CACGTCCACCCGCAACAGCAGCAACTCGGACACCGGCCGGGCGCTGGGTGCGGCGCGATAGAGCTCGTTGAAGGCTTCGTACGAGGTc                                                 >  2:1532281/1‑150 (MQ=255)
                                                                                                       tGCCCGGCAGGCCTGCCAGGATCAGGTCCAGTTCCTCGATGCAGGCATGCAGCGAGCGCGG‑CACGTCCACCCGCAACAGCAGCAACTCGGACACCGGCCGGGCGCTGGGTGCGGCGCGATAGAGCTCGTTGAAGGCTTCGTACGAGGTCa                                                >  2:2387305/1‑150 (MQ=255)
                                                                                                       tGCCCGGCAGGCCTGCCAGGATCAGGTCCAGTTCCTCGATGCAGGCATGCAGCGACCGCGG‑CACGTCCACCCGCAACAACAGCAACTCGGACACCGGCCGGGCCCTGGGTGCGGCGCGATAGATCTCGGTGAAGGCTTCGTACgctggcc                                                >  1:531451/1‑144 (MQ=255)
                                                                                                         cctggCAGGCCTGCCGGGATCAGGTCCAGTTCCTCGATGCAGGCATGCAGCGAGCGCGG‑CACGTCCACCCGCAACAGAAACCACACGGCTAACGTCCGGTCGCTGGGTGCGGCGCGATAGAGCTCGTTGAAGGCTTCGTACGAGGTCAgc                                              >  1:281366/4‑150 (MQ=255)
                                                                                                          ccGGCAGGCCTGCCAGGATCAGGTCCAGTTCCTCGATGCAGGCATGCAGCGAGCGCGG‑CACGTCCACCCGCAACAGCAGCAACTCGGACACCGGCCGGGCGCTGGGTGCGGCGCGATAGAGCTCGTTGAAGGCTTCGTACGAGGTCAgcg                                             <  1:1094107/150‑1 (MQ=255)
                                                                                                                   cTGCCAGGATCAGGTCCAGTTCCTCGATGCAGGCATGCAGCGAGCGCGG‑CACGTCCACCCGCAACAGCAGCAACTCGGACACCGGCCGGGCGCTGGGTGCGGCGCGATAGAGCTCGTTGAAGGCTTCGTACGAGGTCAGCGCGCgcagca                                    >  2:1078963/1‑150 (MQ=255)
                                                                                                                    tGCCAGGATCAGGTCCAGTTCCTCGATGCAGGCATGCAGCGAGCGCGG‑CACGTCCACCCGCAACAGCAGCAACTCGGACACCGGCCGGGCGCTGGGTGCGGCGCGATAGAGCTCGTTGAAGGCTTCGTACGAGGTCAGCGCGCGcagcag                                   >  1:227337/1‑150 (MQ=255)
                                                                                                                     gCCAGGATCAGGTCCAGTTCCTCGATGCAGGCATGCAGCGAGCGCGG‑CACGTCCACCCGCAACAGCAGCAACTCGGACACCGGCCGGGCGCTGGGTGCGGCGCGATAGAGCTCGTTGAAGGCTTCGTACGAg                                                    <  2:2216704/132‑1 (MQ=255)
                                                                                                                     gCCAGGATCAGGTCCAGTTCCTCGATGCAGGCATGCAGCGAGCGCGG‑CACGTCCACCCGCAACAGCAGCAACTCGGACACCGGCCGGGCGCTGGGTGCGGCGCGATAGAGCTCGTTGAAGGCTTCGTACGAg                                                    >  1:2216704/1‑132 (MQ=255)
                                                                                                                     gCCAGGATCAGGTCCAGTTCCTCGATGCAGGCATGCAGCGAGCGCGG‑CACGTCCACCCGCAACAGCAGCAACTCGGACACCGGCCGGGCGCTGGGTGCGGCGCGATAGAGCTCGTTGAAGGCTTCGTACGAGGTCAGCGCGCGCAGCAgg                                  <  1:475602/150‑1 (MQ=255)
                                                                                                                     gCCAGGATCAGGTCCAGTTCCTCGATGCAGGCATGCAGCGAGCGCGG‑CACGTCCACCCGCAACAGCAGCAACTCGGACACCGGCCGGGCGCTGGGTGCGGCGCGATAGAGCTCGTTGAAGGCTTCGTACGAGGTCAGCGCGCGCAGCAgg                                  >  2:865/1‑150 (MQ=255)
                                                                                                                      ccAGGATCAGGTCCAGTTCCTCGATGCAGGCATGCAGCGAGCGCGG‑CACGTCCACCCGCAACAGCAGCAACTCGGACACCGGCCGGGCGCTGGGTGCGGCGCGATAGAGCTCGTTGAAGGCTTCGTACGAGGTCAGCGCGCGCAGCAggg                                 >  2:691606/1‑150 (MQ=255)
                                                                                                                       cAGGATCAGGTCCAGTTCCTCGATGCAGGCATGCAGCGAGCGCGG‑CACGTCCACCCGCAACAGCAGCAACTCGGACACCGGCCGGGCGCTGGGTGCGGCGCGATAGAGCTCGTTGAAGGCTTCGTACGAGGTCAGCACGCGCAGCAGGgc                                >  2:1292255/1‑150 (MQ=255)
                                                                                                                          gATCAGGTCCAGTTCCTCGATGCAGGCATGCAGCGAGCGCGG‑CACGTCCACCCGCAACAGCAGCAACTCGGACACCGGCCGGGCGCTGGGTGCGGCGCGATAGAGCTCGTTGa                                                                  <  2:1413284/113‑1 (MQ=255)
                                                                                                                          gATCAGGTCCAGTTCCTCGATGCAGGCATGCAGCGAGCGCGG‑CACGTCCACCCGCAACAGCAGCAACTCGGACACCGGCCGGGCGCTGGGTGCGGCGCGATAGAGCTCGTTGa                                                                  >  1:1413284/1‑113 (MQ=255)
                                                                                                                               ggTCCAGTTCCTCGATGCAGGCATGCAGCGAGCGCGG‑CACGTCCACCCGCAACAGCAGCAACTCGGACACCGGCCGGGCGCTGGGTGCGGCGCGATAGAGCTCGTTGAAGGCTTCGTACGAGGTCAGCGCGCGCAGCAGGGCGCTCCAtt                        >  1:345885/1‑150 (MQ=255)
                                                                                                                                gTCCAGTTCCTCGATGCAGGCATGCAGCGAGCGCGG‑CACGTCCACCCGCAACAGCAGCAACTCGGACACCGGCCGGGCGCTGGGTGCGGCGCGATCTTGCACCTTGAAGGCTTCGTACGAGGTCAGCGCGCGCAGCAGGGCGCTCCATTg                       <  1:1995672/150‑1 (MQ=255)
                                                                                                                                gTCCAGTTCCTCGATGCAGGCATGCAGCGAGCGCGG‑CACGTCCACCCGCAACAGCAGCAACTCGGACACCGGCCGGGCGCTGGGTGCGGCGCGATAGAGCTCGTTGAAGGCTTCGTACGAGGTCAGCGCGCGCAGCAGGGCGCTCCATTg                       >  1:1599375/1‑150 (MQ=255)
                                                                                                                                gTCCAGTTCCTCGATGCAGGCATGCAGCGAGCGCGG‑CACGTCCACCCGCAACAGCAGCAACTCGGACACCGGCCGGGCGCTGGGTGCGGCGCGATAGAGCTCGTTGAAGGCTTCGTACGAGGTCAGCGCGCGCAGCAGGGCGCTCCATTg                       <  1:1652528/150‑1 (MQ=255)
                                                                                                                                 tCCAGTTCCTCGATGCAGGCATGCAGCGAGCGCGG‑CACGTCCACCCGCAACAGCAGCAACTCGGACACCGGCCGGGCGCTGGGTGCGGCGCGATAGAGCTCGTTGAAGGCTTCGTACGAGGTCAGCGCGCGCAGCAGGGCGCTCCATTgg                      >  2:1137901/1‑150 (MQ=255)
                                                                                                                                  ccAGTTCCTCGATGCAGGCATGCAGCGAGCGCGG‑CACGTCCACCCGCAACAGCAGCAACTCGGACACCGGCCGGGCGCTGGGTGCGGCGCGATAGAGCTCGTTGAAGGCTTCGTACGAGGTCAGCGCGCGCAGCAGGGCGCTCCATTGGt                     <  1:824340/150‑1 (MQ=255)
                                                                                                                                     gTTCCTCGATGCAGGCATGCAGCGAGCGCGG‑CACGTCCACCCGCAACAGCAGCAACTCGGACACCGGCCGGGCGCTGGGTGCGGCGCGATAGAGCTCGTTGAAGGCTTCGTACGAGGTCAGCGCGCGCAGCAGGGCGCTCCATTGGtaat                  <  1:1044467/150‑1 (MQ=255)
                                                                                                                                      ttCCTCGATGCAGGCATGCAGCGAGCGCGG‑CACGTCCACCCGCAACAGCAGCAACTCGGACACCGGCCGGGCGCTGGGTGCGGCGCGATAGAGCTCGtt                                                                    >  1:1541478/1‑99 (MQ=255)
                                                                                                                                      ttCCTCGATGCAGGCATGCAGCGAGCGCGG‑CACGTCCACCCGCAACAGCAGCAACTCGGACACCGGCCGGGCGCTGGGTGCGGCGCGATAGAGCTCGtt                                                                    <  2:1541478/99‑1 (MQ=255)
                                                                                                                                      ttCCTCGATGCAGGCATGCAGCGAGCGCGG‑CACGTCCACCCGCAACAGCAGCAACTCGGACACCGGCCGGGCGCTGGGTGCGGCGCGATAGAGCTCGTTGAAGGCTTCGTACGAGGTCAGCGCGCGCAGCAGGgc                                >  2:2011297/1‑135 (MQ=255)
                                                                                                                                      ttCCTCGATGCAGGCATGCAGCGAGCGCGG‑CACGTCCACCCGCAACAGCAGCAACTCGGACACCGGCCGGGCGCTGGGTGCGGCGCGATAGAGCTCGTTGAAGGCTTCGTACGAGGTCAGCGCGCGCAGCAGGgc                                <  1:2011297/135‑1 (MQ=255)
                                                                                                                                         cTCGATGCAGGCATGCAGCGAGCGCGG‑CACGTGCACCCGCAACAGCAGCAACTCGGACACCGGCCGGGCGCTGGGTGCGGCGCGATAGAGCTCGTTGAAGGCTTCg                                                          >  2:1345386/1‑106 (MQ=255)
                                                                                                                                         cTCGATGCAGGCATGCAGCGAGCGCGG‑CACGTCCACCCGCAACAGCAGCAACTCGGACACCGGCCGGGCGCTGGGTGCGGCGCGATAGAGCTCGTTGAAGGCTTCg                                                          <  1:1345386/106‑1 (MQ=255)
                                                                                                                                         cTCGATGCAGGCATGCAGCGAGCGCGG‑CACGTCCACCCGCAACAGCAGCAACTCGGACACCGGCCGGGCGCTGGGTGCGGCGCGATAGAGCTCGTTGAAGGCTTCGTACGAGGTCAGCGCGCGCAGCAGGGCGCTCCATTGGTAATAAcc              >  2:889079/1‑150 (MQ=255)
                                                                                                                                         cTCGATGCAGGCATGCAGCGAGCGCGG‑CACGTCCACCCGCAACAGCAGCAACTCGGACACCGGCCGGGCGCTGGGTGCGGCGCGATAGAGCTCGTTGAAGGCTTCGTACGAGGTCAGCGCGCGCAGCAGGGCGCTCCATTGGTAATAAcc              >  1:778709/1‑150 (MQ=255)
                                                                                                                                         cTCGATGCAGGCATGCAGCGAGCGCGG‑CACGTCCACCCGCAACAGCAGCAACTCGGACACCGGCCGGGCGCTGGGTGCGGCGCGATAGAGCTCGTTGAAGGCTTCGTACGAGGTCAGCGCGCGCAGCAGGGCGCTCCATTGGTAATAAcc              >  2:316069/1‑150 (MQ=255)
                                                                                                                                         cTCGATGCAGGCATGCAGCGAGCGCGG‑CACGTCCACCCGCAACAGCAGCAACTCGGACACCGGCCGGGCGCTGGGTGCGGCGCGATAGAGCTCGTTGAAGGCTTCGTACGAGGTCAGCGCGCGCAGCAGGGCGCTCCATTGGTAATAAcc              >  2:388124/1‑150 (MQ=255)
                                                                                                                                         cTCGATGCAGGCATGCAGCGAGCGCGG‑CACGTCCACCCGCAACAGCAGCAACTCGGACACCGGCCGGGCGCTGGGTGCGGCGCGATAGAGCTCGTTGAAGGCTTCGTACGAGGTCAGCGCGCGCAGCAGGGCGCTCCATTGGTAATAAcc              >  2:526770/1‑150 (MQ=255)
                                                                                                                                         cTCGATGCAGGCATGCAGCGAGCGCGG‑CACGTCCACCCGCAACAGCAGCAACTCGGACACCGGCCGGGCGCTGGGTGCGGCGCGATAGAGCTCGTTGAAGGCTTCGTACGAGGTAAGCGCGCGCAGCAGGGCGCTCCATTGGTAATAAcc              >  2:1694838/1‑150 (MQ=255)
                                                                                                                                          tCGATGTAGGCATGCAGCGAGCGCGG‑CACGTCCACCCGCAACAGCAGCAACTCGGACACCGGCCGGGCGCTGGGTGCGGCGCGATAGAGCTCGTTGAAGGCTTCGTACGAGGTCAGCGCGCGCAGCAGGGCGCTCCATTGGTAATAACcg             >  2:577855/1‑150 (MQ=255)
                                                                                                                                          tCGATGCAGGCATGCAGCGAGCGCGGACAGAACCA‑CCGCAACAGCAGCACCTTAGACACCGGCCGGGCGCTGGGTGCGGCGCGATAGAGCTCGTTGAAGGCTTCGTACGAGGTCAGCGCGCGCAGCAGGGCGCTCCATTGGTAATAACcg             >  1:299870/1‑150 (MQ=255)
                                                                                                                                          tCGATGCAGGCATGCAGCGAGCGCGG‑CACGTCCACCCGCAACAGTATCAACTCGGACACCGGCCGGGCGCTGGGTGCGGCGCGATAGAGCTCGTTGAAGGCTTCGTACGAGGTCAGCGCGCGCAGCAGGGCGCTCCATTGGTAATAACcg             >  1:2399433/1‑150 (MQ=255)
                                                                                                                                               gCAGGCATGCAGCGAGCGCGG‑CACGTCCACCCGCAACAGCAGCAACTCGGACACCGGCCGGGCGCTGGGTGCGGCGCGATAGAGCTCGTTGAAGGCTTCGTACGAGGTCAGCGCGCGCAGCAGGGCGCTCCATTGGTAATAACCgcgggc        >  2:2119017/1‑150 (MQ=255)
                                                                                                                                                  ggCATGCAGCGAGCGCGG‑CACGTCCACCCGCAACAGCAGCAACTCGGACACCGGCCGGGCGCTGGGTGCGGCGCGATAGAGCTCGTTGAAGGCTTCGTACGAGGTCAGCGCGCGCAGCAGGGCGCTCCATTGGTAATAACCGCGGGCGGa     <  2:778709/150‑1 (MQ=255)
                                                                                                                                                  ggCATGCAGCGAGCGCGG‑CACGTCCACCCGCAACAGCAGCAACTCGGACACCGGCCGGGCGCTGGGTGCGGCGCGATAGAGCTCGTTGAAGGCTTCGTACGAGGTCAGCGCGCGCAGCAGGGCGCTCCATTGGTAATAACCGCGGGCGGa     <  2:1984117/150‑1 (MQ=255)
                                                                                                                                                  ggCATGCAGCGAGCGCGG‑CACGTCCACCCGCAACAGCAGCAACTCGGACACCGGCCGGGCGCTGGGTGCGGCGCGATAGAGCTCGTTGAAGGCTTCGTACGAGGTCAGCGCGCGCAGCAGGGCGCTCCATTGGTAATAACCGCGGGCGGa     <  1:316069/150‑1 (MQ=255)
                                                                                                                                                   gCATGAAGCGAGCGCGG‑CACGTCCACCCGCAACAGCAGCAACTCGGACACCGGCCGGGCGCTGGGTGCGGCGCGATAGAGCTCGTTGAAGGCTTCGTACGAGGTCAGCGCGCGCAGCAGGGCGCTCCATTGGTAATAACCGCGGGCGGAg    <  2:1079430/150‑1 (MQ=255)
                                                                                                                                                     aTGCAGCGAGCGCGG‑CACGTCCACCCGCAACAGCAGCAGCTCGGACACCGGCCGGGCGCTGGGTGCGGCGCGATAGAGCTAGTTGGAGTCTTGGTACGAGGTCAGCGCGCGCAGCAGGGCGCTCCATTGGTAATAACCGCGGGCGGAgtc  <  2:202708/150‑1 (MQ=255)
                                                                                                                                                     aTGCAGCGAGCGCGG‑CACGTCCACCCGCAACAGCAGCAACTCGGACACCGGCCGGGCGCTGGGTGCGGCGCGATAGAGCTCGTTGAAGGCTTCGTACGAGGTCAGCGCGCGCAGCAGGGCGCTCCATTGGTAATAACCGCGGGCGGAgtc  <  2:2132501/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                      
CTTAGCCAGGGATGCAACCCTGCATCGAGGATTTCGTCGATTGCGGTGTAGCGCAGGCGCGCGTTCAGTTCGGCGGCCATGCGCTGGGCGGCACGCCCGGTACTGCCCGGCAGGCCTGCCAGGATCAGGTCCAGTTCCTCGATGCAGGCATGCAGCGAGCGCGG‑CACGTCCACCCGCAACAGCAGCAACTCGGACACCGGCCGGGCGCTGGGTGCGGCGCGATAGAGCTCGTTGAAGGCTTCGTACGAGGTCAGCGCGCGCAGCAGGGCGCTCCATTGGTAATAACCGCGGGCGGAGTC  >  NC_002947/3077012‑3077310

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