Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,201,162 G→C 7.4% T134S (ACC→AGC)  PP_4582 ← membrane protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,201,1620GC 7.4% 239.5 / 6.5 76T134S (ACC→AGC) PP_4582membrane protein
Reads supporting (aligned to +/- strand):  ref base G (34/36);  new base C (3/3);  total (37/39)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GTAGGGACGTTGGCACTCATGCCCAGCAGTTCACACATAGGCGTGTCTCCGTTTACAGACGGGGCTCGACCCGGCCTGCGCCACGCGCAGCGGCTGGCGGCAAGTCGTCACGGTAACGGTCGCTGCGCTCGGCAGCCA‑GCGGCGCCTGGGTGCCGAGCTGGTCATCCTCGGCCCGCTGACGCCTGGCGGCGGCTTCGGCCTGCTCACGGCGTGCGTGCGTGCGCTTTTCCAGCGGCCAGCGCAGCAGCACGAAGAGGAAGTACAGGCCGAAGGCGATCATGCCGTACATGGC  >  NC_002947/5201013‑5201304
                                                                                                                                                      |                                                                                                                                              
gTAGGGACGTTGGCACTCATGCCCAGCAGTTCACACATAGGCGTGTCTCCGTTTACAGACGGGGCTCGACCCGGCCTGCGCCACGCGCAGCGGCTGGCGGCAAGTCGTCACGGTAACGGTCGCTGCGCTCGGCAGCCA‑GCGGCGCCTggg                                                                                                                                                >  1:1228644/1‑150 (MQ=255)
      ccGTTGGCACTCCTGCCCAGCAGTTCACACATAGGCGTGTCTCCGATAACAGACGGGGCTCGACCCGGCCTGCGCCACGCGCAGCGGCTGGCGGCAAGTCGTCACGGTAACGGTCGCTGCGCTCGGCACCCAGGCGCCG‑CTGGCTGCCGa                                                                                                                                          <  1:616113/149‑1 (MQ=255)
        gTTGGCACTCATGCCCAGCAGTTCACACATAGGCGTGTCTCCGTTTACAGACGGGGCTCGACCCGGCCTGCGCCACGCGCAGCGGCTGGCGGCAAGTCGTCACGGTAACGGTCGCTGCGCTCGGCAGCCA‑GCGGCGCCTGGGTGCCGAGc                                                                                                                                        >  1:1199533/1‑150 (MQ=255)
           ggCACTCATGCCCAGCAGTTCACACATAGGCGTGTCTCCGTTTACAGACGGGGCTCGACCCGGCCTGCGCCACGCGCAGCGGCTGGCGGCAAGTCGTCACGGTAACGGTCGCTGCGCTCGGCAGCCA‑GCGGCGCCTGGGTGCCGAGCTgg                                                                                                                                     <  2:784593/150‑1 (MQ=255)
            gCACTCATGCCCAGCAGTTCACACATAGGCGTGTCTCCGTTTACAGACGGGGCTCGACCCGGCCTGCGCCACGCGCAGCGGCTGGCGGCAAGTCGTCACGGTAACGGTCGCTGCGCTCGGCAGCCA‑GCGGCGCCTGGGTGCCGAGCTGGt                                                                                                                                    >  2:1535004/1‑150 (MQ=255)
             cACTCATGCCCAGCAGTTCACACATAGGCGTGTCTCCGTTTACAGACGGGGCTCGACCCGGCCTGCGCCACGCGCAGCGGCTGGCGGCAAGTCGTCACGGTAACGGTCGCTGCGCTCGGCAGCCA‑GCGGCGCCTGGGTGCCGAGCTGGTc                                                                                                                                   >  1:1517206/1‑150 (MQ=255)
              aCTCATGCCCAGCAGTTCACACATAGGCGTGTCTCCGTTTACAGACGGGGCTCGACCCGGCCTGCGCCACGCGCAGCGGCTGGCGGCAAGTCGTCACGGTAACGGTCGCTGCGCTCGGCAGCCA‑GCGGCGCCTGGGTGCCGAGCTGGTCa                                                                                                                                  >  1:1421538/1‑150 (MQ=255)
                tCATGCCCAGCAGTTCACACATAGGCGTGTCTCCGTTTACAGACGGGGCTCGACCCGGCCTGCGCCACGCGCAGCGGCTGGCGGCAAGTCGTCACGGTAACGGTCGCTGCGCTCGGCAGCCA‑GCGGCGCCTGGGTGCCGAGCTGGTCATc                                                                                                                                >  1:1646074/1‑150 (MQ=255)
                  aTGCCCAGCAGTTCACACATAGGCGTGTCTCCGTTTACAGACGGGGCTCGACCCGGCCTGCGCCACGCGCAGCGGCTGGCGGCAAGTCGTCACGGTAACGGTCGCTGCGCTCGGCAGCCA‑GCGGCGCCTGGGTGCCGAGCTGGTCATCCt                                                                                                                              >  2:1014550/1‑150 (MQ=255)
                     cccAGCAGTTCACACATAGGCGTGTCTCCGTTTACAGACGGGGCTCGACCCGGCCTGCGCCACGCGCAGCGGCTGGCGGCAAGTCGTCACGGTAACGGTCGCTGCGCTCGGCAGCCA‑GCGGCGCCTGGGTGCCGAGCTGGTCATCCTCgg                                                                                                                           <  1:1014550/150‑1 (MQ=255)
                      ccAGCAGTTCACACATAGGCGTGTCTCCGTTTACAGACGGGGCTCGACCCGGCCTGCGCCACGCGCAGCGGCTGGCGGCAAGTCGTCACGGTAACGGTCGCTGCGCTCGGCAGCCA‑GCGGCGCCTGGGTGCCGAGCTGGTCATCCTCGGc                                                                                                                          <  2:46803/150‑1 (MQ=255)
                          cagTTCACACATAGGCGCCCCTCCGTTTACAGACGGGGCTCGACCCGGCCTGCGCCACGCGCAGCGGCTGGCGGCAAGTCGTCACGGTAACGGTCGCTGCGCTCGGCAGCCA‑GCGGCGCCTGGGTGCCGAGCTGGTCATCCTCGGCCCGc                                                                                                                      <  1:83421/150‑1 (MQ=255)
                               cacacaTAGGCGTGTCTCCGTTTACAGACGGGGCTCGACCCGGCCTGCGCCACGCGCAGCGGCTGGCGGCAAGTCGTCACGGTAACGGTCGCTGCGCTCGGCAGCCA‑GCGGCGCCTGGGTGCCGAGCTGGTCATCCTCGGCCCGCTGACg                                                                                                                 >  1:24310/1‑150 (MQ=255)
                                acacaTAGGCGTGTCTCCGTTTACAGACGGGGCTCGACCCGGCCTGCGCCACGCGCAGCGGCTGGCGGCAAGTCGTCACGGTAACGGTCGCTGCGCTCGGCAGCCA‑GCGGCGCCTGGGTGCCGAGCTGGTCATCCTCGGCCCGCTGACGc                                                                                                                >  2:789590/1‑150 (MQ=255)
                                acacaTAGGCGTGTCTCCGTTTACAGACGGGGCTCGACCCGGCCTGCGCCACGCGCAGCGGCTGGCGGCAAGTCGTCACGGTAACGGTCGCTGCGCTCGGCAGCCA‑GCGGCGCCTGGGTGCCGAGCTGGTCATCCTCGGCCCGCTGACGc                                                                                                                >  1:287189/1‑150 (MQ=255)
                                  acaTAGGCGTGTCTCCGTTTACAGACGGGGCTCGACCCGGCCTGCGCCACGCGCAGCGGCTGGCGGCAAGTCGTCACGGTAACGGTCGCTGCGCTCGGCAGCCA‑GCGGCGCCTGGGTGCCGAGCTGGTCATCCTCGGCCCGCTGACGCCt                                                                                                              >  1:1266116/1‑150 (MQ=255)
                                   caTAGGCGTGTCTCCGTTTACAGACGGGGCTCGACCCGGCCTGCGCCACGCGCAGCGGCTGGCGGCAAGTCGTCACGGTAACGGTCGCTGCGCTCGGCAGCCA‑GCGGCGCCTGGGTGCCGAGCTGGTCATCCTCGGCCGGCTGACGCCTg                                                                                                             <  2:1517206/150‑1 (MQ=255)
                                    aTAGGCGTGTCTCCGTTTACAGACGGGGCTCGACCCGGCCTGCGCCACGCGCAGCGGCTGGCGGCAAGTCGTCACGGTAACGGTCGCTGCGCTCGGCAGCCA‑GCGGCGCCTGGGTGCCGAGCTGGTCATCCTCGGCCCGCTGACGCCTgg                                                                                                            >  2:708421/1‑150 (MQ=255)
                                    aTAGGCGTGTCTCCGTTTACAGACGGGGCTCGACCCGGCCTGCGCCACGCGCAGCGGCTGGCGGCAAGTCGTCACGGTAACGGTCGCTGCGCTCGGCAGCCA‑GCGGCGCCTGGGTGCCGAGCTGGTCATCCTCGGCCCGCTGACGCCTgg                                                                                                            >  2:248909/1‑150 (MQ=255)
                                          gtgtCTCCGTTTACAGACGGGGCTCGACCCGGCCTGCGCCACGCGCAGCGGCTGGCGGCAAGTCGTCACGGTAACGGTCGCTGCGCTCGGCAGCCA‑GCGGCGCCTGGGTGCCGAGCTGGTCATCCTCGGCCCGCTGACGCCTggcggcgg                                                                                                      >  2:850570/1‑150 (MQ=255)
                                               tcCGTTTACAGACGGGGCTCGACCCGGCCTGCGCCACGCGCAGCGGCTGGCGGCAAGTCGTCACGGTAACTGTCGCTGCGCTCGGCAGCCA‑GCGGCGCCTTGCTGCCGATCTTGTCATCCTCGGCCCGCTGACGCCTGGCGGCGGCTTCg                                                                                                 <  1:850570/150‑1 (MQ=255)
                                               tcCGTTTACAGACGGGGCTCGACCCGGCCTGCGCCACGCGCAGCGGCTGGCGGCAAGTCGTCACGGTAACGGTCGCTGCGCTCGGCAGCCA‑GCGGCGCCTGGGTGCCGAGCTGGTCATCCTCGGc                                                                                                                          >  2:689145/1‑125 (MQ=255)
                                               tcCGTTTACAGACGGGGCTCGACCCGGCCTGCGCCACGCGCAGCGGCTGGCGGCAAGTCGTCACGGTAACGGTCGCTGCGCTCGGCAGCCA‑GCGGCGCCTGGGTGCCGAGCTGGTCATCCTCGGc                                                                                                                          <  1:689145/125‑1 (MQ=255)
                                               tcCGTTTACAGACGGGGCTCGACCCGGCCTGCGCCACGCGCAGCGGCTGGCGGCAAGTCGTCACGGTAACGGTCGCTGCGCTCGGCAGCCA‑GCGGCGCCTGGGTGCCGAGCTGGTCATCCTCAGCCCGCGGACGCCTGGCGGCGGCTTCg                                                                                                 >  2:428883/1‑150 (MQ=255)
                                                ccGTTTACAGACGGGGCTCGACCCGGCCTGCGCCACGCGCAGCGGCTGGCGGCAAGTCGTCACGGTAACGGTCGCTGCGCTCGGCAGCCA‑GCGGCGCCTGGGTGCCGAGCTGGTCATCCTCGGCCCGCTGACGCCTGGCGGCGGCTTCgg                                                                                                >  2:819512/1‑150 (MQ=255)
                                                ccGTTTACAGACGGGGCTCGACCCGGCCTGCGCCACGCGCAGCGGCTGGCGGCAAGTCGTCACGGTAACGGTCGCTGCGCTCGGCAGCCA‑GCGGCGCCTGGGTGCCGAGCTGGTCATCCTCGGCCCGCTGACGCCTGGCGGCGGCTTCgg                                                                                                >  2:1364055/1‑150 (MQ=255)
                                                ccGTTTACAGACGGGGCTCGACCCGGCCTGCGCCACGCGCAGCGGCTGGCGGCAAGTCGTCACGGTAACGGTCGCTGCGCTCGGCAGCCA‑GCGGCGCCTGGGTGCCGAGCTGGTCATCCTAGGCCCGCTGACGCCTGGCGGCGGCTTCgg                                                                                                <  1:1535004/150‑1 (MQ=255)
                                                  gTTTACAGACGGGGCTCGACCCGGCCTGCGCCACGCGCAGCGGCTGGCGGCAAGTGGTCACGGTACCGGTCGCTGCGCTCGGCAGCCA‑GCGGCGCCTGGGGGCCGACCTGGTCATCCTCGGCCCGCGGGCGCCGGGAGGCGGCTTCGGcc                                                                                              >  2:762275/1‑150 (MQ=255)
                                                  gTTTACAGACGGGGCTCGACCCGGCCTGCGCCACGCGCAGCGGCTGGCGGCAAGTCGTCACGGTAACGGTCGCTGCGCTCGGCAGCCA‑GCGGCGCCTGGGTGCCGAGCTGGTCATCCTCGGCCCGCTGACGCCTGGCGGCGGCTTCGGcc                                                                                              <  1:492326/150‑1 (MQ=255)
                                                  gTTTACAGACGGGGCTCGACCCGGCCTGCGCCACGCGCAGCGGCTGGCGGCAAGTCGTCACGGTAACGGTCGCTGCGCTCGGCAGCCA‑GCGGCGCCTGGGTGCCGAGCTGGTCATCCTCGGCCCGCTGACGCCTGGCGGCGGCTTCGGcc                                                                                              >  1:1128714/1‑150 (MQ=255)
                                                   tttACAGACGGGGCTCGACCCGGCCTGCGCCACGCGCAGCGGCTGGCGGCAAGTCGTCACGGTAACGGTCGCTGCGCTCGGCCGCCA‑GCGGCGCCTGGGTGCCGAGCTGGTCATCCTCGGCCCGCTGACGCCTGGCGGCGGCTTCGGCCt                                                                                             >  1:667617/1‑150 (MQ=255)
                                                    ttACAGACGGGGCTCGACCCGGCCTGCGCCACGCGCAGCGGCTGGCGGCAAGTCGTCACGGTAACGGTCGCTGCGCTCGGCAGCCA‑GCGGCGCCTGGGTGCCGAGCTGGTCATCCTCGGCCCGCTGACGCCTGGCGGCGGCTTCGGCCTg                                                                                            <  2:1421538/150‑1 (MQ=255)
                                                         ccgatctcTCGACCCGGCCTGCGCCACGCGCAGCGGCTGGCGGCAAGTCGTCACGGTAACGGTCGCTGCGCTCGGCAGCCA‑GCGGCGCCTGGGTGCCGAGCTGGTCATCCTCGGCCCGCTGACGCCTGGCGGCGGCTTCGGCCTGCTCAc                                                                                       <  2:850927/143‑1 (MQ=255)
                                                          aCGGGGCTCGACCCGGCCTGCGCCACGCGCAGCGGCTGGCGGCAAGTCGTCACGGTAACGGTCGCTGCGCTCGGCAGCCA‑GCGGCGCCTGGGTGCCGAGCTGGTCATCCTCGGCCCGCTGACGCCTGGCGGCGGCTTCGGCCTGCTCACg                                                                                      <  1:1368157/150‑1 (MQ=255)
                                                          aCGGGGCTCGACCCGGCCTGCGCCACGCGCAGCGGCTGGCGGCAAGTCGTCACGGTAACGGTCGCTGCGCTCGGCAGCCA‑GCGGCGCCTGGGTGCCGAGCTGGTCATCCTCGGCCCGCTGACGCCTGGCGGCGGCTTCGGCCTGCTCACg                                                                                      <  2:1134125/150‑1 (MQ=255)
                                                            ggggCTCGACCCGGCCTGCGCCACGCGCAGCGGCTGGCGGCAAGTCGTCACGGTAACGGTCGCTGCGCTCGGCAGCCA‑GCGGCGCCTGGGTGCCGAGCTGGTCATCCTCGGCCCGCTGACGCCTGGCGGCGGCTTCGGCCTAATCACggg                                                                                    >  2:37708/1‑149 (MQ=255)
                                                              ggCTCGACCCGGCCTGCGCCACGCGCAGCGGCTGGCGGCAAGTCGTCACGGTAACGGTCGCTGCGCTCGGCAGCCA‑GCGGCGCCTGGGTGCCGAGCTGGTCATCCTCGGCCCGCTGACGCCTGGCGGCGGCTTCGGCCTGCTCACGgcgt                                                                                  <  2:24310/150‑1 (MQ=255)
                                                                cTCGACCCGGCCTGCGCCACGCGCAGCGGCTGGCGGCAAGTCGTCACGGTAACGGTCGCTGCGCTCGGCAGCCA‑GCGGCGCCTGGGTGCCGAGCTGGTCATCCTCGGCCCGCTGACGCCTGGCGGCGGCTTCGGCCTGCTCACGgcgtgc                                                                                <  1:1292858/150‑1 (MQ=255)
                                                                 tCGACCCGGCCTGCGCCACGCGCAGCGGCTGGCGGCAAGTCGTCACGGTAACGGTCGCTGCGCTCGGCAGCCA‑GCGGCGCCTGGGTGCCGAGCTGGTCATCCTCGGCCCGCTGACGCCTGGCGGCGGCTTCGGCCTGCTCACGgcgtgcg                                                                               >  1:1269635/1‑150 (MQ=255)
                                                                 tCGACCCGGCCTGCGCCACGCGCAGCGGCTGGCGGCAAGTAGTCAAGGGAACGGTCGCTGCGCTCGGCAGCCA‑GCGGCGCCTGGGTGCCGAGCTGGTCATCCTCGGCCCGCTGACGCCTGGCGGCGGCTTCGGCCTGCTCACgagatcgg                                                                               >  1:850927/1‑143 (MQ=255)
                                                                    aCCCGGCCTGCGCCACGCGCAGCGGCTGGCGGCAAGTCGTCACGGTAACGGTCGCTGCGCTCGGCAGCCA‑GCGGCGCCTGGGTGCCGAGCTGGTCATCCTCGGCCCGCTGACGCCTGGCGGCGGCTTCGGCCTGCTCACGgcgtgcgtgc                                                                            <  2:1228644/150‑1 (MQ=255)
                                                                             gcgcCACGCGCAGCGGCTGGCGGCAAGTCGTCACGGTAACGGTCGCTGCGCTCGGCAGCCA‑GCGGCGCCTGGGTGCCGAGCTGGTCATCCTCGGCCCGCTTACGCCTGGCGGCGGCTTCGGCCTGCTCACGGCGTGCGTGCGTGCGCttt                                                                   <  1:248909/150‑1 (MQ=255)
                                                                             gcgcCACGCGCAGCGGCTGGCGGCAAGTCGTCACGGTAACGGTCGCTGCGCTCGGCAGCCA‑GCGGCGCCTGGGTGCCGAGCTGGTCATCCTCGGCCCGCTGACGCCTGGCGGCGGCTTCGGCCTGCTCACGGCGTGCGTGCGTGCGCttt                                                                   <  2:287189/150‑1 (MQ=255)
                                                                                   cgcgCAGCGGCTGGCGGCAAGTCGTCACGGTAACGGTCGCTGCGCTCGGCAGCCA‑GCGGCGCCTGGGTGCCGAGCTGGTCATCCTCGGCCCGCTGACGCCTGGCGGCGGCTTCGGCCTGCTCACGGCGTGCGTGCGTGCGCTTTTCCAGc                                                             <  2:1199533/150‑1 (MQ=255)
                                                                                     cgcAGCGGCTGGCGGCAAGTCGTCACGGTAACGGTCGCTGCGCTCGGCAGCCA‑GCGGCGCCTGGGTGCCGAGCTGGTCATCCTCGGCCCGCTGACGCCTGGCGGCGGCTTCGGCCTGCTCACGGCGTGCGTGCGTGCGCTTTTCCAGCgg                                                           <  1:422799/150‑1 (MQ=255)
                                                                                            gCTGGCGGCAAGTCGTCACGGTAACGGTCGCTGCGCTCGGCAGCCA‑GCGGCGCCTGGGTGCCGAGCTGGTCATCCTCGGCCCGCTGCCGCCTGGCGGCGGCTTCGGCCTGCTCACGGCGTGCGTGCGTGCGCTTTTCCAGCGGCCAgcgc                                                    >  2:389474/1‑150 (MQ=255)
                                                                                            gCTGGCGGCAAGTCGTCACGGTAACGGTCGCTGCGCTCGGCAGCCA‑GCGGCGCCTGGGTGCCGAGCTGGTCATCCTCGGCCCGCTGACGCCTGGCGGCGGCTTCGGCCTGCTCACGGCGTGCGTGCGTGCGCTTTTCCAGCGGCCAgcgc                                                    <  1:389474/150‑1 (MQ=255)
                                                                                                gcggcAAGTCGTCACGGTAACGGTCGCTGCGCTCGGCAGCCA‑GCGGCGCCTGGGTGCCGAGCTGGTCATCCTCGGCCCGCTGACGCCTGGCGGCGGCTTCGGc                                                                                               <  2:382034/103‑1 (MQ=255)
                                                                                                gcggcAAGTCGTCACGGTAACGGTCGCTGCGCTCGGCAGCCA‑GCGGCGCCTGGGTGCCGAGCTGGTCATCCTCGGCCCGCTGACGCCTGGCGGCGGCTTCGGc                                                                                               >  1:382034/1‑103 (MQ=255)
                                                                                                   gcAAGTCGTCACGGTAACGGTCGCTGCGCTCGGCAGCCA‑GCGGCGCCTGGGTGCCGAGCTGGTCATCCTCGGCCAGCTGCACCCAGTCGGCGGCTTCGGCCTGCTCACGGCGTGCGTGCGTGCGCTTTTCCAGCCGCCAGCgcagcagca                                             >  1:1518542/1‑150 (MQ=255)
                                                                                                     aaGTCGTCACGGTAACGGTCGCTGCGCTCGGCAGCCA‑GCGGCGCCTGGGTGCCGAGCTGGTCATCCTCGGCCCGCTGACGCCTGGCGGCGGCTTCGGCCTGCTCACGGCGTGCGTGCGTGCGCTTTTCCAGCGGCCAGCGCAGCAGCACg                                           >  2:1450984/1‑150 (MQ=255)
                                                                                                       gtcgtcACGGTAACGGTCGCTGCGCTCGGCAGCCA‑GCGGCGCCTGGGTACCGAGCTGGTCATCCTCGGCCCGCTGACGCCTGGCGGCGGCTTCGGCCTGCTCACGGCGTGCGTGCGTGCGCTTTTCCAGCGGCCAGCGCAGCGGCACGaa                                         >  1:631689/1‑150 (MQ=255)
                                                                                                        tcgtcACGGTAACGGTCGCTGCGCTCGGCAGCCA‑GCGGCGCCTGGGTGCCGAGCTGGTCATCCTCGGCCCGCTGACGCCTGGCGGCGGCTTCGGCCTGCTCACGGCGTGCGTGCGTGCGCTTTTCCAGCGGCCAGCGCAGCAGCACGAag                                        <  1:255125/150‑1 (MQ=255)
                                                                                                         cgtcACGGTAACGGTCGCTGCGCTCGGCAGCCA‑GCGGCGCCTGGGTGCCGAGCTGGTCATCCTCGGCCCGCTGACGCCTGGCGGCGGCTTCGGCCTGCTCACGGCGTGCGTGCGTGCGCTTTTCCAGCGGCCAGCGCAGCAGCACGAaga                                       <  2:1266116/150‑1 (MQ=255)
                                                                                                         cgtcACGGTAACGGTCGCTGCGCTCGGCAGCCA‑GCGGCGCCTGGGTGCCGAGCTGGTCATCCTCGGCCCGCTGACGCCTGGCGGCGGCTTCGGCCTGCTCACGGCGTGCGTGCGTGCGCTTTTCCAGCGGCCAGCGCAGCAGCACGAaga                                       <  1:819512/150‑1 (MQ=255)
                                                                                                         cgtcACGGTAACGGTCGCTGCGCTCGGCAGCCA‑GCGGCGCCTGGGTGCCGAGCTGGTCATCCTCGGCCCGCTGACGCCTGGCGGCGGCTTCGGCCTGCTCACGGCGTGCGTGCGTGCGCTTTTCCAGCGGCCAGCGCAGCAGCACGAaga                                       <  1:708421/150‑1 (MQ=255)
                                                                                                          gtcACGGTAACGGTCGCTGCGCTCGGCAGCCA‑GCGGCGCCTGGGTGCCGAGCTGGTCATCCTCGGCCCGCTGACGCCTGGCGGCGGCTTCGGCCTGCTCACGgcgtgcgtgcgt                                                                          >  2:1189316/1‑114 (MQ=255)
                                                                                                          gtcACGGTAACGGTCGCTGCGCTCGGCAGCCA‑GCGGCGCCTGGGTGCCGAGCTGGTCATCCTCGGCCCGCTGACGCCTGGCGGCGGCTTCGGCCTGCTCACGgcgtgcgtgcgt                                                                          <  1:1189316/114‑1 (MQ=255)
                                                                                                          gtcACGGTAACGGTCGCTGCGCTCGGCAGCCA‑GCGGCGCCTGGGTGCCGAGCTGGTCATCCTCGGCCCGCTGACGCCTGGCGGCGGCTTCGGCCTGCTCACGGCGTGCGTGCGTGCGCTTTTCCAGCGGCCAGCGCAGCAGCACGAagag                                      <  1:789590/150‑1 (MQ=255)
                                                                                                           tcACGGTAACGGTCGCTGCGCTCGGCAGCCA‑GCGGCGCCTGGGTGCCGAGCTGGTCATCCTCGGCCCGCt                                                                                                                     <  2:1492633/70‑1 (MQ=255)
                                                                                                           tcACGGTAACGGTCGCTGCGCTCGGCAGCCA‑GCGGCGCCTGGGTGCCGAGCTGGTCATCCTCGGCCCGCt                                                                                                                     >  1:1492633/1‑70 (MQ=255)
                                                                                                             aCGGTAACGGTCGCTGCGCTCGGCAGCCA‑GCGGCGCCTGGGTGCCGAGCTGGTCATCCTCGGCCCGCTGACGCCTGGCGGCGGCTTCGGCCTGCTCACGGCGTGCGTGCGTGCGCTTTTCCAGCGGCCAGCGCAGCAGCACGAAGAGGaa                                   <  2:1128714/150‑1 (MQ=255)
                                                                                                                   aCGGTCGCTGCGCTCGGCAGCCA‑GCGGCGCCTGGGTGCCGAGCTGGTCATCCTCGGCCCGCTGACGCCTGGCGGCGGCTTCGGCCTGCTCACGGCGTGCGTGCGTGCGCTTTTCCAGCGGCCAGCGCAGCAGCACGAAGAGGAAGTACAg                             >  2:1415192/1‑150 (MQ=255)
                                                                                                                      gTCGCTGCGCTCGGCAGCCA‑GCGGCGCCTGGGTGCCGAGCTGGTCATCCTCGGCCCGCTGACGCCTGGCGGCGGCTTCGGCCTGCTCACGGCGTGCGTGCGTGCGCTTTTCCAGCGGCCAGCGCAGCAGCACGAAGAGGAAGTACAGGcc                          >  2:1154171/1‑150 (MQ=255)
                                                                                                                           ccagcTCGGCACCCAGGCGCCG‑CTGGCTGCCGAGCTGGTCATCCTCGGCCCGCTGACGCCTGGCGGCGGCTTCGGCCTGCTCACGGCGTGCGTGCGTGCGCTTTTCCAGCGGCCAGCGCAGCAGCACGAAGAGGAAGTACAGGCCGAAgg                     >  1:1281552/4‑150 (MQ=255)
                                                                                                                            cagcTCGGCACCCAGGCGCCG‑CTGGCTGCCGAGCTGGTCATCCTCGGCCCGCTGACGCCTGGCGGCGGCTTCGGCCTGCTCACGGCGTGCGTGCGTGCGCTTTTCCAGCGGCCAGCGCAGCAGCACGAAGAGGAAATACGGGCAGAAGGc                    >  2:1325493/3‑150 (MQ=255)
                                                                                                                            cagcTCGGCACCCA‑GGCGCCGCTGGCTGCCGAGCTGGTCATCCTCGGCCCGCTGACGCCTGGCGGCGGCTTCGGCCTGCTCACGGCGTGCGTGCGTGCGCTTTTCCAGCGGCCAGCGCAGCAGCACGAAGAGGAAGTACAGGCCGAAGGc                    >  1:3587/3‑150 (MQ=255)
                                                                                                                             cgcTCGGCAGCCA‑GCGGCGCCTGGGTGCCGAGCTGGTCATCCTCGGCCCGCTGACGCCTGGCGGCGGCTTCGGCCTGCTCACGGCGTGCGTGCGTGCGCGTTTCCAGCGGCCAGCGCAGCAGCACGAAGAGGAAGTACAGGCCGAAGGCg                   <  1:37708/150‑1 (MQ=255)
                                                                                                                              gcTCGGCAGCCA‑GCGGCGCCTGGGTGCCGAGCTGGTCCTTCTCGGCCCGCTGACGCCTGGCGGCGGCTTCGGCCTGCTCACGGCGTGCGTGCGTGCGCTTTTCCAGCGGCCAGCGCAGCAGCACGAAGAGGAAGTACAGGCCGAAGGCGa                  <  1:428883/150‑1 (MQ=255)
                                                                                                                              gcTCGGCAGCCA‑GCGGCGCCTGGGTGCCGAGCTGGTCATCCTCGGCCCGCTGACGCCTGGCGGCGGCTTCGGCCTGCTCACGgcgtgc                                                                                >  2:441577/1‑88 (MQ=255)
                                                                                                                              gcTCGGCAGCCA‑GCGGCGCCTGGGTGCCGAGCTGGTCATCCTCGGCCCGCTGACGCCTGGCGGCGGCTTCGGCCTGCTCACGgcgtgc                                                                                <  1:441577/88‑1 (MQ=255)
                                                                                                                                tCGGCAGCCA‑GCGGCGCCTGGGTGCCGAGCTGGTCATCCTCGGCCCGCTGACGCCTGGCGGCGGCTTCGGCCTGCTCACGGCGTGCGTGCGTGCGCTTTTCCAGCGGCCAGCGCAGCAGCACGAAGAGGAAGTACAGGCCGAAGGCGATc                <  2:653345/150‑1 (MQ=255)
                                                                                                                                   gCAGCCA‑GCGGCGCCTGGGTGCCGAGCTGGTCATCCTCGGCCCGCTGACGCCTGGCGGCGGCTTCGGCCTGCTCACGGCGTGCGTGCGTGCGCTTTTCCAGCGGCCAGCGCAGCAGCACGAAGAGGAAGTACAGGCCGAAGGCGATCATg             >  2:306630/1‑150 (MQ=255)
                                                                                                                                     agcca‑gcGGCGCCTGGGTGCCGAGCTGGTCATCCTCGGCCCGCTGACGCCTGGCGGCGGCTTCGGCCTGCTCACGGCGTGCGTGCGTGCGCTTTTCCAGCGGCCAGCGCAGCAGCACGAAGAGGAAGTACAGGCCGAAGGCGATCATGcc           <  2:1518542/150‑1 (MQ=255)
                                                                                                                                           ggcgccgcTGGCTGCCGAGCTGGTCATCCTCGGCCCGCTGACGCCTGGCGGCGGCTTCGGCCTGCTCACGGCGTGCGTGCGTGCGCTTTTCCAGCGGCCAGCGCAGCAGCACGAAGAGGAAGTACAGGCCGAAGGCGATCATGCCGTACa      <  1:1325493/143‑1 (MQ=255)
                                                                                                                                              gcgcCTGGGTGCCGAGCTGGTCATCCTCGGCCCGCTGACGCCTGGCGGCGGCTTCGGCCTGCTCACGGCGTGCGTGCGTGCGCTTTTCCAGCGGCCAGCGCAGCAGCACGAAGAGGAAGTACAGGCCGAAGGCGATCATGCCGTACATgg   <  2:1646074/150‑1 (MQ=255)
                                                                                                                                               cgcCTGGGTGCCGAGCTGGTCATCCTCGGCCCGCTGACGCCTGGCGGCGGCTTCGGCCTGCTCACGGCGTGCGTGCGTGCGCTTTTCCAGCGGCCAGCGCAGCAGCACGAAGAGGAAGTACAGGCCGAAGGCGATCATGCCGTACATGGc  <  1:1084559/150‑1 (MQ=255)
                                                                                                                                                      |                                                                                                                                              
GTAGGGACGTTGGCACTCATGCCCAGCAGTTCACACATAGGCGTGTCTCCGTTTACAGACGGGGCTCGACCCGGCCTGCGCCACGCGCAGCGGCTGGCGGCAAGTCGTCACGGTAACGGTCGCTGCGCTCGGCAGCCA‑GCGGCGCCTGGGTGCCGAGCTGGTCATCCTCGGCCCGCTGACGCCTGGCGGCGGCTTCGGCCTGCTCACGGCGTGCGTGCGTGCGCTTTTCCAGCGGCCAGCGCAGCAGCACGAAGAGGAAGTACAGGCCGAAGGCGATCATGCCGTACATGGC  >  NC_002947/5201013‑5201304

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: ‑