Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 863,383 A→C 16.6% S172R (AGC→CGC)  hemH → ferrochelatase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_002947863,3830AC16.6% 1368.2 / 83.7 614S172R (AGC→CGC) hemHferrochelatase
Reads supporting (aligned to +/- strand):  ref base A (386/107);  new base C (12/106);  total (398/216)
Fisher's exact test for biased strand distribution p-value = 3.08e-44
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 680 total aligned reads displayed.

CACCGCTTTATCCACAGTTCGCCGACAGCACGGTGACCACTGTGATCGAGCAGGCAAAACAGACGGTCAGCGAGCACCAGTTACCGTTGCAGATGCGCGTGCTGCAGCCGTTCTATGAGCACCCGGCGTACATCGAGGCGCTGGCAGCCAGCGCCCGGCCGTACCTGGAGCAGGGTTACGACCACCTGCTACTGAGCTTCCATGGTTTGCCGGAACGGCACCTGAAGAAACTGTTTCCGAAGGGGGTCAAGCACGACCTGCGCGCTGCCGATTGCTGCCACGGTGCAACCGCCGAAGTA  >  NC_002947/863234‑863532
                                                                                                                                                     |                                                                                                                                                     
 aCCGCTTTATCCACAGTTCGCCGACAGCACGGTGACCACTGTGATCGAGCAGGCAAAACAGACGGTCAGCGAGCACCAGTTACCGTTGCAGATGCGCGTGCTGCAGCCGTTCTATGAGCACCCGGCGTACATCGAGGCGCTGGcagccag                                                                                                                                                      >  1:8623373/1‑150 (MQ=255)
 aCCGCTTTATCCACAGTTCGCCGACAGCACGGTGACCACTGTGATCGAGCAGGCAAAACAGACGGTCAGCGAGCACCAGTTACCGTTGCAGATGCGCGTGCTGCAGCCGTTCTATGAGCACCCGGCGTACATCGAGGCGCTGGcagccag                                                                                                                                                      >  1:5666723/1‑150 (MQ=255)
  ccGCTTTATCCACAGTTCGCCGACAGCACGGTGACCACTGTGATCGAGCAGGCAAAACAGACGGTCAGCGAGCACCAGTTACCGTTGCAGATGCGCGTGCTGCAGCCGTTCTATGAGCACCCGGCGTACATCGAGGCGCTGGcagccagc                                                                                                                                                     <  2:2648020/150‑1 (MQ=255)
    gCTTTATCCACAGTTCGCCGACAGCACGGTGACCACTGTGATCGAGCAGGCAAAACAGACGGTCAGCGAGCACCAGTTACCGTTGCAGATGCGCGTGCTGCAGCCGTTCTATGAGCACCCGGCGTACATCGAGGCGCTGGCAGCCAgcgc                                                                                                                                                   >  1:5138773/1‑150 (MQ=255)
     cTTTATCCACAGTTCGCCGACAGCACGGTGACCACTGTGATCGAGCAGGCAAAACAGACGGTCAGCGAGCACCAGTTACCGTTGCAGATGCGCGTGCTGCAGCCGTTCTATGAGCACCCGGCGTACATCGAGGCGCTGGCTGCCAGCGcc                                                                                                                                                  <  1:183206/150‑1 (MQ=255)
     cTTTATCCACAGTTCGCCGACAGCACGGTGACCACTGTGATCGAGCAGGCAAAACAGACGGTCAGCGAGCACCAGTTACCGTTGCAGATGCGCGTGCTGCAGCCGTTCTATGAGCACCCGGCGTACATCGAGGCGCTGGCTGCCAGCGcc                                                                                                                                                  <  1:4813264/150‑1 (MQ=255)
      tttATCCACAGTTCGCCGACAGCACGGTGACCACTGTGATCGAGCAGGCAAAACAGACGGTCAGCGAGCACCAGTTACCGTTGCAGATGCGCGTGCTGCAGCCGTTCTATGAGCACCCGGCGTACATCGAGGCGCTGGCTGCCAGCGcct                                                                                                                                                 <  1:4917512/150‑2 (MQ=255)
      tttATCCACAGTTCGCCGACAGCACGGTGACCACTGTGATCGAGCAGGCAAAACAGACGGTCAGCGAGCACCAGTTACCGTTGCAGATGCGCGTGCTGCAGCCGTTCTATGAGCACCCGGCGTACATCGAGGCGCTGGCAGCCAGCGccc                                                                                                                                                 <  2:6363639/150‑1 (MQ=255)
       ttATCCACAGTTCGCCGACAGCACGGTGACCACTGTGATCGAGCAGGCAAAACAGACGGTCAGCGAGCACCAGTTCCCGTTGCAGATGCGCGTGCTGCAGCCGTTCTATGAGCACCCGGCGTACATCGAGGCGCTGGCTGCCAGCGcctc                                                                                                                                                <  2:6768892/150‑3 (MQ=255)
       ttATCCACAGTTCGCCGACAGCACGGTGACCACTGTGATCGAGCAGGCAAAACAGACGGTCAGCGAGCACCAGTTACCGTTGCAGATGCGCGTGCTGCAGCCGTTCTATGAGCACCCGGCGTACATCGAGGCGCTGGCTGCCAGCGcctc                                                                                                                                                <  2:2135103/150‑3 (MQ=255)
        tATCCACAGTTCGCCGACAGCACGGTGACCACTGTGATCGAGCAGGCAAAACAGACGGTCAGCGAGCACCAGTTACCGTTGCAGATGCGCGTGCTGCAGCCGTTCTATGAGCACCCGGCGTACATCGGGGCGCTGGCTGCCAGCGcctcg                                                                                                                                               <  1:7144648/150‑4 (MQ=255)
        tATCCACAGTTCGCCGACAGCACGGTGACCACTGTGATCGAGCAGGCAAAACAGACGGTCAGCGAGCACCAGTTACCGTTGCAGATGCGCGTGCTGCAGCCGTTCTATGAGCACCCGGCGTACATCGAGGCGCTGGCTGCCAGCGcctcg                                                                                                                                               <  1:8323465/150‑4 (MQ=255)
         ctCCCCAGTTCGCCGACAGCACGGTGACCACTGTGATCGAGCAGGCAAACCAGACGGTCAGCGAGCACCAGTTACCGTTGCAGATGCGCGTGCTGCAGCCGTTCTATGAGCACCCGGCGTACATCGAGGCGCTGGCTGCCAGCGcctcga                                                                                                                                              <  2:1931339/149‑5 (MQ=255)
         aTCCACTGTTCGCCGACAGCACGGTGACCACTGTGATGGAGCAGGCAAAACAGACGGTCAGCGAGCACCAGTTACCGTTGCAGATGCGCGTGCTGCAGCCGTTATATGAGCACCCGGCGTACATCGAGGCGCTGGCAGCCAGCGCCCGGc                                                                                                                                              <  2:1250180/150‑1 (MQ=255)
          tCCACAGTTCGCCGACAGCACGGTGACCACTGTGATCGAGCAGGCAAAACAGACGGTCAGCGAGCACCAGTTACCGTTGCAGATGCGCGTGCTGCAGCCGTTCTATGAGCACCCGGCGTACATCGAGGCGCTGGCAGCCCGCGCCCGGcc                                                                                                                                             >  1:6375693/1‑150 (MQ=255)
          tCCACAGTTCGCCGACAGCACGGTGACCACTGTGATCGAGCAGGCAAAACAGACGGTCAGCGAGCACCAGTTACCGTTGCAGATGCGCGTGCTGCAGCCGTTCTATGAGCACCCGGCGTACATCGAGGCGCTGGCAGCCAGCGCCCGGcc                                                                                                                                             >  2:5285572/1‑150 (MQ=255)
          tCCACAGTTCGCCGACAGCACGGTGACCACTGTGATCGAGCAGGCAAAACAGACGGTCAGCGAGCACCAGTTACCGTTGCAGATGCGCGTGCTGCAGCCGTTCTATGAGCACCCGGCGTACATCGAGGCGCTGGCAGCCAGCGCCCGGcc                                                                                                                                             >  2:7064959/1‑150 (MQ=255)
           ccACAGTTCGCCGACAGCACGGTGACCACTGTGATCGAGCAGGCAAAACAGACGGTCAGCGAGCACCAGTTACCGTTGCAGATGCGCGTGCTGCAGCCGTTCTATGAGCACCCGGCGTACATCGAGGCGCTGGCAGCCCGCGCCCGGCCg                                                                                                                                            >  1:3229531/1‑150 (MQ=255)
           ccACAGTTCGCCGACAGCACGGTGACCACTGTGATCGAGCAGGCAAAACAGACGGTCAGCGAGCACCAGTTACCGTTGCAGATGCGCGTGCTGCAGCCGTTCTATGAGCACCCGGCGTACATCGAGGCGCTGGCAGCCAGCGCCCGGCCg                                                                                                                                            >  1:2251271/1‑150 (MQ=255)
            cacaGTTCGCCGACAGCACGGTGACCACTGTGATCGAGCAGGCAAAACAGACGGTCAGCGAGCACCAGTTACCGTTGCAGATGCGCGTGCTGCAGCCGTTCTATGAGCACCCGGCGTACATCGAGGCGCTGGCAGCCAGCGCCCGGCCgg                                                                                                                                           >  1:3551464/1‑149 (MQ=255)
            cacaGTTCGCCGACAGCACGGTGACCACTGTGATCGAGCAGGCAAAACAGACGGTCAGCGAGCACCAGTTACCGTTGCAGATGCGCGTGCTGCAGCCGTTCTATGAGCACCCGGCGTACATCGAGGCGCTGGCAGCCAGCGCCCGGCCGt                                                                                                                                           >  2:4665861/1‑150 (MQ=255)
            cacaGTTCGCCGACAGCACGGTGACCACTGTGATCGAGCAGGCAAAACAGACGGTCAGCGAGCACCAGTTACCGTTGCAGATGCGCGTGCTGCAGCCGTTCTATGAGCACCCGGCGTACATCGAGGCGCTGGCAGCCAGCGCCCGGCCGt                                                                                                                                           >  2:31543/1‑150 (MQ=255)
             acaGTTCGCCGACAGCACGGTGACCACTGTGATCGAGCAGGCAAAACAGACGGTCAGCGAGCACCAGTTACCGTTGCAGATGCGCGTGCTGCAGCCGTTCTATGAGCACCCGGCGTACATCGAGGCGCTGGCAGCCAGCGCCCGGCCGTa                                                                                                                                          >  1:7476995/1‑150 (MQ=255)
             acaGTTCGCCGACAGCACGGTGACCACTGTGATCGAGCAGGCAAAACAGACGGTCAGCGAGCACCAGTTACCGTTGCAGATGCGCGTGCTGCAGCCGTTCTATGAGCACCCGGCGTACATCGAGGCGCTGGCAGCCAGCGCCCGGCCGTa                                                                                                                                          >  2:3907716/1‑150 (MQ=255)
              caGTTCGCCGACAGCACGGTGACCACTGTGATCGAGCAGGCAAAACAGACGGTCAGCGAGCACCAGTTACCGTTGCAGATGCGCGTGCTGCAGCCGTTCTATGAGCACCCGGCGTACATCGAGGCGCTGGCAGCCAGCGCCCGGCCGGAc                                                                                                                                         >  1:6275384/1‑150 (MQ=255)
                gtgcGCCGACAGCACGGTGACCACTGTGATCGAGCCGGCAAACCAGACGGTCAGCGAGCACCTGTTCCGGTTGCTGTGGGGCGTGCTGCAGCGTTTCTATGAGCACCCGGCGTACATCGCGGCGCTGGCCCCCAGCGCCCGGCCGTACCt                                                                                                                                       <  2:6007284/147‑1 (MQ=255)
                     ccGACAGCACGGTGACCACTGTGATCGAGCAGGCAAAACAGACGGTCAGCGAGCACCAGTTACCGTTGCAGATGCGCGTGCTGCAGCCGTTCTATGAGCACCCGGCGTACATCGAGGCGCTGGCAGCCAGCGCCCGGCCGTACCTGGAGc                                                                                                                                  >  1:3043981/1‑150 (MQ=255)
                     ccGACAGCACGGTGACCACTGTGATCGAGCAGGCAAAACAGACGGTCAGCGAGCACCAGTTACCGTTGCAGATGCGCGTGCTGCAGCCGTTCTATGAGCACCCGGCGTACATCGAGGCGCTGGCAGCCAGCGCCCGGCCGTACCTGGAGc                                                                                                                                  >  1:8801364/1‑150 (MQ=255)
                          aGCACGGTGACCACTGTGATCGAGCAGGCAAAACAGACGGTCAGCGAGCACCAGTTACCGTTGCAGATGCGCGTGCTGCAGCCGTTCTATGAGCACCCGGCGTACATCGAGGCGCTGGCAGCCAGCGCCCGGCCGTACCTGGAGCAGGGt                                                                                                                             >  2:2993652/1‑150 (MQ=255)
                            cACGGTGACCACTGTGATCGAGCAGGCAAAACAGACGGTCAGCGAGCACCAGTTACCGTTGCAGATGCGCGTGCTGCAGCCGTTCTATGAGCACCCGGCGTACATCGAGGCGCTGGCAGCCAGCGCCCGGCCGTACCTGGAGCAGGGTTa                                                                                                                           >  2:6835718/1‑150 (MQ=255)
                            cACGGTGACCACTGTGATCGAGCAGGCAAAACAGACGGTCAGCGAGCACCAGTTACCGTTGCAGATGCGCGTGCTGCAGCCGTTCTATGAGCACCCGGCGTACATCGAGGCGCTGGCAGCCAGCGCCCGGCCGTACCTGGAGCAGGGTTa                                                                                                                           >  2:85537/1‑150 (MQ=255)
                             aCGGTGACCACTGTGATCGAGCAGGCAAAACAGACGGTCAGCGAGCACCAGTTACCGTTGCAGATGCGCGTGCTGCAGCCGTTCTATGAGCACCCGGCGTACATCGAGGCGCTGGCAGCCAGCGCCCGGCCGTACCTGGAGCAGGGTTAc                                                                                                                          >  2:2781974/1‑150 (MQ=255)
                             aCGGTGACCACTGTGATCGAGCAGGCAAAACAGACGGTCAGCGAGCACCAGTTACCGTTGCAGATGCGCGTGCTGCAGCCGTTCTATGAGCACCCGGCGTACATCGAGGCGCTGGCAGCCAGCGCCCGGCCGTACCTGGAGCAGGGTTAc                                                                                                                          >  2:1596100/1‑150 (MQ=255)
                                        tgtgATCGAGCAGGCAAAACAGACGGTCAGCGAGCACCAGTTACCGTTGCAGATGCGCGTGCTGCAGCCGTTCTATGAGCACCCGGCGTACATCGAGGCGCTGGCAGCCAGCGCCCGGCCGTACCTGGAGCAGGGTTACGACCACCTGCt                                                                                                               >  2:6281258/1‑150 (MQ=255)
                                        tgtgATCGAGCAGGCAAAACAGACGGTCAGCGAGCACCAGTTACCGTTGCAGATGCGCGTGCTGCAGCCGTTCTATGAGCACCCGGCGTACATCGAGGCGCTGGCAGCCAGCGCCCGGCCGTACCTGGAGCAGGGTTACGACCACCTGCt                                                                                                               >  2:378087/1‑150 (MQ=255)
                                        tgtgATCGAGCAGGCAAAACAGACGGTCAGCGAGCACCAGTTACCGTTGCAGATGCGCGTGCTGCAGCCGTTCTATGAGCACCCGGCGTACATCGAGGCGCTGGCAGCCAGCGCCCGGCCGTACCTGGAGCAGGGTTACGACCACCTGCt                                                                                                               >  2:8639251/1‑150 (MQ=255)
                                        tgtgATCGAGCAGGCAAAACAGACGGTCAGCGAGCACCAGTTACCGTTGCAGATGCGCGTGCTGCAGCCGTTCTATGAGCACCCGGCGTACATCGAGGCGCTGGCAGCCAGCGCCCGGCCGTACCTGGAGCAGGGTTACGACCACCTGCt                                                                                                               >  2:8790164/1‑150 (MQ=255)
                                        tgtgATCGAGCAGGCAAAACAGACGGTCAGCGAGCACCAGTTACCGTTGCAGATGCGCGTGCTGCAGCCGTTCTATGAGCACCCGGCGTACATCGAGGCGCTGGCAGCCAGCGCCCGGCCGTACCTGGAGCAGGGTTACGACCACCTGCt                                                                                                               >  1:7834770/1‑150 (MQ=255)
                                            aTCGAGCAGGCAAAACAGACGGTCAGCGAGCACCAGTTACCGTTGCAGATGCGCGTGCTGCAGCCGTTCTATGAGCACCCGGCGTACATCGAGGCGCTGGCAGCCAGCGCCCGGCCGTACCTGGAGCAGGGTTACGACCACCTGCTACTg                                                                                                           >  1:7586742/1‑150 (MQ=255)
                                             tCGAGCAGGCAAAACAGACGGTCAGCGAGCACCAGTTACCGTTGCAGATGCGCGTGCTGCAGCCGTTCTATGAGCACCCGGCGTACATCGAGGCGCTGGCAGCCAGCGCCCGGCCGTACCTGGAGCAGGGTTACGACCACCTGCTACTGa                                                                                                          >  2:2532957/1‑150 (MQ=255)
                                             tCGAGCAGGCAAAACAGACGGTCAGCGAGCACCAGTTACCGTTGCAGATGCGCGTGCTGCAGCCGTTCTATGAGCACCCGGCGTACATCGAGGCGCTGGCAGCCAGCGCCCGGCCGTACCTGGAGCAGGGTTACGACCACCTGCTACTGa                                                                                                          >  2:4297218/1‑150 (MQ=255)
                                             ccGAGCTGGCAAAACAGCCGGTCCGCGACCACCCGTTCCCGTCGCAGAGGCGCGTGCTGCAGCCTTTCTATGAGCCCCCGGCGTACACCGAGGCGGTGGCAGCCCGCGCCCGGCCGTACCTGGAGCAGGGTTACGACCACCTGCTACTGa                                                                                                          <  2:4562222/149‑1 (MQ=255)
                                                     gCAAAACAGACGGTCAGCGAGCACCAGTTACCGTTGCAGATGCGCGTGCTGCAGCCGTTCTATGAGCACCCGGCGTACATCGAGGCGCTGGCAGCCAGCGCCCGGCCGTACCTGGAGCAGGGTTACGACCACCTGCAACTGAGCTTCCAt                                                                                                  >  1:99390/1‑150 (MQ=255)
                                                     gCAAAACAGACGGTCAGCGAGCACCAGTTACCGTTGCAGATGCGCGTGCTGCAGCCGTTCTATGAGCACCCGGCGTACATCGAGGCGCTGGCAGCCAGCGCCCGGCCGGACCTGGAGCAGGGTTACGACCACCGGCGACTGAGCTTCCAt                                                                                                  >  1:10176/1‑150 (MQ=255)
                                                         aaccgccGGTCCGCTAGCACCAGTTCCCTTTGCAGTGGCGCGTGCTGCAGCCTTTATATGAGCCCCCGGCGTACATCGCGGCGCTGGCAGCCAGCGCCCGGCCGTACCTGGAGCAGGGTTACGACCACCTGCTACTGAGCTTCCATGGtt                                                                                              <  2:1375887/144‑1 (MQ=255)
                                                              aCGGTCAGCGAGCACCAGTTACCGTTGCAGATGCGCGTGCTGCAGCCGTTCTATGAGCACCCGGCGTACATCGAGGCGCTGGCAGCCAGCGCCCGGCCGTACCTGGAGCAGGGTTACGACCACCTGCTACTGAGCTTCCATGGTTTGCCg                                                                                         >  2:7522706/1‑150 (MQ=255)
                                                                   ctgCGACCACCAGTTACCGTTGCAGATGCGCGTGCTGCCGCCTTTCTATGAGCACCCGGCGTACATCGAGGCGCTGGCAGCCAGCGCCCGGCCGTACCTGGAGCAGGGTTACGACCACCTGCTACTGAGCTTCCATGGTTTGCCGGAACg                                                                                    <  2:5026732/148‑1 (MQ=255)
                                                                       gAGCACCTGTTACCGTTGCAGTGGCGCGGGCTGCAGCCGTTCTATGAGCCCCCGGCGTACATCGCGGCGCTGGCACCCCGCGCCCGGCCGTACCTGGAGCAGGGTTACGACCACCTGCTACTGAGCTTCCATGGTTTGCCGGAACGGCAc                                                                                <  2:542758/150‑1 (MQ=255)
                                                                       gAGCACCAGTTACCGTTGCAGATGCGCGTGCTGCAGCCGTTCTATGAGCACCCGGCGTACATCGAGGCGCTGGCAGCCAGTGCCCGGCCGTACCTGGAGCAGGGTTACGACCACCTGCTACTGAGCTTCCATGGTTTGCCGGAACGGCAc                                                                                <  2:3504593/150‑1 (MQ=255)
                                                                          cACCAGTTACCGTTGCAGATGCGCGTGCTGCAGCCGTTCTATGAGCACCCGGCGTACATCGAGGCGCTGGCAGCCAGCGCCCGGCCGTACCTGGAGCAGGGTTACGACCACCTGCTACTGAGCTTCCATGGTTTGCCGGAACGGCACCTg                                                                             >  1:3785951/1‑150 (MQ=255)
                                                                            cctttttccGTTGCTGTTGCGCGTGGTGCCGCCTTTCTATGCGCACCGGGCTTACATCGCGGCGCTGGCCCCCCCCCCCCGGCCGTACCTGGAGCAGGGTTACGACCACCTGCTACTGAGCTTCCATGGTTTGCCGGAACGGCACCTgaa                                                                           <  1:3416244/143‑1 (MQ=255)
                                                                             cAGTTACCGTTGCAGCGGCGCGTGCTGCCGCCGTTCTATCAGCACCCGGCGTCCAGCGCGGCGCGGCCCGCCCGCGCCCGGCCGTACCTGGAGCAGGGTTACGACCACCTGCTACTGAGCTTCCATGGTTTGCCGGAACGGCACCTgaag                                                                          <  1:4302532/150‑1 (MQ=255)
                                                                             cAGTTACCGTTGCAGATGCGCGTGCTGCAGCCGTTCTATGAGCACCCGGCGTACATCGAGGCGCTGGCAGCCAGCGCCCGGCCGTACCTGGAGCAGGGTTACGACCACCTGCTACTGAGCTTCCATGGTTTGCCGGAACGGCACCTgaag                                                                          >  2:8355300/1‑150 (MQ=255)
                                                                             cAGTTACCGTTGCAGATGCGCGTGCTGCAGCCGTTCTATGAGCACCCGGCGTACATCGAGGCGCTGGCAGCCAGCGCCCGGCCGTACCTGGAGCAGGGTTACGACCACCTGCTACTGAGCTTCCATGGTTTGCCGGAACGGCACCTgaag                                                                          >  2:699006/1‑150 (MQ=255)
                                                                               gTTACCGTTGCAGATGCGCGTGCTGCAGCCGTTCTATGAGCACCCGGCGTACATCGAGGCGCTGGCAGCCAGCGCCCGGCCGTACCTGGAGCAGGGTTACGACCACCTGCTACTGAGCTTCCATGGTTTGCCGGAACGGCACCTgaagaa                                                                        >  1:7274355/1‑150 (MQ=255)
                                                                                tcacggTTGCTGAGGGGCGTGCGGCCGCGGGTCTATGCGCCCCCGGCTAACCTGGCGGCGCGGGCCGCCCGCGCCCGGCCGTACCTGGAGCAGGGTTACGACCACCTGCTACTGAGCTTCCATGGTTTGCCGGAACGGCACCTGAAGaaa                                                                       <  2:6089523/145‑1 (MQ=255)
                                                                                  cccGTTGCTGCTGCGCGTGCTGCCGCCTTTCTAGTCGCACCCGGCGTACAGCGTGGCGCGGGCCCCCCGCGCCCGGCCGTACCTGGAGCAGGGTTACGACCACCTGCTACTGAGCTTCCATGGTTTGCCGGAACGGCACCTGAAGAAACt                                                                     <  2:4392697/149‑1 (MQ=255)
                                                                                  aCCGTTGCAGATGCGCGTGCTGCAGCCGTTCTATGAGCACCCGGCGTACATCGAGGCGCTGGCAGCCAGCGCCCGGCCGTACCTGGAGCAGGGTTACGACCACCAGCTACTGAGCTTCCATGGTTTGCCGGAACGGCACCTGAAGAAACt                                                                     >  1:4426358/1‑150 (MQ=255)
                                                                                   ccGTTGCAGATGCGCGTGCTGCAGCCGTTCTATGAGCACCCGGCGTACATCGAGGCGCTGGCAGCCAGCGCCCGGCCGTACCTGGAGCAGGGTTACGACCACCTGCTACTGAGCTTCCATGGTTTGCCGGAACGGCACCTGAAGAAACTg                                                                    >  1:1510134/1‑150 (MQ=255)
                                                                                    cGTTGCAGCTGCGGTGGCGGCCGCCGTTATTTGCGCACCCGGCTTACAGCGCGGCGCTGCCCGCCCGCGCCCGGCCGTACCTGGAGCAGGGTTACGACCACCTGCTACTGAGCTTCCATGGTTTGCCGGAACGGCACCTGAAGAAACTGt                                                                   <  2:1458511/150‑1 (MQ=255)
                                                                                    cGTTGCAGATGCGCGTGCTGCAGCCGTTCTATGAGCACCCGGCGTACATCGAGGCGCTGGCAGCCAGCGCCCGGCCGTACCTGGAGCAGGGTTACGACCACCTGCTACTGAGCTTCCATGGTTTGCCGGAACGGCACCTGAAGAAACTGt                                                                   >  1:456029/1‑150 (MQ=255)
                                                                                     gTTGCAGTTGGGGGTGCTGCAGCCGTTCTTTGAGCACCGGGCTTACATCGCGGCGCTGGCAGCCAGCGCCCGGCCGTACCTGGAGCAGGGTTACGACCACCTGCTACTGAGCTTCCATGGTTTGCCGGAACGGCACCTGAAGAAACTGtt                                                                  <  2:5756052/150‑1 (MQ=255)
                                                                                       tGCAGTTGGGCGTGCAGCAGCCGTTCTATGAGCACCCGGCGTACATCGTGGCGCTGGCAGCCAGCGCCCGGCCGTACCTGGAGCAGGGTTACGACCACCTGCTACTGAGCTTCCATGGTTTGCCGGAACGGCACCTGAAGAAACTGTTTc                                                                <  1:8056771/150‑1 (MQ=255)
                                                                                          aGATGCGCGTGCTGCAGCCGTTCTATGAGCACCCGGCGTACATCGAGGCGCTGGCAGCCAGCGCCCGGCCGTACCTGGAGCAGGGTTACGACCACCTGCTACTGAGCTTCCATGGTTTGCCGGAACGGCACCTGAAGAAACTGTTTCCGa                                                             >  2:7708429/1‑150 (MQ=255)
                                                                                               cgcgTGCGGCAGCCGTTCTTCGAGCACCCGGCGAACAGCGCGGCGCCGGCAGCCCGCCCCCGGCCGTACCTGGAGCAGGGTTACGACCACCTGCTACTGAGCTTCCATGGTTTGCCGGAACGGCACCTGAAGAAACTGTTTCCGAAgggg                                                        <  1:5365741/150‑1 (MQ=255)
                                                                                                 cgTGGTGCCGCCTTTCTATGAGCACCGGGCTTTATTCGCGGGTGTGGCTGCCAGCGCCCGGCCGTACCTGGAGCAGGGTTACGACCACCTGCTACTGAGCTTCCATTGTTTGCCGGAACGGCACCTGAAGAAACTGTTTCCGAAGGGGGt                                                      <  2:5093330/150‑1 (MQ=255)
                                                                                                      tgcAGCCGTTCTATGAGCACCCGGCGTACATCGAGGCGCTGGCAGCCAGCGCCCGGCCGTACCTGGAGCAGGGTTACGACCACCTGCTACTGAGCTTCCATGGTTTGCCGGAACGGCACCTGAAGAAACTGTTTCCGAAGGGGGTCAAGc                                                 >  2:3059379/1‑150 (MQ=255)
                                                                                                         aGCCGTTCTATGAGCACCCGGCGTACATCGAGGCGCTGGCAGCCAGCGCCCGGCCGTACCTGGAGCCGGGTTACGACCACCTGCTACTGAGCTTCCATGGTTTGCCGGAACGGCACCTGAAGAAACTGTTTCCGAAGGGGGTCAAGCACg                                              >  2:1721841/1‑150 (MQ=255)
                                                                                                           cgtttCTTTGAGCCCCGGGCGTACATCGCGGCGCTGGCAGCCCGCGCCCGGCCGTACCTGGAGCAGGGTTACGACCACCTGCTACTGAGCTTCCATGGTTTGCCGGAACGGCACCTGAAGAAACTGTTTCCGAAGGGGGTCAAGCACGAc                                            <  2:3441056/147‑1 (MQ=255)
                                                                                                            cGTTCTATGAGCACCCGGCGTACATCGAGGCGCTGGCAGCCAGCGCCCGGCCGTACCTGGAGCAGGGTTACGACCACCTGCTACTGAGCTTCCATGGTTTGCCGGAACGGCACCTGAAGAAACTGTTTCCGAAGGGGGTCAAGCACGAcc                                           >  2:5547762/1‑150 (MQ=255)
                                                                                                             gTTCTATGAGCACCCGGCGTACATCGAGGCGCTGGCAGCCCGCGCCCGGCCGTACCTGGAGCAGGGTTACGACCACCTGCTACTGAGCTTCCATGGTTTGCCGGAACGGCACCAGAAGAAACTGTTTCCGAAGGGGGTCAAGCACGACCt                                          >  1:8106340/1‑150 (MQ=255)
                                                                                                                 tATGAGCACCCGGCGTACATCGAGGCGCTGGCAGCCAGCGCCCGGCCGTACCTGGAGCAGGGTTACGACCACCTGCTACTGAGCTTCCATGGTTTGCCGGAACGGCACCTGAAGAAACTGTTTCCGAAGGGGGTCAAGCACGACCTgcgc                                      >  2:4084461/1‑150 (MQ=255)
                                                                                                                 tATGAGCACCCGGCGTACATCGAGGCGCTGGCAGCCAGCGCCCGGCCGTACCTGGAGCAGGGTTACGACCACCTGCTACTGAGCTTCCATGGTTTGCCGGAACGGCACCTGAAGAAACTGTTTCCGAAGGGGGTCAAGCACGACCTgcgc                                      >  1:555454/1‑150 (MQ=255)
                                                                                                                   tGAGCACCCGGCGTACATCGAGGCGCTGGCAGCCAGCGCCCGGCCGTACCTGGAGCAGGGTTACGACCACCTGCTACTGAGCTTCCATGGTTTGCCGGAACGGCACCTGAAGAAACTGTTTCCGAAGGGGGTCAAGCACGACCTgcgagc                                    >  2:6458975/1‑150 (MQ=255)
                                                                                                                    gAGCACCCGGCGTACATCGAGGCGCTGGCAGCCAGCGCCCGGCCGTACCTGGAGCAGGGTTACGACCACCTGCTACTGAGCTTCCATGGTTTGCCGGAACGGCACCTGAAGAAACTGTTTCCGAAGGGGGTCAAGCACGACCTGCGCGCt                                   >  1:7737699/1‑150 (MQ=255)
                                                                                                                        aCCCGGCGTACATCGAGGCGCTGGCAGCCAGCGCCCGGCCGTACCTGGAGCAGGGTTACGACCACCTGCTACTGAGCTTCCATGGTTTGCCGGAACGGCACCTGAAGAAACTGTTTCCGAAGGGGGTCAAGCACGACCTGCGCGCTGCCg                               >  1:3603968/1‑150 (MQ=255)
                                                                                                                          ccGGCGTACCCGGGGGCGCGGGCCGCCCGCGCCCGGCCGTACCTGGAGCAGGGTTACGACCACCTGCTACTGAGCTTCCATGGTTTGCCGGAACGGCACCTGAAGAAACTGTTTCCGAAGGGGGTCAAGCACGACCTGCGCGCTGCCGAt                             <  1:8737458/150‑1 (MQ=255)
                                                                                                                           cGGCGTACATCGCGGCGCGGGCAGCCAGCGCCCGGCCGTACCTGGAGCAGGGTTACGACCACCTGCTACTGAGCTTCCATGGTTTGCCGGAACGGCACCTGAAGAAACTGTTTCCGAAGGGGGTCAAGCACGACCTGCGCGCTGCCGAtt                            <  2:2293146/150‑1 (MQ=255)
                                                                                                                           cGGCGTACATCGAGGCGCTGGCAGCCAGCGCCCGGCCGTACCTGGCGCAGGGTTACGACCACCTGCTACTGAGCTTCCATGGTTTGCCGGAACGGCACCTGCAGAAACTGTTTCCGAAGGGGGTCAAGCACGACCTGCGCGCTGCCGAtt                            >  1:5889508/1‑150 (MQ=255)
                                                                                                                            ggCTTACAGCGTGGCGCTGGCAGCCAGCGCCCGGCCGTACCTGGAGCAGGGTTACGACCACCTGCTACTGAGCTTCCATGGTTTGCCGGAACGGCACCTGAAGAAACTGTTTCCGAAGGGGGTCAAGCACGACCTGCGCGCTGCCGATtg                           <  2:6628484/150‑1 (MQ=255)
                                                                                                                               gTACATCGAGGCGCTGGCAGCCAGCGCCCGGCCGTACCTGGAGCAGGGTTACGACCACCTGCTACTGAGCTTCCATGGTTTGCCGGAACGGCACCTGAAGAAACTGTTTCCGAAGGGGGTCAAGCACGACCTGCGCGCTGCCGATtgctg                        >  1:6503048/1‑150 (MQ=255)
                                                                                                                               gTACATCGAGGCGCTGGCAGCCAGCGCCCGGCCGTACCTGGAGCAGGGTTACGACCACCTGCTACTGAGCTTCCATGGTTTGCCGGAACGGCACCTGAAGAAACTGTTTCCGAAGGGGGTCAAGCACGACCTGCGCGCTGCCGATtgctg                        >  1:2099254/1‑150 (MQ=255)
                                                                                                                               gTACATCGAGGCGCTGGCAGCCAGCGCCCGGCCGTACCTGGAGCAGGGTTACGACCACCTGCTACTGAGCTTCCATGGTTTGCCGGAACGGCACCGGAAGAAACTGTTTCCGAAGGGGGTCAAGCACGGCCTGCGCGCTGCGGAGtgctg                        >  2:8410122/1‑150 (MQ=255)
                                                                                                                                   aTCGAGGCGCTGGCAGCCAGCGCCCGGCCGTACCTGGAGCAGGGTTACGACCACCTGCTACTGAGCTTCCATGGTTTGCCGGAACGGCACCTGACGAAACTGTTTCCGAAGGGGGTCAAGCACGACCTGCGCGCTGCCGATTGCTGCCAc                    >  2:82014/1‑150 (MQ=255)
                                                                                                                                     cgcggCGCTGGCAGCCAGCGCCCGGCCGTACCTGGAGCAGGGTTACGACCACCTGCTACTGAGCTTCCATGGTTTGCCGGAACGGCACCTGAAGAAACTGTTTCCGAAGGGGGTCAAGCACGACCTGCGCGCTGCCGATTGCTGCCACgg                  <  2:7308651/147‑1 (MQ=255)
                                                                                                                                      gcggCGGTGCCCGCCCGCCCCCGGCCGTACCTGGAGCAGGGTTACGACCACCTGCTACTGAGCTTCCATGGTTTGCCGGAACGGCACCTGAAGAAACTGTTTCCGAAGGGGGTCAAGCACGACCTGCGCGCTGCCGATTGCTGCCACGGt                 <  2:365126/148‑1 (MQ=255)
                                                                                                                                      gAGGCGCTGGCAGCCAGCGCCCGGCCGTACCTGGAGCAGGGTTACGACCACCTGCTACTGAGCTTCCATGGTTTGCCGGAACGGCACCTGACGAAACTGTTTCCGAAGGGGGTCAAGCACGACCTGCGCGCTGCCGATTGCTGCCACGGt                 >  2:4870137/1‑150 (MQ=255)
                                                                                                                                      gAGGCGCTGGCAGCCAGCGCCCGGCCGTACCTGGAGCAGGGTTACGACCACCTGCTACTGAGCTTCCATGGTTTGCCGGAACGGCACCTGAAGAAACTGTTTCCGAAGGGGGTCAAGCACGACCTGCGCGCTGCCGATTGCTGCCACGGt                 >  1:6032533/1‑150 (MQ=255)
                                                                                                                                      gAGGCGCTGGCAGCCAGCGCCCGGCCGTACCTGGAGCAGGGTTACGACCACCTGCTACTGAGCTTCCATGGTTTGCCGGAACGGCACCTGAAGAAACTGTTTCCGAAGGGGGTCAAGCACGACCTGCGCGCTGCCGATTGCTGCCACGGt                 >  1:843462/1‑150 (MQ=255)
                                                                                                                                       gggCGCTGGCTGCCAGCGCCCGGCCGTACCTGGAGCAGGGTTACGACCACCTGCTACTGAGCTTCCATGGTTTGCCGGAACGGCACCTGAAGAAACTGTTTCCGAAGGGGGTCAAGCACGACCTGCGCGCTGCCGATTGCTGCCACGGTg                >  2:8057900/2‑150 (MQ=255)
                                                                                                                                       gggCGCTGGCTGCCAGCGCCCGGCCGTACCTGGAGCAGGGTTACGACCACCTGCTACTGAGCTTCCATGGTTTGCCGGAACGGCACCTGAAGAAACTGTTTCCGAAGGGGGTCAAGCACGACCTGCGCGCTGCCGATTGCTGCCACGGTg                >  1:2758477/2‑150 (MQ=255)
                                                                                                                                       aGGCGCTGGCAGCCAGCGCCCGGCCGTACCTGGAGCAGGGTTACGACCACCTGCTACTGAGCTTCCATGGTTTGCCGGAACGGCACCTGAAGAAACTGTTTCCGAAGGGGGTCAAGCACGACCTGCGCGCTGCCGATTGCTGCCACGGTg                >  1:6915931/1‑150 (MQ=255)
                                                                                                                                       aGGCGCTGGCAGCCAGCGCCCGGCCGTACCTGGAGCAGGGTTACGACCACCTGCTACTGAGCTTCCATGGTTTGCCGGAACGGCACCTGAAGAAACTGTTTCCGAAGGGGGTCAAGCACGACCTGCGCGCTGCCGATTGCTGCCACGGTg                >  1:1283805/1‑150 (MQ=255)
                                                                                                                                        ggCGCTGGCAGCCAGCGCCCGGCCGTACCTGGAGCAGGGTTACGACCACCTGCTACTGAGCTTCCATGGTTTGCCGGAACGGCACCTGAAGAAACTGTTTCCGAAGGGGGTCAAGCACGACCTGCGCGCTGCCGATTGCTGCCACGGTGc               >  2:3346930/1‑150 (MQ=255)
                                                                                                                                           gctggctgCCAGCGCCTGGCCGTACCTGGAGCAGGGTTACGACCACCTGCTACTGAGCTTCCATGGTTTGCCGGAACGGCACCTGAAGAAACTGTTTCCGAAGGGGGTCAAGCACGACCTGCGCGCTGCCGATTGCTGCCATGGTGCAAc            <  1:4068774/150‑1 (MQ=255)
                                                                                                                                             tGGCAGCCAGCGCCCGGCCGTACCTGGAGCAGGGTTACGACCACCTGCTACTGAGCTTCCATGGTTTGCCGGAACGGCACCTGAAGAAACTGTTTCCGAAGGGGGTCAAGCACGACCTGCGCGCTGCCGATTGCTGCCACGGTGCAAccg          <  2:7834770/150‑1 (MQ=255)
                                                                                                                                               gCAGCCAGCGCCCGGCCGTACCTGGAGCAGGGTTACGACCACCTGCTACTGAGCTTCCATGGTTTGCCGGAACGGCACCTGACGAACCTGTTTCCGAAGGGGGTCAAGCACGACCTGCGCGCTGCCGAGTGCTGCCACGGTGCACccgcc        >  1:2457438/1‑150 (MQ=255)
                                                                                                                                                  gcccgcGCCCGGCCGTACCTGGAGCAGGGTTACGACCACCTGCTACTGAGCTTCCATGGTTTGCCGGAACGGCACCTGAAGAAACTGTTTCCGAAGGGGGTCAAGCACGACCTGCGCGCTGCCGATTGCTGCCACGGTGCAACCGCCGaa     <  2:6181427/150‑1 (MQ=255)
                                                                                                                                                  gccagcGCCCGGCCGTACCTGGAGCAGGGTTACGACCACCTGCTACTGAGCTTCCATGGTTTGCCGGAACGGCACCTGAAGAAACTGTTTCCGAAGGGGGTCAAGCACGACCTGCGCGCTGCCGATTGCTGCCACGGTGCAACCGCCGaa     <  2:4455611/150‑1 (MQ=255)
                                                                                                                                                    cagcGCCCGGCCGTACCTGGAGCAGGGTTACGACCACCTGCTACTGAGCTTCCATGGTTTGCCGGAAAGGCACCTGAAGAAACTGTTTCCGAAGGGGGTCAAGAACGACCTGCGCGCTGCCGATTGCTGCCACGGTGCAACCGCCGaaga   >  2:3960840/1‑149 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
CACCGCTTTATCCACAGTTCGCCGACAGCACGGTGACCACTGTGATCGAGCAGGCAAAACAGACGGTCAGCGAGCACCAGTTACCGTTGCAGATGCGCGTGCTGCAGCCGTTCTATGAGCACCCGGCGTACATCGAGGCGCTGGCAGCCAGCGCCCGGCCGTACCTGGAGCAGGGTTACGACCACCTGCTACTGAGCTTCCATGGTTTGCCGGAACGGCACCTGAAGAAACTGTTTCCGAAGGGGGTCAAGCACGACCTGCGCGCTGCCGATTGCTGCCACGGTGCAACCGCCGAAGTA  >  NC_002947/863234‑863532

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 11 ≤ ATCG/ATCG < 25 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