Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,595,751 C→A 6.0% G110C (GGT→TGT)  thiC ← phosphomethylpyrimidine synthase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,595,7510CA 6.0% 215.3 / 2.7 67G110C (GGT→TGT) thiCphosphomethylpyrimidine synthase
Reads supporting (aligned to +/- strand):  ref base C (33/30);  new base A (2/2);  total (35/32)
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.34e-01

GCCACGTATTCCATTTCGGCGGTGATGATGCCCTGACGGGCATAGTGCATCTGCGAGACGTTGGCGCCGGCCTTGGCCCGGCGCGGGTTGCGCACATGGGCAAAACGCAGCTTGGCCAGCTCTGCATCGTTCAGGCGCTGCTGACCGAAGTCGGAGCTCAGGCCTTCCAGGCGCTCGGTGTCACCGCGGGCTTCGATCCACGCCGAACGCACATCGGGCAGGCCTTTGCGCACGTCGATGATCACGTCGGGGTCGGTGTACGGGCCAGACGTGTCATACACCAGCACCGGGGCG  >  NC_002947/5595606‑5595899
                                                                                                                                                 |                                                                                                                                                    
tccACGTATTCCATTTCGGCGGTGATGATGCCCTGACGGGCATAGTGCATCTGCGAGACGTTGGCGCCGGCCTTGGCCCGGCGCGGGTTGCGCACATGGGCAAAACGCAGCTTGGCCAGCTCTGCATCGTTCAGGCGCTGCTGACCGaa                                                                                                                                                   <  7:327672/148‑1 (MQ=255)
 ccACGTATTCCATTTCGGCGGTGATGATGCCCTGACGGGCATAGTGCATCTGCGAGACGTTGGCGCCGGCCTTGGCCCGGCGCGGGTTGCGCACATGGGCAAAACGCAGCTTGGCCAGCTCTGCATCGTTCAGGCGCTGCTGACCGAAg                                                                                                                                                  >  3:296931/1‑149 (MQ=255)
  cACGTATTCCATTTCGGCGGTGATGATGCCCTGACGGGCATAGTGCATCTGCGAGACGTTGGCGCCGGCCTTGGCCCGGCGCGGGTTGCGCACATGGGCAAAACGCAGCTTGGCCAGCTCTGCATCGTTCAGGCGCTGCTGACCGaa                                                                                                                                                   >  8:327672/1‑147 (MQ=255)
     gTATTCCATTTCGGCGGTGATGATGCCCTGACGGGCATAGTGCATCTGCGAGACGTTGGCGCCGGCCTTGGCCCGGCGCGGGTTGCGCACATGGGCAAAACGCAGCTTGGCCAGCTCTGCATCGTTCAGGCGCTGCTGACCGAAGTCgg                                                                                                                                              <  2:235095/149‑1 (MQ=255)
     gTATTCCATTTCGGCGGTGATGATGCCCTGACGGGCATAGTGCATCTGCGAGACGTTGGCGCCGGCCTTGGCCCGGCGCGGGTTGCGCACATGGGCAAAACGCAGCTTGGCCAGCTCTGCATCGTTCAGGCGCTGCTGACCGAAGTCg                                                                                                                                               <  8:88078/148‑1 (MQ=255)
     gTATTCCATTTCGGCGGTGATGATGCCCTGACGGGCATAGTGCATCTGCGAGACGTTGGCGCCGGCCTTGGCCCGGCGCGGGTTGCGCACATGGGCAAAACGCAGCTTGGCCAGCTCTGCATCGTTCAGGCGCTGCTGACCGAAGTCg                                                                                                                                               >  7:88078/1‑148 (MQ=255)
                     gTGATGATGCCCTGACGGGCATAGTGCATCTGCGAGACGTTGGCGCCGGCCTTGGCCCGGCGCGGGTTGCGCACATGGGCAAAACGCAGCTTGGCCAGCTCTGCATCGTTCAGGCGCTGCTGACCGAAGTCGGAGCTCAGGCCTTCCAg                                                                                                                              >  5:221473/1‑149 (MQ=255)
                      tgatgaTGCCCTGACGGGCATAGTGCATCTGCGAGACGTTGGCGCCGGCCTTGGCCCGGCGCGGGTTGCGCACATGGGGAACACGCAGCTTGGCCAGCTCTGCATCGTTCAGGCGCTGCTGAGCGAAGTCGGAGCTCAGGCCTTCCAgg                                                                                                                             >  2:21204/1‑149 (MQ=255)
                          gatgCCCTGACGGGCATAGTGCATCTGCGAGACGTTGGCGCCGGCCTTGGCCCGGCGCGGGTTGCGCACATGGGCAAAACGCAGCTTGGCCAGCTCTGCATCGTTCAGGCGCTGCTGACCGAAGTCGGAg                                                                                                                                            >  4:321453/1‑130 (MQ=255)
                          gatgCCCTGACGGGCATAGTGCATCTGCGAGACGTTGGCGCCGGCCTTGGCCCGGCGCGGGTTGCGCACATGGGCAAAACGCAGCTTGGCCAGCTCTGCATCGTTCAGGCGCTGCTGACCGAAGTCGGAg                                                                                                                                            <  3:321453/130‑1 (MQ=255)
                          gatgCCCTGACGGGCATAGTGCATCTGCGAGACGTTGGCGCCGGCCTTGGCCCGGCGCGGGTTGCGCACATGGGCAAAACGCAGCTTGGCCAGCTCTGCATCGTTCAGGCGCTGCTGACCGAAGTCGGAGCTCAGGCCTTCCAGGCGCt                                                                                                                         >  3:312043/1‑149 (MQ=255)
                            tgCCCTGACGGGCATAGTGCATCTGCGAGACGTTGGCGCCGGCCTTGGCCCGGCGCGGGTTGCGCACATGGGCAAAACGCAGCTTGGCCAGCTCTGCATCGTTCAGGCGCTGCTGACCGAAGTCGGAGCTCa                                                                                                                                        >  4:300026/1‑132 (MQ=255)
                            tgCCCTGACGGGCATAGTGCATCTGCGAGACGTTGGCGCCGGCCTTGGCCCGGCGCGGGTTGCGCACATGGGCAAAACGCAGCTTGGCCAGCTCTGCATCGTTCAGGCGCTGCTGACCGAAGTCGGAGCTCa                                                                                                                                        <  3:300026/132‑1 (MQ=255)
                             gCCCTGACGGGCATAGTGCATCTGCGAGACGTTGGCGCCGGCCTTGGCCCGGCGCGGGTTGCGCACATGGGCAAAACGCAGCTTGGCCAGCTCTGCATCGTTCAGGCGCTGCTGACCGAAGTCGGAGCTCAGGCCTTCCAGGCGCTCgg                                                                                                                      >  7:197043/1‑149 (MQ=255)
                             gCCCTGACGGGCATAGTGCATCTGCGAGACGTTGGCGCCGGCCTTGGCCCGGCGCGGGTTGCGCACATGGGCAAAACGCAGCTTGGCCAGCTCTGCATCGTTCAGGCGCTGCTGACCGAAGTCGGAGCTCAGGCCTTCCAGGCGCTCg                                                                                                                       <  2:341621/148‑1 (MQ=255)
                              cccTGACGGGCATAGTGCATCTGCGAGACGTTGGCGCCGGCCTTGGCCCGGCGCGGGTTGCGCACATGGGCAAAACGCAGCTTGGCCAGCTCTGCATCGTTTAGGCGCTGCTGACCGAAGTCGGAGCTCAGGCCTTCCAGGCGCTCGgt                                                                                                                     <  1:77225/149‑1 (MQ=255)
                                 tGACGGGCATAGTGCATCTGCGAGACGTTGGCGCCGGCCTTGGCCCGGCGCGGGTTGCGCACATGGGCAAAACGCAGCTTGGCCAGCTCTGCATCGTTCAGGCGCTGCTGACCGAAGTCGGAGCTCAGGCCTTCCAGGCGCTCGgtgt                                                                                                                   >  7:115613/1‑148 (MQ=255)
                                     gggCATAGTGCATCTGCGAGACGTTGGCGCCGGCCTTGGCCCGGCGCGGGTTGCGCACATGGGCAAAACGCAGCTTGGCCAGCTCTGCATCGTTCAGGCGCTGCTGACCGAAGTCg                                                                                                                                               <  5:171921/116‑1 (MQ=255)
                                     gggCATAGTGCATCTGCGAGACGTTGGCGCCGGCCTTGGCCCGGCGCGGGTTGCGCACATGGGCAAAACGCAGCTTGGCCAGCTCTGCATCGTTCAGGCGCTGCTGACCGAAGTCg                                                                                                                                               >  6:171921/1‑116 (MQ=255)
                                            gTGCATCTGCGAGACGTTGGCGCCGGCCTTGGCCCGGCGCGGGTTGCGCACATGGGCAAAACGCAGCTTGGCCAGCTCTGCATCGTTCAGGCGCTGCTGACCGAAGTCGGAGCTCAGGCCTTCCAGGCGCTCGGTGTCACCGCGGGCtt                                                                                                       <  8:197043/149‑1 (MQ=255)
                                                  cTGCGAGACGTTGGCGCCGGCCTTGGCCCGGCGCGGGTTGCGCACATGGGCAAAACGCAGCTTGGCCAGCTCTGCATCGTTCAGGCGCTGCTGACCGAAGTCGGAGCTCAGGCCTTCCAGGCGCTCGGTGTCACCGCGGGCTTCGATcc                                                                                                 >  3:188898/1‑149 (MQ=255)
                                                           gTTGGCGCCGGCCTTGGCCCGGCGCGGGTTGCGCACATGGGCAAAACGCAGCTTGGCCAGCTCTGCATCGTTCAGGCGCTGCTGACCGAAGTCGGAGCTCAGGCCTTCCAGGCGCTCGGTGTCACCGCGGGCTTCGATCCAc                                                                                               <  8:200230/142‑1 (MQ=255)
                                                           gTTGGCGCCGGCCTTGGCCCGGCGCGGGTTGCGCACATGGGCAAAACGCAGCTTGGCCAGCTCTGCATCGTTCAGGCGCTGCTGACCGAAGTCGGAGCTCAGGCCTTCCAGGCGCTCGGTGTCACCGCGGGCTTCGATCCAc                                                                                               >  7:200230/1‑142 (MQ=255)
                                                           gTTGGCGCCGGCCTTGGCCCGGCGCGGGTTGCGCACATGGGCAAAACGCAGCTTGGCCAGCTCTGCATCGTTCAGGCGCTGCTGACCGAAGTCGGAGCTCAGGCCTTCCAGGCGCTCGGTGTCACCGCGGGCTTCGATCCACGCCGAAc                                                                                        >  4:58365/1‑149 (MQ=255)
                                                           gTTGGCGCCGGCCTTGGCCCGGCGCGGGTTGCGCACATGGGCAAAACGCAGCTTGGCCAGCTCTGCATCGTTCAGGCGCTGCTGACCGAAGTCGGAGCTCAGGCCTTCCAGGCGCTCGGTGTCACCGCGGGCTTCGATCCACGCCGAAc                                                                                        <  4:312043/149‑1 (MQ=255)
                                                           gTTGGCGCCGGCCTTGGCCCGGCGCGGGTTGCGCACATGGGCAAAACGCAGCTTGGCCAGCTCTGCATCGTTCAGGCGCTGCTGACCGAAGTCGGAGCTCAGGCCTTCCAGGCGCTCGGTGTCACCGCGGGCTTCGATCCACGCCGAAc                                                                                        >  8:141937/1‑149 (MQ=255)
                                                              ggCGCCGGCCTTGGCCCGGCGCGGGTTGCGCACATGGGCAAAACGCAGCTTGGCCAGCTCTGCATCGTTCAGGCGCTGCTGACCGAAGTCGGAGCTCAGGCCTTCCAGGCGCTCGGTGTCACCGCGGGCTTCGATCCACGCCGAAc                                                                                        >  6:359022/1‑146 (MQ=255)
                                                              ggCGCCGGCCTTGGCCCGGCGCGGGTTGCGCACATGGGCAAAACGCAGCTTGGCCAGCTCTGCATCGTTCAGGCGCTGCTGACCGAAGTCGGAGCTCAGGCCTTCCAGGCGCTCGGTGTCACCGCGGGCTTCGATCCACGCCGAAc                                                                                        <  5:359022/146‑1 (MQ=255)
                                                                cgcCGGCCTTGGCCCGGCGCGGGTTGCGCACATGGGCAAAACGCAGCTTGGCCAGCTCTGCATCGTTCAGGCGCTGCTGACCGAAGTCGGAGCTCAGGCCTTc                                                                                                                                 >  8:177916/1‑103 (MQ=255)
                                                                 gcCGGCCTTGGCCCGGCGCGGGTTGCGCACATGGGCAAAACGCAGCTTGGCCAGCTCTGCATCGTTCAGGCGCTGCTGACCGAAGTCGGAGCTCAGGCCTTCCAGGCGCTCGGTGTCACCGCGGGCTTCGATCCACGCCGAACGCACAt                                                                                  >  7:278114/1‑149 (MQ=255)
                                                                     gCCTTGGCCCGGCGCGGGTTGCGCACATGGGCAAAACGCAGCTTGGCCAGCTCTGCATCGTTCAGGCGCTGCTGACCGAAGTCGGAGCTCAGGCCTTCCAGGCGCTCGGTGTCACCGCGGGCTTCGATCCACGCCGAACGCACATCggg                                                                              <  8:162952/149‑1 (MQ=255)
                                                                      ccTTGGCCCGGCGCGGGTTGCGCACATGGGCAAAACGCAGCTTGGCCAGCTCTGCATCGTTCAGGCGCTGCTGACCGAAGTCGGAGCTCAGGCCTTCCAGGCGCTCGGTGTCACCGCGGGCTTCGATCCACGCCGAACGCACATCGGGc                                                                             <  4:188898/149‑1 (MQ=255)
                                                                      ccTTGGCCCGGCGCGGGTTGCGCACATGGGCAAAACGCAGCTTGGCCAGCTCTGCATCGTTCAGGCGCTGCTGACCGAAGTCGGAGCTCAGGCCTTCCAGGCGCTCGGTGTCACCGCGGGCTTCGATCCACGCCGAACGCACATCGGGc                                                                             <  8:278114/149‑1 (MQ=255)
                                                                       cTTGGCCCGGCGCGGGTTGCGCACATGGGCAAAACGCAGCTTGGCCAGCTCTGCATCGTTCAGGCGCTGCTGACCGAAGTCGGAGCTCAGGCCTTCCAGGCGCTCGGTGTCACCGCGGGCTTCGATCCACGCCGAACGCACATCggg                                                                              >  3:254583/1‑147 (MQ=255)
                                                                       cTTGGCCCGGCGCGGGTTGCGCACATGGGCAAAACGCAGCTTGGCCAGCTCTGCATCGTTCAGGCGCTGCTGACCGAAGTCGGAGCTCAGGCCTTCCAGGCGCTCGGTGTCACCGCGGGCTTCGACCCACGCCGAACGCACATCggg                                                                              >  3:22454/1‑147 (MQ=255)
                                                                       cTTGGCCCGGCGCGGGTTGCGCACACGGGCAAAACGCAGCTTGGCCAGCTCTGCATCGTTCAGGCGCTGCTGACCGAAGTCGGAGCTCAGGCCTTCCAGGCGCTCGGTGTCACCGCGGGCTTCGATCCACGCCGAACGCACATCggg                                                                              >  4:49808/1‑147 (MQ=255)
                                                                       cTTGGCCCGGCGCGGCTTGCGCACATGGGCAAAACGCAGCTTGGCCAGCTCTGCATCGTTCAGGCGCTGCTGACCGAAGTCGGAGCTCAGGCCTTCCAGGCGCTCGGTGTCACCGCGGGCTTCGATCCACGCCGAACGCACATCggg                                                                              >  1:33629/1‑147 (MQ=255)
                                                                          ggCCCGGCGCGGGTTGCGCACATGGGCAAAACGCAGCTTGGCCAGCTCTGCATCGTTCAGGCGCTGCTGACCGAAGTCGGAGCTCAGGCCTTCCAGGCGCTCGGTGTCACCGCGGGCTTCGATCCACGCCGAACGCACATCGGGCAgg                                                                          <  4:296931/148‑1 (MQ=255)
                                                                          ggCCCGGCGCGGGTTGCGCACATGGGCAAAACGCAGCTTGGCCAGCTCTGCATCGTTCAGGCGCTGCTGACCGAAGTCGGAGCTCAGGCCTTCCAGGCGCTCGGTGTCACCGCGGGCTTCGATCCACGCCGAACGCACATCGGGCAg                                                                           <  3:58365/147‑1 (MQ=255)
                                                                          ggCCCGGCGCGGGTTGCGCACATGGGCAAAACGCAGCTTGGCCAGCTCTGCATCGTTCAGGCGCTGCTGACCGAAGTCGGAGCTCAGGCCTTCCAGGCGCTCGGTGTCACCGCGGGCTTCGATCCACGCCGAACGCACATCGGGCAg                                                                           <  7:141937/147‑1 (MQ=255)
                                                                           gCCCGGCGCGGGTTGCGCACATGGGCAAAACGCAGCTTGGCCAGCTCTGCATCGTTCAGGCGCTGCTGACCGAAGTCGGAGCTCAg                                                                                                                                       <  4:316175/86‑1 (MQ=255)
                                                                           gCCCGGCGCGGGTTGCGCACATGGGCAAAACGCAGCTTGGCCAGCTCTGCATCGTTCAGGCGCTGCTGACCGAAGTCGGAGCTCAg                                                                                                                                       >  3:316175/1‑86 (MQ=255)
                                                                            cccGGCGCGGGTTGCGCACATGGGCAAAACGCAGCTTGGCCAGCTCTGCATCGTTCAGGCGCTGCTGACCGAAGTCGGAGCTCAGGCCTTCCAGGCGCTCGGTGTCACCGCGGGCTTCGATCCACGCCGAACGCACATCGGGCAGGCCt                                                                       >  7:307156/1‑149 (MQ=255)
                                                                                          cgcACATGGGCAAAACGCAGCTTGGCCAGCTCTGCATCGTTCAGGCGCTGCTGACCGAAGTCGGAGCTCAGGCCTTCCAGGCGCTCGGTGTCACCGCGGGCTTCGATCCACGCCGAACGCACATCGGGCAGGCCTTTGCGCACGTCga                                                          <  4:165312/148‑1 (MQ=255)
                                                                                           gcACATGGGCAAAACGCAGCTTGGCCAGCTCTGCATCGTTCAGGCGCTGCTGACAGAAGTCGGAGCTCAGGCCTTCCAGGCGCTCGGTGTCACCGCGGGCTTCGATCCACGCCGAACGCACATCGGGCAGGCCTTTGCGCACGTCgat                                                         >  3:218677/1‑148 (MQ=255)
                                                                                           gcACATGGGCAAAACGCAGCTTGGCCAGCTCTGCATCGTTCAGGCGCTGCTGACAGAAGTCGGAGCTCAGGCCTTCCAGGCGCTCGGTGTCACCGCGGGCTTCGATCCACGCCGAACGCACATCGGGCAGGCCTTTGCGCACGTCgat                                                         >  3:223198/1‑148 (MQ=255)
                                                                                            cacaTGGGCAAAACGCAGCTTGGCCAGCTCTGCATCGTTCAGGCGCTGCTGACCGAAGTCGGAGCTCAGGCCTTCCAGGCGCTCGGTGTCACCGCGGGCTTCGATCCACGCCGAACGCACATCGGGCAGGCCTTTGCGCACGTCgatga                                                       <  8:115613/149‑1 (MQ=255)
                                                                                             acaTGGGCAAAACGCAGCTTGGCCAGCTCTGCATCGTTCAGGCGCTGCTGACCGAAGTCGGAGCTCAGGCCTTCCAGGCGCTCGGTGTCACCGCGGGCTTCGATCCACGCCGAACGCACATCGGGCAGGCCTTTGCGCACGTCgatgat                                                      <  4:254583/149‑1 (MQ=255)
                                                                                                tGGGCAAAACGCAGCTTGGCCAGCTCTGCATCGTTCAGGCGCTGCTGACAGAAGTCGGAGCTCAGGCCTTCCAGGCGCTCGGTGTCACCGCGGGCTTCGATCCACGCCGAACGCACATCGGGCAGGCCTTTGCGCACGTCGATGATCAc                                                   <  4:223198/149‑1 (MQ=255)
                                                                                                tGGGCAAAACGCAGCTTGGCCAGCTCTGCATCGTTCAGGCGCTGCTGACAGAAGTCGGAGCTCAGGCCTTCCAGGCGCTCGGTGTCACCGCGGGCTTCGATCCACGCCGAACGCACATCGGGCAGGCCTTTGCGCACGTCGATGATCAc                                                   <  4:218677/149‑1 (MQ=255)
                                                                                                  ggCAAAACGCAGCTTGGCCAGCTCTGCATCGTTCAGGCGCTGCTGACCGAAGTCGGAGCTCAGGCCTTCCAGGCGCTCGGTGTCACCGCGGGCTTCGATCCACGCCGAACGCACATCGGGCAGGCCTTTGCGCACGTCGATGATCACGt                                                 >  1:333156/1‑149 (MQ=255)
                                                                                                         cGCAGCTTGGCCAGCTCTGCATCGTTCAGGCGCTGCTGACCGAAGTCGGAGCTCAGGCCTTCCAGGCGCTCGGTGTCACCGCGGGCTTCGATCCACGCCGAACGCACATCGGGCAGGCCTTTGCGCACGTCGATAATCACGTCtggccc                                          >  1:200420/1‑146 (MQ=255)
                                                                                                         cGCAGCTTGGCCAGCTCTGCATCGTTCAGGCGCTGCTGACCGAAGTCGGAGCTCAGGCCTTCCAGGCGCTCGGTGTCACCGCGGGCTTCGAGCCACGCCGAACGCACATCGGGCAGGCCTTTGCGCACGTCGATGATCACGTCGGGGTc                                          <  4:22454/149‑1 (MQ=255)
                                                                                                          gCAGCTTGGCCAGCTCTGCATCGTTCAGGCGCTGCTGACCGAAGTCGGAGCTCAGGCCTTCCAGGCGCTCGGTGTCACCGCGGGCTTCGATCCACGCCGAACGCACATCGGGCAg                                                                           <  5:104531/115‑1 (MQ=255)
                                                                                                          gCAGCTTGGCCAGCTCTGCATCGTTCAGGCGCTGCTGACCGAAGTCGGAGCTCAGGCCTTCCAGGCGCTCGGTGTCACCGCGGGCTTCGATCCACGCCGAACGCACATCGGGCAg                                                                           >  6:104531/1‑115 (MQ=255)
                                                                                                            aGCTTGGCCAGCTCTGCATCGTTCAGGCGCTGCTGACCGAAGTCGGAGCTCAGGCCTTCCAGGCGCTCGGTGTCACCGCGGGCTTCGATCCACGCCGAACGCACATCGGGCAGGCCTTTGCGCACGTCGATGATCACGTCGGGGTCGgt                                       >  7:136992/1‑149 (MQ=255)
                                                                                                                      gCTCTGCATCGTTCAGGCGCTGCTGACCGAAGTCGGAGCTCAGGCCTTCCAGGCGCTCGGTGTCACCGCGGGCTTCGATCCACGCCGAACGCACATCGGGCAGGCCTTTGCGCACGTCGATGATCACGTCGGGGTCGGTGTACGGGc                               >  5:61901/1‑147 (MQ=255)
                                                                                                                          tGCATCGTTCAGGCGCTGCTGACCGAAGTCGGAGCTCAGGCCTTCCAGGCGCTCGGTGTCACCGCGGGCTTCGATCCACGCCGAACGCACATCGGGCAGGCCTTTGCGCACGTCGATGATCACGTCgggg                                            >  3:14392/1‑130 (MQ=255)
                                                                                                                          tGCATAGTTCAGGCGCTGCTGACCGAAGTCGGAGCTCAGGCCTTCCAGGCGCTCGGTGTCACCGCGGGCTTCGATCCACGCCGAACGCACATCGGGCAGGCCTTTGCGCACGTCGATGATCACGTCgggg                                            <  4:14392/130‑1 (MQ=255)
                                                                                                                            cATCGTTCAGGCGCTGCTGACCGAAGTCGGAGCTCAGGCCTTCCAGGCGCTCGGTGTCACCGCGGGCTTCGATCCACGCCGAACGCACATCGGGCAGGCCTTTGCGCACGTCGATGATCACGTCGGGGTCGGTGTACGGGCCAGAc                          <  2:333156/146‑1 (MQ=255)
                                                                                                                            cATCGTTCAGGCGCTGCTGACCGAAGTCGGAGCTCAGGCCTTCCAGGCGCTCGGTGTCACCGCGGGCTTCGATCCACGCCGAACGCACATCGGGCAGGCCTTTGCGCACGTCGATGATCACGTCGGGGTCGGTGTACGGGCCAGAc                          <  2:33629/146‑1 (MQ=255)
                                                                                                                            cATCGTTCAGGCGCTGCTGACCGAAGTCGGAGCTCAGGCCTTCCAGGCGCTCGGTGTCACCGCGGGCTTCGATCCACGCCGAACGCACATCGGGCAGGCCTTTGCGCACGTCGATGATCACGTCGGGGTCGGTGTACGGGCCAGACgtg                       >  3:215249/1‑149 (MQ=255)
                                                                                                                                 ttCAGGCGCTGCTGACCGAAGTCGGAGCTCAGGCCTTCCAGGCGCTCGGTGTCACCGCGGGCTTCGATCCACGCCGAACGCACATCGGGCAGGCCTTTGCGCACGTCGATGATCACGTCGGGGTCGGTGTACGGGCCAGACGTGTCATa                  >  4:54459/1‑149 (MQ=255)
                                                                                                                                     ggCGCTGCTGACCGAAGTCGGAGCTCAGGCCTTCCAGGCGCTCGGTGTCACCGCGGGCTTCGATCCACGCCGAACGCACATCGGGCAGGCCTTTGCGCACGTCGATGATCACGTCGGGGTCGGTGTACGGGCCAGACGTGTCATACAcc              <  6:61901/149‑1 (MQ=255)
                                                                                                                                        gctgctGACCGAAGTCGGAGCTCAGGCCTTCCAGGCGCTCGGTGTCACCGCGGGCTTCGATCCACGCCGAACGCACATCGGGCAGGCCTTTGCGCACGTCGATGATCACGTCGGGGTCGGTGTACGGGCCAGACGTGTCATACACCAGc           >  7:140475/1‑149 (MQ=255)
                                                                                                                                            ctgACCGAAGTCGGAGCTCAGGCCTTCCAGGCGCTCGGTTTCACCGCGGGCTTCGATCCACGCCGAACGCACATCGGGCAGGCCTTTGCGCACGTCGATGATCACGTCGGGGTCGGTGTACGGGCCAGACGTGTCATACACCAGCACCg       <  8:136992/149‑1 (MQ=255)
                                                                                                                                            ctgACCGAAGTCGGAGCTCAGGCCTTCCAGGCGCTCGGTGTCACCGCGGGCTTCGATCCACGCCGAACGCACATCGGGCAGGCCTTTGCGCACGTCGATGATCACGTCGGGGTCGGTGTACGGGCCAGACGTGTCATACACCAGCACCg       <  8:307156/149‑1 (MQ=255)
                                                                                                                                                 cGAAGTCGGAGCTCAGGCCTTCCAGGCGCTCGGTGTCACCGCGGGCTTCGATCCACGCCGAACGCACATCGGGCAGGCCTTTGCGCACGTCGATGATCACGTCGGGGTCGGTGTACGGGCCAGACGTGTCATACACCAGCACCGGGGCg  >  5:199642/1‑149 (MQ=255)
                                                                                                                                                 |                                                                                                                                                    
GCCACGTATTCCATTTCGGCGGTGATGATGCCCTGACGGGCATAGTGCATCTGCGAGACGTTGGCGCCGGCCTTGGCCCGGCGCGGGTTGCGCACATGGGCAAAACGCAGCTTGGCCAGCTCTGCATCGTTCAGGCGCTGCTGACCGAAGTCGGAGCTCAGGCCTTCCAGGCGCTCGGTGTCACCGCGGGCTTCGATCCACGCCGAACGCACATCGGGCAGGCCTTTGCGCACGTCGATGATCACGTCGGGGTCGGTGTACGGGCCAGACGTGTCATACACCAGCACCGGGGCG  >  NC_002947/5595606‑5595899

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