Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,843,873 C→A 7.8% L241M (CTG→ATG) ‡ PP_4259 → iron‑sulfur cluster‑binding protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,843,8730CA 7.8% 203.8 / 6.0 64L241M (CTG→ATG) ‡PP_4259iron‑sulfur cluster‑binding protein
Reads supporting (aligned to +/- strand):  ref base C (34/25);  new base A (3/2);  total (37/27)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.72e-01

CACGGGCCAGGCCGATGGCTGGGCCTACTTCTGGGTCGGCTTCTTCACCCTCGCCACTTACGGCAACGCCGGCTGGCTACGCGAACAGGTGTGCGTGTACATGTGCCCCTATGCCCGCTTCCAGAGCGTGATGTTCGACAAGGACACGCTGATCGTGTCCTACGACCCACGCCGTGGCGAAACCCGCGGCCCGCGCAAGAAAGATACCGACTACAAGGCCCAGGGCCTGGGTGACTGCATCGACTGCACCATGTGCGTGCAGGTCTGCCCGACCGGCATCGACATCCGTGACGGCCT  >  NC_002947/4843725‑4844021
                                                                                                                                                    |                                                                                                                                                      
cACGGGCCAGGCCGATGGCTGGGCCTACTTCTGGGTCGGCTTCTTCACCCTCGCCACTTACGGCAACGCCGGCTGGCTACGCGAACAGGTGTGCGTGTACATGTGCCCCTATGCCCGCTTCCAGAGCGTGATGTTCGACAAGGACACGc                                                                                                                                                        >  5:268995/1‑149 (MQ=255)
 aCGGGCCAGGCCGATGGCTGGGCCTACTTCTGGGTCGGCTTCTTCACCCTCGCCACTTACGGCAACGCCGGCTGGCTACGCGAACAGGTGTGCGTGTACATGTGCCCCTATGCCCGCTTCCAGAGCGTGATGTTCGACAAGGACACGc                                                                                                                                                        <  2:229877/148‑1 (MQ=255)
     gCCAGGCCGATGGCTGGGCCTACTTCTGGGTCGGCTTCTTCACCCTCGCCACTTACGGCAACGCCGGCTGGCTACGCGAACAGGTGTGCGTGTACATGTGCCCCTATGCCCGCTTCCAGAGCGTGATGTTCGACAAGGACACGCTGATc                                                                                                                                                   >  7:239141/1‑149 (MQ=255)
      ccAGGCCGATGGCTGGGCCTACTTCTGGGTCGGCTTCTTCACCCTCGCCACTTACGGCAACGCCGGCTGGCTACGCGAACAGGTGTGCGTGTACATGTGCCCCTATGCCCGCTTCCAGAGCGTGATGTTCGACAAGGACACGCTGATCg                                                                                                                                                  >  7:249603/1‑149 (MQ=255)
              aTGGCTGGGCCTACTTCTGGGTCGGCTTCTTCACCCTCGCCACTTACGGCAACGCCGGCTGGCTACGCGAACAGGTGTGCGTGTACATGTGCCCCTATGCCCGCTTCCAGAGCGTGATGTTCGACAAGGACACGCTGATCGTGTCCTAc                                                                                                                                         <  8:156948/149‑1 (MQ=255)
               tGGCTGGGCCTACTTCTGGGTCGGCTTCTTCACCCTCGCCACTTACGGCAACGCCGGCTGGCTACGCGAACAGGTGTGCGTGTACATGTGCCCCTATGCCCGCTTCCAGAGCGTGATGTTCGACAAGGACACGATCAGCGTGTCCttgt                                                                                                                                        <  2:219679/149‑4 (MQ=255)
                 gCTTGGCCTACTTCTGGGTCGGCTTCTTCACCCTCGCCACTTACGGCAACGCCGGCTGGCTACGCGAACAGGTGTGCGTGTACATGTGCCCCTATGCCCGCTTCCAGAGCGTTATGTTCGACAAGTAGACGCTGATCGTGTCCTACGa                                                                                                                                       <  4:200206/148‑1 (MQ=255)
                  cTGGGCCTACTTCTGGGTCGGCTTCTTCACCCTCGCCACTTACGGCAACGCCGGCTGGCTACGCGAACAGGTGTGCGTGTACATGTGCCCCTATGCCCGCTTCCAGAGCGTGATGTTCGACAAGGACACGCTGATCGTGTCCTACGAcc                                                                                                                                     >  2:288641/1‑149 (MQ=255)
                        cTACTTCTGGGTCGGCTTCTTCACCCTCGCCACTTACGGCAACGCCGGCTGGCTACGCGAACAGGTGTGCGTGTACATGTGCCCCTATGCCCGCTTCCAGAGCGTGATGTTCGACAAGGACACGCTGATCGTGTCCTACGACCCAc                                                                                                                                  >  4:174937/1‑146 (MQ=255)
                        cTACTTCTGGGTCGGCTTCTTCACCCTCGCCACTTACGGCAACGCCGGCTGGCTACGCGAACAGGTGTGCGTGTACATGTGCCCCTATGCCCGCTTCCAGAGCGTGATGTTCGACAAGGACACGCTGATCGTGTCCTACGACCCAc                                                                                                                                  <  3:174937/146‑1 (MQ=255)
                         tACTTCTGGGTCGGCTTGTTCACCCTCGCCACTTACGGCAACGCCGGCTGGCTACGCGAACAGGTGTGCGTGTTCATGTGCCCCTATGCCCGCTTCCAGAGCGTGATGTTCGACAAGGACACGCTGATCGTGTCCTACGACCCACGCCg                                                                                                                              >  8:39761/1‑149 (MQ=255)
                              cTGGGTCGGCTTCTTCACCCTCGCCACTTACGGCAACGCCGGCTGGCTACGCGAACAGGTGTGCGTGTACATGTGCCCCTATGCCCGCTTCCAGAGCGTGATGTTCGACAAGGACACGCTGATCGTGTCCTACGACCCACGCCGTgg                                                                                                                           <  1:288641/147‑1 (MQ=255)
                                      gCTTCTTCACCCTCGCCACTTACGGCAACGCCGGCTGGCTACGCGAACAGGTGTGCGTGTACATGTGCCCCTATGCCCGCTTCCAGAGCGTGATGTTCGACAAGGACACGCTGATCGTGTCCTACGACCCACGCCGTGGCGAAACCcgc                                                                                                                 >  7:170881/1‑149 (MQ=255)
                                        ttcttcACCCTCGCCACTTACGGCAACGCCGGCTGGCTACGCGAACAGGTGTGCGTGTACATGTGCCCCTATGCCCGCTTCCAGAGCGTGATGTTCGACAAGGACACGCTGATCGTGTCCTACGACCCACGCCGTGGCGAAACCCGCgg                                                                                                               >  7:234615/1‑149 (MQ=255)
                                          cttcACCCTCGCCACTTACGGCAACGCCGGCTGGCTACGCGAACAGGTGTGCGTGTACATGTGCCCCTATGCCCGCTTCCAGAGCGTGATGTTCGACAAGGACACGCTGATCGTGTCCTACGACCCACGCCGTGGCGAAACCCGCgg                                                                                                               <  1:50649/147‑1 (MQ=255)
                                          cttcACCCTCGCCACTTACGGCAACGCCGGCTGGCTACGCGAACAGGTGTGCGTGTACATGTGCCCCTATGCCCGCTTCCAGAGCGTGATGTTCGACAAGGACACGCTGATCGTGTCCTACGACCCACGCCGTGGCGAAACCCGCGGcc                                                                                                             >  7:223957/1‑149 (MQ=255)
                                          cttcACCCTCGCCACTTACGGCAACGCCGGCTGGCTACGCGAACAGGTGTGCGTGTACATGTGCCCCTATGCCCGCTTCCAGAGCGTGATGTTCGACAAGGACACGCTGATCGTGTCCTACGACCCACGCCGTGGCGAAACCCGCGGc                                                                                                              <  7:61636/148‑1 (MQ=255)
                                              aCCCTCGCCACTTACGGCAACGCCGGCTGGCTACGCGAACAGGTGTGCGTGTACATGTGCCCCTATGCCCGCTTCCAGAGCGTGATGTTCGACAAGGACACGCTGATCGTGTCCTAc                                                                                                                                         >  6:203187/1‑117 (MQ=255)
                                              aCCCTCGCCACTTACGGCAACGCCGGCTGGCTACGCGAACAGGTGTGCGTGTACATGTGCCCCTATGCCCGCTTCCAGAGCGTGATGTTCGACAAGGACACGCTGATCGTGTCCTAc                                                                                                                                         <  5:203187/117‑1 (MQ=255)
                                               cccTCGCCACTTACGGCAACGCCGGCTGGCTACGCGAACAGGTGTGCGTGTACATGTGCCCCTATGCCCGCTTCCAGAGCGTGATGTTCGACAAGGACACGCTGATCGTGTCCTACGACCCACGCCGTGGCGAAACCCGCGGCCcgcg                                                                                                         >  3:371374/1‑148 (MQ=255)
                                                 cTCGCCACTTACGGCAACGCCGGCTGGCTACGCGAACAGGTGTGCGTGTACATGTGCCCCTATGCCCGCTTCCAGAGCGTGATGTTCGACAAGGACACGCTGATCGTGTCCTACGACCCACGCCGTGGCGAAACCCGCGGCCCgcgc                                                                                                        >  3:28941/1‑147 (MQ=255)
                                                 cTCGCCACTTACGGCAACGCCGGCTGGCTACGCGAACAGGTGTGCGTGTACATGTGCCCCTATGCCCGCTTCCAGAGCGTGATGTTCGACAAGGACACGCTGATCGTGTCCTACGACCCACGCCGTGGCGAAACCCGCGGCCCGCGCaa                                                                                                      <  6:268995/149‑1 (MQ=255)
                                                    gCCACTTACGGCAACGCCGGCTGGCTACGCGAACAGGTGTGCGTGTACATGTGCCCCTATGCCCGCTTCCAGAGCGTGATGTTCGACAAGGACACGCTGATCGTGTCCTACGACCCACGCCGTGGCGAAACCCGCGGCCCGCGCaaga                                                                                                    >  6:354176/1‑148 (MQ=255)
                                                           aCGGCAACGCCGGCTGGCTACGCGAACAGGTGTGCGTGTACATGTGCCCCTATGCCCGCTTCCAGAGCGTGATGTTCGACAAGGACACGCTGATCGTGTCCTTACGACCCACGCCGTGGCGAAACCCGCGGCCCGCGCAAGAAAGAt                                                                                               <  1:80130/147‑1 (MQ=255)
                                                                 gcGCCGGCTGGCTACGCGAACAGGTGTGCGTGTACATGTGCCCCTATGCCCGCTTCCAGAGCGTGATGTTCGACAAGGACACGCTGATCGTGTCCTACGACCCACGCCGTGGCGAAACCCGCGGCCCGCGCAAGAAAGATACCGACTa                                                                                       >  3:363450/2‑148 (MQ=255)
                                                                   gCCGGCTGGCTACGCGAACTGGTGTGCGTGTACATGTGCCCCTATGCCCGCTTCCAGAGCGTGATGTTCGACAAGGACACGCTGATCGTGTCCTACGACCCACGCCGTGGCGAAACCCGCGGCCCGCGCAAGAAAGATACCGACTACaa                                                                                   <  8:249603/149‑1 (MQ=255)
                                                                           gctACGCGAACAGGTGTGCGTGTACATGTGCCCCTATGCCCGCTTCCAGAGCGTGATGTTCGACAAGGACACGCTGATCGTGTCCTACGACCCACGCCGTGGCGAAACCCGCGGCCCGCGCAAGAAAGATACCGACTACAAGGCCCAgg                                                                           >  5:98159/1‑149 (MQ=255)
                                                                             tACGCGAACAGGTGTGCGTGTACATGTGCCCCTATGCCCGCTTCCAGAGCGTGATGTTCGACAAGGACACGCTGATCGTGTCCTACGACCCACGCCGTGGCGAAACCCGCGGCCCGCGCAAGAAAGATACCGACTACAAGGCCCAGGGc                                                                         >  2:252771/1‑149 (MQ=255)
                                                                                 cgAACAGGTGTGCGTGTACATGTGCCCCTATGCCCGCTTCCAGAGCGTGATGTTCGACAAGGACACGCTGATCGTGTCCTACGACCCACGCCGTGGCGAAACCCGCGGCCCGCGCAAGAAAGATACCGACTACAAGGCCCAGGGCCTgg                                                                     <  8:234615/149‑1 (MQ=255)
                                                                                  gAACAGGTGTGCGTGTACATGTGCCCCTATGCCCGCTTCCAGAGCGTGATGTTCGACAAGGACACGCTGATCGTGTCCTACGACCCACGCCGTGGCGAAACCCGCGGCCCGCGCAAGAAAGATACCGACTACAAGGCCCAGGGCCTggg                                                                    <  8:303832/149‑1 (MQ=255)
                                                                                  gAACAGGTGTGCGTGTACATGTGCCCCTATGCCCGCTTCCAGAGCGTGATGTTCGACAAGGACACGCTGATCGTGTCCTACGACCCACGCCGTGGCGAAACCCGCGGCCCGCGCAAGAAAGATACCGACTACAAGGCCCAGGGCCTggg                                                                    >  5:357015/1‑149 (MQ=255)
                                                                                    aCAGGTGTGCGTGTACATGTGCCCCTATGCCCGCTTCCAGAGCGTGATGTTCGACAAGGACACGCTGATCGTGTCCTACGACCCACGCCGTGGCGAAACCCGCGGCCCGCGCAAGAAAGATACCGACTACAAGGCCCAGGGCCTGGGt                                                                   >  7:112615/1‑148 (MQ=255)
                                                                                    aCAGGTGTGCGTGTACATGTGCCCCTATGCCCGCTTCCAGAGCGTGATGTTCGACAAGGACACGCTGATCGTGTCCTACGACCCACGCCGTGGCGAAACCCGCGGCCCGCGCAAGAAAGATACCGACTACAAGGCCCAGGGCCTGGGt                                                                   >  3:261438/1‑148 (MQ=255)
                                                                                       tgggtgCGTGTACATGTGCCCCTATTCCCGCTTCCAGAGAGTGATGTTCGACAAGGACACGCTGATCGTGTCCTACGACCCACGCCGTGGCGAAACCCGCGGCCCGCGCAAGATAGATACCGACTACAAGGCCCAGGGCCTGGGTGACt                                                               <  3:224897/146‑1 (MQ=255)
                                                                                           tgCGTGTACATGTGCCCCTATGCCCGCTTCCAGAGCGTGAGGTTCGACAAGGACACGCTGCTCGTGTCCTACGCCCCACGCCGTGGCGAACCCCGCGGCACGCGCAAGAAAGATACCGACTACAAGGCCCCGGGCCTGGGTGACTGCAt                                                           >  5:16397/1‑149 (MQ=255)
                                                                                              gtgtACATGTGCCCCTATGCCCGCTTCCAGAGCGTGATGTTCGACAAGGACACGCTGATCGTGTCCTACGACCCACGCCGTGGCGAAACCCGCGGCCCGCGCAAGAAAGATACCGACTACAAGGCCCAGGGCCTGGGTGACTGCATCGa                                                        <  5:192857/149‑1 (MQ=255)
                                                                                              gtgtACATGTGCCCCTATGCCCGCTTCCAGAGCGTGATGTTCGACAAGGACACGCTGATCGTGTCCTACGACACACGCCGTGGCGAAACCCGCGGCCCGCGCAAGAAAGATACCGACTACAAGGCCCAGGGCCTGGGTGACTGCATCGa                                                        >  5:380402/1‑149 (MQ=255)
                                                                                                   cATGTGCCCCTATGCCCGCTTCCAGAGCGTGATGTTCGACAAGGACACGCTGATCGTGTCCTACGACCCACGCCGTGGCGAAACCCGCGGCCCGCGCAAGAAAGATACCGACTACAAGGCCCAGGGCCTGGGTGACTGCATCGACTGc                                                    <  8:223957/148‑1 (MQ=255)
                                                                                                   cATGTGCCCCTATGCCCGCTTCCAGAGCGTGATGTTCGACAAGGACACGCTGATCGTGTCCTACGACCCACGCCGTGGCGAAACCCGCGGCCCGCGCAAGAAAGATACCGACTACAAGGCCCAGGGCCTGGGTGACTGCATCGACTGCa                                                   <  1:252771/149‑1 (MQ=255)
                                                                                                      gtgCCCCTATGCCCGCTTCCAGAGCGTGATGTTCGACAAGGACACGCTGATCGTGTCCTACGACCCACGCCGTGGCGAAACCCGCGGCCCGCGCAAGAAAGATACCGACTACAAGGCCCAGGGCCTGGGTGACTGCATCGACTGCACCa                                                >  5:296880/1‑149 (MQ=255)
                                                                                                      gtgCCCCTATGCCCGCTTCCAGAGCGTGATGTTCGACAAGGACACGCTGATCGTGTCCTACGACCCACGCCGTGGCGAAACCCGCGGCCCGCGCAAGAAAGATACCGACTACAAGGCCCAGGGCCTGGGTGACTGCATCGACTGCACCa                                                >  6:313034/1‑149 (MQ=255)
                                                                                                         ccccTATGCCCGCTTCCAGAGCGTGATGTTCGACAAGGACACGCTGATCGTGTCCTACGACCCACGCCGTGGCGAAACCCGCGGCCCGCGCAAGAAAGATACCGACTACAAGGCCCAGGGCCTGGGTGACTGCATCGACTGCACCAtg                                              <  5:187185/148‑1 (MQ=255)
                                                                                                                 cccGCTTCCAGAGCGTGATGTTCGACAAGGACACGCTGATCGTGTCCTACGACCCACGCCGTGGCGAAACCCGCGGCCCGCGCAAGAAAGATACCGACTACAAGGCCCAGGGCCTGGGTGACTGCATCGACTGCACCATGTGCGTGCAg                                     <  4:361668/149‑1 (MQ=255)
                                                                                                                  ccGCTTCCAGAGCGTGATGTTCGACAAGGACACGCTGATCGTGTCCTACGACCCACGCCGTGGCGAAACCCGCGGCCCGCGCAAGAAAGATACCGACTACAAGGCCCAGGGCCTGGGTGACTGCATCGACTGCACCATGTGCGTGCAgg                                    >  2:79179/1‑149 (MQ=255)
                                                                                                                   cggttgcAGAGCGTGATGTTCGATAAGGACACGCTGATCGTGTACTACGACCCACGCGGTGGCGTAAGCCGCGGCCCGTGCAATAATGATACCGACTACAAGGCCCAGGGCCTGGGTGTGTGCATGGACTGCAGGATGTGCGTGCAGGt                                   <  6:98159/143‑1 (MQ=255)
                                                                                                                   cGCTTCCAGAGCGTGATGTTCGACAAGGACACGCTGATCGTGTCCTCCGACCCACGCCGTGGCGAAACCCGCGGCCCGCGCAAGAAAGATACCGACTACAAGGCCCAGGGCCTGGGTGACTGCATCGACTGCACCATGTGCGTGCAGGt                                   >  5:299017/1‑149 (MQ=255)
                                                                                                                   cGCTTCCAGAGCGTGATGTTCGACAAGGACACGCTGATCGTGTCCTACGACCCACGCCGTGGCGAAACCCGCGGCCCGCGCAAGAAAGATACCGACTACAAGGCCCAGGGCCTGGGTGACTGCATCGACTGCGCCATGTGCGTGCAGGt                                   >  8:326939/1‑149 (MQ=255)
                                                                                                                          agagCGTGATGTTCGACAAGGACACGCTGATCGTGTCCTACGACCCACGCCGTGGCGAAACCCGCGGCCCGCGCAAGAAAGATACCGACTACAAGGCCCAGGGCCTGGGTGACTGCATCGACTGCACCATGTGCGTGCAGGTCTGccc                             >  1:94406/1‑148 (MQ=255)
                                                                                                                          agagCGTGATGTTCGACAAGGACACGCTGATCGTGTCCTACGACCCACGCCGTGGCGAAACCCGCGGCCCGCGCAAGAAAGATACCGACTACAAGGCCAAGGGCCTGGGTGACTGCATCGACTGCACCATGTGCGTGCAGGTCTGccc                             >  1:20957/1‑148 (MQ=255)
                                                                                                                              cGTGATGTTCGACAAGGACACGCTGATCGTGTCCTACGACCCACGCCGTGGCGAAACCCGCGGCCCGCGCAAGAAAGATACCGACTACAAGGCCCAGGGCCTGGGTGACTGCATCGACTGCACCATGTGCGTGCAGGTCTGCCCGACCg                        <  5:354176/149‑1 (MQ=255)
                                                                                                                                 gATGTTCGACAAGGACACGCTGATCGTGTCCTACGACCCACGCCGTGGCGAAACCCGCGGCCCGCGCAAGAAAGATACCGACTACAAGGCCCAGGGCCTGGGTGACTGCATCGACTGCACCATGTGCGTGCAGGTCTGCCCGACCGGCa                     >  4:55531/1‑149 (MQ=255)
                                                                                                                                  aTGTTCGACAAGGACACGCTGATCGTGTCCTACGACCCACGCCGTGGCGAAACCCGCGGCCCGCGCAAGAAAGATACCGACTACAAGGCCCAGGGCCTGGGTGACTGCATCGACTGCACCATGTGCGTGCAGGTCTGCCCGACCGGCAt                    <  6:357015/149‑1 (MQ=255)
                                                                                                                                    gTTCGACAAGGACACGCTGATCGTGTCCTACGACCCACGCCGTGGCGAAACCCGCGGCCCGCGCAAGAAAGATACCGACTACAAGGCCCAGGGCCTGGGTGACTGCATCGACTGCACCATGTGCGTGCAGGTCTGCCCGACCGGCATCg                  >  6:313989/1‑149 (MQ=255)
                                                                                                                                     ttCGACAAGGACACGCTGATCGTGTCCTACGACCCACGCCGTGGCGAAACCCGCGGCCCGCGCAAGAAAGATACCGACTACAAGGCCCAGGGCCTGGGTGACTGCATCGACTGCACCATGTGCGTGCAGGTCTGCCCGACCGGCATCGa                 <  1:79179/149‑1 (MQ=255)
                                                                                                                                       cGACAAGGACACGCTGATCGTGTCCTACGACCCACGCCGTGGCGAAACCCGCGGCCCGCGCAAGAAAGATACCGACTACAAGGCCCAGGGCCTGGGTGACTGCATCGACTGCACCATGTGCGTGCAGGTCTGCCCGACCGGCATCGACa               >  7:53165/1‑149 (MQ=255)
                                                                                                                                          cAAGGACACGCTGATCGTGTCCTACGACCCACGCCGTGGCGAAACCCGCGGCCCGCGCAAGAAAGATACCGACTACAAGGCCCAGGGCCTGGGTGACTGCATCGACTGCACCATGTGCGTGCAGGTCTGCCCGACCGGCATCGACATcc            <  4:28941/149‑1 (MQ=255)
                                                                                                                                          cAAGGACACGCTGATCGTGTCCTACGACCCACGCCGTGGCGAAACCCGCGGCCCGCGCAAGAAAGATACCGAAAAGCAAGGCCCAGGGCCTGGGTGACTGCATCGACTGCACCATGTGCGTGCAGGTCTGCCCGACCGGCATCGACATc             >  1:122011/1‑149 (MQ=255)
                                                                                                                                           taGGACACGATCAGCGTGTCCTACGACCCACGCCGTGGCGAAACCCGCGGCCCGCGCAAGAAAGATACCGACTACAAGGCCCAGGGCCTGGGTGACTGCATCGACTGCACCATGTGCGTGCAGGTCTGCCCGACCGGCATCGACATCCg           >  1:5177/2‑149 (MQ=255)
                                                                                                                                            aGGACACGATCAGCGTGTCCTACGACCCACGCCGTGGCGAAACCCGCGGCCCGCGCAAGAAAGATACCGACTACAAGGCCCAGGGCCTGGGTGACTGCATCGACTGCACCATGTGCGTGCAGGTCTGCCCGACCGGc                      <  4:115284/137‑1 (MQ=255)
                                                                                                                                            aGGACACGATCAGCGTGTCCTACGACCCACGCCGTGGCGAAACCCGCGGCCCGCGCAAGAAAGATACCGACTACAAGGCCCAGGGCCTGGGTGACTGCATCGACTGCACCATGTGCGTGCAGGTCTGCCCGACCGGc                      >  3:115284/1‑137 (MQ=255)
                                                                                                                                            aGGACACGATCAGCGTGTCCTACGACCCACGCCGTGGCGAAACCCGCGGCCCGCGCAAGAAAGATACCGACTACAAGGCCCAGGGCCTGGGTGACTGCATCGACTGCACCATGTGCGTGCAGGTCTGCCCGACCGGCATCGACATCCGt          >  5:189159/1‑149 (MQ=255)
                                                                                                                                             ggccacGCTGATCGTGTCCTACGACCCACGCCGTGGCGAAACCCGCGGCCCGCGCAAGAAAGATACCGACTACAAGGCCCAGGGCCTGGGTGACTGCATCGACTGCACCATGTGCGTGCAGGTCTGCCCGACCGGCATCGACATCCGt          >  6:305828/4‑148 (MQ=255)
                                                                                                                                                   gCTGATCGTGTCCTACGACCCACGCCGTGGCGAAACCCGCGGCCCGCGCAAGAAAGATACCGACTACAAGGCCCAGGGCCTGGGTGACTGCATCGACTGCACCATGTGCGTGCAGGTCTGCCCGACCGGCATCGACATCCGTGACGGcc   <  6:299017/149‑1 (MQ=255)
                                                                                                                                                   gCTGATCGTGTCCTACGACCCACGCCGTGGCGAAACCCGCGGCCCGCGCAAGAAAGATACCGACTACAAGGCCCAGGGCCTGGGTGACTGCATCGACTGCACCATGTGCGTGCAGGTCTGCCCGACCGGCATCGACATCCGTGACGGcc   >  2:24282/1‑149 (MQ=255)
                                                                                                                                                   gCTGATCGTGTCCTACGACCCACGCCGTGGCGAAACCCGCGGCCCGCGCAAGAAAGATACCGACTACAAGGCCCAGGGCCTGGGTGACTGCATCGACTGCACCATGTGCGTGCAGGTCTGCCCGACCGGCATCGACATCCGTGACGGc    <  8:226306/148‑1 (MQ=255)
                                                                                                                                                   gCTGATCGTGTCCTACGACCCACGCCGTGGCGAAACCCGCGGCCCGCGCAAGAAAGATACCGACTACAAGGCCCAGGGCCTGGGTGACTGCATCGACTGCACCATGTGCGTGCAGGTCTGCCCGACCGGCATCGACATCCGTGACGGc    >  7:226306/1‑148 (MQ=255)
                                                                                                                                                    cTGATCGTGTCCTACGACCCACGCCGTGGCGAAACCCGCGGCCCGCGCAAGAAAGATACCGACTACAAGGCCCAGGGCCTGGGTGACTGCATCGACTGCACCATGTGCGTGCAGGTCTGCCCGACCGGCATCGACATCCGTGACGGCCt  <  1:195809/149‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                      
CACGGGCCAGGCCGATGGCTGGGCCTACTTCTGGGTCGGCTTCTTCACCCTCGCCACTTACGGCAACGCCGGCTGGCTACGCGAACAGGTGTGCGTGTACATGTGCCCCTATGCCCGCTTCCAGAGCGTGATGTTCGACAAGGACACGCTGATCGTGTCCTACGACCCACGCCGTGGCGAAACCCGCGGCCCGCGCAAGAAAGATACCGACTACAAGGCCCAGGGCCTGGGTGACTGCATCGACTGCACCATGTGCGTGCAGGTCTGCCCGACCGGCATCGACATCCGTGACGGCCT  >  NC_002947/4843725‑4844021

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: