Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 2,651,276 G→C 13.9% intergenic (+41/+29) PP_2323 → / ← aroF‑I hypothetical protein/tyr‑sensitive phospho‑2‑dehydro‑3‑deoxyheptonate aldolase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029472,651,2760GC13.9% 183.9 / 18.4 71intergenic (+41/+29)PP_2323/aroF‑Ihypothetical protein/tyr‑sensitive phospho‑2‑dehydro‑3‑deoxyheptonate aldolase
Reads supporting (aligned to +/- strand):  ref base G (34/27);  new base C (4/6);  total (38/33)
Fisher's exact test for biased strand distribution p-value = 4.97e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GCTCTGGAGTACGCCGAGCAGATGGGGTACACGGTGATCCCGTCTTGCTCATATGTGGAGCGTTACATGGAGCGCCAGCAGCGGCATTCGAGCAAGGCCTGAGCTGACGACAGGCACCAGACACAACAAAAAACGCCGGGCAGTGCCCGGCGTTTTTCATTCCCGCAATGTTCAGCCGCGCTTACGCTGCGGCAGCACATCCTTGAGCTTGGCATGCATGCTGCGCAGGGTTTTCTCGGTAGCCGACCAGTCGATGCAGGCGTCCGTGATCGACACGCCATA  >  NC_002947/2651134‑2651415
                                                                                                                                              |                                                                                                                                           
