Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 2,526,412 A→T 7.3% K219M (AAG→ATG) ‡ PP_2215 → acetyl‑CoA acetyltransferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029472,526,4120AT 7.3% 196.5 / 5.4 68K219M (AAG→ATG) ‡PP_2215acetyl‑CoA acetyltransferase
Reads supporting (aligned to +/- strand):  ref base A (31/32);  new base T (2/3);  total (33/35)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.93e-01

CAATGCCTTCAGCCGCGAGGCCCAGGACCAGTTCGCCATCGCCTCGCTGACCCGAGCGCAGGAAGCCATCAGCAGCGGCCGTTTTGCCGCCGAGATCGTGCCGGTGGAAGTCACCGAGGGCAAGGAAAAGCGCGTCATCAAGGATGACGAGCAGCCGCCCAAGGCGCGTCTGGACAAGATTGCGCAGCTCAAACCGGCGTTTCGTGAAGGCGGCACCGTGACGGCGGCCAACGCCAGTTCGATTTCCGACGGCGCTGCGGCGCTGGTACTGATGCGCCGCTCCGAGGC  >  NC_002947/2526272‑2526559
                                                                                                                                            |                                                                                                                                                    
cAATGCCTTCAGCCGCGAGGCCCAGGACCAGTTCGCCATCGCCTCGCTGACCCGAGCGCAGGAAGCCATCAGCAGCGGCCGTTTTGCCGCCGAGATCGTGCCGGTGGAAGTCACCGAGGGCAAGGAAAAGCGCGTCATCCTTGATGAc                                                                                                                                               <  8:200863/148‑1 (MQ=255)
cAATGCCTTCAGCCGCGAGGCCCAGGACCAGTTCGCCATCGCCTCGCTGACCCGAGCGCAGGAAGCCATCAGCAGCGGCCGTTTTGCCGCCGAGATCGTGCCGGTGGAAGTCACCGAGGGCAAGGAAAAGCGCGTCATCCTTGATGAc                                                                                                                                               >  7:200863/1‑148 (MQ=255)
 aaTGCCTTCAGCCGCGAGGCCCAGGACCAGTTCGCCATCGCCTCGCTGACCCGAGCGCAGGAAGCCATCAGCAGCGGCCGTTTTGCCGCCGAGATCGTGCCGGTGGAAGTCACCGAGGGCAAGGAAAAGCGCGTCATCCTTGATGACgc                                                                                                                                             <  7:299552/149‑2 (MQ=255)
       ttCAGCCGCGAGGCCCAGGACCAGTTCGCCATCGCCTCGCTGACCCGAGCGCAGGAAGCCATCAGCAGCGGCCGTTTTGCCGCCGAGATCGTGCCGGTGGAAGTCACCGAGGGCAAGGAAAAGCGCGTCATCAAGGATGACGAGCAgcc                                                                                                                                       >  5:177371/1‑149 (MQ=255)
       ttCAGCCGCGAGGCCCAGGACCAGTTCGCCATCGCCTCGCTGACCCGAGCGCAGGAAGCCATCAGCAGCGGCCGTTATGCCGCCGAGATCGTGCCGGTGGAAGTCACCGAGGGCAAGGAAAAGCGCGTCATCAAGGATGACGCGCAgcc                                                                                                                                       >  6:43864/1‑149 (MQ=255)
         cAGCCGCGAGGCCCAGGACCAGTTCGCCATCGCCTCGCTGACCCGAGCGCAGGAAGCCATCAGCAGCGGCCGTTTTGCCGCCGAGATCGTGCCGGTGGAAGTCACCGAGGGCAAGGAAAAGCGCGTCATCAAg                                                                                                                                                     >  4:268148/1‑133 (MQ=255)
         cAGCCGCGAGGCCCAGGACCAGTTCGCCATCGCCTCGCTGACCCGAGCGCAGGAAGCCATCAGCAGCGGCCGTTTTGCCGCCGAGATCGTGCCGGTGGAAGTCACCGAGGGCAAGGAAAAGCGCGTCATCAAg                                                                                                                                                     <  3:268148/133‑1 (MQ=255)
                gAGGCCCAGGACCAGTTCGCCATCGCCTCGCTGACCCGAGCGCAGGAAGCCATCAGCAGCGGCCGTTTTGCCGCCGAGATCGTGCCGGTGGAAGTCACCGAGGGCAAGGAAAAGCGCGTCATCAAGGATGAc                                                                                                                                               >  5:282761/1‑132 (MQ=255)
                gAGGCCCAGGACCAGTTCGCCATCGCCTCGCTGACCCGAGCGCAGGAAGCCATCAGCAGCGGCCGTTTTGCCGCCGAGATCGTGCCGGTGGAAGTCACCGAGGGCAAGGAAAAGCGCGTCATCAAGGATGAc                                                                                                                                               <  6:282761/132‑1 (MQ=255)
                   gCCCAGGACCAGTTCGCCATCGCCTCGCTGACCCGAGCGCAGGAAGCCATCAGCAGCGGCCGTTTTGCCGCCGAGATCGTGCCGGTGGAAGTCACCGAGGGCAAGGAAAAGCGCGTCATCAAGGATGACGAGCAGCCGCCCAAGgcgc                                                                                                                            <  4:35471/148‑1 (MQ=255)
                   gCCCAGGACCAGTTCGCCATCGCCTCGCTGACCCGAGCGCAGGAAGCCATCAGCAGCGGCCGTTTTGCCGCCGAGATCGTGCCGGTGGAAGTCACCGAGGGCAAGGAAAAGCGCGTCATCAAGGATGACGAGCAGCCGCCCAAGGcgcg                                                                                                                           >  5:205245/1‑149 (MQ=255)
                     ccAGGACCAGTTCGCCATCGCCTCGCTGACCCGAGCGCAGGAAGCCATCAGCAGCGGCCGTTTTGCCGCCGAGATCGTGCCGGTGGAAGTCACCGAGGGCAAGGAAAAGCGCGTCATCAAGGATGACGagca                                                                                                                                          >  7:184521/1‑132 (MQ=255)
                     ccAGGACCAGTTCGCCATCGCCTCGCTGACCCGAGCGCAGGAAGCCATCAGCAGCGGCCGTTTTGCCGCCGAGATCGTGCCGGTGGAAGTCACCGAGGGCAAGGAAAAGCGCGTCATCAAGGATGACGagca                                                                                                                                          <  8:184521/132‑1 (MQ=255)
                      cAGGACCAGTTCGCCATCGCCTCGCTGACCCGAGCGCAGGAAGCCATCAGCAGCGGCCGTTTTGCCGCCGAGATCGTGCCGGTGGAAGTCACCGAGGGCAAGGAAAAGCGCGTCATCAAGGATGACGAGCAGCCGCCCAAGGCGCGTCt                                                                                                                        >  7:115949/1‑149 (MQ=255)
                            cAGTTCGCCATCGCCTCGCTGACCCGAGCGCAGGAAGCCATCAGCAGCGGCCGTTTTGCCGCCGAGATCGTGCCGGTGGAAGTCACCGAGGGCAAGGAAAAGCGCGTCATCAAGGATGACGAGCAGCCGCCCAAGGCGCGTCTGGAc                                                                                                                    >  3:6279/1‑147 (MQ=255)
                               ttCGCCATCGCCTCGCTGACCCGAGCGCAGGAAGCCATCAGCAGCGGCCGTTTTGCCGCCGAGATCGTGCCGGTGGAAGTCACCGAGGGCAAGGAAAAGCGCGTCATCAAGGATGACGAGCAGCCGCCCAAGGCGCGTCTGGACAAGAt                                                                                                               >  1:16699/1‑149 (MQ=255)
                                   ccATCGCCTCGCTGACCCGAGCGCAGGATGCCATCAGCAGCGGCCGTTTTGCCGCCGAGATCGTGCCGGTGGAAGTCACCGAGGGCAAGGAAAAGCGCGTCATCAAGGATGACGAGCAGCCGCCCAAGGCGCGTCTGGACAAGATTgcg                                                                                                           >  1:196811/1‑149 (MQ=255)
                                   ccATCGCCTCGCTGACCCGAGCGCAGGAAGCCATCAGCAGCGGCCGTTTTGCCGCCGAGATCGTGCCGGTGGAAGTCACCGAGGGCAAGGAAAAGCGCGTCATCAAGGATGACGAGCAGCCGCCCAAGGCGCGTCTGGACAAGAtt                                                                                                              >  2:289995/1‑146 (MQ=255)
                                   ccATCGCCTCGCTGACCCGAGCGCAGGAAGCCATCAGCAGCGGCCGTTTTGCCGCCGAGATCGTGCCGGTGGAAGTCACCGAGGGCAAGGAAAAGCGCGTCATCAAGGATGACGAGCAGCCGCCCAAGGCGCGTCTGGACAAGAtt                                                                                                              <  1:289995/146‑1 (MQ=255)
                                         ccTCGCTGACCCGAGCGCAGGAAGCCATCAGCAGCGGCCGTTTTGCCGCCGAGATCGTGCCGGTGGAAGTCACCGAGGGCAAGGAAAAGCGCGTCATCAAgg                                                                                                                                                    <  7:322467/102‑1 (MQ=255)
                                         ccTCGCTGACCCGAGCGCAGGAAGCCATCAGCAGCGGCCGTTTTGCCGCCGAGATCGTGCCGGTGGAAGTCACCGAGGGCAAGGAAAAGCGCGTCATCAAgg                                                                                                                                                    >  8:322467/1‑102 (MQ=255)
                                         ccTCGCTGACCCGAGCGCAGGAAGCCATCAGCAGCGGCCGTTTTGCCGCCGAGATCGTGCCGGTGGAAGTCACCGAGGGCAAGGAAAAGCGCGTCATCAAGGATGACGAGCAGCCGCCCAAGGCGCGTCTGGACAAGATTGCGCAGCt                                                                                                      >  7:52781/1‑148 (MQ=255)
                                                 aCCCGAGCGCAGGAAGCCATCAGCAGCGGCCGTTTTGCCGCCGAGATCGTGCCGGTGGAAGTCACCGAGGGCAAGGAAAAGCGCGTCATCAAGGATGACGAGCAGCCGCCCAAGGCGCGTCTGGACAAGATTGCGCAGCTCAAACCgg                                                                                             <  8:115949/148‑1 (MQ=255)
                                                   ccGAGCGCAGGAAGCCATCAGCAGCGGCCGTTTTGCCGCCGAGATCGTGCCGGTGGAAGTCACCGAGGGCAAGGAAAAGCGCGTCATCAAGGATGACGAGCAGCCGCCCAAGGCGCGTCTGGACAAGATTGCGCAGCTCAAACCGGCGt                                                                                          >  1:84820/1‑149 (MQ=255)
                                                   ccGAGCGCAGGAAGCCATCAGCAGCGGCCGTTTTGCCGCCGAGATCGTGCCGGTGGAAGTCACCGAGGGCAAGGAAAAGCGCGTCATCAAGGATGACGAGCAGCCGCCCAAGGCGCGTCTGGACAAGATTGCGCAGCTCAAACCGGCGt                                                                                          <  6:177371/149‑1 (MQ=255)
                                                   ccGAGCGCAGGAAGCCATCAGCAGCGGCCGTTTTGCCGCCGAGATCGTGCCGGTGGAAGTCACCGAGGGCAAGGAAAAGCGCGTCATCAAGGATGACGAGCAGCCGCCCAAGGCGCGTCTGGACAAGATTGCGCAGCTCAAACCGGCGt                                                                                          >  6:98779/1‑149 (MQ=255)
                                                         gcAGGAAGCCATCAGCAGCGGCCGTTTTGCCGCCGAGATCGTGCCGGTGGAAGTCACCGAGGGCAAGGAAAAGCGCGTCATCAAGGATGACGAGCAGCCGCCCAAGGCGCGTCTGGACAAGATTGCGCAGCTCAAACCGGCGTTTCGTg                                                                                    <  4:6279/149‑1 (MQ=255)
                                                          cAGGAAGCCATCAGCAGCGGCCGTTTTGCCGCCGAGATCGTGCCGGTGGAAGTCACCGAGGGCAAGGAAAAGCGCGTCATCAAGGATGAc                                                                                                                                               <  4:185806/90‑1 (MQ=255)
                                                          cAGGAAGCCATCAGCAGCGGCCGTTTTGCCGCCGAGATCGTGCCGGTGGAAGTCACCGAGGGCAAGGAAAAGCGCGTCATCAAGGATGAc                                                                                                                                               >  3:185806/1‑90 (MQ=255)
                                                           aGGAAGCCATCAGCAGCGGCCGTTTTGCAGCCGAGATCGTGCCGGTGGAAGTCCCCGAGGGCAATGAAAAGCGCGTCATCAAGGATGACGTGCAGCCGCCCAAGGCGCGTCTGGACAAGATTGGGCAGCTCAATCCGGCGTTTCGTGaa                                                                                  <  7:300379/149‑1 (MQ=255)
                                                               aGCCATCAGCAGCGGCCGTTTTGCCGCCGAGATCGTGCCGGTGGAAGTCACCGAGGGCAAGGAAAAGCGCGTCATCAAGGATGACGAGCAGCCGCCCAAGGCGCGTCTGGACAAGATTGCGCAGCTCAAACCGGCGTTTCGTGAAggcg                                                                              <  7:193877/149‑1 (MQ=255)
                                                                  cATCAGCAGCGGCCGTTTTGCCGCCGAGATCGTGCCGGTGGAAGTCACCGAGGGCAAGGAAAAGCGCGTCATCAAGGATGACGAGCAGCCGCCCAAGGCGCGTCTGGACAAGATTGCGCAGCTCAAACCGGCGTTTCGTGAAGGCGGCa                                                                           <  8:52781/149‑1 (MQ=255)
                                                                  cATCAGCAGCGGCCGTTTTGCCGCCGAGATCGTGCCGGTGGAAGTCACCGAGGGCAAGGAAAAGCGCGTCATCAAGGATGACGAGCAGCCGCCCAAGGCGCGTCTGGACAAGATTGCGCAGCTCAAACCGGCGTTTCGTGAAGGCGGCa                                                                           >  5:129796/1‑149 (MQ=255)
                                                                   aTCAGCAGCGGCCGTTTTGCCGCCGAGATCGTGCCGGTGGAAGTCACCGAGGGCAAGGAAAAGCGCGTCATCAAGGATGACGAGCAGCCGCCCAAGGCGCGTCTGGACAAGATTGCGCAGCTCAAACCGGCGTTTCGTGAAGGCGGCAc                                                                          >  7:222938/1‑149 (MQ=255)
                                                                     cagcagCGGCCGTTTTGCCGCCGAGATCGTGCCGGTGGAAGTCACCGAGGGCAAGGAAAAGCGCGTCATCAAGGATGACGAGCAGCCGCCCAAGGCGCGTCTGGACAAGATTGCGCAGCTCAAACCGGCGTTTCGTGAAGGCGGCACCg                                                                        >  4:233765/1‑149 (MQ=255)
                                                                        cagcGGCCGTTTTGCCGCCGAGATCGTGCCGGTGGAAGTCACCGAGGGCAAGGAAAAGCGCGTCATCAAGGATGACGAGCAGCCGCCCAAGGCGCGTCTGGACAAGATTGCGCAGCTCAAACCGGCGTTTCGTGAAGGCGGCACCGTGa                                                                     <  3:233765/149‑1 (MQ=255)
                                                                        cagcGGCCGTTTTGCCGCCGAGATCGTGCCGGTGGAAGTCACCGAGGGCAAGGAAAAGCGCGTCATCAAGGATGACGAGCAGCCGCCCAAGGCGCGTCTGGACAAGATTGCGCAGCTCAAACCGGCGTTTCGTGAAGGCGGCACCGTGa                                                                     <  2:16699/149‑1 (MQ=255)
                                                                          gcGGCCGTTTTGCCGCCGAGATCGTGCCGGTGGAAGTCACCGAGGGCAAGGAAAAGCGCGTCATCAAGGATGACGAGCAGCCGCCCAAGGCGCGTCTGGACAAGATTGCGCAGCTCAAACCGGCGTTTCGTGAAGGCGGCACCGTGAcg                                                                   <  6:162455/149‑1 (MQ=255)
                                                                            ggCCGTTTTGCCGCCGAGATCGTGCCGGTGGAAGTCACCGAGGGCAAGGAAAAGCGCGTCATCAAGGATGACGAGCAGCCGCCCAAGGCGCGTCTGGACAAGATTGCGCAGCTCAAACCGGCGTTTCGTGAAGGCGGCACCGTGAcggc                                                                 <  5:98779/149‑1 (MQ=255)
                                                                             gCCGTTTTGCCGCCGAGATCGTGCCGGTGGAAGTCACCGAGGGCAAGGAAAAGCGCGTCATCAAGGATGACGagca                                                                                                                                          <  5:104834/76‑1 (MQ=255)
                                                                             gCCGTTTTGCCGCCGAGATCGTGCCGGTGGAAGTCACCGAGGGCAAGGAAAAGCGCGTCATCAAGGATGACGagca                                                                                                                                          >  6:104834/1‑76 (MQ=255)
                                                                             gCCGTTTTGCCGCCGAGATCGTGCCGGTGGAAGTCACCGAGGGCAAGGAAAAGCGCGTCATCAAGGATGACGAGCAGCCGCCCAAGGCGCGTCTGGACAAGATTGCGCAGCTCAAACCGGCGTTTCGTGAAGGCGGCACCGTGAcggcg                                                                <  2:196811/149‑1 (MQ=255)
                                                                               cGTTTTGCCGCCGAGATCGTGCCGGTGGAAGTCACCGAGGGCAAGGAAAAGCGCGTCATCAAGGATGACGAGCAGCCGCCCAAGGCGCGTCTGGACAAGATTGCGCAGCTCAAACCGGCGTTTCGTGAAggcg                                                                              <  1:14792/133‑1 (MQ=255)
                                                                               cGTTTTGCCGCCGAGATCGTGCCGGTGGAAGTCACCGAGGGCAAGGAAAAGCGCGTCATCAAGGATGACGAGCAGCCGCCCAAGGCGCGTCTGGACAAGATTGCGCAGCTCAAACCGGCGTTTCGTGAAggcg                                                                              >  2:14792/1‑133 (MQ=255)
                                                                                      ccgccgAGATCGTGCCGGTGGAAGTCACCGAGGGCAAGGAAAAGCGCGTCATCAAGGATGACGAGCAGCCGCCCAAGGCGCGTCTGGACAAGATTGCGCAGCTCAAACCGGCGTTTCGTGAAGGCGGCACCGTGACGGCGGCCAACGc                                                        <  4:126538/148‑1 (MQ=255)
                                                                                       cgccgAGATCGTGCCGGAGGAAGTCACCGAGGGCAAGGAAAAGCGCGTCATCAAGGATGACGAGCAGCCGCCCAAGGCGCGTCTGGACAAGATTGCGCAGCTCAAACCGGCGTTTCGTGAAGGCGGCACCGTGACGGCGGCCAACGc                                                        >  6:6226/1‑147 (MQ=255)
                                                                                          cgAGATCGTGCCGGTGGAAGTCACCGAGGGCAAGGAAAAGCGCGTCATCAAGGATGACGAGCAGCCGCCCAAGGCGCGTCTGGACAAGATTGCGCAGCTCAAACCGGCGTTTCGTGAAGGCGGCACCGTGACGGCGGCCAACGCCAGtt                                                   <  6:205245/149‑1 (MQ=255)
                                                                                           gagaTCGTGCCGGTGGAAGTCACCGAGGGCAAGGAAAAGCGCGTCATCAAGGATGACGAGCAGCCGCCCAAGGCGCGTCTGGACAAGATTGCGCAGCTCAAACCGGCGTTTCGTGAAGGCGGCACCGTGACGGCGGCCAACGCCAGTTc                                                  >  4:320734/1‑149 (MQ=255)
                                                                                                  tGCCGGTGGAAGTCACCGAGGGCAAGGAAAAGCGCGTCATCAAGGATGACGAGCAGCCGCCCAAGGCGCGTCTGGACAAGATTGCGCAGCTCAAACCGGCGTTTCGTGAAGGCGGCACCGTGACGGCGGCCAACGCCAGTTCGATTTcc                                           <  2:194827/149‑1 (MQ=255)
                                                                                                       gTGGAAGTCACCGAGGGCAAGGAAAAGCGCGTCATCAAGGATGACGAGCAGCCGCCCAAGGCGCGTCTGGACAAGATTGCGCAGCTCAAACCGGCGTTTCGTGAAGGCGGCACCGTGACGGCGGCCAACGCCAGTTCGATTTCCGAcgg                                      >  6:227775/1‑149 (MQ=255)
                                                                                                                   gAGGGCAAGGAAAAGCGCGTCATCAAGGATGACGAGCAGCCGCCCAAGGCGCGTCTGGACAAGATTGCGCAGCTCAAACCGGCGTTTCGTGAAGGCGGCACCGTGACGGCGGCCAACGCCAGTTCGATTTCCGAcggcgctgcggcgct                          <  5:227775/149‑1 (MQ=255)
                                                                                                                    aGGGCAAGGAAAAGCGCGTCATCAAGGATGACGAGCAGCCGCCCAAGGCGCGTCTGGACAAGATTGCGCAGCTCAAACCGGCGTTTCGTGAAGGCGGCACCGTGACGGCGGCCAACGCCAGTTCGATTTCCGAcggcgctgcggcgctg                         <  8:222938/149‑1 (MQ=255)
                                                                                                                         aaGGAAAAGCGCGTCATCAAGGATGACGAGCAGCCGCCCAAGGCGCGTCTGGACAAGATTGCGCAGCTCAAACCGGCGTTTCGTGAAGGCGGCACCGTGACGGCGGCCAACGCCAGTTCGATTTCCGACGGCGCTGCGGCGCTGGt                       <  8:136352/146‑1 (MQ=255)
                                                                                                                         aaGGAAAAGCGCGTCATCAAGGATGACGAGCAGCCGCCCAAGGCGCGTCTGGACAAGATTGCGCAGCTCAAACCGGCGTTTCGTGAAGGCGGCACCGTGACGGCGGCCAACGCCAGTTCGATTTCCGACGGCGCTGCGGCGCTGGTaag                    >  7:136352/1‑147 (MQ=255)
                                                                                                                         aaGGAAAAGCGCGTCATCAAGGACGACGAGCAGCCGCCCAAGGCGCGTCTGGACAAGATTGCGCAGCTCAAACCGGCGTTTCATGAAGGCGGCACCGTGACGGCGGCCAACGCCAGTTCGATTTCCGACGGCGCTGCGGCGCTGGTACt                    <  3:46987/149‑1 (MQ=255)
                                                                                                                           ggAAAAGCGCGTCATCAAGGATGACGAGCAGCCGCCCAAGGCGCGTCTGGACAAGATTGCGCAGCTCAAACCGGCGTTTCGTGAAGGCGGCACCGTGACGGCGGCCAACGCCAGTTCGATTTCCGACGGCGCTGCGGCGCTGGTACTg                   <  7:29379/148‑1 (MQ=255)
                                                                                                                           ggAAAAGCGCGTCATCAAGGATGACGAGCAGCCGCCCAAGGCGCGTCTGGACAAGATTGCGCAGCTCAAACCGGCGTTTCGTGAAGGCGGCACCGTGACGGCGGCCAACGCCAGTTCGATTTCCGACGGCGCTGCGGCGCTGGTACTg                   >  8:29379/1‑148 (MQ=255)
                                                                                                                           ggAAAAGCGCGTCATCAAGGATGACGAGCAGCCGCCCAAGGCGCGTCTGGACAAGATTGCGCAGCTCAAACCGGCGTTTCGTGAAGGCGGCACCGTGACGGCGGCCAACGCCAGTTCGATTTCCGACGGCGCTGCGGCGCTGGTACTGa                  <  2:84820/149‑1 (MQ=255)
                                                                                                                              aaaGCGCGTCATCAAGGATGACGAGCAGCCGCCCAAGGCGCGTCTGGACAAGATTGCGCAGCTCAAACCGGCGTTTCGTGAAGGCGGCACCGTGACGGCGGCCAACGCCAGTTCGATTTCCGACGGCGCTGCGGCGCTGGTACTGATgc               <  6:129796/149‑1 (MQ=255)
                                                                                                                                 gctcgTCATCCTTGATGACGCGCAGCCGCCCAAGGCGCGTCTGGACAAGATTGCGCAGCTCAAACCGGCGTTTCGTGAAGGCGGCACCGTGACGGCGGCCAACGc                                                        <  2:290592/102‑1 (MQ=255)
                                                                                                                                  ctcgTCATCCTTGATGACGAGCAGCCGCCCAAGGCGCGTCTGGACAAGATTGCGCAGCTCAAACCGGCGTTTCGTGAAGGCGGCACCGTGACGGCGGCCAACGCCAGTTCGATTTCCGACGGCGCTGCGGCGCTGGTACTGATGcgccg           >  6:353837/3‑149 (MQ=255)
                                                                                                                                     gTCATCAAGGATGACGAGCAGCCGCCCAAGGCGCGTCTGGACAAGATTGCGCAGCTCAAACCGGCGTTTCGTGAAGGCGGCACCGTGAcggcgg                                                               <  2:16537/94‑1 (MQ=255)
                                                                                                                                     gTCATCAAGGATGACGAGCAGCCGCCCAAGGCGCGTCTGGACAAGATTGCGCAGCTCAAACCGGCGTTTCGTGAAGGCGGCACCGTGAcggcgg                                                               >  1:16537/1‑94 (MQ=255)
                                                                                                                                         tcaAGGATGACGAGCAGCCGCCCAAGGCGCGTCTGGACAAGATTGCGCAGCTCGAAACCGGCGTTTCGTGAAGGCGGCACCGTGACGGCGGCCAACGCCAGTTCGATTTCCGACGGCGCTGCGGCGCTGGTACTGATGCGCCGCTcc       >  1:89640/1‑147 (MQ=255)
                                                                                                                                         tcaAGGATGACGAGCAGCCGCCCAAGGCGCGTCTGGACAAGATTGCGCAGCTCAAACCGGCGTTTCGTGAAGGCGGCACCGTGACGGCGGCCAACGCCAGTTCGATTTCCGACGGCGCTGCGGCGCTGGTACTGATGCGCCGCTCCGAg    >  7:202034/1‑149 (MQ=255)
                                                                                                                                          caAGGATGACGAGCAGCCGCCCAAGGCGCGTCTGGACAAGATTGCGCAGCTCAAACCGGCGTTTCGTGAAGGCGGCACCGTGACGGCGGCCAACGCCAGTTCGATTTCCGAcggcgctgcggcgc                           >  2:11272/1‑125 (MQ=255)
                                                                                                                                          caAGGATGACGAGCAGCCGCCCAAGGCGCGTCTGGACAAGATTGCGCAGCTCAAACCGGCGTTTCGTGAAGGCGGCACCGTGACGGCGGCCAACGCCAGTTCGATTTCCGAcggcgctgcggcgc                           <  1:11272/125‑1 (MQ=255)
                                                                                                                                          caAGGATGACGAGCAGCCGCCCAAGGCGCGTCTGGACAAGATTGCGCAGCTCAAACCGGCGTTTCGTGAAGGCGGCACCGTGACGGCGGCCAACGCCAGTTCGATTTCCGACGGCGCTGCGGCGCTGGTACTGATGCGCCGCTCCGAgg   >  1:228268/1‑149 (MQ=255)
                                                                                                                                           aaGGATGACGAGCAGCCGCCCAAGGCGCGTCTGGACAAGATTGCGCAGCTCAAACCGGCGTTTCGTGAAGGCGGCACCGTGACGGCGGCCAACGCCAGTTCGATTTCCGACGGCGCTGCGGCGCTGGTACTGATGCGCCGCTCCGAGGc  >  3:155684/1‑149 (MQ=255)
                                                                                                                                            |                                                                                                                                                    
CAATGCCTTCAGCCGCGAGGCCCAGGACCAGTTCGCCATCGCCTCGCTGACCCGAGCGCAGGAAGCCATCAGCAGCGGCCGTTTTGCCGCCGAGATCGTGCCGGTGGAAGTCACCGAGGGCAAGGAAAAGCGCGTCATCAAGGATGACGAGCAGCCGCCCAAGGCGCGTCTGGACAAGATTGCGCAGCTCAAACCGGCGTTTCGTGAAGGCGGCACCGTGACGGCGGCCAACGCCAGTTCGATTTCCGACGGCGCTGCGGCGCTGGTACTGATGCGCCGCTCCGAGGC  >  NC_002947/2526272‑2526559

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: