Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 892,015 C→G 5.4% R567P (CGG→CCG) ‡ pta ← phosphate acetyltransferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_002947892,0150CG 5.4% 241.3 / 2.9 74R567P (CGG→CCG) ‡ptaphosphate acetyltransferase
Reads supporting (aligned to +/- strand):  ref base C (39/31);  new base G (2/2);  total (41/33)
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.17e-01

GCCGGGTTGGCCGCTGCGTCATACTGCAACGGCCCATCGATCAGTAGCTCCTGCGCCGCGTCCTGTGCCAGGCGTGTCGCTTCGCGCACTTTCTCGACCTCTTCATCGCTGGCCGAATCGCTTGAATAGCTGATCATCGCCACCCGTGGCGCGATGCCGAAGGCCTGCGCCGAGGCGGCACTCTGTTGGGCAATCTCCGCCAGCTCGGCAGCGCTTGGGTGCGGGTTCATCACGCAGTCGCCGTACACCAGCACCTGCTCGGGGAACAGCATGAAGAACACCGACGAGACCAG  >  NC_002947/891871‑892163
                                                                                                                                                |                                                                                                                                                    
gCCGGGTTGGCCGCTGCGTCATACTGCAACGGCCCATCGATCAGTAGCTCCTGCGCCGCGTCCTGTGCCAGGCGTGTCGCTTCGCGCACTTTCTCGACCTCTTCATCGCTGGCCGAATCGCTTGAATAGCTGATCATCGCCACCCGTgg                                                                                                                                                  <  3:99095/149‑1 (MQ=255)
 ccGGGTTGGCCGCTGCGTCATACTGCAACGGCCCATCGATCAGTAGCTCCTGCGCCGCGTCCTGTGCCAGGCGTGTCGCTTCGCGCACTTTCTCGACCTCTTCATCGCTGGCCGAATCGCTTGAATAGCTGATCATCGCCACGGGTGgc                                                                                                                                                 >  6:330675/1‑149 (MQ=255)
   gggTTGGCCGCTGCGTCATACTGCAACGGCCCATCGATCAGTAGCTCCTGCGCCGCGTCCTGTGCCAGGCGTGTCGCTTCGCGCACTTTCTCGACCTCTTCATCGCTGGCCGAATCGCTTGAATAGCTGATCATCGCCACCCGTGgcgc                                                                                                                                               <  1:227652/149‑1 (MQ=255)
    ggTTGGCCGCTGCGTCATACTGCAACGGCCCATCGATCAGTAGCTCCTGCGCCGCGTCCTGTGCCAGGCGTGTCGCTTCGCGCACTTTCTCGACCTCTTCATCGCTGGCCGAATCGCTTGAATAGCTGATCATCGCCACGGGTGgcgat                                                                                                                                              <  5:330675/149‑3 (MQ=255)
         gCCGCTGCGTCATACTGCAACGGCCCATCGATCAGTAGCTCCTGCGCCGCGTCCTGTGCCAGGCGTGTCGCTTCGCGCACTTTCTCGACCTCTTCATCGCTGGCCGAATCGCTTGAATAGCTGATCATCGCCACCCGTGGCGCGATGcc                                                                                                                                         >  6:112656/1‑149 (MQ=255)
           cGCTGCGTCATACTGCAACGGCCCATCGATCAGTAGCTCCTGCGCCGCGTCCTGTGCCAGGCGTGTCGCTTCGCGCACTTTCTCGACCTCTTCATCGCTGGCCGAATCGCTTGAATAGCTGATCATCGCCACCCGTGGCGCGATGCCg                                                                                                                                        >  5:117955/1‑148 (MQ=255)
                 gTCATACTGCAACGGCCCATCGATCAGTAGCTCCTGCGCCGCGTCCTGTGCCAGGCGTGTCGCTTCGCGCACTTTCTCGACCTCTTCATCGCTGGCCGAATCGCTTGAATAGCTGATCATCGCCACCCGTGGCGCGATGCCGAAGGCCt                                                                                                                                 >  3:328672/1‑149 (MQ=255)
                    aTACTGCAACGGCCCATCGATCAGTAGCTCCTGCGCCGCGTCCTGTGCCAGGCGTGTCGCTTCGCGCACTTTCTCGACCTCTTCATCGCTGGCCGAATCGCTTGAATAGCTGATCATCGCCACCCGTGGCGCGATGCCGAAgg                                                                                                                                    <  3:270602/143‑1 (MQ=255)
                    aTACTGCAACGGCCCATCGATCAGTAGCTCCTGCGCCGCGTCCTGTGCCAGGCGTGTCGCTTCGCGCACTTTCTCGACCTCTTCATCGCTGGCCGAATCGCTTGAATAGCTGATCATCGCCACCCGTGGCGCGATGCCGAAgg                                                                                                                                    >  4:270602/1‑143 (MQ=255)
                        tGCAACGGCCCATCGATCAGTAGCTCCTGCGCCGCGTCCTGTGCCAGGCGTGTCGCTTCGCGCACTTTCTCGACCTCTTCATCGCTGGCCGAATCGCTTGAATAGCTGATCATCGCCACCCGTGGCGCGATGCCGAAGGCCTGCGcc                                                                                                                            >  7:79162/1‑147 (MQ=255)
                        tGCAACGGCCCATCGATCAGTAGCTCCTGCGCCGCGTCCTGTGCCAGGCGTGTCGCTTCGCGCACTTTCTCGACCTCTTCATCGCTGGCCGAATCGCTTGAATAGCTGATCATCGCCACCCGTGGCGCGATGCCGAAGGCCTGCGcc                                                                                                                            >  1:3061/1‑147 (MQ=255)
                          cAACGGCCCATCGATCAGTAGCTCCTGCGCCGCGTCCTGTGCCAGGCGTGTCGCTTCGCGCACTTTCTCGACCTCTTCATCGCTGGCCGAATCGCTTGAATAGCTGATCATCGCCACCCGTGGCGCGATGCCGAAGGCCTGCGCCGAg                                                                                                                         <  4:138104/148‑1 (MQ=255)
                             cGGCCCATCGATCAGTAGCTCCTGCGCCGCGTCCTGTGCCAGGCGTGTCGCTTCGCGCACTTTCTCGACCTCTTCATCGCTGGCCGAATCGCTTGAATAGCTGATCATCGCCACCCGTGGCGCGATGCCGAAGGCCTGCGCCGAggcgg                                                                                                                     <  5:112656/149‑1 (MQ=255)
                                 ccATCGATCAGTAGCTCCTGCGCCGCGTCCTGTGCCAGGCGTGTCGCTTCGCGCACTTTCTCGACCTCTTCATCGCTGGCCGAATCGCTTGAATAGCTGATCATCGCCACCCGTGGCGCGATGCCGAAGGCCTGCGCCGAGGCGGCAct                                                                                                                 <  7:111676/149‑1 (MQ=255)
                                 ccATCGATCAGTAGCTCCTGCGCCGCGTCCTGTGCCAGGCGTGTCGCTTCGCGCACTTTCTCGACCTCTTCATCGCTGGCCGAATCGCTTGAATAGCTGATCATCGCCACCCGTGGCGCGATGCCGAAGGCCTGCGCCGAGGCGGCAct                                                                                                                 >  5:252983/1‑149 (MQ=255)
                                 ccATCGATCAGTAGCTCCTGCGCCGCGTCCTGTGCCAGGCGTGTCGCTTCGCGCACTTTCTCGACCTCTTCATCGCTGGCCGAATCGCTTGAATAGCTGATCATCGCCACCCGTGGCGCGATGCCGAAGGCCTGCGCCGAGGCGGCAct                                                                                                                 >  1:253274/1‑149 (MQ=255)
                                 ccATCGATCAGTAGCTCCTGCGCCGCGTCCTGTGCCAGGCGTGTCGCTTCGCGCACTTTCTCGACCTCTTCATCGCTGGCCGAATCGCTTGAATAGCTGATCATCGCCACCCGTGGCGCGATGCCGAAGGCCTGCGCCGAGGCGGCAc                                                                                                                  >  3:142250/1‑148 (MQ=255)
                                    tCGATCAGTAGCTCCTGCGCCGCGTCCTGTGCCAGGCGTGTCGCTTCGCGCACTTTCTCGACCTCTTCATCGCTGGCCGAATCGCTTGAATAGCTGATCATCGCCACCCGTGGCGCGATGCCGAAGGCCTGCGCCGAGGCGGCACTCTg                                                                                                              >  6:118199/1‑149 (MQ=255)
                                      gATCAGTAGCTCCTGCGCCGCGTCCTGTGCCAGGCGTGTCGCTTCGCGCACTTTCTCGACCTCTTCATCGCTGGCCGAATCGCTTGAATAGCTGATCATCGCCACCCGTGGCGCGATGCCGAAGGCCTGCGCCGAGGCGGCACTCTGtt                                                                                                            <  4:328672/149‑1 (MQ=255)
                                       aTCAGTAGCTCCTGCGCCGCGTCCTGTGCCAGGCGTGTCGCTTCGCGCACTTTCTCGACCTCTTCATCGCTGGCCGAATCGCTTGAATAGCTGATCATCGCCACGGGTGgcgatg                                                                                                                                             >  7:109200/1‑112 (MQ=255)
                                       aTCAGTAGCTCCTGCGCCGCGTCCTGTGCCAGGCGTGTCGCTTCGCGCACTTTCTCGACCTCTTCATCGCTGGCCGAATCGCTTGAATAGCTGATCATCGCCACGGGTGgcgatg                                                                                                                                             <  8:109200/115‑4 (MQ=255)
                                       aTCAGTAGCTCCTGCGCCGCGTCCTGTGCCAGGCGTGTCGCTTCGCGCACTTTCTCGACCTCTTCATCGCTGGCCGAATCGCTTGAATAGCTGATCATCGCCACCCGTGGCGCGATGCCGAAGGCCTGCGCCGAGGCGGCACTCTGTTg                                                                                                           >  3:157773/1‑149 (MQ=255)
                                        tCAGTAGCTCCTGCGCCGCGTCCTGTGCCAGGCGTGTCGCTTCGCGCACTTTCTCGACCTCTTCATCGCTGGCCGAATCGCTTGAATAGCTGATCATCGCCACCCGTGGTGCGATGCCGAAGGCCTGCGCCGAGGCGTCACTCTGTTgg                                                                                                          <  3:240259/149‑1 (MQ=255)
                                        tCAGTAGCTCCTGCGCCGCGTCCTGTGCCAGGCGTGTCGCTTCGCGCACTTTCTCGACCTCTTCATCGCTGGCCGAATCGCTTGAATAGCTGATCATCGCCACCCGTGGCGCGATGCCGAAGGc                                                                                                                                   <  5:283105/124‑1 (MQ=255)
                                        tCAGTAGCTCCTGCGCCGCGTCCTGTGCCAGGCGTGTCGCTTCGCGCACTTTCTCGACCTCTTCATCGCTGGCCGAATCGCTTGAATAGCTGATCATCGCCACCCGTGGCGCGATGCCGAAGGc                                                                                                                                   >  6:283105/1‑124 (MQ=255)
                                        tCAGTAGCTCCTGCGCCGCGTCCTGTGCCAGGCGTGTCGCTTCGCGCACTTTCTCGACCTCTTCATCGCTGGCCGAATCGCTTGAATAGCTGATCATCGCCACCCGTGGCGCGATGCCGAAGGCCTGCGCCGAGGCGGCACTCTGTTgg                                                                                                          <  6:117955/149‑1 (MQ=255)
                                          aGTAGCTCCTGCGCCGCGTCCTGTGCCAGGCGTGTCGCTTCGCGCACTTTCTCGACCTCTTCATCGCTGGCCGAATCGCTTGAATAGCTGATCATCGCCACCCGTGGCGCGATGCCGAAGGCCTGCGCCGAGGCGGCACTCTGTTGGGc                                                                                                        <  4:35210/149‑1 (MQ=255)
                                           gTAGCTCCTGCGCCGCGTCCTGTGCCAGGCGTGTCGCTTCGCGCACTTTCTCGACCTCTTCATCGCTGGCCGAATCGCTTGAATAGCTGATCATCGCCACCCGTGGCGCGATGCCGAAGGCCTGCGCCGAGGCGGCACTCTGTTGGGCa                                                                                                       <  5:320745/149‑1 (MQ=255)
                                              gCTCCTGCGCCGCGTCCTGTGCCAGGCGTGTCGCTTCGCGCACTTTCTCGACCTCTTCATCGCTGGCCGAATCGCTTGAATAGCTGATCATCGCCACCCGTGGCGCGATGCCGAAGGCCTGCGCCGAGGCGGCACTCTGTTGGGCAAtc                                                                                                    >  8:205013/1‑149 (MQ=255)
                                                  cTGCGCCGCGTCCTGTGCCAGGCGTGTCGCTTCGCGCACTTTCTCGACCTCTTCATCGCTGGCCGAATCGCTTGAATAGCTGATCATCGCCACCCGTGGCGCGATGCCGAAGGCCTGCGCCGAGGCGGCACTCTGTTGGGCAATCTCCg                                                                                                >  4:325770/1‑149 (MQ=255)
                                                   tGCGCCGCGTCCTGTGCCAGGCGTGTCGCTTCGCGCACTTTCTCGACCTCTTCATCGCTGGCCGAATCGCTTGAATAGCTGATCATCGCCACCCGTGGCGCGATGCCGAAGGCCTGCGCCGAGGCGGCACTCTGTTGGGCAATCTCCGc                                                                                               >  7:250969/1‑149 (MQ=255)
                                                          cgTCCTGTGCCAGGCGTGTCGCTTCGCGCACTTTCTCGACCTCTTCATCGCTGGCCGAATCGCTTGAATAGCTGATCATCGCCACCCGTGGCGCGATGCCGAAGGCCTGCGCCGAGGCGGCACTCTGTTGGGCAATCTCCGCCAGCTCg                                                                                        <  4:42237/149‑1 (MQ=255)
                                                                   ccAGGCGTGTCGCTTCGCGCACTTTCTCGACCTCTTCATCGCTGGCCGAATCGCTTGAATAGCTGATCATCGCCACCCGTGGCGCGATGCCGAAGGCCTGCGCCGAGGCGGCAct                                                                                                                 >  5:238348/1‑115 (MQ=255)
                                                                   ccAGGCGTGTCGCTTCGCGCACTTTCTCGACCTCTTCATCGCTGGCCGAATCGCTTGAATAGCTGATCATCGCCACCCGTGGCGCGATGCCGAAGGCCTGCGCCGAGGCGGCAct                                                                                                                 <  6:238348/115‑1 (MQ=255)
                                                                   ccAGGCGTGTCGCTTCGCGCACTTTCTCGACCTCTTCATCGCTGGCCGAATCGCTTGAATAGCTGATCATCGCCACCCGTGGCGCGATGCCGAAGGCCTGCGCCGAGGCGGCACTCTGTTGGGCAATCTCCGCCAGCTCGGCAGCGCtt                                                                               >  1:116626/1‑149 (MQ=255)
                                                                     aGGCGTGTCGCTTCGCGCACTTTCTCGACCTCTTCATCGCTGGCCGAATCGCTTGAATAGCTGATCATCGCCACCCGTGGCGCGATGCCGAAGGCCTGCGCCGa                                                                                                                          >  7:137160/1‑104 (MQ=255)
                                                                     aGGCGTGTCGCTTCGCGCACTTTCTCGACCTCTTCATCGCTGGCCGAATCGCTTGAATAGCTGATCATCGCCACCCGTGGCGCGATGCCGAAGGCCTGCGCCGa                                                                                                                          <  8:137160/104‑1 (MQ=255)
                                                                         gtgtCGCTTCGCGCACTTTCTCGACCTCTTCATCGCTGGCCGAATCGCTTGAATAGCTGATCATCGCCACCCGTGGCGCGATGCCGAAGGCCTGCGCCGAGGCGGCACTCTGTTGGGCAATCTCCGCCAGCTCGGCAg                                                                                    <  1:95139/138‑1 (MQ=255)
                                                                         gtgtCGCTTCGCGCACTTTCTCGACCTCTTCATCGCTGGCCGAATCGCTTGAATAGCTGATCATCGCCACCCGTGGCGCGATGCCGAAGGCCTGCGCCGAGGCGGCACTCTGTTGGGCAATCTCCGCCAGCTCGGCAg                                                                                    >  2:95139/1‑138 (MQ=255)
                                                                         gtgtCGCTTCGCGCACTTTCTCGACCTCTTCATCGCTGGCCGAATCGCTTGAATAGCTGATCATCGCCACCCGTGGCGCGATGCCGAAGGCCTGCGCCGAGGCGGCACTCTGTTGGGCAATCTCCGCCAGCTCGGCAGCGCTTGGGTGc                                                                         <  5:271002/149‑1 (MQ=255)
                                                                           gtCGCTTCGCGCACTTTCTCGACCTCTTCATCGCTGGCCGAATCGCTTGAATAGCTGATCATCGCCACCCGTGGCGCGATGCCGAAGGCCTGCGCCGAGGCGGCACTCTGTTGGGCAATCTCCGCCAGCTCGGCAGCGCTTGGGTGCgg                                                                       <  2:3061/149‑1 (MQ=255)
                                                                            tCGCTTCGCGCACTTTCTCGACCTCTTCATCGCTGGCCGAATCGCTTGAATAGCTGATCATCGCCACCCGTGGCGCGATGCCGAAGGCCTGCGCCGAGGCGGCACTCTGTTGGGCAATCTCCGCCAGCTCGGCAGCGCTTGGGTGCgg                                                                       <  7:205013/148‑1 (MQ=255)
                                                                              gCTTCGCGCACTTTCTCGACCTCTTCATCGCTGGCCGAATCGCTTGAATAGCTGATCATCGCCACCCGTGGCGCGATGCCGAAGGCCTGCGCCGAGGCGGCACTCTGTTGGGCAATCTCCGCCAGCTCGGCAGCGCTTGGGTGCGGGt                                                                     >  5:37228/1‑148 (MQ=255)
                                                                               cTTCGCGCACTTTCTCGACCTCTTCATCGCTGGCCGAATCGCTTGAATAGCTGATCATCGCCACCCGTGGCGCGATGCCGAAGGCCTGCGCCGAGGCGGCACTCTGTTGGGCAATCTCCGCCAGCTCGGCAGCGCTTGGGTGCGGGtt                                                                    >  6:274819/1‑148 (MQ=255)
                                                                                           tctcGACCTCTTCATCGCTGGCCGAATCGCTTGAATAGCTGATCATCGCCACCCGTGGCGCGATGCCGAAGGCCTGCGCCGAGGCGGCACTCTGTTGGGCAATCTCCGCCAGCTCGGCAGCGCTTGGGTGCGGGTTCATCACGCAGTCg                                                       >  3:59089/1‑149 (MQ=255)
                                                                                               gACCTCTTCATCGCTGGCCGAATCGCTTGAATAGCTGATCATCGCCACCCGTGGCGCGATGCCGAAGGCCTGCGCCGAGGCGGCACTCTGTTGGGCAATCTCCGCCAGCTCGGCAGCGCTTGGGTGCGGGTTCATCACGCAGTCGcc                                                     <  2:35068/147‑1 (MQ=255)
                                                                                               gACCTCTTCATCGCTGGCCGAATCGCTTGAATAGCTGATCATCGCCACCCGTGGCGCGATGCCGAAGGCCTGCGCCGAGGCGGCACTCTGTTGGGCAATCTCCGCCAGCTCGGCAGCGCTTGGGTGCGGGTTCATCACGCAGTCGCCg                                                    >  2:23372/1‑148 (MQ=255)
                                                                                                     ttCATCGCTGGCCGAATCGCTTGAATAGCTGATCATCGCCACCCGTGGCGCGATGCCGAAGGCCTGCGCCGAGGCGGCACTCTGTTGGGCAATCTCCGCCAGCTCGGCAGCGCTTGGGTGCGGGTTCATCACGCAGTc                                                        >  2:90222/1‑138 (MQ=255)
                                                                                                     ttCATCGCTGGCCGAATCGCTTGAATAGCTGATCATCGCCACCCGTGGCGCGATGCCGAAGGCCTGCGCCGAGGCGGCACTCTGTTGGGCAATCTCCGCCAGCTCGGCAGCGCTTGGGTGCGGGTTCATCACGCAGTc                                                        <  1:90222/138‑1 (MQ=255)
                                                                                                        aTCGCTGGCCGAATCGCTTGAATAGCTGATCATCGCCACCCGTGGCGCGATGCCGAAGGCCTGCGCCGAGGCGGCACTCTGTTGGGCAATCTCCGCCAGCTCGGCAGCGCTTGGGTGCgg                                                                       >  2:266364/1‑120 (MQ=255)
                                                                                                        aTCGCTGGCCGAATCGCTTGAATAGCTGATCATCGCCACCCGTGGCGCGATGCCGAAGGCCTGCGCCGAGGCGGCACTCTGTTGGGCAATCTCCGCCAGCTCGGCAGCGCTTGGGTGCGGGTTCATCACGCAGTCGCCGTACACCAGc                                           <  5:274819/148‑1 (MQ=255)
                                                                                                        aTCGCTGGCCGAATCGCTTGAATAGCTGATCATCGCCACCCGTGGCGCGATGCCGAAGGCCTGCGCCGAGGCGGCACTCTGTTGGGCAATCTCCGCCAGCTCGGCAGCGCTTGGGTGCGGGTTCATCACGCAGTCGCCGTACACCAGCa                                          <  8:250969/149‑1 (MQ=255)
                                                                                                         tCGCTGGCCGAATCGCTTGAATAGCTGATCATCGCCACCCGTGGCGCGATGCCGAAGGCCTGCGCCGAGGCGGCACTCTGTTGGGCAATCTCCGCCAGCTCGGCAGCGCTTGGGTGCGGGTTCATCACGCAGTCGCCGTACACCAGCAc                                         <  8:79162/149‑1 (MQ=255)
                                                                                                           gCTGGCCGAATCGCTTGAATAGCTGATCATCGCCACCCGTGGCGCGATGCCGAAGGCCTGCGCCGAGGCGGCACTCTGTTGGGCAATCTCCGCCAGCTCGGCAGCGCTTGGGTGCGGGTTCATCACGCAGTCGCCGTACACCAGCACCt                                       >  3:172283/1‑149 (MQ=255)
                                                                                                             tGGCCGAATCGCTTGAATAGCTGATCATCGCCACCCGTGGCGCGATGCCGAAGGCCTGCGCCGAGGCGGCACTCTGTTGGGCAATCTCCGCCAGCTCGGCAGCGCTTGGGTGCggg                                                                      <  8:240033/116‑1 (MQ=255)
                                                                                                             tGGCCGAATCGCTTGAATAGCTGATCATCGCCACCCGTGGCGCGATGCCGAAGGCCTGCGCCGAGGCGGCACTCTGTTGGGCAATCTCCGCCAGCTCGGCAGCGCTTGGGTGCggg                                                                      >  7:240033/1‑116 (MQ=255)
                                                                                                               gCCGAATCGCTTGAATAGCTGATCATCGCCACCCGTGGCGCGATGCCGAAGGCCTGCGCCGAGGCGGCACTCTGTTGGGCAATCTCCGCCAGCTCGGCAGCGCTTGGGTGCGGGTTCATCACGCAGTCGCCGTACACCAGCACCTGCTc                                   >  6:174374/1‑149 (MQ=255)
                                                                                                               gCCGAATCGCTTGAATAGCTGATCATCGCCACCCGTGGCGCGATGCCGAAGGCCTGCGCCGAGGCGGCACTCTGTTGGGCAATCTCCGCCAGCTCGGCAGCGCTTGGGTGCGGGTTCATCACGCAGTCGCCGTACACCAGCACCTGCTc                                   >  2:39639/1‑149 (MQ=255)
                                                                                                                  gAATCGCTTGAATAGCTGATCATCGCCACCCGTGGCGCGATGCCGAAGGCCTGCGCCGAGGCGGCACTCTGTTGGGCAATCTCCGCCAGCTCGGCAGCGCTTGGGTGCGGGTTCATCACGCAGTCGCCGTACACCAGCACCTGCTCggg                                <  4:157773/149‑1 (MQ=255)
                                                                                                                        cTTGAATAGCTGATCATCGCCACCCGTGGCGCGATGCCGAAGGCCTGCGCCGAGGCGGCACTCTGTTGGGCAATCTCCGCCAGCTCGGCAGCGCTTGGGTGCGGGTTCATCACGCAGTCGCCGTACACCAGCACCTGCTCGGGGAACAg                          >  5:68681/1‑149 (MQ=255)
                                                                                                                           gAATAGCTGATCATCGCCACCCGTGGCGCGATGCCGAAGGCCTGCGCCGAGGCGGCACTCTGTTGGGCAATCTCCGCCAGCTCGGCAGCGCTTGGGTGCGGGTTCATCACGCAGTCGCCGTACACCAGCACCTGCTCGGGGAACAGc                         <  1:23372/147‑1 (MQ=255)
                                                                                                                           gAATAGCTGATCATCGCCACCCGTGGCGCGATGCCGAAGGCCTGCGCCGAGGCGGCACTCTGTTGGGCAATCTCCGCCAGCTCGGCAGCGCTTGGGTGCGGGTTCATCACGCAGTCGCCGTACACCAGCACCTGCTCGGGGAACAGCAt                       >  4:219866/1‑149 (MQ=255)
                                                                                                                             aTAGCTGATCATCGCCACCCGTGGCGCGATGCCGAAGGCCTGCGCCGAGGCGGCACTCTGTTGGGCAATCTCCGCCAGCTCGGCAGCGCTTGGGTGCGGGTTCATCACGCAGTCGCCGTACACCAGCACCTGCTCGGGGAACAGCATg                      >  6:247042/1‑148 (MQ=255)
                                                                                                                              tAGCTGATCATCGCCACCCGTGGCGCGATGCCGAAGGCCTGCGCCGAGGCGGCACTCTGTTGGGCAATCTCCGCCAGCTCGGCAGCGCTTGGGTGCGGGTTCATCACGCAGTCGCCGTACACCAGCACCTGCTCGGGGAACAGCATgaa                    <  5:247042/149‑1 (MQ=255)
                                                                                                                                gCTGATCATCGCCACCCGTGGCGCGATGCCGAAGGCCTGCGCCGAGGCGGCACTCTGTTGGGCAATCTCCGCCAGCTCGGCAGCGCTTGGGTGCGGGTTCATCACGCAGTCGCCGTACACCAGCACCTGCTCGGGGAACAGCATgaag                   <  4:142250/148‑1 (MQ=255)
                                                                                                                                 cTGATCATCGCCACCCGTGGCGCGATGCCGAAGGCCTGCGCCGAGGCGGCACTCTGTTGGGCAATCTCCGCCAGCTCGGCAGCGCTTGGGTGCGGGTTCATCACGCAGTCGCCGTACACCAGCACCTGCTCGGGGAACAGCATgaagaa                 >  6:221959/1‑149 (MQ=255)
                                                                                                                                   gATCATCGCCACCCGTGGCGCGATGCCGAAGGCCTGCGCCGAGGCGGCACTCTGTTGGGCAATCTCCGCCAGCTCGGCAGCGCTTGGGTGCGGGTTCATCACGCAGTCGCCGTACACCAGCACCTGCTCGGGGAACAGCATGAAGAaca               >  3:39378/1‑149 (MQ=255)
                                                                                                                                    atcatcGCCACCCGTGGCGCGATGCCGAAGGCCTGCGCCGAGGCGGCACTCTGTTGGGCAATCTCCGCCAGCTCGGCAGCGCTTGGGTGCGGGTTCATCACGCAGTCGCCGTACACCAGCACCTGCTCGGGGAACAGCATGAAGAacac              >  5:268864/1‑149 (MQ=255)
                                                                                                                                    atcatcGCCACCCGTGGCGCGATGCCGAAGGCCTGCGCCGAGGCGGCACTCTGTTGGGCAATCTCCGCCAGCTCGGCAGCGCTTGGGTGCGGGTTCATCACGCAGTCGCCGTACACCAGCACCTGCTCGGGGAACAGCATGAAGAacac              <  2:116626/149‑1 (MQ=255)
                                                                                                                                     tcatcGCCACCCGTGGCGCGATGCCGAAGGCCTGCGCCGAGGCGGCACTCTGTTGGGCAATCTCCGCCAGCTCGGCAGCGCTTGGGTGCGGGTTCATCACGCAGTCGCCGTACACCAGCACCTGCTCGGGGAACAGCATGAAGAACAcc             <  5:174374/149‑1 (MQ=255)
                                                                                                                                       atcGCCACCCGTGGCGCGATGCCGAAGGCCTGCGCCGAGGCGGCACTCTGTTGGGCAATCTCCGCCAGCTCGGCAGCGCTTGGGTGCGGGTTCATCACGCAGTCGCCGTACACCAGCACCTGCTCGGGGAACAGCATGAAGAACACcga           >  7:338986/1‑149 (MQ=255)
                                                                                                                                         cGCCACCCGTGGCGCGATGCCGAAGGCCTGCGCCGAGGCGGCACTCTGTTGGGCAATCTCCGCCAGCTCGGCAGCGCTTGGGTGCGGGTTCATCACGCAGTCGCCGTACACCAGCACCTGCTCGGGGAACAGCATGAAGAACACcgac          >  7:243905/1‑148 (MQ=255)
                                                                                                                                           ccACCCGTGGCGCGATGCCGAAGGCCTGCGCCGAGGCGGCACTCTGTTGGGCAATCTCCGCCAGCTCGGCAGCGCTTGGGTGCGGGTTCATCACGCAGTCGCCGTACACCAGCACCTGCTCGGGGAACAGCATGAAGAACACcgac          >  2:39119/1‑146 (MQ=255)
                                                                                                                                             aCCCGTGGCGCGATGCCGAAGGCCTGCGCCGAGGCGGCACTCTGTTGGGCAATCTCCGCCAGCTCGGCAGCGCTTGGGTGCGGGTTCATCACGCAGTCGCCGTACACCAGCACCTGCTCGGGGAACAGCATGAAGAACACCGACGAGAc     >  3:163947/1‑149 (MQ=255)
                                                                                                                                              cccGTGGCGCGATGCCGAAGGCCTGCGCCGAGGCGGCACTCTGTTGGGCAATCTCCGCCAGCTCGGCAGCGCTTGGGTGCGGGTTCATCACGCAGTCGCCGTACACCAGCACCTGCTCGGGGAACAGCATGAAGAACACCGACGAGAcc    <  1:39119/149‑1 (MQ=255)
                                                                                                                                                cGTGGCGCGATGCCGAAGGCCTGCGCCGAGGCGGCACTCTGTTGGGCAATCTCCGCCAGCTCGGCAGCGCTTGGGTGCGGGTTCATCACGCAGTCGCCGTACACCAGCACCTGCTCGGGGAACAGCATGAAGAACACCGACGAGACCAg  <  1:39639/149‑1 (MQ=255)
                                                                                                                                                |                                                                                                                                                    
GCCGGGTTGGCCGCTGCGTCATACTGCAACGGCCCATCGATCAGTAGCTCCTGCGCCGCGTCCTGTGCCAGGCGTGTCGCTTCGCGCACTTTCTCGACCTCTTCATCGCTGGCCGAATCGCTTGAATAGCTGATCATCGCCACCCGTGGCGCGATGCCGAAGGCCTGCGCCGAGGCGGCACTCTGTTGGGCAATCTCCGCCAGCTCGGCAGCGCTTGGGTGCGGGTTCATCACGCAGTCGCCGTACACCAGCACCTGCTCGGGGAACAGCATGAAGAACACCGACGAGACCAG  >  NC_002947/891871‑892163

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