gCTCTGGAGTACGCCGAGCAGATGGGGTACACGGTGATCCCGTCTTGCTCATATGTGGAGCGTTACATGGAGCGCCAGCAGCGGCATTCGAGCAAGGCCTGAGCTGACGACAGGCACCAGACACAACAAAAAACGCCGGGCAGTGCCCg                                                                                                                                       <  3:42890/149‑1 (MQ=255)
      gAGTACGCCGAGCAGATGGGGTACACGGTGATCCCGTCTTGCTCATATGTGGAGCGTTACATGGAGCGCCAGCAGCGGCATTCGAGCAAGGCCTGAGCTGACGACAGGCACCAGACACAACAAAAAACGCCGGGCAGTGCCCGGCGttt                                                                                                                                 >  2:157609/1‑149 (MQ=255)
      gAGTACGCCGAGCAGATGGGGTACACGGTGATCCCGTCTTGCTCATATGTGGAGCGTTACATGGAGCGCCAGCAGCGGCATTCGAGCAAGGCCTGAGCTGACAACAGGCACCAGACACAACAAAAAACGCCGGGCACTGCCCGGCGttt                                                                                                                                 <  5:308198/149‑1 (MQ=255)
       aGTACGCCGAGCAGATGGGGTACACGGTGATCCCGTCTTGCTCATATGTGGAGCGTTACATGGAGCGCCAGCAGCGGCATTCGAGCAAGGCCTGAGCTGACGACAGGCACCAGACACAACAAAAAACGCCGGGCACTGCCCGGCGtttt                                                                                                                                <  8:184434/149‑1 (MQ=255)
        gTACGCCGAGCAGATGGGGTACACGGTGATCCCGTCTTGCTCATATGTGGAGCGTTACATGGAGCCCCAGCAGCGGCATTCGAGCAAGGCCTGAGCTGACGACAGGCACCAGACACAACAAAAAACGCCGGGCACTGCCCGGCGttttt                                                                                                                               <  4:370709/149‑1 (MQ=255)
         tACGCCGAGCAGATGGGGTACACGGTGATCCCGTCTTGCTCATATGTGGAGCGTTACATGGAGCGCCAGCAGCGGCATTCGAGCAAGGCCTGAGCTGACGACAGGCACCAGACACAACAAAAAACGCCGGGCAGTGCCCGGCGTTTTTc                                                                                                                              <  2:11366/149‑1 (MQ=255)
           cGCCGAGCAGATGGGGTACACGGTGATCCCGTCTTGCTCATATGTGGAGCGTTACATGGAGCGCCAGCAGCGGCATTCGAGCAAGGCCTGAGCTGACGACAGGCACCAGACACAACAAAAAACGCCGGGCACTGCCCGGCGttttttg                                                                                                                             <  6:184131/148‑3 (MQ=255)
            gCCGAGCAGATGGGGTACACGGTGATCCCGTCTTGCTCATATGTGGAGCGTTACATGGAGCGCCAGCAGCGGCATTCGAGCAAGGCCTGAGCTGACGACAGGCACCAGACACAACAAAAAACGCCGGGCAGTGCCCGGCGTTTTTCAtt                                                                                                                           >  5:84421/1‑149 (MQ=255)
                   aGATGGGGTACACGGTGATCCCGTCTTGCTCATATGTGGAGCGTTACATGGAGCGCCAGCAGCGGCATTCGAGCAAGGCCTGAGCTGACGACAGGCACCAGACACAACAAAAAACGCCGGGCAGTGCCCGGCGTTTTTCATTCCCGCaa                                                                                                                    <  1:157609/149‑1 (MQ=255)
                                    aTCCCGTCTTGCTCATATGTGGAGCGTTACATGGAGCGCCAGCAGCGGCATTCGAGCAAGGCCTGAGCTGACGACAGGCACCAGACACAACAAAAAACGCCGGGCAGTGCCCGGCGTTTTTCATTCCCGCAATGTTCAGCCGCGCTTa                                                                                                    <  6:84421/148‑1 (MQ=255)
                                        cGTCTTGCTCATATGTGGAGCGTTACATGGAGCGCCAGCAGCGGCATTCGAGCAAGGCCTGAGCTGACGACAGGCACCAGACACAACAAAAAACGCCGGGCAGTGCCCGGCGTTTTTCATTCCCGCAATGTTCAGCCGCGCTTACGCTg                                                                                               >  5:257905/1‑149 (MQ=255)
                                              gCTCATATGTGGAGCGTTACATGGAGCGCCAGCAGCGGCATTCGAGCAAGGCCTGAGCTGACGACAGGCACCAGACACAACAAAAAACGCCGGGCAGTGCCCGGCGTTTTTCATTCCCGCAATGTTCAGCCGCGCTTACGCTGCGgcag                                                                                         >  6:71972/1‑149 (MQ=255)
                                              gCTCATATGTGGAGCGTTACATGGAGCGCCAGCAGCGGCATTCGAGCAAGGCCTGAGCTGACGACAGGCACCAGACACAACAAAAAACGCCGGGCAGTGCCCGGCGTTTTTCATTCCCGCAATGTTCAGCCGCGCTTACGCTGCGgcag                                                                                         >  1:32405/1‑149 (MQ=255)
                                               cTCATATGTGGAGCGTTACATGGAGCGCCAGCAGCGGCATTCGAGCAAGGCCTGAGCTGACGACAGGCACCAGACACAACAAAAAACGCCGGGCAGTGCCCGGCGTTTTTCATTCCCGCAATGTTCAGCCGCGCTTACGCTGCGgcag                                                                                         <  8:349429/148‑1 (MQ=255)
                                                tCATATGTGGAGCGTTACATGGAGCGCCAGCAGCGGCATTCGAGCAAGGCCTGAGCTGACGACAGGCACCAGACACAACAAAAAACGCCGGGCAGTGCCCGGCGTTTTTCATTCCCGCAATGTTCAGCCGCGCTTACGCTGCGgcagca                                                                                       >  5:200714/1‑149 (MQ=255)
                                                  atatGTGGAGCGTTACATGGAGCGCCAGCAGCGGCATTCGAGCAAGGCCTGAGCTGACGACAGGCACCAGACACAACAAAAAACGCCGGGCAGTGCCCGGCGTTTTTCATTCCCGCAATGTTCAGCCGCGCTTACGCTGCGGCAGcaca                                                                                     >  2:278852/1‑149 (MQ=255)
                                                  atatGTGGAGCGTTACATGGAGCGCCAGCAGCGGCATTCGAGCAAGGCCTGAGCTGACGACAGGCACCAGACACAACAAAAAACGCCGGGCAGTGCCCGGCGTTTTTCATTCCCGCAATGTTCAGCCGCGCTTACGCTGCGGCAGcaca                                                                                     <  4:327027/149‑1 (MQ=255)
                                                  atatGTGGAGCGTTACATGGAGCGCCAGCAGCGGCATTCGAGCAAGGCCTGAGCTGACGACAGGCACCAGACACAACAAAAAACGCCGGGCAGTGCCCGGCGTTTTTCATTCCCGCAATGTTCAGCCGCGCTTACGCTGCGGCAGcaca                                                                                     >  3:327027/1‑149 (MQ=255)
                                                  atatGTGGAGCGTTACATGGAGCGCCAGCAGCGGCATTCGAGCAAGGCCTGAGCTGACGACAGGCACCAGACACAACAAAAAACGCCGGGCAGTGCCCGGCGTTTTTCATTCCCGCAATGTTCAGCCGCGCTTACGCTGCGGCAGcac                                                                                      >  4:130116/1‑148 (MQ=255)
                                                  atatGTGGAGCGTTACATGGAGCGCCAGCAGCGGCATTCGAGCAAGGCCTGAGCTGACGACAGGAACCAGACACAACAAAAAACGACGGGAAGTGCCCGGCGTTTTTCATTCCCGCAATGTTCAGCCGCGCTTACGCTGCGGCAGcac                                                                                      >  7:255960/1‑148 (MQ=255)
                                                        ggAGCGTTACATGGAGCGCCAGCAGCGGCATTCGAGCAAGGCCTGAGCTGACGACAGGCACCAGACACAACAAAAAACGCCGGGCAGTGCCCGGCGTTTTTCATTCCCGCAATGTTCAGCCGCGCTTACGCTGCGGCAGCACATCCtt                                                                                >  1:185832/1‑148 (MQ=255)
                                                            cGTTACATGGAGCGCCAGCAGCGGCATTCGAGCAAGGCCTGAGCTGACGACAGGCACCAGACACAACAAAAAACGCCGGGCAGTGCCCGGCGTTTTTCATTCCCGCAATGTTCAGCCGCGCTTACGCTGCGGCAGCACATCCTTGAGCt                                                                           >  1:320932/1‑149 (MQ=255)
                                                              ttACATGGAGCGCCAGCAGCGGCATTCGAGCAAGGCCTGAGCTGACGACAGGCACCAGACACAACAAAAAACGCCGGGCAGTGCCCGGCGTTTTTCATTCCCGCAATGTTCAGCCGCGCTTACGCTGCGGCAGCACATCCTTGAGCTTg                                                                         >  5:124459/1‑149 (MQ=255)
                                                                aCATGGAGCGCCAGCAGCGGCATTCGAGCAAGGCCTGAGCTGACGACAGGCACCAGACACAACAAAAAACGCCGGGCAGTGCCCGGCGTTTTTCATTCCCGCAATGTTCAGCCGCGCTTACGCTGc                                                                                              >  1:339154/1‑126 (MQ=255)
                                                                aCATGGAGCGCCAGCAGCGGCATTCGAGCAAGGCCTGAGCTGACGACAGGCACCAGACACAACAAAAAACGCCGGGCAGTGCCCGGCGTTTTTCATTCCCGCAATGTTCAGCCGCGCTTACGCTGc                                                                                              <  2:339154/126‑1 (MQ=255)
                                                                  aTGGAGCGCCAGCAGCGGCATTCGAGCAAGGCCTGAGCTGACGACAGGCACCAGACACAACAAAAAACGCCGGGCAGTGCCCGGCGTTTTTCATTCCCGCAATGTTCAGCCGCGCTTACGCTGCGGCAGCACATCCTTGAGCTTGgcat                                                                     >  6:137925/1‑149 (MQ=255)
                                                                        cgcCAGCAGCGGCATTCGAGCAAGGCCTGAGCTGACGACAGGCACCAGACACAACAAAAAACGCCGGGCAGTGCCCGGCGTTTTTCAGTCCCGCAATGTTCAGCCGCGCTTACGCTGCGGCAGCACATCCTTGAGCTTGGCatgcatgc                                                               >  5:131428/1‑149 (MQ=255)
                                                                           cagcagCGGCATTCGAGCAAGGCCTGAGCTGACGACAGGCACCAGACACAACAAAAAACGCCGGGCAGTGCCCGGCGTTTTTCATTCCCGCAATGTTCAGCCGCGCTTACGCTGCGGCAGCACATCCTTGAGCTTGGCATGCAtgctgc                                                            >  1:374394/1‑149 (MQ=255)
                                                                            agcagcGGCATTCGAGCAAGGCCTGAGCTGACGACAGGCACCAGACACAACAAAAAACGCCGGGCAGTGCCCGGCGTTTTTCATTCCCGCAATGTTCAGCCGCGCTTACGCTGCGGCAGCACATCCTTGAGCTTGGCATGCATGCTgcg                                                           <  6:124459/149‑1 (MQ=255)
                                                                               agcGGCATTCGAGCAAGGCCTGAGCTGACGACAGGCACCAGACACAACAAAAAACGCCGGGCAGTGCCCGGCGTTTTTCATTCCCGCAATGTTCAGCCGCGCTTACGCTGCGGCAGCACATCCTTGAGCTTGGCATGCATGCTGCGCAg                                                        <  5:137925/149‑1 (MQ=255)
                                                                                gcGGCATTCGAGCAAGGCCTGAGCTGACGACAGGCACCAGACACAACAAAAAACGCCGGGCAGTGCCCGGCGTTTTTCATTCCCGCAATGTTCAGCCGCGCTTACGCTGCGGCAGCACATCCTTGAGCTTGGCATGCATGCTGCGCAgg                                                       <  2:32405/149‑1 (MQ=255)
                                                                                      ttCGAGCAAGGCCTGAGCAGACGACAGGCACCAGCCACAACAAAAAACGCCGGGCAGTGCCCGGCGTTTTTAATTCCCGCAATGTTCAGCCGCGCTTACGCTGCGCCAGCACATCCTTGAGCTTGCCATGCATGCTGCGCAgcgattaa                                                 >  5:98660/1‑142 (MQ=255)
                                                                                       tCGAGCAAGGCCTGAGCTGACGACAGGCACCAGACACAACAAAAAACGCCGGGCAGTGCCCGGCGTTTTTCATTCCCGCAATGTTCAGCCGCGCTTACGCTGCGGCAGCACATCCTTGAGCTTGGCATGCATGCTGAGCAGGGTTTtct                                                >  8:63385/1‑149 (MQ=255)
                                                                                         gAGCAAGGCCTGAGCTGACGACAGGCACCAGACACAACAAAAAACGCCGGGCAGTGCCCGGCGTTTTTCATTCCCGCAATGTTCAGCCGCGCTTACGCTGCGGCAGCACATCCTTGAGCTTGGcatgcatg                                                                >  7:146823/1‑131 (MQ=255)
                                                                                         gAGCAAGGCCTGAGCTGACGACAGGCACCAGACACAACAAAAAACGCCGGGCAGTGCCCGGCGTTTTTCATTCCCGCAATGTTCAGCCGCGCTTACGCTGCGGCAGCACATCCTTGAGCTTGGcatgcatg                                                                <  8:146823/131‑1 (MQ=255)
                                                                                            cAAGGCCTGAGCTGACGACAGGCACCAGACACAACAAAAAACGCCGGGCAGTGCCCGGCGTTTTTCATTCCCGCAATGTTCAGCCGCGCTTACGCTGCGGCAGCACATCCTTGAGCTTGGCATGCATGCTGCGCAGGGTTTTCTCGGTa                                           <  2:185832/149‑1 (MQ=255)
                                                                                               ggCCTGAGCTGACGACAGGCACCAGACACAACAAAAAACGCCTGGCAGTGCCCGGCGTTTTTCATTCCCGCAATGTTCAGCCGCGCTTACGCTGCGGCAGCACATCCTTGAGCTTGGCATGCATGCTGCGCAGGGTTTTCTCGGTAGcc                                        <  6:257905/149‑1 (MQ=255)
                                                                                                 ccTGAGCTGACGACAGGCACCAGACACAACAAAAAACGCCGGGCAGTGCCCGGCGTTTTTCATTCCCGCAATGTTCAGCCGCGCTTACGCTGCGGCAGCACATCCTTGAGCTTGGCATGCATGCTGCGCAGGGTTTTCTCGGTAGCCGa                                      >  7:354942/1‑149 (MQ=255)
                                                                                                  cTGAGCTGACGACAGGCACCAGACACAACAAAAAACGCCGGGCAGTGCCCGGCGTTTTTCATTCCCGCAATGTTCAGCCGCGCTTACGCTGCGGCAGCACATCCTTGAGCTTGGCATGCATGCTGCGCAGGGTTTTCTCGGTAGCCGAc                                     <  6:200714/149‑1 (MQ=255)
                                                                                                     aGCTGACGACAGGCACCAGACACAACAAAACACGCCGGGCAGTGCCCGGCGTTTTTCATTCCCGCAATGTTCAGCCGCGCTTACGCTGCGGCAGCACATCCTTGa                                                                              <  7:73271/105‑1 (MQ=255)
                                                                                                     aGCTGACGACAGGCACCAGACACAACAAAAAACGCCGGGCAGTGCCCGGCGTTTTTCATTCCCGCAATGTTCAGCCGCGCTTACGCTGCGGCAGCACATCCTTGa                                                                              >  8:73271/1‑105 (MQ=255)
                                                                                                     aGCTGACGACAGGCACCAGACACAACAAAAAACGCCGGGCAGGGCCCGGCGTTTTTCATTCCCGCAATGTTCAGCCGCGCTTACGCTGCGGCAGCACATCCTTGAGCTTGGCATGCATGCTGCGCAGGGTTTTCTCGGTAGCCGACCa                                   <  5:172602/148‑1 (MQ=255)
                                                                                                     aGCTGACGACACGCACCAGACACAACAAAAAACGCCGGGCAGTGCCCGGCGTTTTTCATTCCCGCAATCTTCAGCCGCGCTTACGCTGCGGCAGCACATCCTTGAGCTTGGCATGCATGCTGCGCAGGGTTTTCTCGGTAGCCGACCa                                   <  5:178079/148‑1 (MQ=255)
                                                                                                          acgacAGGCACCAGACACAACAAAAAACGCCGGGCAGTGCCCGGCGTTTTTCATTCCCGCAATGTTCAGCCGCGCTTACGCTGCGGCAGCACATCCTTGAGCTTGGCATGCATGCTGCGCAGGGTTTTCTCGGTAGCCGACCAGTCg                               <  8:354942/147‑1 (MQ=255)
                                                                                                           cgacAGGCACCAGACACAACAAAAAACGCCGGGCAGTGCCCGGCGTTTTTCATTCCCGCAATGTTCAGCCGCGCTTACGCTGCg                                                                                             >  6:368065/1‑84 (MQ=255)
                                                                                                           cgacAGGCACCAGACACAACAAAAAACGCCGGGCAGTGCCCGGCGTTTTTCATTCCCGCAATGTTCAGCCGCGCTTACGCTGCg                                                                                             <  5:368065/84‑1 (MQ=255)
                                                                                                            gacAGGCACCAGACACAACAAAAAACGCCGGGCAGTGCCCGGCGTTTTTCATTCCCGCAATGTTCAGCCGCGCTTACGCTGCGGCAGCACATCCTTGAGCTTGGCATGCATGCTGCGCAGGGTTTTCTCGGTAGCCGACCAGTCGATg                            >  5:204917/1‑148 (MQ=255)
                                                                                                                ggCACCAGACACAACAAAAAACGCCGGGCAGTGCCCGGCGTTTTTCATTCCCGCAATGTTCAGCCGCGCTTACGCTGCGGCAGCACATCCTTGAGCTTGGCATGCATGCTGCGCAGGGTTTTCTCGGTAGCCGACCAGTCGATGCAGGc                       <  2:374394/149‑1 (MQ=255)
                                                                                                                  cACCAGACACAACAAAAAACGCCGGGCAGTGCCCGGCGTTTTTCATTCCCGCAATGTTCAGCCGCGCTTACGCTGCGGCAGCACATCCTTGAGCTTGGCATGCATGCTGCGCAGGGTTTTCTCGGTAGCCGACCAGTCGATGCAGGCGt                     <  3:130116/149‑1 (MQ=255)
                                                                                                                  cACCAGACACAACAAAAAACGCCGGGCAGTGCCCGGAGTTTTTCATTCCCGCAATGTTCAGCCGCGCTTACGCTGCGGCAGCACTTCGTTGAGATTGGCATTCATGCTGCGCAGTGTTTTCTCGGTAGCCGAACAGTCGATGCAGGCGt                     <  6:98660/149‑1 (MQ=255)
                                                                                                                       gACACAACAAAAAACGCCGGGCAGTGCCCGGCGTTTTTCATTCCCGCAATGTTCAGCCGCGCTTACGCTGCGGCAGCACATCCTTGa                                                                              >  7:73642/1‑87 (MQ=255)
                                                                                                                       gACACAACAAAAAACGCCGGGCAGTGCCCGGCGTTTTTCATTCCCGCAATGTTCAGCCGCGCTTACGCTGCGGCAGCACATCCTTGa                                                                              <  8:73642/87‑1 (MQ=255)
                                                                                                                          acaacaAAAAACGCCGGGCAGTGCCCGGCGTTTTTCATTCCCGCAATGTTCAGCCGCGCTTATGCTGCGGCAGCACATCCTTGAGCTTGGCATGCATGCTGCGCAGGGTTTTCTCGGTAGCCGACCAGTCGATGCAGGCGTCCGTGATc             >  3:246443/1‑149 (MQ=255)
                                                                                                                           caacaaAAAACGCCGGGCAGTGCCCGGCGTTTTTCATTCCCGCAATGTTCAGCCGCGCTTACGCTGCGGCAGCACATCCTTGAGCTTGGCATGCATGCTGCGCAGGGTTTTCTCg                                              <  8:252604/115‑1 (MQ=255)
                                                                                                                           caacaaAAAACGCCGGGCAGTGCCCGGCGTTTTTCATTCCCGCAATGTTCAGCCGCGCTTACGCTGCGGCAGCACATCCTTGAGCTTGGCATGCATGCTGCGCAGGGTTTTCTCg                                              >  7:252604/1‑115 (MQ=255)
                                                                                                                              caaAAAACGCCGGGCAGTGCCCGGCGTTTTTCATTCCCGCAATGTTCAGCCGCGCTTACGCTGCGGCAGCACATCCTTGAGCTTGGCATGCATGCTGCGCAGGGTTTTCTCGGTAGCCGACCAGTCGATGc                           <  5:365762/131‑1 (MQ=255)
                                                                                                                              caaAAAACGCCGGGCAGTGCCCGGCGTTTTTCATTCCCGCAATGTTCAGCCGCGCTTACGCTGCGGCAGCACATCCTTGAGCTTGGCATGCATGCTGCGCAGGGTTTTCTCGGTAGCCGACCAGTCGATGc                           >  6:365762/1‑131 (MQ=255)
                                                                                                                              caaAAAACGCCGGGCAGTGCCCGGCGTTTTTCATTCCCGCAATGTTCAGCCGCGCTTACGCTGCGGCAGCACATCCTTGAGCTTGGCATGCATGCTGCGCAGGGTTTTCTCGGTAGCCGACCAGTCGATGCAg                         >  8:342521/1‑133 (MQ=255)
                                                                                                                              caaAAAACGCCGGGCAGTGCCCGGCGTTTTTCATTCCCGCAATGTTCAGCCGCGCTTACGCTGCGGCAGCACATCCTTGAGCTTGGCATGCATGCTGCGCAGGGTTTTCTCGGTAGCCGACCAGTCGATGCAg                         <  7:342521/133‑1 (MQ=255)
                                                                                                                               aaaaaaCGCCGGGCAGTGCCCGGCGTTTTTCATTCCCGCAATGTTCAGCCGCGCTTACGCTGCGGCAGCACATCCTTGAGCTTGGCATGCATGCTGCGCAGGGTTTTCTCGGTAGCCGACCAGTCGATGCAGGCGTCCGTGATCGacac        >  1:124028/1‑149 (MQ=255)
                                                                                                                                 aaaaCGCCGGGCAGTGCCCGGCGTTTTTCATTCCCGCAATGTTCAGCCGCGCTTACGCTGCGGCAGCACATCCTTGAGCTTGGCATGCATGCTGCGCAGGGTTTTCTCGGTAGCCGACCAGTCGATGCAGGCGTCCGTGATCGac          >  3:317516/1‑145 (MQ=255)
                                                                                                                                 aaaaCGCCGGGCAGTGCCCGGCGTTTTTCATTCCCGCAATGTTCAGCCGCGCTTACGCTGCGGCAGCACATCCTTGAGCTTGGCATGCATGCTGCGCAGGGTTTTCTCGGTAGCCGACCAGTCGATGCAGGCGTCCGTGATCGac          <  4:317516/145‑1 (MQ=255)
                                                                                                                                  aaaCGCCGGGCAGTGCCCGGCGTTTTTCATTCCCGCAATGTTCAGCCGCGCTTACGCTGCGGCAGCACATCCTTGAGCTTGGCATGCATGCTGCGCAGGGTTTTCTCGGTAGCCGACCAGTCGATGCAGGCGTCCGTGATCGACACGc      >  5:310880/1‑148 (MQ=255)
                                                                                                                                  aaaCGCCGGGCACTGCCCGGCGATTTTCATTCCCGCAATGTTCAGCCGCGCTTACGCTGCGGCAGCACATCCTTGAGCTTGGCATGCATGCTGCGCAGGGTTTTCTCGGTAGCCGACCAGTCGATGCAGGCGTCCGTGATCGACACGc      >  4:264447/1‑148 (MQ=255)
                                                                                                                                   aaCGCCGGGCAGTGCCCGGCGTTTTTCATTCCCGCAATGTTCAGCCGCGCTTACGCTGCGGCAGCACATCCTTGAGCTTGGCATGCATGCTGCGCAGGGTTTTCTCGGTAGCCGACCAGTCGATGCAGGCGTCCGTGATCGACACGCCa    >  4:230899/1‑149 (MQ=255)
                                                                                                                                   aaCGCCGGGCACTGCCCGGCGTTTTTAATTCCCGCAATGTTCAGACGCGCTTACGCTCCAGCAGCACATCCTTGAGCTTCGCATGCATGCTGCGCAGGCTTTTCTCGGTAGACGACCAGTCGATGCAGGCGTCCGTGATCGACACGcc     >  6:335099/1‑148 (MQ=255)
                                                                                                                                    aCGCCGGGCACTGCCCGGCGTTTTTCATTCCCGCAATGTTCAGGCGCGCTTACGCTGCGGCAGCACATCCTTGAGCTTGGCATGCATGCTGCGCAGGGTTTTCTCGGTAGCCGACCAGTCGATGCAGGCGTCCGTGATc             <  1:308834/139‑1 (MQ=255)
                                                                                                                                    aCGCCGGGCACTGCCCGGCGTTTTTCATTCCCGCAATGTTCAGCCGCGCTTACGCTGCGGCAGCACATCCTTGAGCTTGGCATGCATGCTGCGCAGGGTTTtct                                                >  2:343972/1‑104 (MQ=255)
                                                                                                                                    aCGCCGGGCACTGCCCGGCGTTTTTCATTCCCGCAATGTTCAGCCGCGCTTACGCTGCGGCAGCACATCCTTGAGCTTGGCATGCATGCTGCGCAGGGTTTtct                                                <  1:343972/104‑1 (MQ=255)
                                                                                                                                    aCGCCGGGCACTGCCCGGCGTTTTTCATTCCCGCAATGTTCAGCCGCGCTTACGCTGCGGCAGCACATCCTTGAGCTTGGCATGCATGCTGCGCAGGGTTTTCTCGGTAGCCGACCAGTCGATGCAGGCGTCCGTGATc             >  2:308834/1‑139 (MQ=255)
                                                                                                                                     cGCCGGGCAGTGCCCGGCGTTTTTCATTCCCGCAATGTTCAGCCGCGCTTACGCTGCGGCAGCACATCCTTGAGCTTGGCATGCATGCTGCGCAGGGTTTTCTCGGTAGCCGACCAGTCGATGCAGGCGTCCGTGATCGACACGCCATa  >  7:39470/1‑149 (MQ=255)
                                                                                                                                              |                                                                                                                                           
GCTCTGGAGTACGCCGAGCAGATGGGGTACACGGTGATCCCGTCTTGCTCATATGTGGAGCGTTACATGGAGCGCCAGCAGCGGCATTCGAGCAAGGCCTGAGCTGACGACAGGCACCAGACACAACAAAAAACGCCGGGCAGTGCCCGGCGTTTTTCATTCCCGCAATGTTCAGCCGCGCTTACGCTGCGGCAGCACATCCTTGAGCTTGGCATGCATGCTGCGCAGGGTTTTCTCGGTAGCCGACCAGTCGATGCAGGCGTCCGTGATCGACACGCCATA  >  NC_002947/2651134‑2651415

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: