Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 6,122,012 T→A 5.1% K46N (AAA→AAT) ‡ PP_5371 ← rubredoxin/rubredoxin reductase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029476,122,0120TA 5.1% 302.2 / 3.9 97K46N (AAA→AAT) ‡PP_5371rubredoxin/rubredoxin reductase
Reads supporting (aligned to +/- strand):  ref base T (47/45);  new base A (2/3);  total (49/48)
Fisher's exact test for biased strand distribution p-value = 6.78e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.13e-01

CCTGTGCCAGGCCGTAGCCGGCATGGCCGCTACCGATGATGACCACCGGTTGGGCGGGACGAGCGGCAACCAGTGGTTGCGCCTGGATAACCGGCGCGGTAGCCACTGCTGCGATCGGGGTGCTGAGTTCGATCATTTCAAAGTCGGCTTTGCCGACCTTGCACTCAGGGCATAGCCAGTCTTCGGGCACGTCCTCCCAGCGGGTACCGGGGGCAATGCCGTCGTCGGGCCAACCTTTGGCTTCGTCGTAGATCAGACCACAGATGATGCAGAGCCAGGTTTTCATTATGCTGCC  >  NC_002947/6121864‑6122158
                                                                                                                                                    |                                                                                                                                                  
ccTGTGCCAGGCCGTAGCCGGCATGGCCGCTACCGATGATGACCACCGGTTGGGCGGGACGAGCGGCAACCAGTGGTTGCGCCTGGATAACCGGCGCGGTAGCCACTGCTGCGATCGGGGTGCTGAGTTCGATCATTTCAAAGTCGGCt                                                                                                                                                    >  2:180515/1‑149 (MQ=255)
   gtgCCAGGCCGTAGCCGGCATGGCCGCTACCGATGATGACCACCGGTTGGGCGGGACGAGCGGCAACCAGTGGTTGCGCCTGGATAACCGGCGCGGTAGCCACTGCTGCGATCGGGGTGCTGAGTTCGATCATTTCAAAGTCGGCTTTg                                                                                                                                                 <  4:106207/149‑1 (MQ=255)
    tgCCAGGCCGTAGCCGGCATGGCCGCTACCGATGATGACCACCGGTTGGGCGGGACGAGCGGCAACCAGTGGTTGCGCCTGGATAACCGGCGCGGTAGCCACTGCTGCGATCGGGGTGCTGAGTTCGATCATTTCAAAGTCGGCTTTGc                                                                                                                                                <  2:64667/149‑1 (MQ=255)
      ccAGGCCGTAGCCGGCATGGCCGCTACCGATGATGACCACCGGTTGGGCGGGACGAGCGGCAACCAGTGGTTGCGCCTGGATAACCGGCGCGGTAGCCACTGCTGCGATCGGGGTGCTGAGTTCGATCATTTCAAAGTCGGCTTTGcc                                                                                                                                               >  1:311345/1‑148 (MQ=255)
      ccAGGACGTAGCCGGCATGGCCGCTACCGATGATGACCACCGGTTGGGCGGGACGAGCGGCAACCAGTGGTTGCGCCTGGATAACCGGCGCGGTAGCCACTGCTGCGATCGGGGTGCTGAGTTCGATCATTTCAAAGTCGGCTTTGcc                                                                                                                                               >  7:48060/1‑148 (MQ=255)
       cAGGCCGTAGCCGGCATGGCCGCTACCGATGATGACCACCGGTTGGGCGGGACGAGCGGCAACCAGTGGTTGCGCCTGGATAACCGGCGCGGTAGCCACTGCTGCGATCGGGGTGCTGAGTTCGATCATTTCAAAGTCGGCTTTGCCGa                                                                                                                                             >  7:291483/1‑149 (MQ=255)
           cggTAGCGGCCATGCCGGCTACCGATGATGACCACCGGTTGGGCGGGACGAGCGGCAACCAGTGGTTGCGCCTGGATAACCGGCGCGGTAGCCACTGCTGCGATCGGGGTGCTGAGTTCGATCATTTCAAAGTCGGCTTTGCCGACCtt                                                                                                                                         >  5:226521/3‑149 (MQ=255)
            cGTAGCCGGCATGGCCGCTACCGATGATGACCACCGGTTGGGCGGGACGAGCGGCAACCAGTGGTTGCGCCTGGATAACCGGCGCGGTAGCCACTGCTGCGATCGGGGTGCTGAGTTCGATCATTTCAAAGTCGGCTTTGCCGACCTTg                                                                                                                                        <  3:134714/149‑1 (MQ=255)
             gTAGCCGGCATGGCCGCTACCGATGATGACCACCGGTTGGGCGGGACGAGCGGCAACCAGTGGTTGCGCCTGGATAACCGGCGCGGTAGCCACTGCTGCGATCGGGGTGCTGAGTTCGATCATTTCAAAGTCGGCTTTGCCGACCTTGc                                                                                                                                       <  4:350075/149‑1 (MQ=255)
                gCCGGCATGGCCGCTACCGATGATGACCACCGGTTGGGCGGGACGAGCGGCAACCAGTGGTTGCGCCTGGATAACCGGCGCGGTAGCCACTGCTGCGATCGGGGTGCTGAGTTCGATCATTTCAAAGTCGGCAAAGCCGACTTTgaaa                                                                                                                                     <  5:167697/148‑4 (MQ=255)
                      aTGGCCGCTACCGATGATGACCACCGGTTGGGCGGGACGAGCGGCAACCAGTGGTTGCGCCTGGATAACCGGCGCGGTAGCCACTGCTGCGATCGGGGTGCTGAGTTCGATCATTTCAAAGTCGGCTTTGCCGACCTTGCACTCAGGGc                                                                                                                              <  5:134378/149‑1 (MQ=255)
                       tGGCCGCTACCGATGATGACCACCGGTTGGGCGGGACGAGCGGCAACCAGTGGTTGCGCCTGGATAACCGGCGCGGTAGCCACTGCTGCGATCGGGGTGCTGAGTTCGATCATTTCAAAGTCGGCTTTGCCGACCTTGCACTCAGGGCa                                                                                                                             >  5:184595/1‑149 (MQ=255)
                         gCCGCTACCGATGATGACCACCGGTTGGGCGGGACGAGCGGCAACCAGTGGTTGCGCCTGGATAACCGGCGCGGTAGCCACTGCTGCGATCGGGGTGCTGAGTTCGATCATTTCAAAGTCGGCTTTGCCGACCTTGCACTCAGGGCATa                                                                                                                           <  6:184595/149‑1 (MQ=255)
                             cTACCGATGATGACCACCGGTTGGGCGGGACGAGCGGCAACCAGTGGTTGCGCCTGGATAACCGGCGCGGTAGCCACTGCTGCGATCGGGGTGCTGAGTTCGATCATTTCAAAGTCGGCTTTGCCGACCTTGCACTCAGGGCATAGCCa                                                                                                                       <  2:311345/149‑1 (MQ=255)
                              tACCGATGATGACCACCCGATCGGCGGCACGAGCCGCAACCAGTGGTTGCGCCTGGAAAACCGGCGCGGTAGCCACTGCTGCGATCGGCGTGCTTAGTTCGATCATTTCAAAGTCAGCTTAGAATACCTAGCACTCAGGGCAAAGCAAg                                                                                                                      >  8:251328/1‑149 (MQ=255)
                                ccGATGATGACCACCGGTTGGGCGGGACGAGCGGCAACCAGTGGTTGCGCCTGGATAACCGGCGCGGTAGCCACTGCTGCGATCGGGGTGCTGAGTTCGATCATTTCAAAGTCGGCTTTGCCGACCTTGCACTCAGGGCATAGCCAGTc                                                                                                                    >  1:252694/1‑149 (MQ=255)
                                 cGATGATGACCACCGGTTGGGCGGGACGAGCGGCAACCAGTGGTTCCCCAGGGAAAACCGCCGCGTTAGCCACTCATGCGTGCGGGGTGGTGAGTAGGAACATTTCAAAGACGGGTTTCCCCACATTGCCGTCAGGGCATTGcaggaat                                                                                                                   >  4:96050/1‑143 (MQ=255)
                                   atgatgACCACCGGTTGGGCGGGACGAGCGGCAACCAGTGGTTGCGCCTGGATAACCGGCGCGGTAGCCACTGCTGCGATCGGGGTGCTGAGTTCGATCATTTCAAAGTCGGCTTTGCCGACCTTGCACTCAGGGCATAGCCAGTCTTc                                                                                                                 <  2:127702/149‑1 (MQ=255)
                                   atgatgACCACCGGTTGGGCGGGACGAGCGGCAACCAGTGGTTGCGCCTGGATAACCGGCGCGGTAGCCACTGCTGCGATCGGGGTGCTGAGTTCGATCATTTCAAAGTCGGCTTTGCCGACCTTGCACTCAGGGCATAGCCAGTCTTc                                                                                                                 <  7:10507/149‑1 (MQ=255)
                                    tgatgaCCACCGGTTGGGCGGGACGAGCGGCAACCAGTGGTTGCGCCTGGATAACCGGCGCGGTAGCCACTGCTGCGATCGGGGTGCTGAGTTCGATCATTTCAAAGTCGGCTTTGCCGACCTTGCACTCAGGGCATAGCCAGTCTTCg                                                                                                                >  6:331677/1‑149 (MQ=255)
                                      atgaCCACCGGTTGGGCGGGACGAGCGGCAACCAGTGGTTGCGCCTGGATAACCGGCGCGGTAGCCACTGCTGCGATCGGGGTGCTGAGTTCGATCATTTCAAAGTCGGCTTTGCCGACCTTGCACTCAGGGCATAGCCAGTCTTCgg                                                                                                               <  5:380752/148‑1 (MQ=255)
                                        gaCCACCGGTTGGGCGGGTCGTTCGGCAACCAGTGGTTGGGCCTGTATAACCGGCGCGGTTTCCACTGCTGCGTTCGGGGTGCTGAGTTCGATCATTTCAAATTCGGCTTTGCCGACCTTGCACTCAGGGCATAGCCAGTCTTCGGGc                                                                                                             <  7:57975/148‑1 (MQ=255)
                                         accaccGGTTGGGCGGGACGAGCGGCAACCAGTGGTTGCGCCTGGATAACCGGCGCGGTAGCCACTGCTGCGATCGGGGTGCTGAGTTCGATCATTTCAAAGTCGGCTTTGCCGACCTTGCACTCAGGGCATAGCCAGTCTTCGGGCAc                                                                                                           <  6:226521/149‑1 (MQ=255)
                                             ccGGTTGGGCGGGACGAGCGGCAACCAGTGGTTGCGCCTGGATAACCGGCGCGGTAGCCACTGCTGCGATCGGGGTGCTGAGTTCGATCATTTCAAAGTCGGCTTTGCCGACCTTGCACTCAGGGCATAGCCAGTCTTCGGGCACGtcc                                                                                                       <  2:252694/149‑1 (MQ=255)
                                                   gggCGGGACGAGCGGCAACCAGTGGTTGAGCCTGGATAACCGGCGCGGTAGCCACTGCTGCGATCGGGGTGCTGAGTTCGATCATTTCAAAGTCGGCTTTGCCGACCTTGCACTCAGGGAATAGCCAGTCTTCGGGCACGTCCTCCCac                                                                                                 >  4:77009/1‑148 (MQ=255)
                                                    ggCGGGACGAGCGGCAACCAGTGGTTGCGCCTGGATAACCGGCGCGGTAGCCACTGCTGCGATCGGGGTGCTGAGTTCGATCATTTCAAAGTCGGCTTTACCGACCTTGCACt                                                                                                                                    >  1:63176/1‑113 (MQ=255)
                                                    ggCGGGACGAGCGGCAACCAGTGGTTGCGCCTGGATAACCGGCGCGGTAGCCACTGCTGCGATCGGGGTGCTGAGTTCGATCATTTCAAAGTCGGCTTTACCGACCTTGCACt                                                                                                                                    <  2:63176/113‑1 (MQ=255)
                                                      cgtgACGAGCGGCAACCAGTGGTTGCGGCTGGATAACCTGCGCGGTAGCCACTGCTGCGATCGGGGTGCTGAGTTCGATCATTTCTAAGTGGGCTTTGCCGACCTTGCACTCAGGGCATAGCCAGTCTTCGGGCACGTCCTCCCAGCgg                                                                                              <  7:251328/146‑1 (MQ=255)
                                                         gACGAGCGGCAACCAGTGGTTGCGCCTGGATAACCGGCGCGGTAGCCACTGCTGCGATCGGGGTGCTGAGTTCGATCATTTCAAAGTCGGCTTTGCCGACCTTGCACTCAGGGCATAGCCGGTCTTCGGGCACGTCCTCCCAGCGGGTa                                                                                           >  5:228340/1‑149 (MQ=255)
                                                         gACGAGCGGCAACCAGTGGTTGCGCCTGGATAACCGGCGCGGTAGCCACTGCTGCGATCGGGGTGCTGAGTTCGATCATTTCAAAGTCGGCTTTGCCGACCTTGCACTCAGGGCATAGCCAGTCTTCGGGCACGTCCTCCCAGCGGGTa                                                                                           >  5:267310/1‑149 (MQ=255)
                                                          aCGAGCGGCAACCAGTGGTTGCGCCTGGATAACCGGCGCGGTAGCCACTGCTGCGATCGGGGTGCTGAGTTCGATCATTTCAAAGTCGGCTTTGCCGACCTTGCACTCAGGGCATAGCCAGTCTTCGGGCACGTCCTCCCAGc                                                                                                >  5:293277/1‑143 (MQ=255)
                                                          aCGAGCGGCAACCAGTGGTTGCGCCGGGATAACCGGCGCGGTAGCCACTGCTGCGATCGGGGTGCTGAGTTCGATCATTTCAAAGTCGGCTTTGCCGACCTTGCACTCAGGGCATAGCCAGTCTTCGGGCACGTCCTCCCAGc                                                                                                <  6:293277/143‑1 (MQ=255)
                                                               gcgCAACCACTGGTTGCGCCTGGATAACCGGCGCGGTAGCCACTGCTGCGATCGGGGTGCTGAGTTCGATCATTTCAAAGTCGGCTTTGCCGACCTTGCACTCAGGGCATAGCCAGTCTTCGGGCACGTCCTCCCAGCGGGTACCgggg                                                                                     >  7:105294/3‑149 (MQ=255)
                                                               cGGCAACCAGTGGTTGCGCCTGGATAACCGGCGCGGTAGCCACTGCTGCGATCGGGGTGCTGAGTTCGATCATTTCAAAGTCGGCTTTGCCGACCTTGCACTCAGGGCATAGCCAGTCTTCGGGCACGTCCTCCCAGCGGGTACCgggg                                                                                     >  1:82877/1‑149 (MQ=255)
                                                               cGGCAACCAGTGGTTGCGCCTGGATAACCGGCGCGGTAGCCACTGCTGCGATCGGGGTGCTGAGTTCGATCATTTCAAAGTCGGCTTTGCCGACCTTGCACTCAGGGCATAGCCAGTCTTCGGGCACGTCCTCCCAGCGGGTACCgggg                                                                                     >  7:73483/1‑149 (MQ=255)
                                                                ggCAACCAGTGGTTGCGCCTGGATAACCGGCGCGGTAGCCACTGCTGCGATCGGGGTGCTGAGTTCGATCATTTCAAAGTCGGCTTTGCCGACCTTGCACTCAGGGCATAGCCAGTCTTCGGGCACGTCCTCCCAGCGGGTACCggggg                                                                                    <  8:129425/149‑1 (MQ=255)
                                                                ggCAACCAGTGGTTGCGCCTGGATAACCGGCGCGGTAGCCACTGCTGCGATCGGGGTGCTGAGTTCGATCATTTCAAAGTCGGCTTTGCCGACCTTGCACTCAGGGCATAGCCAGTCTTCGGGCACGTCCTCCCAGCGGGTACCggggg                                                                                    <  1:180515/149‑1 (MQ=255)
                                                                  cAACCAGTGGTTGCGCCTGGATAACCGGCGCGGTAGCCACTGCTGCGATCGGGGTGCTGAGTTCGATCATTTCAAAGTCGGCTTTGCCGACCTTGCACTCAGGGCATAGCCAGTCTTCGGGCACGTCCTCCCAGCGGGTACCGGGGGc                                                                                   >  4:30893/1‑148 (MQ=255)
                                                                   aaCCAGTGGTTGCGCCTGGATAACCGGCGCGGTAGCCACTGCTGCGATCGGGGTGCTGAGTTCGATCATTTCAAAGTCGGCTTTGGCGACCTTGCACTCAGGGCATAGCCAGTCTTCGGGCACGTCCTCCCAGCGGGTACCGGGGGCaa                                                                                 >  1:162183/1‑149 (MQ=255)
                                                                   aaCCAGTGGTTGCGCCTGGATAACCGGCGCGGTAGCCACTGCTGCGATCGGGGTGCTGAGTTCGATCATTTCAAAGTCGGCTTTGCCGACCTTGCACTCAGGGCATAGCCAGTCTTCGGGCACGTCCTCCCAGc                                                                                                <  6:262872/134‑1 (MQ=255)
                                                                   aaCCAGTGGTTGCGCCTGGATAACCGGCGCGGTAGCCACTGCTGCGATCGGGGTGCTGAGTTCGATCATTTCAAAGTCGGCTTTGCCGACCTTGCACTCAGGGCATAGCCAGTCTTCGGGCACGTCCTCCCAGc                                                                                                >  5:262872/1‑134 (MQ=255)
                                                                     ccAGTGGTTGCGCCTGGATAACCGGCGCGGTAGCCACTGCTGCGATCGGGGTGCTGAGTTCGATCATTTCAAAGTCGGCTTTGCCGACCTTGCACTCAGGGCATAGCCAGTCTTCGGGCACGTCCTCCCAGCGGGTACCAGGGGCAATg                                                                               <  6:228340/149‑1 (MQ=255)
                                                                            ttGCGCCTGGATAACCGGCGCGGTAGCCACTGCTGCGATCGGGGTGCTGAGTTCGATCATTTCAAAGTCGGCTTTGCCGACCTTGCACTCAGGGCATAGCCAGTCTTCGGGCACGTCCTCCCAGCGGGTACCGGGGGCAATGCcgtcgt                                                                        <  7:128032/149‑1 (MQ=255)
                                                                                gcCTGGATAACCGGCGCGGTAGCCACTGCTGCGATCGGGGTGCTGAGTTCGATCATTTCAAAGTCGGCTTTGCCGACCTTGCACTCAGGGCATAGCCAGTCTTCGGGCACGTCCTCCCAGCGGGTACCGGGGGCAATGCCGTCGTCggg                                                                    >  3:196710/1‑149 (MQ=255)
                                                                                     gATAACCGGCGCGGTAGCCACTGCTGCGATCGGGGTGCTGAGTTCGATCATTTCAAAGTCGGCTTTGCCGACCTTGCACTCAGGGCATAGCCAGTCTTCGGGCACGTCCTCCCAGCGGGTACCGGGGGCAATGCcg                                                                            <  2:233206/136‑1 (MQ=255)
                                                                                     gATAACCGGCGCGGTAGCCACTGCTGCGATCGGGGTGCTGAGTTCGATCATTTCAAAGTCGGCTTTGCCGACCTTGCACTCAGGGCATAGCCAGTCTTCGGGCACGTCCTCCCAGCGGGTACCGGGGGCAATGCcg                                                                            >  1:233206/1‑136 (MQ=255)
                                                                                        aaCCGGCGCGGTAGCCACTGCTGCGATCGGGGTGCTGAGTTCGATCATTTCAAAGTCGGCTTTGCCGACCTTGCACTCAGGGCATAGCCAGTCTTCGGGCACGTCCTCCCAGCGGGTACCGGGGGCAATGCCGTCGTCGGGCCAACCtt                                                            >  6:288211/1‑149 (MQ=255)
                                                                                         aCCGGGGGGGTAGCCACTGCTTCGATCGGGGTGCTTAGTTCTATCATTTCAAAGTCGGCTTTGCCGACCTTGCACTCAGTGCATAGCCAGTCTTCCGGCACGTCCTCCCAGCGGGTACCGGGGGCAATGCCGTCGTCGGGCCAACCttt                                                           <  3:77009/149‑1 (MQ=255)
                                                                                         aCCGGCGCGGTAGCCACTGCTGCGATCGGGGTGCTGAGTTCGATCATTTCAAAGTCGGCTTTGCCGACCTTGCACTCAGGGCATAGCCAGTCTTCGGGCACGTCCTCCCAGCGGGTACCGGGGGCAATGCCGTCGTCGGGCCAACCttt                                                           <  3:30893/149‑1 (MQ=255)
                                                                                         aCCGGCGCGGTAGCCACTGCTGCGATCGGGGTGCTGAGTTCGATCATTTCAAAGTCGGCTTTGCCGACCTTGCACTCAGGGCATAGCCAGTCTTCGGGCACGTCCTCCCAGCGGGTACCGGGGGCAATGCCGTCGTCGGGCCAACCttt                                                           <  8:291483/149‑1 (MQ=255)
                                                                                          ccGGCGCGGTAGCCACTGCTGCGATCGGGGTGCTGAGTTCGATCATTTCAAAGTCGGCTTTGCCGACCTTg                                                                                                                                        <  1:193751/71‑1 (MQ=255)
                                                                                          ccGGCGCGGTAGCCACTGCTGCGATCGGGGTGCTGAGTTCGATCATTTCAAAGTCGGCTTTGCCGACCTTg                                                                                                                                        >  2:193751/1‑71 (MQ=255)
                                                                                          ccGGCGCGGTAGCCACTGCTGCGATCGGGGTGCTGAGTTCGATCATTTCAAAGTCGGCTTTGCCGACCTTGCACTCAGGGCATAGCCAGTCTTCGGGCACGTCCTCCCAGCGGGTACCGGGGGCAATGCCGTCGTCGGGCCAACCTTTg                                                          <  8:48060/149‑1 (MQ=255)
                                                                                          ccGGCGCGGTAGCCACTGCTGCGATCGGGGTGCTGAGTTCGATCATTTCAAAGTCGGCTTTGCCGACCTTGCACTCAGGGCATAGCCAGTCTTCGGGCACGTCCTCCCAGCGGGTACCGGGGGCAATGCCCTCGTCGGGCCAACCTTTg                                                          >  2:105030/1‑149 (MQ=255)
                                                                                          ccGGCGCGGTAGCCACTGCTGCGATCGCGGTGCTGAGTTCGATCATTTCAAAGTCGGCTTTGCCGACCTTGCACTCAGGGCATAGCCAGTCTTCGGGCACGTCCTCCCAGCGGGTACCGGGGGCAATGCCGTCGTCGGGCCAACCTTTg                                                          >  3:116509/1‑149 (MQ=255)
                                                                                               gcgGTAGCCACTGCTGCGATCGGGGTGCTGAGTTCGATCATTTCAAAGTCGGCTTTGGCGACCTTGCACTCAGGGCATAGCCAGTCTTCGGGCACGTCCTCCCAGCGGGTACCGGGGGCAATGCCGTCGTCGGGCCAACCTTTGGCTtc                                                     <  2:162183/149‑1 (MQ=255)
                                                                                                       cACTGCTGCGATCGGGGTGCTGAGTTCGATCATTTCAAAGTCGGCTTTGCCGACCTTGCACTCAGGGCATAGCCAGTCTTCGGGCACGTCCTCCCAGCGGGTACCGGGGGCAATGCCGTCGTCGGGCCAACCTTTGGCTTCGTCGTAGa                                             <  5:331677/149‑1 (MQ=255)
                                                                                                        aCTGCTGCGATCGGGGTGCTGAGTTCGATCATTTCAAAGTCGGCTTTGCCGACCTTGCACTCAGGGCATAGCCAGTCTTCGGGCACGTCCTCCCAGCGGGTACCGGGGGCAATGCCGTCGTCGGGCCAACCTTTGGCTtcgtcg                                                 >  1:50050/1‑144 (MQ=255)
                                                                                                          tgctgcGATCGGGGTGCTGAGTTCGATCATTTCAAAGTCGGCTTTGCCGACCTTGCACTCAGGGCATAGCCAGTCTTCGGGCACGTCCTCCCAGCGGGTACCGGGGGCAATGCCGTCGTCGGGCCAACCTTTGGCTTCGTCGTAGATCa                                          >  1:159996/1‑149 (MQ=255)
                                                                                                            ctgcGATCGGGGTGCTGAGTTCGATCATTTCAAAGTCGGCTTTGCCGACCTTGCACTCAGGGCATAGCCAGTCTTCGGGCACGTCCTCCCAGCGGGTACCGGGGGCAATGCCGTCGTCGGGCCAACCTTTGGCTTCGTCGTAGATCa                                          <  8:73483/147‑1 (MQ=255)
                                                                                                            ctgcGATCGGGGTGCTGAGTTCGATCATTTCAAAGTCGGCTTTGCCGACCTTGCACTCAGGGCATAGCCAGTCTTCGGGCACGTCCTCCCAGCGGGTACCGGGGGCAATGCCGTCGTCGGGCCAACCTTTGGCTTCGTCGTAGATCAGa                                        >  1:57593/1‑149 (MQ=255)
                                                                                                             tgcGATCGCGGTGCTGAGTTCGATCATTTCAAAGTCGGCTTTGCCGACCTTGCACTCAGGGCATAGCCAGTCTTCGGGCACGTCCTCCCAGCGGGTACCCGGGGCAATGCCGTCGTCGGGCCAACCTTTGGCTTCGTCGTAGATCAGAc                                       >  1:293475/1‑149 (MQ=255)
                                                                                                               cGATCGGGGTGCTGAGTTCGATCATTTCAAAGTCGGCTTTGCCGACCTTGCACTCAGGGCATAGCCAGTCTTCGGGCACGTCCTCCCAGCGGGTACCGGGGGCAATGCCGTCGTCGGGCCAACCTTTGGCTTCGTCGTAGATCAGAc                                       <  2:82877/147‑1 (MQ=255)
                                                                                                                gATCGGGGTGCTGAGTTCGATCATTTCAAATTCGGCTTTGCCGACGTTGCACTCAGGGCATAGCCAGTCTTCGGGCACGTGCTCCCAGCGGGTACCGGGGGCAATGCCGTCGTCGGGCCAACCTTTGGCTTCGTCGTAGATCAGACca                                     <  2:293475/148‑1 (MQ=255)
                                                                                                                     gggTGCTGAGTTCGATCATTTCAAAGTCGGCTTTGCCGACCTTGCACTCAGGGCATAGCCAGTCTTCGGGCACGTCCTCCCAGCGGGTACCGGGGGCAATGCCGTCGTCGGGCCAACCTTTGGCTTCGTCGTAGATCAGACCACAgat                                >  4:183852/1‑148 (MQ=255)
                                                                                                                     gggTGCTGAGTTCGATCATTTCAAAGTCGGCTTTGCCGACCTTGCACTCAGGGCATAGCCAGTCTTCGGGCACGTCCTCCCAGCGGGTACCGGGGGCAATGCCGTCGTCGGGCCAACCTTTGGCTTCGTCGTAGAGCAGAACACAgat                                >  6:168030/1‑148 (MQ=255)
                                                                                                                     gggTGCTGAGTTCGATCATTTCAAAGTCGGCTTTGCCGACCTTGCACTCAGGGCATAGCCAGTCTTCGGGCACGTCCTCCCAGCGGGTACCCGGGGCAATGCCGTCGTCGGGCCAACCTTTGGCTTCGTCGTa                                               >  6:261097/1‑133 (MQ=255)
                                                                                                                     gggTGCTGAGTTCGATCATTTCAAAGTCGGCTTTGCCGACCTTGCACTCAGGGCATAGCCAGTCTTCGGGCACGTCCTCCCAGCGGGTACCCGGGGCAATGCCGTCGTCGGGCCAACCTTTGGCTTCGTCGTa                                               <  5:261097/133‑1 (MQ=255)
                                                                                                                        tGCTGAGTTCGATCATTTCAAAGTCGGCTTTGCCGACCTTGCACTCAGGGCATAGCCAGTCTTCGGGCACGTCCTCCCAGCGGGTACCGGGGGCAATGCCGTCGTCGGGCCAACCTTTGGCTTCGTCGTAGATCAGACCACAGatgatg                            <  5:84044/149‑1 (MQ=255)
                                                                                                                          cTGAGGTCGATCATTTCAAAGTCGGCTTTGCCGACCTTGCACTCAGGGCATAGCCAGTCTTCGGGCACGTCCTCCCAGCGGGTACCGGGGGCAATGCCGTCGTCGGGCCAACCTTTGGCTTCGTCGTAGATCAGACCACAGATGATGc                           <  2:245634/148‑1 (MQ=255)
                                                                                                                           tGAGTCCGATCATTTCAAAGTCGGCTTTGCCGACCTTGCACTCAGGGCATAGCCAGTCTTCGGGCACGTCCTCCCAGCGGGTACCGGGGGCAATGCCGTCGTCGGGCCAACCTTTGGCTTCGTCGTAGATCAGACCACAGATGATGCag                         >  6:32312/1‑149 (MQ=255)
                                                                                                                              gTTCGATCATTTCAAAGTCGGCTTTGCCGACCTTGCACTCAGGGCATAGCCAGTCTTCGGGCACGTCCTCCCAGCGGGTACCGGGGGCAATGCCGTCGTCGGGCCAACCTTTGGCTTCGTCGTAGATCAGACCACAGATGATGCAGAGc                      <  4:196710/149‑1 (MQ=255)
                                                                                                                              gTTCGATCATTTCAAAGTCGGCTTTGCCGACCTTGCACTCAGGGCATAGCCAGTCTTCGGGCACGTCCTCCCAGCGGGTACCGGGGGCAATGCCGTCGTCGGGCCAACCTTTGGCTTCGTCGTAGATCAGACCACAGATGATGCAGAGc                      <  6:267310/149‑1 (MQ=255)
                                                                                                                               ttCGATCATTTCAAAGTCGGCTTTGCCGACCTTGCACTCAGGGCATAGCCAGTCTTCGGGCACGTCCTCCCAGCGGGTACCGGGGGCAATGCCGTCGTCGGGCCAACCTTTGGCTTCGTCGTAGATCAGACCACAGATGATGCAGAGc                      >  5:156460/1‑148 (MQ=255)
                                                                                                                                 cGATCATTTCAAAGTCGGCTTTGCCGACCTTGCACTCAGGGCATAGCCAGTCTTCGGGCACGTCCTCCCAGCGGGTACCGGGGGCAATGCCGTCGTCGGGCCAACCTTTGGCTTCGTCGTAGATCAGACCACAGATGATGCAGAGCCAg                   >  8:111644/1‑149 (MQ=255)
                                                                                                                                  gATCATTTAAAAGTCGGCTTTGCCGACCTTGCACTCAGGGCATAGCCAGTCTTCGGGCACGTCCTCCCAGCGGGTACCGGGGGCAATGCCGTCGTCGGGCCAACCTTTGTCTTCGTCGTAGATAAGACCACAGATGATGCAGAGCCAgg                  >  4:101412/1‑149 (MQ=255)
                                                                                                                                   ctgagtgcAAGGTCGGCAAAGCCGACTTTGCACTCAGGGCATAGCGAGTCTTCGGGCACGTCCTCCCAGCGGGTACCGGGGGCAATGCCGTCGTCGGGCCAACCTTTGGCTTCGTCGTAGATCAGACCACAGATGATGCAGAGCCAGGt                 >  1:266875/8‑149 (MQ=255)
                                                                                                                                   ctgagtgcAAGGTCGGCAAAGCCGACTTTGCACTCAGGGCATAGCCAGTCTTCGGGCACGTCCTCCCAGCGGGTACCGGGGGCAATGCCGTCGTCGGGCCAACCTTTGGCTTCGTCGTAGATCAGACCACAGATGATGCAGAGCCAGGt                 >  1:343709/8‑149 (MQ=255)
                                                                                                                                   aTCATTTCAAAGTCGGCTTTGCCGACCTTGCACTCAGGGCATAGCCAGTCTTCGGGCACGTCCTCCCAGCGGGTACCGGGGGCAATGCCGTCGTCGGGCCAACCTTTGGCTTCGTCGTAGATCAGACCACAGATGATGCAGAGCCAGGt                 >  3:365755/1‑149 (MQ=255)
                                                                                                                                     cATTTCAAAGTCGGCTTTGCCGACCTTGCACTCAGGGCATAGCCAGTCTTCGGGCACGTCCTCCCAGCGGGTACCGGGGGCAATGCCGTCGTCGGGCCAACCTTTGGCTTCGTCGTAGATCAGACCACAGatgatg                            <  4:325970/136‑1 (MQ=255)
                                                                                                                                     cATTTCAAAGTCGGCTTTGCCGACCTTGCACTCAGGGCATAGCCAGTCTTCGGGCACGTCCTCCCAGCGGGTACCGGGGGCAATGCCGTCGTCGGGCCAACCTTTGGCTTCGTCGTAGATCAGACCACAGatgatg                            >  3:325970/1‑136 (MQ=255)
                                                                                                                                       tttCAAAGTCGGCTTTGCCGACCTTGCACTCAGGGCATAGCCAGTCTTCGGGCACGTCCTCCCAGCGGGTACCGGGGGCAATGCCGTCGTCGGGCCAACCTTTGGCTTCGTCGTAGATCAGACCACAGATGATGCAGAGCCAGGtttt              >  6:103238/1‑148 (MQ=255)
                                                                                                                                       gtgcAAGGTCGGCAAAGCCGACTTTGCACTCAGGGCATAGCCAGTCTTCGGGCACGTCCTCCCAGCGGGTACCGGGGGCAATGCCGTCGTCGGGCCAACCTTTGGCTTCGTCGTAGATCAGACCACAGATGATGCAGAGCCAGGTTTTc             <  2:266875/146‑1 (MQ=25)
                                                                                                                                       gtgcAAGGTCGGCAAAGCCGACTTTGCACTCAGGGCATAGCCAGTCTTCGGGCACGTCCTCCCAGCGGGTACCGGGGGCAATGCCGTCGTCGGGCCAACCTTTGGCTTCGTCGTAGATCAGACCACAGATGATGCAGAGCCAGGTTTTc             <  2:343709/146‑1 (MQ=25)
                                                                                                                                         tCAAAGTCGGCTTTGCCGACCTTGCACTCAGGGCATAGCCAGTCTTCGGGCACGTCCTCCCAGCGGGTACCGGGGGCAATGCCGTCGTCGGGCCAACCTTTGGCTTCGTCGTAGATCAGACCACAGATGATGCAGAGCCAGGTTTTCa            <  2:50050/148‑1 (MQ=255)
                                                                                                                                         tCAAAGTCGGCTTTGCCGACCTTGCACTCAGGGCATAGCCAGTCTTCGGGCACGTCCTCCCAGCGGGTACCGGGGGCAATGCCGTCGTCGGGCCAACCTTTGGCTTCGTCGTAGATCAGACCACAGATGATGCAGAGCCAGGTTTTCa            <  4:365755/148‑1 (MQ=255)
                                                                                                                                         tCAAAGTCGGCTTTGCCGACCTTGCACTCAGGGCATAGCCAGTCTTCGGGCACGTCCTCCCAGCGGGTACCGGGGGCAATGCCGTCGTCGGGCCAACCTTTGGCTTCGTCGTAGATCAGACCACAGATGATGCAGAGCCAGGTTTTCAt           >  3:22787/1‑149 (MQ=255)
                                                                                                                                         tCAAAGTCGGCTTTGCCGACCTTGCACTCAGGGCATAGCCAGTCTTCGGGCACGTCCTCCCAGCGGGTACCGGGGGCAATGCCGTCGTCGGGCCAACCTTTGGCTTCGTCGTAGATCAGACCACAGATGATGCAGAGCCAGGTTTTCAt           >  7:185550/1‑149 (MQ=255)
                                                                                                                                          cAAAGTCGGCTTTGCCGACCTTGCACTCAGGGCATAGCCAGTCTTCGGGCACGTCCTCCCAGCGGGTACCGGGGGCAATGCCGTCGTCGGGCCAACCTTTGGCTtcgtcg                                                 >  5:3873/1‑110 (MQ=255)
                                                                                                                                          cAAAGTCGGCTTTGCCGACCTTGCACTCAGGGCATAGCCAGTCTTCGGGCACGTCCTCCCAGCGGGTACCGGGGGCAATGCCGTCGTCGGGCCAACCTTTGGCTtcgtcg                                                 <  6:3873/110‑1 (MQ=255)
                                                                                                                                           aaaGTCGGCTTTGCCGACCTTGCACTCAGGGCATAGCCAGTCTTCGGGCACGTCCTCCCAGCGGGTACCGGGGGCAATGCCGTCGTCGGGCCAACCTTTGGCTTCGTCGTAGATCAGACCACAGATGATGCAGAGCCAGGTTTTCATTa         <  3:96050/149‑1 (MQ=255)
                                                                                                                                            aaGTCGGCTTTGCCGACCTTGCACTCAGGGCATAGCCAGTCTTCGGGCACGTCCTCCCAGCGGGTACCGGGGGCAATGCCGTCGTCGGGCCAACCTTTGGCTTCGTCGTAGATCAGACCACAGATGATGCAGAGCCAGGTTTTCATTAt        >  1:331621/1‑149 (MQ=255)
                                                                                                                                             aGTCGGCTTTGCCGACCTTGCACTCAGGGCATAGCCAGTCTTCGGGCACGTCCTCCCAGCGGGTACCGGGGGCAATGCCGTCGTCGGGCCAACCTTTGGCTTCGTCGTAGATCAGACCACAGATGATGCAGAGCCAGGTTTTCATTAtg       <  6:156460/149‑1 (MQ=255)
                                                                                                                                              gTCGGCTTTGCCGACCTTGCACTCAGGGCATAGCCAGTCTTCGGGCACGTCCTCCCAGCGGGTACCGGGGGCAATGCCGTCGTCGGGCCAACCTTTGGCTTCGTCGTAGATCAGACCACAGATGATGCAGAGCCAGGTTTTCATTAtgc      >  6:226837/1‑149 (MQ=255)
                                                                                                                                               tCGGCTTTGCCGACCTTGCACTCAGGGCATAGCCAGTCTTCGGGCACGTCCTCCCAGCGGGTACCGGGGGCAATGCCGTCGTCGGGCCAACCTTTGGCTTCGTCGTAGATCAGACCACAGATGATGCAGAGCCAGGTTTTCATTAtgct     <  8:105294/149‑1 (MQ=255)
                                                                                                                                               tCGGCTTTGCCGACCTTGCACTCAGGGCATAGCCAGTCTTCGGGCACGTCCTCCCAGCGGGTACCGGGGGCAATGCCGTCGTCGGGCCAACCTTTGGCTTCGTCGTAGATCAGACCACAGATGATGCAGAGCCAGGTTTTCATTAtgct     <  8:185550/149‑1 (MQ=255)
                                                                                                                                                cGGCTTTGCCGACCTTGCACTCAGGGCATAGCCAGTCTTCGGGCACGTCCTCCCAGCGGGTACCGGGGGCAATGCCGTCGTCGGGCCAACCTTTGGCTTCGTCGTAGATCAGACCACAGATGATGCAGAGCCAGGTTTTCATTAtgctg    >  1:188717/1‑149 (MQ=255)
                                                                                                                                                  gCTTTGCCGACCTTGCACTCAGGGCATAGCCAGTCTTCGGGCACGTCCTCCCAGCGGGTACCGGGGGCAATGCCGTCGTCGGGCCAACCTTTGGCTTCGTCGTAGATCAGACCACAGATGATGCAGAGCCAGGTTTTCATTATGCTGcc  >  5:150502/1‑149 (MQ=255)
                                                                                                                                                    |                                                                                                                                                  
CCTGTGCCAGGCCGTAGCCGGCATGGCCGCTACCGATGATGACCACCGGTTGGGCGGGACGAGCGGCAACCAGTGGTTGCGCCTGGATAACCGGCGCGGTAGCCACTGCTGCGATCGGGGTGCTGAGTTCGATCATTTCAAAGTCGGCTTTGCCGACCTTGCACTCAGGGCATAGCCAGTCTTCGGGCACGTCCTCCCAGCGGGTACCGGGGGCAATGCCGTCGTCGGGCCAACCTTTGGCTTCGTCGTAGATCAGACCACAGATGATGCAGAGCCAGGTTTTCATTATGCTGCC  >  NC_002947/6121864‑6122158

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