Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 401,889 C→T 5.8% intergenic (+416/+13) PP_0334 → / ← PP_0335 transposase/glutathione S‑transferase family protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_002947401,8890CT 5.8% 207.0 / 2.1 69intergenic (+416/+13)PP_0334/PP_0335transposase/glutathione S‑transferase family protein
Reads supporting (aligned to +/- strand):  ref base C (31/34);  new base T (2/2);  total (33/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 = 7.05e-01

GAGTCCATACATAGGAGCGGGCTTGCCCGCGAACACACATAATGGGGTGTCAGGCACCGGCTTCGCCGGTGTTCGCGGGCACGCCCGCTCCCAAAAGCCCACCACAGATAGCGAACACTGTGGAAAACCTTGTAGGAGCGGGTTTACCCGCGAAGGGGCCCTCAGGCCACTGGCGCTGTAGCGACCATCGATGGCCGCTTGACCACCTCGGCATACCACGCCGCCAACCCAGGCTGGCGCTCACGCCAGCCAAACCCGGGCTGACGCAAATCCAGATAGCCCAGCGCACAGGCCAC  >  NC_002947/401741‑402036
                                                                                                                                                    |                                                                                                                                                    
gAGTCCATACATAGGAGCGGGCTTGCCCGCGAACACACATAATGGGGTGTCAGGCACCGGCTTCGCCGGTGTTCGCGGGCACGCCCGCTCCCAAAAGCCCACCACAGATAGCGAACACTGTGGAAAACCTTGTAGGAGCGGGTTTAccc                                                                                                                                                      <  2:42469/149‑1 (MQ=255)
         cataGGAGCGGGCTTGCCCGCGAACACACATAATGGGGTGTCAGGCACCGGCTTCGCCGGTGTTCGCGGGCACGCCCGCTCCCAAAAGCCCACCACAGATAGCGAACACTGTGGAAAACCTTGTAGGAGCGGGTTTACCCGCGAAgggg                                                                                                                                             <  8:91515/149‑1 (MQ=255)
         cataGGAGCGGGCTTGCCCGCGAACACACATAATGGGGTGTCAGGCACCGGCTTCGCCGGTGTTCGCGGGCACGCCCGCTCCCAAAAGCCCACCACAGATAGCGAACACTGTGGAAAACCTTGTAGGAGCGGGTTTACCCGCGAAggg                                                                                                                                              <  1:148373/148‑1 (MQ=255)
         cataGGAGCGGGCTTGCCCGCGAACACACATAATGGGGTGTCAGGCACCGGCTTCGCCGGTGTTCGCGGGCACGCCCGCTCCCAAAAGCCCACCACAGATAGCGAACACTGTGGAAAACCTTGTAGGAGCGGGTTTACCCGCGAAggg                                                                                                                                              <  2:327978/148‑1 (MQ=255)
            aGGAGCGGGCTTGCCCGCGAACACACATAATGGGGTGTCAGGCACCGGCTTCGCCGGTGTTCGCGGGCACGCCCGCTCCCAAAAGCCCACCACAGATAGCGAACACTGTGGAAAACCTTGTAGGAGCGGGTTTACCCGCGAAGGGGccc                                                                                                                                          >  3:39996/1‑149 (MQ=255)
               aGCGGGCTTGCCCGCGAACACACATAATGGGGTGTCAGGCACCGGCTTCGCCGGTGTTCGCGGGCACGCCCGCTCCAAAAAGCCCACCACAGATAGCGAACACTGTGGAAAACCTTGTAGTAGCGGGTTTACCCGCGAAGGGGCCCTCa                                                                                                                                       <  4:244270/149‑1 (MQ=255)
                gCGGGCTTGCCCGCGAACACACATAATGGTGTGTCAGGCACCGGCTTCGCCGGTGTTCGCGGGCACGCCCGCTCCCAAAAGCCCACCACAGATAGCGAACACTGTGGAAAACCTTGTAGGAGCGGGTTTACCCGCGAAGGGGCCCTCAg                                                                                                                                      <  4:145196/149‑1 (MQ=255)
                          ccGCGAACACACATAATGGGGTGTCAGGCACCGGCTTCGCCGGTGTTCGCGGGCACGCCCGCTCCCAAAAGCCCACCACAGATAGCGAACACTGTGGAAAACCTTGTAGGAGCGGGTTTACCCGCGAAGGGGCCCTCAGGCCACTGgcg                                                                                                                            <  3:153495/149‑1 (MQ=255)
                               aaCACACATAATGGGGTGTCAGGCACCGGCTTCGCCGGTGTTCGCGGGCACGCCCGCTCCCAAAAGCCCACCACAGATAGCGAACACTGTGGAAAACCTTGTAGGAGCGGGTTTACCCGCGAAGGGGCCCTCAGGCCACTGGCGCTGTa                                                                                                                       <  6:13446/149‑1 (MQ=255)
                                  acacaTAATGGGGTGTCAGGCACCGGCTTCGCCGGTGTTCGCGGGCACGCCCGCTCCCAAAAGCCCACCACAGATAGCGAACACTGTGGAAAACCTTGTAGGAGCGGGTTTACCCGCGaa                                                                                                                                                 >  5:35416/1‑120 (MQ=255)
                                  acacaTAATGGGGTGTCAGGCACCGGCTTCGCCGGTGTTCGCGGGCACGCCCGCTCCCAAAAGCCCACCACAGATAGCGAACACTGTGGAAAACCTTGTAGGAGCGGGTTTACCCGCGaa                                                                                                                                                 <  6:35416/120‑1 (MQ=255)
                                  acacaTAATGGGGTGTCAGGCACCGGCTTCGCCGGTGTTCGCGGGCACGCCCGCTCCCAAAAGCCCACCACAGATAGCGAACACTGTGGAAAACCTTGTAGGAGCGGGTTTACCCGCGAAGGGGCCCTCAGGCCACTGGCGCTGTAGc                                                                                                                     >  1:20474/1‑148 (MQ=255)
                                  acacaTAATGGGGTGTCAGGCACCGGCTTCGCCGGTGTTCGCGGGCACGCCCGCTCCCAAAAGCCCACCACAGATAGCGAACACTGTGGAAAACCTTGTAGGAGCGGGTTTACCCGCGAAGGGGCCCTCAGGCCACTGGCGCTGTAGc                                                                                                                     >  4:124181/1‑148 (MQ=255)
                                   cacaTAATGGGGTGTCAGGCACCGGCTTCGCCGGTGTTCGCGGGCACGCCCGCTCCCAAAAGCCCACCACAGATAGCGAACACTGTGGAAAACCTTGTAGGAGCGGGTTTACCCGCGAAGGGGCCCTCAGGCCACTGGCGCTGTAGCGa                                                                                                                   <  4:320387/149‑1 (MQ=255)
                                    acaTAATGGGGTGTCAGGCACCGGCTTCGCCGGTGTTCGCGGGCACGCCCGCTCCCAAAAGCCCACCACAGATAGCGAACACTGTGGAAAACCTTGTAGGAGCGGGTTTACCCGCGAAGGGGCCCTCAGGCCACTGGCGCTGTAGCGAc                                                                                                                  >  8:164784/1‑149 (MQ=255)
                                     caTAATGGGGTGTCAGGCACCGGCTTCGCCGGTGTTCGCGGGCACGCCCGCTCCCAAAAGCCCACCACAGATAGCGAACACTGTGGAAAACCTTGTAGGAGCGGGTTTACCCGCGAAGGGGCCCTCAGGCCACTGGCGCTGTAGCGAcc                                                                                                                 <  8:140791/149‑1 (MQ=255)
                                     caTAATGGGGTGTCAGGCACCGGCTTCGCCGGTGTTCGCGGGCACGCCCGCTCCCAAAAGCCCACCACAGATAGCGAACACTGTGGAAAACCTTGTAGGAGCGGGTTTACCCGCGAAGGGGCCCTCAGGCCACTGGCGCTGTAGCGAcc                                                                                                                 <  4:53354/149‑1 (MQ=255)
                                          tGGGGTGTCAGGCACCGGCTTCGCCGGTGTTCGCGGGCACGCCCGCTCCCAAAAGCCCACCACAGATAGCGAACACTGTGGAAAACCTTGTAGGAGCGGGTTTACCCGCGAAGGGGCCCTCAGGCCACTGGCGCTGTAGCGACCATCGa                                                                                                            >  2:249823/1‑149 (MQ=255)
                                            gggTGTCAGGCACCGGCTTCGCCGGTGTTCGCGGGCACGCCCGCTCCCAAAAGCCCACCACAGATAGCGAACACTGTGGAAAACCTTGTAGGAGCGGGTTTACCCGCGAAGGGGCCCTCAGGCCACTGGCGCTGTAGCGACCATCGAt                                                                                                           <  2:62316/148‑1 (MQ=255)
                                                 tCAGGCACCGGCTTCGCCGGTTTTCGCGGGCACGCCCGCTCCCAAAAGCCCACCACAGATAGCGTACACTGTGGAAAACCTTGTAGGAGCGGGTTTACCCGCGAAGGGGCCCTCAGGCCACTGGCGCTGTAGCGACCATCGATGGCCGc                                                                                                     <  8:330399/149‑1 (MQ=255)
                                                   gaacACCGGCGAAGCCGGTGTTCGCGGGCACGCCCGCTCCCAAAAGCCCACCACAGATAGCGAACACTGTGGAAAACCTTGTAGGAGCGGGTTTACCCGCGAAGGGGCCCTCAGGCCACTGGCGCTGTAGCGACCATCGATg                                                                                                          >  5:113228/4‑142 (MQ=255)
                                                   gaacACCGGCGAAGCCGGTGTTCGCGGGCACGCCCGCTCCCAAAAGCCCACCACAGATAGCGAACACTGTGGAAAACCTTGTAGGAGCGGGTTTACCCGCGAAGGGGCCCTCAGGCCACTGGCGCTGTAGCGACCATCGATg                                                                                                          <  6:113228/139‑1 (MQ=255)
                                                     gCACCGGCTTCGCCGGTGTTCGCGGGCACGCCCGCTCCCAAAAGCCCACCACAGATAGCGAACACTGTGGAAAACCTTGTAGGAGCGGGTTTACCCGCGAAGGGGCCCTCAGGCCACTGGCGCTGTAGCGACCATCGATGGCCGCTTGa                                                                                                 >  8:218783/1‑149 (MQ=255)
                                                     gCACCGGCTTCGCCGGTGTTCGCGGGCACGCCCGCTCCCAAAAGCCCACCACAGATAGCGAACACTGTGGAAAACCTTGTAGGAGCGGGTTTACCCGCGAAGGGGCCCTCAGGCCACTGGCGCTGTAGCGACCATCGATGGCCGCTTGa                                                                                                 >  4:315018/1‑149 (MQ=255)
                                                      cACCGGCTTCGCCGGTGTTCGCGGGCACGCCCGCTCCCAAAAGCCCACCACAGATAGCGAACACTGTGGAAAACCTTGTAGGAGCGGGTTTACCCGCGAAGGGGCCCTCAGGCCACTGGCGCTGTAGCGACCATCGATGGCCGCTTGac                                                                                               <  3:124181/149‑1 (MQ=255)
                                                             ttCGCCGGTGTTCGCGGGCACGCCCGCTCCCAAAAGCCCACCACAGATAGCGAACACTGTGGAAAACCTTGTAGGAGCGGGTTTACCCGCGAAGGGGCCCTCAGGCCACTGGCGCTGTAGCGACCATCGATGGCCGCTTGACCACCTCg                                                                                        >  3:18278/1‑149 (MQ=255)
                                                                gCCGGTGTTCGCGGGCACGCCCGCTCCCAAAAGCCCACCACAGATAGCGAACACTGTGGAAAACCTTGTAGGAGCGGGTTTACCCGCGAAGGGGCCCTCAGGCCACTGGCGCTGTAGCGACCATCGATGGCCGCTTGACCACCTCGGCa                                                                                     >  7:79565/1‑149 (MQ=255)
                                                                 ccGGTGTTCGCGGGCACGCCCGCTCCCAAAAGCCCACCACAGATAGCGAACACTGTGGAAAACCTTGTAGGAGCGGGTTTACCCGCGAAGGGGCCCTCAGGCCACTGGCGCTGTAGCGACCATCGATGGCCGCTTGACCACCTCGGCa                                                                                     >  8:185235/1‑148 (MQ=255)
                                                                  cGGTGTTCGCGGGCACGCCCGCTCCCAAAAGCCCACCACAGATAGCGAACACTGTGGAAAACCTTGTAGGAGCGGGTTTACCCGCGAAGGGGCCCTCAGGCCACTGGCGCTGTAGCGACCATCGATGGCCGCTTGACCACCTc                                                                                         <  6:217022/143‑1 (MQ=255)
                                                                  cGGTGTTCGCGGGCACGCCCGCTCCCAAAAGCCCACCACAGATAGCGAACACTGTGGAAAACCTTGTAGGAGCGGGTTTACCCGCGAAGGGGCCCTCAGGCCACTGGCGCTGTAGCGACCATCGATGGCCGCTTGACCACCTc                                                                                         >  5:217022/1‑143 (MQ=255)
                                                                   ggTGTTCGCGGGCACGCCCGCTCCCAAAAGCCCACCACAGATAGCGAACACTGTGGAAAACCTTGTAGGAGCGGGTTTACCCGCGAAGGGGCCCTCAGGCCACTGGCGCTGTAGCGACCATCGATGGCCGCTTGACCACCTCGGCATAc                                                                                  <  8:79565/149‑1 (MQ=255)
                                                                      gtTCGCGGGCACGCCCGCTCCCAAAAGCCCACCACAGATAGCGAACACTGTGGAAAACCTTGTAGGAGCGGGTTTACCCGCGAAGGGGCCCTCAGGCCACTGGCGCTGTAGCGCCCATCGATGGCCGCTTGACCACCTCGGCATACCAc                                                                               >  2:141494/1‑149 (MQ=255)
                                                                       ggagcgGGCGTGCCCGCTCCCAAAAGCCCACCACAGATAGCGAACACTGTGGAAAACCTTGTAGGAGCGGGTTTACCCGCGAAGGGGCCCTCAGGCCACTGGCGCTGTAGCGACCATCGATGGCCGCTTGACCACCTCGGCATACCAcg                                                                              >  4:65491/4‑149 (MQ=255)
                                                                                     cGCTCCCAAAAGCCCACCACAGATAGCGAACACTGTGGAAAACCTTGTAGGAGCGGGTTTACCCGCGAAGGGGCCCTCAGGCCACTGGCGCTGTAGCGACCATCGATGGCCGCTTGACCACCTCGGCATACCACGCCGCCAACCCAGGc                                                                <  7:137105/149‑1 (MQ=255)
                                                                                      gCTCCCAAAAGCCCACCACAGATAGCGAACACTGTGGAAAACCTTGTAGGAGCGGGTTTACCCGCGAAGGGGCCCTCAGGCCACTGGCGCTGTAGCGACCATCGATGGCCGCTTGACCACCTCGGCATACCACGCCGCCAACCCAGGCt                                                               >  5:32388/1‑149 (MQ=255)
                                                                                        tCCCAAAAGCCCACCACAGATAGCGAACACTGTGGAAAACCTTGTAGGAGCGGGTTTACCCGCGAAGGGGCCCTCAGGCCACTGGCGCTGTAGCGACCATCGATGGCCGCTTGACCACCTCGGCATACCACGCCGCCAACCCAGGCTgg                                                             >  3:75789/1‑149 (MQ=255)
                                                                                         cccAAAAGCCCACCACAGATAGCGAACACTGTGGAAAACCTTGTAGGAGCGGGTTTACCCGCGAAGGGGCCCTCAGGCCACTGGCGCTGTAGCGACCATCGATGGCCGCTTGACCACCTCGGCATACCACGCCGCCAACCCAGGCTGgc                                                            >  4:279486/1‑149 (MQ=255)
                                                                                         cccAAAAGCCCACCACAGATAGCGAACACTGTGGAAAACCTTGTAGGAGCGGGTTTACCCGCGAAGGGGCCCTCAGGCCACTGGCGCTGTAGCGACCATCGATGGCCGCTTGACCACCTCGGCATACCACGCCGCCAACCCAGGCTGgc                                                            >  8:190494/1‑149 (MQ=255)
                                                                                          ccAAAAGCCCACCACAGATAGCGAACACTGTGGAAAACCTTGTAGGAGCGGGTTTACCCGCGAAGGGGCCCTCAGGCCACTGGCGCTGTAGCGACCATCGATGGCCGCTTGACCACCTCGGCATACCACGCCGCCAACCCAGGCTGgcg                                                           <  3:202463/149‑1 (MQ=255)
                                                                                           cAAAAGCCCACCACAGATAGCGAACACTGTGGAAAACCTTGTAGGAGCGGGTTTACCCGCGAAGGGGCCCTCAGGCCACTGGCGCTGTAGCGACCATCGATGGCCGCTTGACCCACCTCGGCATACCACGCCGCCAACCCAGGCTGgcg                                                           <  1:249823/149‑1 (MQ=255)
                                                                                            aaaaGCCCACCACAGATAGCGAACACTGTGGAAAACCTTGTAGGAGCGGGTTTACCCGCGAAGGGGCCCTCAGGCCACTGGCGCTGTAGCGACCATCGATGGCCGCTTGACCACCTCGGCATACCACGCCGCCAACCCAGGCTGGCGCt                                                         >  8:54763/1‑149 (MQ=255)
                                                                                            aaaaGCCCACCACAGATAGCGAACACTGTGGAAAACCTTGTAGGAGCGGGTTTACCCGCGAAGGGGCCCTCAGGCCACTGGCGCTGTAGCGACCATCGATGGCCGCTTGACCACCTCGGCATACCACGCCGCCAAACCAGGCTGGCGCt                                                         >  1:286075/1‑149 (MQ=255)
                                                                                               aGCCCACCACAGATAGCGAACACTGTGGAAAACCTTGTAGGAGCGGGTTTACCCGCGAAGGGGCCCTCAGGCCACTGGCGCTGTAGCGACCATCGATGGCCGCTTGACCACCTCGGCATACCACGCCGCCAACCCAGGCTGGCGCTCAc                                                      <  2:20474/149‑1 (MQ=255)
                                                                                                     ccacAGATAGCGAACACTGTGGAAAACCTTGTAGGAGCGGGTTTACCCGCGAAggg                                                                                                                                              <  1:163606/56‑1 (MQ=255)
                                                                                                     ccacAGATAGCGAACACTGTGGAAAACCTTGTAGGAGCGGGTTTACCCGCGAAggg                                                                                                                                              >  2:163606/1‑56 (MQ=255)
                                                                                                     ccacAGATAGCGAACACTGTGGAAAACCTTGTAGGAGCGGGTTTACCCGCGAAGGGGCCCTCAGGCCTCTGGCGCTGTAGCGACCATCGATGGCCGCTTGACCACCTCGGCATACCACGCCGCCAACCCAGGCTGGCGCTCACgccagc                                                >  8:95198/1‑149 (MQ=255)
                                                                                                       acaGATAGCGAACACTGTGGAACACCTTGTAGGAGCGGGTTTACCCGCGAAGGGGCCCTCAGGCCACTGGCGCTGTAGCGACCATCGATGGCCGCTTGACCACCTCGGCATACCACGCCGCCAACCCAGGCTGGCGTGAGCgccagcct                                              <  7:164784/149‑2 (MQ=255)
                                                                                                        caGATAGCGAACACTGTGGAAAACCTTGTAGGAGCGGGTTTACCCGCGAAGGGGCCCTCAGGCCACTGGCGCTGTAGCGACCATCGATGGCCGCTTGACCACCTCGGCATACCACGCCGCCAACCCAGGCTGGCGCTCACGCCAGCCaa                                             <  3:279486/149‑1 (MQ=255)
                                                                                                        caGATAGCGAACACTGTGGAAAACCTTGTAGGAGCGGGTTTACCCGCGAAGGGGCCCTCAGGCCACTGGCGCTGTAGCGACCATCGATGGCCGCTTGACCACCTCGGCATACCACGCCGCCAACCCAGGCTGGCGCTCACGCCAGCCaa                                             <  7:190494/149‑1 (MQ=255)
                                                                                                          gATAGCGAACACTGTGGAAAACCTTGTAGGAGCGGGTTTACCCGCGAAGGGGCCCTCAGGCCACTGGCGCTGTAGCGACCATCGATGGCCGCTTGACCACCTCGGCATACCACGCCGCCAACCCAGGCTGGCGTGAGCgccagcctggg                                           <  7:185235/149‑5 (MQ=255)
                                                                                                              gCGAACACTGTGGAAAACCTTGTAGGAGCGGGTTTACCCGCCAAGGGGCCCTCAGGCCACTGGCGCTGTAGCGACCATCGATGGCCGCTTGACCACCTCGGCATACCACGCCGCCAACCCAGGCTGGCGCTCACGCCAGCCAAACCCgg                                       <  5:323008/149‑1 (MQ=255)
                                                                                                               cGAACACTGTGGAAAACCTTGTAGGAGCGGGTTTACCTGCGAAGGGGCCCTCAGGCCACTGGCGCTGTAGCGACCATTGATGGCCGCTTGACCACCTCGGCATACCACGCCGCCAACCCAGGCTGGCGCTCACGCCAGCCAAACCCggg                                      >  5:30999/1‑149 (MQ=255)
                                                                                                               cGAACACTGTGGAAAACCTTGTAGGAGCGGGTTTACCTGCGAAGGGGCCCTCAGGCCACTGGCGCTGTAGCGACCATTGATGGCCGCTTGACCACCTCGGCATACCACGCCGCCAACCCAGGCTGGCGCTCACGCCAGCCAAACCCggg                                      >  5:25542/1‑149 (MQ=255)
                                                                                                               cGAACACTGTGGAAAACCTTGTAGGAGCGGGTTTACCCGCGAAGGGGCCCTCAGGCCACTGGCGCTGTAGCGACCATCGATGGCCGCTTGACCACCTCGGCATACCACGCCGCCAACCCAGGCTGGCGCTCACGCCAGCCAAACCCggg                                      <  6:294806/149‑1 (MQ=255)
                                                                                                               cGAACACTGTGGAAAACCTTGTAGGAGCGGGTTTACCCGCGAAGGGGCCCTCAGGCCACTGGCGCTGTAGCGACCATCGATGGCCGCTTGACCACCTCGGCATACCACGCCGCCAACCCAGGCTGGCGCTCACGCCAGCCAAACCCggg                                      <  4:75789/149‑1 (MQ=255)
                                                                                                               cGAACACTGTGGAAAACCTTGTAGGAGCGGGTTTACCCGCGAAGGGGCCCTCAGGCCACTGGCGCTGTAGCGACCATCGATGGCCGCTTGACCACCTCGGCATACCACGCCGCCAACCCAGGCTGGCGCTCACGCCAGCCAAACCCggg                                      >  4:96182/1‑149 (MQ=255)
                                                                                                                gAACACTGTGGAAAACCTTGTAGGAGCGGGTTTACCCGCGAAGGGGCCCTCAGGCCACTGGCGCTGTAGCGACCATCGATGGCCGCTTGACCACCTCGGCATACCACGCCGCCACCCCAGGCTGGCGCTCACGCCAGCCAAACCCGGGc                                     >  5:342619/1‑149 (MQ=255)
                                                                                                                        tgGAAAACCTTGTAGGAGCGGGTTTACCCGCGAAGGGGCCCTCAGGCCACTGGCGCTGTAGCGACCATCGATGGCCGCTTGACCACCTCGGCATACCACGCCGCCAACCCAGGCTGGCGCTCACGCCAGCCAAACCCGGGCTGCCGCa                              >  7:7359/1‑148 (MQ=255)
                                                                                                                          gAAAACCTTGTAGGAGCGGGTTTACCCGCGAAGGGGCCCTCAGGCCACTGGCGCTGTAGCGACCATCGATGGCCGCTTGac                                                                                               <  4:120211/81‑1 (MQ=255)
                                                                                                                          gAAAACCTTGTAGGAGCGGGTTTACCCGCGAAGGGGCCCTCAGGCCACTGGCGCTGTAGCGACCATCGATGGCCGCTTGac                                                                                               >  3:120211/1‑81 (MQ=255)
                                                                                                                          gAAAACCTTGTAGGAGCGGGTTTACCCGCGAAGGGGCCCTCAGGCCACTGGCGCTGTAGCGACCATCGATGGCCGCTTGACCACCTCGGCATACCACGCCGCCAACCCAGGCTGGCGCTCACGCCAGCCAAACCCGGGCTGACGCaaa                            <  6:32388/148‑1 (MQ=255)
                                                                                                                          gAAAACCTTGTAGGAGCGGGTTTACCCGCGAAGGGGCCCTCAGGCCACTGGCGCTGTAGCGACCATCGATGGCCGCTTGACCACCTCGGCATACCACGCCGCCAACCCAGGCTGGCGCTCACGCCAGCCAAACCCGGGCTGACGCAAAt                           <  4:18278/149‑1 (MQ=255)
                                                                                                                          gAAAACCTTGTAGGAGCGGGTTTACCCGCGAAGGGGCCCTCAGGCCACTGGCGCTGTAGCGACCATCGATGGCCGCTTGACCACCTCGGCATACCACGCCGCCAACCCAGGCTGGCGCTCACGCCAGCCAAACCCGGGCTGACGCAAAt                           >  5:376982/1‑149 (MQ=255)
                                                                                                                            aaaCCTTGTAGGAGCGGGTTTACCCGCGAAGGGGCCCTCAGGCCACTGGCGCTGTAGCGACCATCGATGGCCGCTTGACCACCTCGGCATACCACGCCGCCAACCCAGGCTGGCGCTCACGCCAGCCAAACCCGGGCTGACGCAAATc                          >  7:351512/1‑148 (MQ=255)
                                                                                                                               ccTTGTAGGAGCGGGTTTACCCGCGAAGGGGCCCTCAGGCCACTGGCGCTGTAGCGACCATCGATGGCCGCTTGACCACCTCGGCATACCACGCCGCCAACCCAGGCTGGCGCTCACGCCAGCCAAACCCGGGCTGACGCAAATCCAGa                      <  1:238811/149‑1 (MQ=255)
                                                                                                                                    tAGGAGCGGGTTTACCTGCGAAGGGGCCCTCAGGCCACTGGCGCTGTAGCGACCATTGATGGCCGCTTGACCACCTCGGCATACCACGCCGCCAACCCAGGCTGGCGCTCACGCCAGCCAAACCCGGGCTGACGCAAATCCAGATAGcc                 <  6:30999/149‑1 (MQ=255)
                                                                                                                                    tAGGAGCGGGTTTACCTGCGAAGGGGCCCTCAGGCCACTGGCGCTGTAGCGACCATTGATGGCCGCTTGACCACCTCGGCATACCACGCCGCCAACCCAGGCTGGCGCTCACGCCAGCCAAACCCGGGCTGACGCAAATCCAGATAGcc                 <  6:25542/149‑1 (MQ=255)
                                                                                                                                               ttACCCGCGAAGGGGCCCTCAGGCCACTGGCGCTGTAGCGACCATCGATGGCCGCTTGACCACCTCGGCATACCACGCCGCCAACCCAGGCTGGCGCTCACGCCAGCCAAACCCTGGCTGACGCACATCCAGATAGCCCAGCGCACAgg      >  3:295458/1‑149 (MQ=255)
                                                                                                                                                tACCCGCGAAGGGGCCCTCAGGCCACTGGCGCTGTAGCGACCATCGATGGCCGCTTGACCACCTCGGCATACCACGCCGCCAACCCAGGCTGGCGCTCACGCCAGCCAAACCCGGGCTGACGCAAATCCAGATAGCCCAGCGCACAgg      <  3:65491/148‑1 (MQ=255)
                                                                                                                                                 aCCCGCGAAGGGGCCCTCAGGCCACTGGCGCTGTAGCGACCATCGATGGCCGCTTGACCACCTCGGCATACCACGCCGCCAACCCAGGCTGGCGCTCACGCCAGCCAAACCCGGGCTGACGCAAATCCAGATAGCCCAGCGCACAGGcc    >  8:252934/1‑149 (MQ=255)
                                                                                                                                                   ccGCGAAGGGGCCCTCAGGCCACTGGCGCTGTAGCGACCATCGATGGCCGCTTGACAACCTCGGCATACCACGCCGCCAACCCAGGCTGGCGCTCACGCCAGCCAAACCCGGGCTGACGCAAATCCAGATAGCCCAGCGCACAGGCCAc  >  6:344039/1‑149 (MQ=255)
                                                                                                                                                   ccGCGAAGCGGCCCTCAGGCCACTGCCGCTGTAGCGACCATCGATGGCCGCATGACCACCTCGGCATACCACGCCGCAAACCCAGGCTAGCGCTCCCGCCAGCAAAACCCGCGCTCACGCAAATCCAGATAGCCAAGCGCACAGGCCAc  >  4:60769/1‑149 (MQ=255)
                                                                                                                                                    |                                                                                                                                                    
GAGTCCATACATAGGAGCGGGCTTGCCCGCGAACACACATAATGGGGTGTCAGGCACCGGCTTCGCCGGTGTTCGCGGGCACGCCCGCTCCCAAAAGCCCACCACAGATAGCGAACACTGTGGAAAACCTTGTAGGAGCGGGTTTACCCGCGAAGGGGCCCTCAGGCCACTGGCGCTGTAGCGACCATCGATGGCCGCTTGACCACCTCGGCATACCACGCCGCCAACCCAGGCTGGCGCTCACGCCAGCCAAACCCGGGCTGACGCAAATCCAGATAGCCCAGCGCACAGGCCAC  >  NC_002947/401741‑402036

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: