Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,002,525 G→C 5.6% G189A (GGT→GCT)  PP_2626 → hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029473,002,5250GC 5.6% 229.1 / 2.2 71G189A (GGT→GCT) PP_2626hypothetical protein
Reads supporting (aligned to +/- strand):  ref base G (33/34);  new base C (2/2);  total (35/36)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.41e-01

GATCGTTTTACCTACACCCCGACCCACGCAGCGGCAGCGCAGGCTGCCTTGAAACGCCTGCAACAACAGGAAACAGTACAGCTGGACACGGCAGCACTGGGCGAGCTTGAGCGCCTGCTAGGCAAGTACTCGGAACTTGCCAACGGTGTTGGCGAGGCT‑GCTGTGCCCAAAGCGCCGGCTGCAACCGAGAAGCCTGCACCGTTGGCCGATGTCATCGCCACACCCGTGTCACGCATGTCGATGGGCAATGAGCGTGAGACACGCCGGGCCGTGCTGTCCATGGCCGAGCTCAT  >  NC_002947/3002380‑3002672
                                                                                                                                                 |                                                                                                                                                    
gATCGTTTTACCTACACCCCGACCCACGCAGCGGCAGCGCAGGCTGCCTTGAAACGCCTGCAACAACAGGAAACAGTACAGCTGGACACGGCAGCACTGGGCGAGCTTGAGCGCCTGCTAGGCAAGTACTCGGAACTTGCCAACGgt                                                                                                                                                     <  5:77025/147‑1 (MQ=255)
           cTACACCCCGACCCACGCAGCGGCAGCGCAGGCTGCCTTGAAACGCCTGCAACAACAGGAAACAGTACAGCTGGACACGGCAGCACTGGGCGAGCTTGAGCGCCTGCTAGGCAAGTACTCGGAACTTGCCAACGGTGTTGGCGAGgct‑g                                                                                                                                       >  5:134202/1‑149 (MQ=255)
             acacCCCGACCCACGCAGCGGCAGCGCAGGCTGCCTTGAAACGCCTGCAACAACAGGAAACAGTACAGCTGGACACGGCAGCACTGGGCGAGCTTGAGCGCCTGCTAGGCAAGTACTCGGAACTTGCCAACGGTGTTGGCGAGgct‑gct                                                                                                                                     >  1:50033/1‑149 (MQ=255)
                 cccGACCCACGCAGCGGCAGCGCAGGCTGCCTTGAAACGCCTGCAACAACAGGAAACAGTACAGCTGGACACGGCAGCACTGGGCGAGCTTGAGCGCCTGCTAGGCAAGTACTCGGAACTTGCCAACGGTGTTGGCGAGGCT‑GCtgt                                                                                                                                   <  5:222679/147‑1 (MQ=255)
                 cccGACCCACGCAGCGGCAGCGCAGGCTGCCTTGAAACGCCTGCAACAACAGGAAACAGTACAGCTGGACACGGCAGCACTGGGCGAGCTTGAGCGCCTGCTAGGCAAGTACTCGGAACTTGCCAACGGTGTTGGCGAGGCT‑GCtgt                                                                                                                                   >  6:222679/1‑147 (MQ=255)
                 cccGACCCACGCAGCGGCAGCGCAGGCTGCCTTGAAACGCCTGCAACAACAGGAAACAGTACAGCTGGACACGGCAGCACTGGGCGAGCTTGAGCGCCTGCTAGGCAAGTACTCGGAACTTGCCAACGGTGTTGGCGAGGCT‑GCTGTGc                                                                                                                                 >  3:212222/1‑149 (MQ=255)
                         aCGCAGCGGCAGCGCAGGCTGCCTTGAAACGCCTGCAACAACAGGAAACAGTACAGCTGGACACGGCAGCACTGGGCGAGCTTGAGCGCCTGCTAGGCAAGTACTCGGAACTTGCCAACGGTGTTGGCGAGGCT‑GCTGTGc                                                                                                                                 <  3:245935/141‑1 (MQ=255)
                         aCGCAGCGGCAGCGCAGGCTGCCTTGAAACGCCTGCAACAACAGGAAACAGTACAGCTGGACACGGCAGCACTGGGCGAGCTTGAGCGCCTGCTAGGCAAGTACTCGGAACTTGCCAACGGTGTTGGCGAGGCT‑GCTGTGc                                                                                                                                 >  4:245935/1‑141 (MQ=255)
                            cagcggcagcgCAGGCTGCCTTGAAACGCCTGCAACAACAGGAAACAGTACAGCTGGACACGGCAGCACTGGGCGAGCTTGAGCGCCTGCTAGGCAAGTACTCGGAACTTGCCAACGGTGTTGGCGAGGCT‑GCTGTGCCCAAAGCGCCg                                                                                                                      <  3:129050/149‑1 (MQ=255)
                               cggcagcgCAGGCTGCCTTGAAACGCCTGCAACAACAGGAAACAGTACAGCTGGACACGGCAGCACTGGGCGAGCTTGAGCGCCTGCTAGGCAAGTACTCGGAACTTGCCAACGGTGTTGGCGAGGCT‑GCTGTGCCCAAAGCGCCGGCt                                                                                                                   <  7:220037/149‑1 (MQ=255)
                                ggcagcgCAGGCTGCCTTGAAACGCCTGCAACAACAGGAAACAGTACAGCTGGACACGGCAGCACTGGGCGAGCTTGAGCGCCTGCTAGGCAAGTACTCGGAACTTGCCAACGGTGTTGGCGAGGCT‑GCTGTGCCCAAAGCGCCGGCTg                                                                                                                  >  7:336771/1‑149 (MQ=255)
                                ggcagcgCAGGCTGCCTTGAAACGCCTGCAACAACAGGAAACAGTACAGCTGGACACGGCAGCACTGGGCGAGCTTGAGCGCCTGCTAGGCAAGTACTCGGAACTTGCCAACGGTGTTGGCGAGGCT‑GCTGTGCCCAAAGCGCCGGCTg                                                                                                                  >  7:167638/1‑149 (MQ=255)
                                                  gAAACGCCTGCAACAACAGGAAACAGTACAGCTGGACACGGCAGCACTGGGCGAGCTTGAGCGCCTGCTAGGCAAGTACTCGGAACTTGCCAACGGTGTTGGCGAGGCT‑GCTGTGCCCAAAGCGCCGGCTGCAACCGAGAAGCCTGCa                                                                                                 <  6:134202/148‑1 (MQ=255)
                                                     aCGCCTGCAACAACAGGAAACAGTACAGCTGGACACGGCAGCACTGGGCGAGCTTGAGCGCCTGCTAGGCAAGTACTCGGAACTTGCCAACACCGTTGGCAAgttcc                                                                                                                                        >  2:138886/1‑103 (MQ=255)
                                                     aCGCCTGCAACAACAGGAAACAGTACAGCTGGACACGGCAGCACTGGGCGAGCTTGAGCGCCTGCTAGGCAAGTACTCGGAACTTGCCAACACCGTTGGCAAgttcc                                                                                                                                        <  1:138886/107‑5 (MQ=255)
                                                     aCGCCTGCAACAACAGGAAACAGTACAGCTGGACACGGCAGCACTGGGCGAGCTTGAGCGCCTGCTAGGCAAGTACTCGGAACTTGCCAACACCGTTGGCAAgttcc                                                                                                                                        <  7:196687/107‑5 (MQ=255)
                                                     aCGCCTGCAACAACAGGAAACAGTACAGCTGGACACGGCAGCACTGGGCGAGCTTGAGCGCCTGCTAGGCAAGTACTCGGAACTTGCCAACACCGTTGGCAAgttcc                                                                                                                                        >  8:196687/1‑103 (MQ=255)
                                                           gCAACAACAGGAAACAGTACAGCTGGACACGGCAGCACTGGGCGAGCTTGAGCGCCTGCTAGGCAAGTACTCGGAACTTGCCAACGGTGTTGGCGAGGCT‑GCTGTGCCCAAAGCGCCGGCTGCAACCGAGAAGCCTGCACCGTTGGcc                                                                                        >  1:274746/1‑148 (MQ=255)
                                                             aacaacAGGAAACAGTACAGCTGGACACGGCAGCACTGGGCGAGCTTGAGCGCCTGCTAGGCAAGTCCTCGGAACTTGCCAACGGTGTTGGCGAGGCT‑GCTGTGCCCAAAGCGCCGGCt                                                                                                                   <  1:93881/119‑1 (MQ=255)
                                                             aacaacAGGAAACAGTACAGCTGGACACGGCAGCACTGGGCGAGCTTGAGCGCCTGCTAGGCAAGTACTCGGAACTTGCCAACGGTGTTGGCGAGGCT‑GCTGTGCCCAAAGCGCCGGCt                                                                                                                   >  2:93881/1‑119 (MQ=255)
                                                                aacaGGAAACAGTACAGCTGGACACGGCAGCACTGGGCGAGCTTGAGCGCCTGCTAGGCAAGTACTCGGAACTTGCCAACGGGGTTGGCGAGGCT‑GCTGTGCCCAAAGCGCCGGCTGCAACCGAGAAGCCTGCACCGTTGGCCGATGt                                                                                   >  2:240410/1‑148 (MQ=255)
                                                                 acaGGAAACAGTACAGCTGGACACGGCAGCACTGGGCGAGCTTGAGCGCCTGCTAGGCAAGTACTCGGAACTTGCCAACGGTGTTGGCGAGGCT‑GCTGTGCCCAAAGCGCCGGCTGCAACCGAGAAGCCTGCACCGTTGGCCGATGTCa                                                                                 >  3:292021/1‑149 (MQ=255)
                                                                  caGGAAACAGTACAGCTGGACACGGCAGCACTGGGCGAGCTTGAGCGCCTGCTAGGCAAGTACTCGGAACTTGCCAACGGTGTTGGCGAGGCT‑GCTGTGCCCAAAGCGCCGGCTGCAACCGAGAAGCCTGCACCGTTGGCCGATGTCAt                                                                                <  1:241862/149‑1 (MQ=255)
                                                                    ggAAACAGTACAGCTGGACACGGCAGCACTGGGCGAGCTTGAGCGCCTGCCAGGCAAGTACTCGGAACTTGCCAACGGTGTTGGCGAGGCT‑GCTGTGCCCAAAGCGCCGGCTGCAAcc                                                                                                             <  2:144525/118‑1 (MQ=255)
                                                                    ggAAACAGTACAGCTGGACACGGCAGCACTGGGCGAGCTTGAGCGCCTGCCAGGCAAGTACTCGGAACTTGCCAACGGTGTTGGCGAGGCT‑GCTGTGCCCAAAGCGCCGGCTGCAAcc                                                                                                             >  1:144525/1‑118 (MQ=255)
                                                                     gAAACAGTACAGCTGGACACGGCAGCACTGGGCGAGCTTGAGCGCCTGCTAGGCAAGTACTCGGAACTTGCCAACGGTGTTGGCGAGGCT‑GCTGTGCCCAAAGCGCCGGCTGCAACCGAGAAGCCTGCACCGTTGGCCGATGTCATCGc                                                                             <  2:50033/149‑1 (MQ=255)
                                                                     gAAACAGTACAGCTGGACACGGCAGCACTGGGCGAGCTTGAGCGCCTGCTAGGCAAGTACTCGGAACTTGCCAACGGTGTTGGCGAGGCT‑GCTGTGCCCAAAGCGCCGGCTGCAACCGAGAAGCCTGCACCGTTGGCCGATGTCATCGc                                                                             >  4:227086/1‑149 (MQ=255)
                                                                       aaCAGTACAGCTGGACACGGCAGCACTGGGCGAGCTTGAGCGCCTGCTAGGCAAGTACTCGGAACTTGCCAACGGTGTTGGCGAGGCT‑GCTGTGCCCAAAGCGCCGGCTGCAACCGAGAAGCCTGCACCGTTGGCCGATGTCATCGCca                                                                           <  4:29129/149‑1 (MQ=255)
                                                                          aGTACAGCTGGACACGGCAGCACTGGGCGAGCTTGAGCGCCTGCTAGGCAAGTACTCGGAACTTGCCAACGGTGTTGGCGAGGCT‑GCTGTGCCCAAAGCGCCGGCTGCAACCGAGAAGCCTGCACCGTTGGCCGATGTCATCGCCacac                                                                        <  6:277285/149‑1 (MQ=255)
                                                                               aGCTTGACACGGCAGCACTGGGCGAGCTTGAGCGCCTGCTAGGCAAGTACTCGGAACTTGCCAACGGTGTTGGGGTGGCTAGCTGTGCCCAAAGCGCCGGCTGCAACCGAGAAGCCTg                                                                                                   <  1:63973/118‑1 (MQ=255)
                                                                               aGCTGGACACGGCAGCACTGGGCGAGCTTGAGCGCCTGCTAGGCAAGTACTCGGAACTTGCCAACGGTGTTGGCGAGGCT‑GCTGTGCCCAAAGCGCCGGCTGCAACCGAGAAGCCTg                                                                                                   >  2:63973/1‑117 (MQ=255)
                                                                                 cTGGACACGGCAGCACTGGGCGAGCTTGAGCGCCTGCTAGGCAAGTACTCGGAACTTGCCAACGGTGTTGGCGAGGCT‑GCTGTGCCCAAAGCGCCGGCTGCAACCGAGAAGCCTGCACCGTTGGCCGATGTCATCGCCACACCCgtgt                                                                  >  6:38345/1‑148 (MQ=255)
                                                                                   ggACACGGCAGCACTGGGCGAGCTTGAGCGCCTGCTAGGCAAGTACTCGGAACTTGCCAACCGTGTTGGCGAGGCT‑GCTGTGCCCAAAGCGCCGGCTGCAACCGAGAAGCCTGCACCGTTGGCCGATGTCATCGCCACACCCGTGTCAc                                                               >  2:251987/1‑149 (MQ=255)
                                                                                    gACACGGCAGCACTGGGCGAGCTTGAGCGCCTGCTAGGCAAGTACTCGGAACTTGCCAACGGTGTTGGCGAGGCT‑GCTGTGCCCAAAGCGCCGGCTGCAACCGAGAAGCCTGCACCGTTGGCCGATGTCATCGCCACACCCGTGTCAc                                                               <  8:167638/148‑1 (MQ=255)
                                                                                         ggCAGCACTGGGCGAGCTTGAGCGCCTGCTAGGCAAGTACTCGGAACTTGCCAACGGTGTTGGCGAGGCT‑GTTGTGCCCAAAGCGCCGGCTGCAACCGAGAAGCCTGCACCGTTGGCCGAt                                                                                     <  6:283283/121‑1 (MQ=255)
                                                                                         ggCAGCACTGGGCGAGCTTGAGCGCCTGCTAGGCAAGTACTCGGAACTTGCCAACGGTGTTGGCGAGGCT‑GTTGTGCCCAAAGCGCCGGCTGCAACCGAGAAGCCTGCACCGTTGGCCGAt                                                                                     >  5:283283/1‑121 (MQ=255)
                                                                                         ggCAGCACTGGGCGAGCTTGAGCGCCTGCTAGGCAAGTACTCGGAACTTGCCAACGGTGTTGGCGAGGCT‑GCTGTGCCCAAAgcgc                                                                                                                        <  6:334283/86‑1 (MQ=255)
                                                                                         ggCAGCACTGGGCGAGCTTGAGCGCCTGCTAGGCAAGTACTCGGAACTTGCCAACGGTGTTGGCGAGGCT‑GCTGTGCCCAAAgcgc                                                                                                                        >  5:334283/1‑86 (MQ=255)
                                                                                          gcagcaCTGGGCGAGCTTGAGCGCCTGCTAGGCAAGTACTCGGAACTTGCCAACGGTGTTGGCGAGGCT‑GCTGTGCCCAAAGCGCCGGCTGCAACCGAGAAGCCTGCACCGTTGGCCGATGTCATCGCCACACCCGTGTCACGCATGTc                                                        >  1:280370/1‑149 (MQ=255)
                                                                                              caCTGGGCGAGCTTGAGCGCCTGCTAGGCAAGTACTCGGAACTTGCCAACGGTGTTGGCGAGGCT‑GCTGTGCCCAAAGCGCCGACTGCAACCGAGAAGCCTGCACCGTTGGCCGATGTCATCGCCACACCCGTGTCACGCATGTCGATg                                                    <  4:212222/149‑1 (MQ=255)
                                                                                               aCTGGGCGAGCTTGAGCGCCTGCTAGGCAAGTACTCGGAACTTGCCAACGGTGTTGGCGAGGCT‑GCTGTGCCCAAAGCGCCGGCTGCAACCGAGAAGCCTGCACCGTTGGCCGATGTCATCGCCACACCCGTGTCACGCATGTCGATgg                                                   <  4:292021/149‑1 (MQ=255)
                                                                                                 tGGGCGAGCTTGAGCGCCTGCTAGGCAAGTACTCGGAACTTGCCAACGGTGTTGGCGAGGCT‑GCTGTGCCCAAAGCGCCGGCTGCAACCGAGAAGCCTGCACCGTTGGCCGATGTCATCGCCACACCCGTGTCACGCATGTCGATgg                                                   <  5:279747/147‑1 (MQ=255)
                                                                                                 tGGGCGAGCTTGAGCGCCTGCTAGGCAAGTACTCGGAACTTGCCAACGGTGTTGGCGAGGCT‑GCTGTGCCCAAAGCGCCGGCTGCAACCGAGAAGCCTGCACCGTTGGCCGATGTCATCGCCACACCCGTGTCACGCATGTCGATgg                                                   >  6:279747/1‑147 (MQ=255)
                                                                                                  gggCGAGCTTGAGCGCCTGCTAGGCAAGTACTCGGAACTTGCCAACGGTGTTGGCGAGGCT‑GCTGTGCCCAAAGCGCCGGCTGCAACCGAGAAGCCTGCACCGTTGGCCGATGTCATCGCCACACCCGTGTCACGCATGTCGATGGGc                                                 >  7:100263/1‑148 (MQ=255)
                                                                                                         cTTGAGCGCCTGCTAGGCAAGTACTCGGAACTTGCCAACGGTGTTGGCGAGGCT‑GCTGTGCCCAAAGCGCCGGCTGCAACCGAGAAGCCTGCACCGTTgg                                                                                          >  3:281010/1‑100 (MQ=255)
                                                                                                         cTTGAGCGCCTGCTAGGCAAGTACTCGGAACTTGCCAACGGTGTTGGCGAGGCT‑GCTGTGCCCAAAGCGCCGGCTGCAACCGAGAAGCCTGCACCGTTgg                                                                                          <  4:281010/100‑1 (MQ=255)
                                                                                                             aGCGCCTGCTAGGCAAGTACTCGGAACTTGCCAACGGTGTTGGCGAGGCT‑GCTGTGCCCAAAGCGCCGGCTGCAACCGAGAAGCCTGCACCGTTGGCCGATGTCATCGCCACACCCGTGTCACGCATGTCGATGGGCAATGAGCGTga                                      <  2:274746/148‑1 (MQ=255)
                                                                                                             aGCGCCTGCTAGGCAAGTACTCGGAACTTGCCAACGGTGTTGGCGAGGCT‑GCTGTGCCCAAAGCGCCGGCTGCAACCGAGAAGCCTGCACCGTTGGCCGATGTCATCGCCACACCCGTGTCACGCATGTCGATGGGCAATGAGCGTga                                      <  5:38345/148‑1 (MQ=255)
                                                                                                                 ccTGCTAGGCAAGTACTCGGAACTTGCCAACGGTGTTGGCGAGGCT‑GCTGTGCCCAAAGCGCCGGCTGCAACCGAGAAGCCTGCACCGTTGGCCGATGTCATCGCcaca                                                                         >  6:209646/1‑109 (MQ=255)
                                                                                                                 ccTGCTAGGCAAGTACTCGGAACTTGCCAACGGTGTTGGCGAGGCT‑GCTGTGCCCAAAGCGCCGGCTGCAACCGAGAAGCCTGCACCGTTGGCCGATGTCATCGCcaca                                                                         <  5:209646/109‑1 (MQ=255)
                                                                                                                  cTGCTAGGCAAGTACTCGGAACTTGCCAACGGTGTTGGCGAGGCT‑GCTGTGCCCAAAGCGCCGGCTGCAACCGAGAAGCCTGCACCGTTGGCCGATGTCATCGCCACACCCGTGTCACGCATGTCGATGGGCAAt                                              <  1:147221/135‑1 (MQ=255)
                                                                                                                  cTGCTAGGCAAGTACTCGGAACTTGCCAACGGTGTTGGCGAGGCT‑GCTGTGCCCAAAGCGCCGGCTGCAACCGAGAAGCCTGCACCGTTGGCCGATGTCATCGCCACACCCGTGTCACGCATGTCGATGGGCAAt                                              >  2:147221/1‑135 (MQ=255)
                                                                                                                         gCAAGTACTCGGAACTTGCCAACGGTGTTGGCGAGGCT‑GCTGTGCCCAAAGCGCCGGCTGCAACCGAGAAGCCTGCACCGTTGGCCGATGTCATCGCCACACCCGTGTCACGCATGTCGATGGGCAATGAGCGTGAGACACGCCGGGcc                         >  7:76182/1‑149 (MQ=255)
                                                                                                                         gCAAGTACTCGGAACTTGCCAACGGTGTTGGCGAGGCT‑GCTGTGCCCAAAGCGCCGGCTGCAACCGAGAAGCCTGCACCGTTGGCCGATGTCATCGCCACACCCGTGTCACGCATGTCGATGGGCAATGAGCGTGAGACACGCCGGGcc                         >  7:266246/1‑149 (MQ=255)
                                                                                                                           aaGTACTCGGAACTTGCCAACGGTGTTGGCGAGGCT‑GCTGTGCCCAAAGCGCCGGCTGCAACCGAGAAGCCTGCACCGTTGGCCGATGTCATCGCCACACCCGTGTCACGCATGTCGATGGGCAATGAGCGTGAGACACGCCGGGCCGt                       >  4:168382/1‑149 (MQ=255)
                                                                                                                           aaGTACTCGGAACTTGCCAACCGTGTTGGCGAGGCT‑GCTGTGCCCAAAGCGCCGGCTGCAACCGAGAAGCCTGCACCGTTGGCCGATGTCATCGCCACACCCGAGTCACGCATGTCGATAGGCAATGAGCGTGAGACACGCCGGGCCGt                       <  1:251987/149‑1 (MQ=255)
                                                                                                                            aGTACTCGGAACTTGCCAACGGTGTTGGCGAGGCT‑GCTGTGCCCAAAGCGCCGGCTGCAACCGAGAAGCCTGCACCGTTGGCCGATGTCATCGCCACACCCGTGTCACGCATGTCGATGGGCAATGAGCGTGAGACACGCCGGGCCGTg                      <  8:336771/149‑1 (MQ=255)
                                                                                                                               aCTCGGAACTTGCCAACGGTGTTGGCGAGGCT‑GCTGTGCCCAAAGCGCCGGCTGCAACCGAGAAGCCTGCACCGTTGGCCGATGTCATCGCCACACCCGTGTCACGCATGTCGATGGGCAATGAGCGTGAGACACGCCGGGCCGTGCt                    <  3:227086/148‑1 (MQ=255)
                                                                                                                                cTCGGAACTTGCCAACGGTGTTGGCGAGGCT‑GCTGGGCCCAAAGCGCCGGCTGCAACCGAGAAGCCTGCACCGTTGGCCGATGTCATCGCCACACCCGTGTCACGCATGTCGATGGGCAATGAGCGTGAGACACGCCGGGCCGTGCTGt                  <  7:327088/149‑1 (MQ=255)
                                                                                                                                 tCGGAACTTGCCAACGGTGTTGGCGAGGCT‑GCTGTGCCCAAAGCGCCGGCTGCAACCGAGAAGCCTGCACCGTTGGCCGATGTCATCGCCACACCCGTGTCACGCATGTCGATGGGCAATGAGCGTGAGACACGCCGGGCCGTGCTGTc                 >  7:155350/1‑149 (MQ=255)
                                                                                                                                 tCGGAACTTGCCAACGGTGTTGGCGAGGCT‑GCTGTGCCCAAAGCGCCGGCTGCAACCGAGAAGCCTGCACCGTTGGCCGATGTCATCGCCACACCCGTGTCACGCATGTCGATGGGCAATGAGCGTGAGACACGCCGGGCCGTGCTGTc                 >  6:313414/1‑149 (MQ=255)
                                                                                                                                 tCGGAACTTGCCAACGGTGTTGGCGAGGCT‑GCTGTGCCCAAAGCGCCGGCTGCAACCGAGAAGCCTGCACCGTTGGCCGATGTCATCGCCACACCCGTGTCACGCATGTCGATGGGCAATGAGCGTGAGACACGCCGGGCCGTGCTGTc                 >  4:327574/1‑149 (MQ=255)
                                                                                                                                   ggAACTTGCCAACGGTGTTGGCGAGGCT‑GCTGTGCCCAAAGCGCCGGCTGCAACCGAGAAGCCTGCACCGTTGGCCGATGTCATCGCCACACCCGTGTCAc                                                               <  3:151880/101‑1 (MQ=255)
                                                                                                                                   ggAACTTGCCAACGGTGTTGGCGAGGCT‑GCTGTGCCCAAAGCGCCGGCTGCAACCGAGAAGCCTGCACCGTTGGCCGATGTCATCGCCACACCCGTGTCAc                                                               >  4:151880/1‑101 (MQ=255)
                                                                                                                                          gCCAACGGTGTTGGCGAGGCT‑GCTGTGCCCAAAGCGCCGGCTGCAACCGAGAAGCCTGCACCGTTGGCCGATGTCATCGCCACACCCGTGTCACGCATGTCGATGGGCAATGAGCGTGAGACACGCCGGGCCGTGCTGTCCATGGcc          <  2:280370/147‑1 (MQ=255)
                                                                                                                                          gCCAACGGTGTTGGCGAGGCT‑GCTGTGCCCAAAGCGCCGGCTGCAACCGAGAAGCCTGCACCGTTGGCCGATGTCATCGCCACACCCGTGTCACGCATGTCGATGGGCAATGAGCGTGAGACACGCCGGGCCGTGCTGTCCATGGcc          >  7:46836/1‑147 (MQ=255)
                                                                                                                                          gCCAACGGTGTTGGCGAGGCT‑GCTGTGCCCAAAGCGCCGGCTGCAACCGAGAAGCCTGCACCGTTGGCCGATGTCATCGCCACACCCGTGTCACGCATGTCGATGGGCAATGAGCGTGAGACACGCCGGGCCGTGCTGTCCATGGcc          <  8:266246/147‑1 (MQ=255)
                                                                                                                                          gCCAACGGTGTTGGCGAGGCT‑GCTGTGCCCAAAGCGCCGGCTGCAACCGAGAAGCCTGCACCGTTGGCCGATGTCATCGCCACACCCGTGTCACGCATGTCGATGGGCAATGAGCGTGAGACACGCCGGGCCGTGCTGTCCATGGcc          <  8:76182/147‑1 (MQ=255)
                                                                                                                                            cAACGGTGTTGGCGAGGCT‑GCTGTGCCCAAAGCGCCGGCTGCAACCGAGAAGCCTGCACCGTTGGCCGATGTCATCGCCACACCCGTGTCACGCATGTCGATGGGCAATGAGCGTGAGACACGCCGGGCCGTGCTGTCCATGGCCGAGc      >  3:72208/1‑149 (MQ=255)
                                                                                                                                              aCGGTGTTGGCGAGGCT‑GCTGTGCCCAAAGCGCCGGCTGCAACCGAGAAGCCTGCACCGTTGGCCGATGTCATCGCCACACCCGTGTCACGCATGTCGATGGGCAATGAGCGTGAGACACGCCGGGCCGTGCTGTCCATGGCCGAGCtc    <  8:100263/149‑1 (MQ=255)
                                                                                                                                              aCGGTGTTGGCGAGGCT‑GCTGTGCCCAAAGCGCCGGCTGCAACCGAGAAGCCTGCACCGTTGGCCGATGTCATCGCCACACCCGTGTCACGCATGTCGATGGGCAATGAGCGTGAGACACGCCGGGCCGTGCTGTCCATGGCCGAGCtc    <  8:46836/149‑1 (MQ=255)
                                                                                                                                               cGGTGTTGGCGAGGCT‑GCTGTGCCCAAAGCGCCGGCTGCAACCGAGAAGCCTGCACCGTTGGCCGATGTCATCGCCACACCCGTGTCACGCATGTCGATGGGCAATGAGCGTGAGACACGCCGGGCCGTGCTGTCCATGGCCGAGCtca   <  4:72208/149‑1 (MQ=255)
                                                                                                                                                ggTGTTGGCGAGGCT‑GCTGTGCCCAAAGCGCCGGCTGCAACCGAGAAGCCTGCACCGTTGGCCGATGTCATCGCCACACCCGTGTCACGCATGTCGATGGGCAATGAGCGTGAGACACGCCGGGCCGTGCTGTCCATGGCCGAGCtcat  <  1:240410/149‑1 (MQ=255)
                                                                                                                                                 |                                                                                                                                                    
GATCGTTTTACCTACACCCCGACCCACGCAGCGGCAGCGCAGGCTGCCTTGAAACGCCTGCAACAACAGGAAACAGTACAGCTGGACACGGCAGCACTGGGCGAGCTTGAGCGCCTGCTAGGCAAGTACTCGGAACTTGCCAACGGTGTTGGCGAGGCT‑GCTGTGCCCAAAGCGCCGGCTGCAACCGAGAAGCCTGCACCGTTGGCCGATGTCATCGCCACACCCGTGTCACGCATGTCGATGGGCAATGAGCGTGAGACACGCCGGGCCGTGCTGTCCATGGCCGAGCTCAT  >  NC_002947/3002380‑3002672

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