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

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_002947863,3830AC19.9% 671.9 / 31.0 303S172R (AGC→CGC) hemHferrochelatase
Reads supporting (aligned to +/- strand):  ref base A (188/45);  new base C (4/64);  total (193/110)
Fisher's exact test for biased strand distribution p-value = 1.38e-30
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 338 total aligned reads displayed.

CACCGCTTTATCCACAGTTCGCCGACAGCACGGTGACCACTGTGATCGAGCAGGCAAAACAGACGGTCAGCGAGCACCAGTTACCGTTGCAGATGCGCGTGCTGCAGCCGTTCTATGAGCACCCGGCGTACATCGAGGCGCTGGCAGCCAGCGCCCGGCCGTACCTGGAGCAGGGTTACGACCACCTGCTACTGAGCTTCCATGGTTTGCCGGAACGGCACCTGAAGAAACTGTTTCCGAAGGGGGTCAAGCACGACCTGCGCGCTGCCGATTGCTGCCACGGTGCAACCGCCGAAGTA  >  NC_002947/863234‑863532
                                                                                                                                                     |                                                                                                                                                     
cACCGCTTTATCCACAGTTCGCCGACAGCACGGTGACCACTGTGATCGAGCAGGCAAAACAGACGGTCAGCGAGCACCAGTTACCGTTGCAGATGCGCGTGCTGCAGCCGTTCTATGAGCACCCGGCGTACATCGAGGCGCTGGcagcca                                                                                                                                                       <  2:1240163/150‑1 (MQ=255)
cACCGCTTTATCCACAGTTCGCCGACAGCACGGTGACCACTGTGATCGAGCAGGCAAAACAGACGGTCAGCGAGCACCAGTTACCGTTGCAGATGCGCGTGCTGCAGCCGTTCTATGAGCACCCGGCGTACATCGAGGCGCTGGcagcca                                                                                                                                                       >  1:5122955/1‑150 (MQ=255)
 aCCGCTTTATCCACAGTTCGCCGACAGCACGGTGACCACTGTGATCGAGCAGGCAAAACAGACGGTCAGCGAGCACCAGTTACCGTTGCAGATGCGCGTGCTGCAGCCGTTCTATGAGCACCCGGCGTACATCGAGGCGCTGGcagccag                                                                                                                                                      >  1:6072930/1‑150 (MQ=255)
 aCCGCTTTATCCACAGTTCGCCGACAGCACGGTGACCACTGTGATCGAGCAGGCAAAACAGACGGTCAGCGAGCACCAGTTACCGTTGCAGATGCGCGTGCTGCAGCCGTTCTATGAGCACCCGGCGTACATCGAGGCGCTGGcagccag                                                                                                                                                      >  1:574055/1‑150 (MQ=255)
 aCCGCTTTATACACAGTTCGCCGACAGCACGGTGACCACTGTGATCGAGAAGGCAAAACAGACGGTCAGCGAGCACCAGTTACCGTTGCAGATGCGCGTGCTGCAGCCGTTCTATGAGCACCCGGCGTACATCGAGGCGCTGGcagccag                                                                                                                                                      >  2:6855138/1‑150 (MQ=255)
    gCTTTATCCACAGTTCGCCGACAGCACGGTGACCACTGTGATCGAGCAGGCAAAACAGACGGTCAGCGAGCACCAGTTACCGTTGCAGATGCGCGTGCTGCAGCCGTTCTATGAGCACCCGGCGTACATCGAGGCGCTGGCCGACAgcgc                                                                                                                                                   >  1:2830108/1‑150 (MQ=255)
      tttATCCACAGTTCGCCGACAGCACGGTGACCACTGTGATCGAGCAGGCAAAACAGACGGTCAGCGAGCACCAGTTACCGTTGCAGATGCGCGTGCTGCAGCCGTTCTATGAGCACCCGGCGTACATCGGGGCGCTGGCTGCCAGCGcct                                                                                                                                                 <  2:4994541/150‑2 (MQ=255)
      tttATCCACAGTTCGCCGACAGCACGGTGACCACTGTGATCGAGCAGGCAAAACAGACGGTCAGCGAGCACCAGTTACCGTTGCAGATGCGCGTGCTGCAGCCGTTCTATGAGCACCCGGCGTACATCGAGGCGCTGGCTGCCAGCGcct                                                                                                                                                 >  1:30443/1‑149 (MQ=255)
      tttATCCACAGTTCGCCGACAGCACGGTGACCACTGTGATCGAGCAGGCAAAACAGACGGTCAGCGAGCACCAGTTACCGTTGCAGATGCGCGTGCTGCAGCCGTTCTATGAGCACCCGGCGTACATCGAGGCGCTGGCAGCCAGCGccc                                                                                                                                                 >  2:6448146/1‑150 (MQ=255)
        taaccACAGTTCGCCGACAGCACGGTGACCACTGTGATCGAGCAGGCAAAACAGACGGTCAGCGAGCACCAGTTACCGTTGCAGATGCGCGTGCTGCAGCCGTTCTATGAGCACCCGGCGTACATCG‑GGCGCTGGCTGCCAGCGcctcga                                                                                                                                              <  1:6242463/147‑5 (MQ=255)
         aTCCACAGTTCGCCGACCGCCGGGGGACCCCTGTGATCGAGCAGGCAAAACAGACGGTCAGCGAGCACCCGTTACCGTTGCTGTTGCGCGTGGTGCAGCCGTTCTATGAGCACCGGGCGTACATCGGGGCGCTGGCAGCCAGCGCCCGGc                                                                                                                                              <  2:5400716/150‑1 (MQ=255)
          tCCACAGTTCGCCGACAGCACGGTGACCACTGTGATCGAGCAGGCAAAACAGACGGTCAGCGAGCACCAGTTACCGTTGCAGATGCGCGTGCTGCAGCCGTTCTATGAGCACCCGGCGTACATCGAGGCGCTGGCAGCCAGCGCCCGGcc                                                                                                                                             >  1:2329824/1‑150 (MQ=255)
          tCCACAGTTCGCCGACAGCACGGTGACCACTGTGATCGAGCAGGCAAAACAGACGGTCAGCGAGCACCAGTTACCGTTGCAGATGCGCGTGCTGCAGCCGTTCTATGAGCACCCGGCGTACATCGAGGCGCTGGCAGCCAGCGCCCGGcc                                                                                                                                             >  2:2912394/1‑150 (MQ=255)
          tCCACAGTTCGCCGACAGCACGGTGACCACTGTGATCGAGCAGGCAAAACAGACGGTCAGCGAGCACCAGTTACCGTTGCAGATGCGCGTGCTGCAGCCGTTCTATGAGCACCCGGCGTACATCGAGGCGCTGGCAGCCAGCGCCCGGcc                                                                                                                                             >  2:2087832/1‑150 (MQ=255)
            cacaGTTCGCCGACAGCACGGTGACCACTGTGATCGAGCAGGCAAAACAGACGGTCAGCGAGCACCAGTTACCGTTGCAGATGCGCGTGCTGCAGCCGTGCTATGAGCACCCGGCGTACATCGAGGCGCTGGCAGCCAGCGCCCGGCCGt                                                                                                                                           >  1:3784042/1‑150 (MQ=255)
                gTTCGCCGACAGCACGGTGACCACTGTGATCTAGCAGGAAAAACAGACGGTCAGCTATCACCAGTTACCGTTGCAGATGCGCGTGCTGCAGCCTTTATATGAGCACCCGGCTTACATCGAGGCGCTGGCACCCACCGCCCGGCCGTACCt                                                                                                                                       <  2:1679377/150‑1 (MQ=255)
                gTTCGCCGACAGCACGGTGACCACTGTGATCGCGCAGGCAAAACAGACGGTCAGCGAGCACCAGTTACGGTTGCAGATGCGCGTGCTGCAGCCGTTCTTTGAGCACCCGGCTTACATCGAGGCGCTGGCCGCCAGCGCCCGGCCGTACCt                                                                                                                                       <  2:4823414/150‑1 (MQ=255)
                gTTCGCCGACAGCACGGTGACCACTGTGATCGAGCAGGCAAAACAGACGGTCAGCGAGCACCAGTTACCGTTGCAGATGCGCGTGCTGCAGCCGTTCTATGAGCACCCGGCGTACATCGAGGCGCTGGCAGCCAGCGCCCGGCCGTACCt                                                                                                                                       >  2:3057423/1‑150 (MQ=255)
                   cgccgcCCGCCCGGTCACCCCTGTGCTCGCGCAGCAAAAACCGCCGGCCCGCGAGCACCTGTGCCGGTTGCCGCGGCGGGGGCTGCCGCCTTTGTATGAGCCCCCGGCGCACAGCGCGGCGGTGGCAGCCCGCGCCCGGCCGTACCTGGa                                                                                                                                    <  2:765469/150‑1 (MQ=255)
                   cGCCGACAGCCGGGTCACCACTTTTATCGAGCAGGCAAAACAGACGGTCAGCGAGCACCAGTTACGTTTGCTTCTGCGGGTGCTGCCGCCTTTATTTGAGCACCCGGCTTACATGGTGGCGCTGGCAGCCAGCGCCCGGCCGTACCTGGa                                                                                                                                    <  2:1490311/150‑1 (MQ=255)
                        aCAGCACGGTGACCACTGTGATCGAGCAGGCAAAACAGACGGTCAGCGAGCACCAGTTACCGTTGCAGATGCGCGTGCTGCAGCCGTTCTATGAGCACCCGGCGTACATCGAGGCGCTGGCAGCCAGCGCCCGGCCGTACCTGGAGCAgg                                                                                                                               >  2:60223/1‑150 (MQ=255)
                            cACGGTGACCACTGTGATCGAGCAGGCAAAACAGACGGTCAGCGAGCACCAGTTACCGTTGCAGATGCGCGTGCTGCAGCCGTTCTATGAGCACCCGGCGTACATCGAGGCGCTGGCAGCCAGCGCCCGGCCGTACCTGGAGCAGGGttc                                                                                                                           >  2:4212104/1‑149 (MQ=255)
                                      aCTGTGATCGAGCAGGCAAAACAGACGGTCAGCGAGCACCAGTTACCGTTGCAGATGCGCGTGCTGCAGCCGTTCTATGAGCACCCGGCGTACATCGAGGCGCTGGCAGCCAGCGCCCGGCCGTACCTGGAGCAGGGTTACGACCACCTg                                                                                                                 >  2:6531792/1‑150 (MQ=255)
                                        tgtgATCGAGCAGGCAAAACAGACGGTCAGCGAGCACCAGTTACCGTTGCAGATGCGCGTGCTGCAGCCGTTCTATGAGCACCCGGCGTACATCGAGGCGCTGGcagcca                                                                                                                                                       >  2:5584955/1‑110 (MQ=255)
                                        tgtgATCGAGCAGGCAAAACAGACGGTCAGCGAGCACCAGTTACCGTTGCAGATGCGCGTGCTGCAGCCGTTCTATGAGCACCCGGCGTACATCGAGGCGCTGGCAGCCAGCGCCCGGCCGTACCTGGAGCAGGGTTACGACCACCTGCt                                                                                                               >  1:1942731/1‑150 (MQ=255)
                                        tgtgATCGAGCAGGCAAAACAGACGGTCAGCGAGCACCAGTTACCGTTGCAGATGCGCGTGCTGAAGCCGTTCTATGAGCACCCGGCGTACATCGAGTCGCTGGCAGCCAGCCCCCGGCCGTACCTGGAGCAGGGTTACGACCACCTGCt                                                                                                               >  1:2643219/1‑150 (MQ=255)
                                             tCGAGCAGGCAAAACAGACGGTCAGCGAGCACCAGTTACCGTTGCAGATGCGCGTGCTGCAGCCGTTCTATGAGCACCCGGCGTACATCGAGGCGGTGGCAGCCAGCGCCCGGCCGTACCTGGAGCAGGGTTACGACCACCTGCTACTGa                                                                                                          >  1:2208154/1‑150 (MQ=255)
                                             tCGAGCACCCAAAACCGCCGGTCCGCGACCACCAGTCACCGGTGCCGAGGCGCGTGCGGCCGCCGTTCTAGGAGCCCCCGGCGTACCCCGCGGCGCGCGCCCCCCGCCCCCGGCCGTACCTGGAGCAGGGTTACGACCACCTGCTACTGa                                                                                                          <  2:3642370/150‑1 (MQ=255)
                                               gAGCAGGCAAAACAGACGGTCAGCGAGCACCAGTTACCGTTGCAGATGCGCGTGCTGCAGCCGTTCTATGAGCACCCGGCGTACATCGAGGCGCTGGCAGCCAGCGCCCGGCCGTACCTGGAGCAGGGTTACGACCACCTGCTACTGAGc                                                                                                        >  2:743900/1‑150 (MQ=255)
                                                 tccggCAAAACAGTCGGTCCGCTAGCACCTGTTACGGTGGCAGAGGCGGGTGGTGCCGCCGTTCTATGACCCCCCGGCTTACTTCGCGGCGGGGGCACCCCGCCCCCGGCCGTACCTGGAGCAGGGTTACGACCACCTGCTACTGAGCtt                                                                                                      <  2:5057823/147‑1 (MQ=255)
                                                 gCAGGCAAAACAGAGGCACCGCGCCCACCCGTTCCCGTGGCAGCGGCGCGTGCTGCCGCCGTTCTCTGAGCCCCCGGCGTACAGCGCGGCGGGGCCCGCCCGCCCCCGGCCGTACCTGGAGCAGGGTTACGACCACCTGCTACTGAGCtt                                                                                                      <  2:4414011/150‑1 (MQ=255)
                                                 gCAGGCAAAACAGACGGTCAGCGAGCACCAGTTACCGTTGCAGATGCGCGTGCTGCAGCCGTTCTATGAGCACCCGGCGTACATCGAGGCGCTGGCAGCCCGCGCCCGGCCGTACCTGGAGCCGGGGTACGACCACCTGCTACTGAGCtt                                                                                                      >  1:2490397/1‑150 (MQ=255)
                                                    ggCAAAACAGACGGTCAGCGAGCACCAGTTACCGTTGCAGATGCGCGTGCTGCCGCCGTTCTATGCGCCCCCGGCGTACATCGAGGCGCTGGCCGCCCGCGCCCGGCCGTACCTGGAGCAGGGTTACGACCACCTGCTACTGAGCTTCCa                                                                                                   <  2:6698503/150‑1 (MQ=255)
                                                     gCAAAACAGACGGTCAGCGAGCACCAGTTACCGTTGCAGATGCGCGTGCTGCAGCCGTTCTATGAGCACCCGGCGTACATCGAGGCGCTGGCAGCCAGCGCCCGGCCGTACCTGGAGCAGGGTTACGACCACCTGCTACTGAGCTTCCAt                                                                                                  >  1:5982109/1‑150 (MQ=255)
                                                       aaaaCAGACGGTCAGCGAGCACCAGTTACCGTTGCAGATGCGCGTGCTGCAGCCGTTCTATGAGCACCCGGCGTACATCGAGGCGCTGGCAGCCAGCGCCCGGCCGTACCTGGAGCAGGGTTACGACCACCTGCTACTGAGCTTCCATgg                                                                                                >  1:4200122/1‑150 (MQ=255)
                                                       aaaaCAGACGGTCAGCGAGCACCAGTTACCGTTGCAGATGCGCGTGCTGCAGCCGTTCTATGAGCACCCGGCGTACATCGAGGCGCTGGCAGCCAGCGCCCGGCCGTACCTGGAGCAGGGTTACGACCACCTGCTACTGAGCTTCCATgg                                                                                                >  1:1630022/1‑150 (MQ=255)
                                                         aaCAGCCGGACAGCCAGCACCATTCCCCGTTGCAGATGCGCGGGCCGCCGCCTTTATTGGAGCACCCGGCGTACAGCGAGGCGCGGGCACCCAGCGCCCGGCCGTACCTGGAGCAGGGTTACGACCACCTGCTACTGAGCTTCCATGGtt                                                                                              <  2:7172520/150‑1 (MQ=255)
                                                             gcgggTCAGCGCGCACCAGTTCCCGTTGCTGTTGCGCGTGCTGCAGCCGTTCTATGAGCCCCGGGCGTACATCGAGGCGCTGGCTGCCAGCGCCCGGCCGTACCTGGAGCAGGGTTACGACCACCTGCTACTGAGCTTCCATGGTTTGcc                                                                                          <  2:5307315/147‑1 (MQ=255)
                                                             gccGGTCCGCGAGCACCTGTGACGGTTGCAGCTGCGCGGGCTGCCGCCGTTCTTTGAGCACCCGGCGTACATCGCGGCGCTGGCCGCCCGCGCCCGGCCGTACCTGGAGCAGGGTTACGACCACCTGCTACTGAGCTTCCATGGTTTGcc                                                                                          <  2:4725467/148‑1 (MQ=255)
                                                                ggTCAGCGAGCACCAGTTACCGTTGCAGATGCGCGTGCTGCAGCCGTTCTATGAGCACCCGGCGTACATCGAGGCGCTGGCAGCCAGCGCCCGGCCGTACCTGGAGCAGGGTTACGACCACCTGCTACTGAGCTTCCATGGTTTGCCGGa                                                                                       >  2:5612288/1‑150 (MQ=255)
                                                                  tCAGCGAGCACCAGTTACCGTTGCAGATGCGCGTGCTGCAGCCGTTCTATGAGCACCCGGCGTACATCGAGGCGCTGGCAGCCAGCGCCCGGCCGTACCTGGAGCAGGGTTACGACCACCTGCTACTGAGCTTCCATGGTTTGCCGGAAc                                                                                     >  2:5748303/1‑150 (MQ=255)
                                                                  tCAGCGAGCACCAGTTACCGTTGCAGATGCGCGTGCTGCAGCCGTTCTATGAGCACCCGGCGTACATCGAGGCGCTGGCAGCCAGCGCCCGGCCGTACCTGGAGCAGGGTTACGACCACCTGCTACTGAGCTTCCATGGTTTGCCGGAAc                                                                                     >  1:3424424/1‑150 (MQ=255)
                                                                   cAGCGAGCACCAGTTACCGTTGCAGATGCGCGTGCTGCAGCCGTTCTATGAGCACCCGGCGTACATCGAGGCGCTGGCAGCCAGCGCCCGGCCGTACCTGGAGCAGGGTTACGACCACCTGCTACTGAGCTTCCATGGTTTGCCGGAACg                                                                                    >  2:1666132/1‑150 (MQ=255)
                                                                         gCACCAGTTACCGTTGCAGATGCGCGTGCTGCAGCCGTTCTATGAGCACCCGGCGTACATCGAGGCGCTGGCAGCCAGCGCCCGGCCGTACCTGGAGCAGGGTTACGACCACCTGCTACTGAGCTTCCATGGTTTGCCGGAACGGCAccg                                                                              >  1:5914045/1‑149 (MQ=255)
                                                                         gCACCAGTTACCGTTGCAGATGCGCGTGCTGCAGCCGTTCTATGAGCACCCGGCGTACATCGAGGCGCTGGCAGCCAGCGCCCGGCCGTACCTGGAGCAGGGTTACGACCACCTGCTACTGAGCTTCCATGGTTTGCCGGAACGGCACCt                                                                              >  1:5489907/1‑150 (MQ=255)
                                                                            cctttgtccGTTGCTGCGGCGCGTGCTGCAGCCCTTCTATGAGCCCCCGGCGTCCAGCGCGGCGCGGGCACCCCCCCCCCGGCCGTACCTGGAGCAGGGTTACGACCACCTGCTACTGAGCTTCCATGGTTTGCCGGAACGGCACCTgaa                                                                           <  1:1321931/143‑1 (MQ=255)
                                                                            cctgtgcccGTTGCCGCTGCGCGGGCTGCCGCCGTTCTCTGCGCCCCCGGCGTCCCGCGCGGCGCGGGCCGCCCGCGCCCGGCCGTACCTGGAGCAGGGTTACGACCACCTGCTACTGAGCTTCCATGGTTTGCCGGAACGGCACCTgaa                                                                           <  1:1047613/143‑1 (MQ=255)
                                                                            ccAGTTACCGTTGCAGATGCGCGTGCTGCAGCCGTTCTATGAGCACCCGGCGTACATCGAGGCGCTGGcagccagc                                                                                                                                                     >  2:2804016/1‑76 (MQ=255)
                                                                            ccAGTTACCGTTGCAGATGCGCGTGCTGCAGCCGTTCTATGAGCACCCGGCGTACATCGAGGCGCTGGCAGCCAGCGCCCGGCCGTACCTGGAGCAGGGTTACGACCACCTGCTACTGAGCTTCCATGGTTTGCCGGACCGGCACCTgaa                                                                           >  2:3003072/1‑150 (MQ=255)
                                                                            ccAGTTACCGTTGCAGATGCGCGTGCTGCAGCCGTTCTATGAGCACCCGGCGTACATCGAGGCGCTGGCAGCCAGCGCCCGGCCGTACCTGGAGCAGGGTTACGACCACCTGCTACTGAGCTTCCATGGTTTGCCGGAACGGCACCGgaa                                                                           >  2:2256256/1‑150 (MQ=255)
                                                                             cAGTTACCGTTGCAGATGCGCGTGCTGCAGCCGTTCTATGAGCACCCGGCGTACATCGAGGCGCTGGCAGCCAGCGCCCGGCCGTACCTGGAGCAGGGTTACGACCACCTGCTACTGAGCTTCCATGGTTTGCCGGAACGGCACCTgaag                                                                          >  2:5225676/1‑150 (MQ=255)
                                                                              ctttACGGTTGCAGCGGGGCGTGCGGCAGCCGTTCTATGAGCCCCCGGCGTACAGCGCGGCGCTGGCAGCCCGCGCCCGGCCGTACCTGGAGCAGGGTTACGACCACCTGCTACTGAGCTTCCATGGTTTGCCGGAACGGCACCTgaaga                                                                         <  2:6101043/148‑1 (MQ=255)
                                                                                ttAGCGTTGCTGTGGCGCGGGGTGCTGCCGTTCACCGCGCCCCGGGCGCACCTGGCGGCGCGGCCCCCCCGCCCCCGGCCGTACCTGGAGCAGGGTTACGACCACCTGCTACTGAGCTTCCATGGTTTGCCGGAACGGCACCTGAAGaaa                                                                       <  1:6474942/150‑1 (MQ=255)
                                                                                  aCCGTTGCAGATGCGCGTGCTGCAGCCGTTCTATGAGCACCCGGCGTACATCGAGGCGCTGGCAGCCAGCGCCCGGCCGTACCTGGAGCAGGGTTACGACCACCTGCTACTGAGCTTCCATGGTTTGCCGGAACGGCACCTGAAGAAACt                                                                     >  1:472082/1‑150 (MQ=255)
                                                                                    cGTTGCAGATGCGCGTGCTGCAGCCGTTCTATGAGCACCCGGCGTACATCGAGGCGCTGGCAGCCAGCGCCCGGCCGTACCTGGAGCAGGGTTACGACCACCTGCTACTGAGCTTCCATGGTTTGCCGGAACGGCACCTGAAGAAACTGt                                                                   >  2:4325688/1‑150 (MQ=255)
                                                                                      ttGCAGATGCGCGTGCTGCAGCCGTTCTATGAGCACCCGGCGTACATCGAGGCGCTGGCAGCCAGCGCCCGGCCGTACCTGGAGCAGGGTTACGACCACCTGCTACTGAGCTTCCATGGTTTGCCGGAACGGCACCTGAAGAAACTGttt                                                                 >  1:7145417/1‑150 (MQ=255)
                                                                                       ggCAGAGGCGCGTGCTGCAGCCGTTCTATGAGCACCCGGCGTACATCGCGGCGCTGGCAGCCCGCGCCCGGCCGTACCTGGAGCAGGGTTACGACCACCTGCTACTGAGCTTCCATGGTTTGCCGGAACGGCACCTGAAGAAACTGTTTc                                                                <  1:6138535/149‑1 (MQ=255)
                                                                                          aGATGCGCGTGCTGCAGCCGTTCTATGAGCACCCGGCGTACATCGAGGCGCTGGCAGCCAGCGCCCGGCCGTACCTGGAGCAGGGTTACGACCACCTGCTACTGAGCTTCCATGGTTTGCCGGAACGGCACCTGAAGAAACTGTTTCCGa                                                             >  2:5115810/1‑150 (MQ=255)
                                                                                           gATGCGCGTGCTGCAGCCGTTCTATGAGCACCCGGCGTACATCGAGGCGCTGGCAGCCAGCGCCCGGCCGTACCTGGAGCAGGGTTACGTCCACCTGCTACTGAGCTTCCATGGTTTGCCGGCACGGCGCCTGAAGAAACTGTTTCCgct                                                            >  1:1839735/1‑148 (MQ=255)
                                                                                           gATGCGCGTGCTGCAGCCGTTCTATGAGCACCCGGCGTACATCGAGGCGCTGGCAGCCAGCGCCCGGCCGTACCTGGAGCAGGGTTACGACCACCTGCTACTGAGCTTCCATGGTTTGCCGGAACGGCACCTGAAGAAACTGTTTCCGaa                                                            >  1:7425900/1‑150 (MQ=255)
                                                                                               ggggTGCTGCCTCCTTTCGATGATCCCCGGGCTTACAGGGCGGCGCGGCCACCCCGCCCCCGGCCGTACCTGGAGCAGGGTTACGACCACCTGCTACTGAGCTTCCATGGTTTGCCGGAACGGCACCTGAAGAAACTGTTTCCGAAgggg                                                        <  1:1538939/147‑1 (MQ=255)
                                                                                                        cAGCCGTTCTATGAGCACCCGGCGTACATCGAGGCGCTGGCAGCCAGCGCCCGGCCGTACCTGGAGCAGGGTTACGACCACCTGCTACTGAGCTTCCATGGTTTGCCGGAACGGCACCTGAAGAAACTGTTTCCGAAGGGGGTCAAGCAc                                               >  1:2424030/1‑150 (MQ=255)
                                                                                                        cAGCCGTTCTATGAGCACCCGGCGTACATCGAGGCGCTGGCAGCCAGCGCCCGGCCGTACCTGGAGCAGGGTTACGACCACCTGCTACTGAGCTTCCATGGTTTGCCGGAACGGCACCTGAAGAAACTGTTTCCGAAGGGGGTCAAGCAc                                               >  1:6799396/1‑150 (MQ=255)
                                                                                                            cGTTCTATGAGCACCCGGCGTACATCGAGGCGCTGGCAGCCAGCGCCCGGCCGTACCTGGAGCAGGGTTACGACCACCTGCTACTGAGCTTCCATGGTTTGCCGGAACGGCACCTGAAGAAACTGTTTCCGAAGGGGGTCAAGCACGAcc                                           >  1:2143968/1‑150 (MQ=255)
                                                                                                             gTTCTATGAGCACCCGGCGTACATCGAGGCGCTGGCAGCCAGCGCCCGGCCGTACCTGGAGCAGGGTTACGACCACCTGCTACTGAGCTTCCATGGTTTGCCGGAACGGCACCTGAAGAAACTGTTTCCGAAGGGGGTCAAGCACGACCt                                          >  2:2293820/1‑150 (MQ=255)
                                                                                                             gTTCTATGAGCACCCGGCGTACATCGAGGCGCTGGCAGCCAGCGCCCGGCCGTACCTGGAGCAGGGTTACGACCACCTGCTACTGAGCTTCCATGGTTTGCCGGAACGGCACCTGAAGAAACTGTTTCCGAAGGGGGTCAAGCACGACCt                                          >  1:475061/1‑150 (MQ=255)
                                                                                                               tCTATGAGCCCCCGGCGTACCTCGCGGCGCGGGCCGCCAGCCCCCGGCCGTACCTGGAGCAGGGTTACGACCACCTGCTACTGAGCTTCCATGGTTTGCCGGAACGGCACCTGAAGAAACTGTTTCCGAAGGGGGTCAAGCACGACCTgc                                        <  2:5247403/150‑1 (MQ=255)
                                                                                                                cTATGAGCCCCCGGCGTACCTCGCGGCGCGGGCAGCCCGCCCCCGGCCGTACCTGGAGCAGGGTTACGACCACCTGCTACTGAGCTTCCATGGTTTGCCGGAACGGCACCTGAAGAAACTGTTTCCGAAGGGGGTCAAGCACGACCTgcg                                       <  2:3240451/150‑1 (MQ=255)
                                                                                                                 tATGAGCACCCGGCGTACATCGAGGCGCTGGCAGCCAGCGCCCGGCCGTACCTGGAGCAGGGTTACGACCACCTGCTACTGAGCTTCCATGGTTTGCCGGAACGGCACCTGAAGAAACTGTTTCCGAAGGGGGTCAAGCACGACCTgcgc                                      >  1:4457191/1‑150 (MQ=255)
                                                                                                                   tGAGCACCCGGCGTACATCGAGGCGCTGGCAGCCAGCGCCCGGCCGTACCTGGAGAAGGGTTACGACCACCTGCTACTGAGCTTCCATGGTTTGCCGGAACGGCACCTGAAGAAACTGTTTCCGAAGGGGGGCAAGAACGACCTgagcgc                                    >  1:3266729/1‑150 (MQ=255)
                                                                                                                      gCACCCGGCGTACATCGCGGCGCTGGCCGCCCGCGCCCGGCCGTACCTGGAGCAGGGTTACGACCACCTGCTACTGAGCTTCCATGGTTTGCCGGAACGGCACCTGAAGAAACTGTTTCCGAAGGGGGTCAAGCACGACCTGCGCGCTGc                                 <  1:3666927/150‑1 (MQ=255)
                                                                                                                      gCACCCGGCGTACATCGAGGCGCTGGCAGCCAGCGCCCGGCCGTACCTGGAGCAGGGTTACGACCACCTGCTACTGAGCTTCCATGGTTTGCCGGAACGGCACCTGAAGAAACTGTTTCCGAAGGGGGTCAAGCACGACCTGCGCGCTGc                                 >  2:4677934/1‑150 (MQ=255)
                                                                                                                       cACCCGGCGTACATCGAGGCGCTGGCAGCCAGCGCCCGGCCGTACCTGGAGCAGGGTTACGACCACCTGCTACTGAGCTTCCATGGTTTGCCGGAACGGCACCTGAAGAAACTGTTTCCGAAGGGGGTCAAGCACGACCTGCGCGCTGcc                                >  2:1370742/1‑150 (MQ=255)
                                                                                                                       cACCCGGCGTACATCGAGGCGCTGGCAGCCAGCGCCCGGCCGTACCTGGAGCAGGGTTACGACCACCTGCTACTGAGCTTCCATGGTTTGCCGGAACGGCACCTGAAGAAACTGTTTCCGAAGGGGGTCAAGCACGACCTGCGCGCTGcc                                >  1:4024026/1‑150 (MQ=255)
                                                                                                                          ccGGCTTACATCGCGGCGCTGGCAGCCAGCGCCCGGCCGTACCTGGAGCAGGGTTACGACCACCTGCTACTGAGCTTCCATGGTTTGCCGGAACGGCACCTGAAGAAACTGTTTCCGAAGGGGGTCAAGCACGACCTGCGCGCTGCCGAt                             <  2:3407901/150‑1 (MQ=255)
                                                                                                                          ccGGCGTACCGCGCGGCGCGGGCCGCCCGCGCCCGGCCGTACCTGGAGCAGGGTTACGACCACCTGCTACTGAGCTTCCATGGTTTGCCGGAACGGCACCTGAAGAAACTGTTTCCGAAGGGGGTCAAGCACGACCTGCGCGCTGCCGAt                             <  2:2494271/150‑1 (MQ=255)
                                                                                                                              cgtagctCGCGGCGGTGGCCGCCCCCCCCCGGCCGTACCTGGAGCAGGGTTACGACCACCTGCTACTGAGCTTCCATGGTTTGCCGGAACGGCACCTGAAGAAACTGTTTCCGAAGGGGGTCAAGCACGACCTGCGCGCTGCCGATtgct                         <  1:6855138/144‑1 (MQ=255)
                                                                                                                              cGTACCTCGCGGCGCGGGCCCCCCGCCCCCGGCCGTACCTGGAGCAGGGTTACGACCACCTGCTACTGAGCTTCCATGGTTTGCCGGAACGGCACCTGAAGAAACTGTTTCCGAAGGGGGTCAAGCACGACCTGCGCGCTGCCGATtgct                         <  1:2087832/150‑1 (MQ=255)
                                                                                                                              cGTACATCGCGGCGCTGGCCGCCCGCGCCCGGCCGTACCTGGAGCAGGGTTACGACCACCTGCTACTGAGCTTCCATGGTTTGCCGGAACGGCACCTGAAGAAACTGTTTCCGAAGGGGGTCAAGCACGACCTGCGCGCTGCCGATtgct                         <  1:3498858/150‑1 (MQ=255)
                                                                                                                              cGTACATCGCGGCGCTGGCAGCCAGCGCCCGGCCGTACCTGGAGCAGGGTTACGACCACCTGCTACTGAGCTTCCATGGTTTGCCGGAACGGCACCTGAAGAAACTGTTTCCGAAGGGGGTCAAGAACGACCTGCGCGCTGCCGATtgct                         <  2:1738827/150‑1 (MQ=255)
                                                                                                                              cGTACATCGAGGCGCTGGCAGCCAGCGCCCGGCCGTACCTGGAGCAGGGTTACGACCACCTGCTACTGAGCTTCCATGGTTTGCCGGAACGGCACCTGAAGAAACTGTTTCCGAAGGGGGTCAAGCACGACCTGCGCGCTGCCGATtgct                         >  2:1853733/1‑150 (MQ=255)
                                                                                                                               gTACCTCGCGGCGCGGGCCGCCCGCGCCCGGCCGTACCTGGAGCAGGGTTACGACCACCTGCTACTGAGCTTCCATGGTTTGCCGGAACGGCACCTGAAGAAACTGTTTCCGAAGGGGGTCAAGCACGACCTGCGCGCTGCCGATtgctg                        <  2:6961918/150‑1 (MQ=255)
                                                                                                                               gTACATCGCGGCGCCGGCCGCCCGCCCCCGGCCGTACCTGGAGCAGGGTTACGACCACCTGCTACTGAGCTTCCATGGTTTGCCGGAACGGCACCTGAAGAAACTGTTTCCGAAGGGGGTCAAGCACGACCTGCGCGCTGCCGATtgctg                        <  1:5321210/150‑1 (MQ=255)
                                                                                                                               gTACACCGCGGCGCGGGCACCCAGCCCCCGGCCGTACCTGGAGCAGGGTTACGACCACCTGCTACTGAGCTTCCATGGTTTGCCGGAACGGCACCTGAAGAAACTGTTTCCGAAGGGGGTCAAGCACGACCTGCGCGCTGCCGATtgctg                        <  1:5619769/150‑1 (MQ=255)
                                                                                                                                 aCATCGAGGCGCTGGCAGCCAGCGCCCGGCCGTACCTGGAGCAGGGTTACGACCACCTGCTACTGAGCTTCCATGGTTTGCCGGAACGGCACCTGAAGAAACAGTTTCCGAAGGGGGTCAATCACGACCTGCGCGCTGCCGATTGCTGcc                      >  1:7037763/1‑150 (MQ=255)
                                                                                                                                    gcgcggCGCGGGCAGCCCGCGCCCGGCCGTACCTGGAGCAGGGTTACGACCACCTGCTACTGAGCTTCCATGGTTTGCCGGAACGGCACCTGAAGAAACTGTTTCCGAAGGGGGTCAAGCACGACCTGCGCGCTGCCGATTGCTGCCACg                   <  2:2023069/146‑1 (MQ=255)
                                                                                                                                    cggAGGCGCTGGCCGCCACCCCCCGGCCGTACCTGGAGCAGGGTTACGACCACCTGCTACTGAGCTTCCATGGTTTGCCGGAACGGCACCTGAAGAAACTGTTTCCGAAGGGGGTCAAGCACGACCTGCGCGCTGCCGATTGCTGCCACg                   <  2:2568714/148‑1 (MQ=255)
                                                                                                                                     cGAGGCGCTGGCAGCCAGCGCCCGGCCGTACCTGGAGCAGGGTTACGACCACCTGCTACTGAGCTTCCATGGTTTGCCGGAACGGCACCTGAAGAAACTGTTTCCGAAGGGGGTCAAGCACGACCTGCGCGCTGCCGATTGCTGCCACgg                  >  1:4835553/1‑150 (MQ=255)
                                                                                                                                      gAGGGGCTGGCCGCCAGCGCCCGGCCGTACCTGGAGCAGGGTTACGACCACCTGCTACTGAGCTTCCATGGTTTGCCGGAACGGCACCTGAAGAAACTGTTTCCGAAGGGGGTCAAGCACGACCTGCGCGCTGCCGATTGCTGCCACGGt                 <  2:4518816/150‑1 (MQ=255)
                                                                                                                                      gAGGCGCTGGCAGCCAGCGCCCGGCCGTACCTGGAGCAGGGTTACGACCACCTGCTACTGAGCTTCCATGGTTTGCCGGAACGGCACCTGAAGAAACTGTTTCCGAAGGGGGTCAAGCACGACCTGCGCGCTGCCGATTGCTGCCACGGt                 >  1:2708215/1‑150 (MQ=255)
                                                                                                                                      caGGCGCTGGCAGCCAGCGCCCGGCCGTACCTGGAGCAGGGTTACGACCACCTGCTACTGAGCTTCCATGGTTTGCCGGAACGGCACCTGAAGAAACTGTTTCCGAAGGGGGTCAAGCACGACCTGCGCGCTGCCGATTGCTGCCACGGt                 >  1:937866/2‑150 (MQ=255)
                                                                                                                                        ggCGCTGGCTGCCAGCGCCCGGCCGTACCTGGAGCAGGGTTACGACCACCTGCTACTGAGCTTCCATGGTTTGCCGGAACGGCACCTGAAGAAACTGTTTCCGAAGGGGGTCAAGCACGACCTGCGCGCTGCCGATTGCTGCCACGGTGc               >  1:5592070/1‑150 (MQ=255)
                                                                                                                                         gcgcTGGCAGCCAGCGCCCGGCCGTACCTGGAGCAGGGTTACGACCACCTGCTACTGAGCTTCCATGGTTTGCCGGAACGGCACCTGAAGAAACTGTTTCCGAAGGGGGTCAAGCACGACCTGCGCGCTGCCGATTGCTGCCACGGTGCa              >  2:5963756/1‑150 (MQ=255)
                                                                                                                                             tGGCAGCCAGCGCCCGGCCGTACCTGGAGCAGGGTTACGACCACCTGCTACTGAGCTTCCATGGTTTGCCGGAACGGCACCTGAAGAAACTGTTTCCGAAGGGGGTCAAGCACGACCTGCGCGCTGCCGATTGCTGCCACGGTGCAAccg          <  2:2435083/150‑1 (MQ=255)
                                                                                                                                              ggCAGCCAGCGCCCGGCCGTACCTGGAGCAGGGTTACGACCACCTGCTACTGAGCTTCCATGGTTTGCCGGAACGGCACCTGAAGAAACTGTTTCCGAAGGGGGTCAAGCACGACCTGCGCGCTGCCGATTGCTGCCACGGTGCAAccgc         >  1:5646679/1‑150 (MQ=255)
                                                                                                                                                  gccagcGCCCGGCCGTACCTGGAGCAGGGTTACGACCACCTGCTACTGAGCTTCCATGGTTTGCCGGAACGGCACCTGAAGAAACTGTTTCCGAAGGGGGTCAAGCACGACCTGCGCGCTGCCGATTGCTGCCACGGTGCAACCGCCGaa     >  2:3554327/1‑150 (MQ=255)
                                                                                                                                                   ccagcGCCCGGCCGTACCTGGAGCAGGGTTACGACCACCTGCTACTGAGCTTCCATGGTTTGCCGGAACGGCACCTGAAGAAACTGTTTCCGAAGGGGGTCAAGCACGACCTGCGCGCTGCCGATTGCTGCCACGGTGCAACCGCCGaag    >  2:4134350/1‑150 (MQ=255)
                                                                                                                                                     agcGCCCGGCCGTACCTGGAGCAGGGTTACGACCACCTGCTACTGAGCTTCCATGGTTTGCCGGAACGGCACCTGAAGAAACTGTTTCCGAAGGGGGTCAAGCACGACCTGCGCGCTGCCGATTGCTGCCACGGTGCAAccgc         >  2:6903296/1‑143 (MQ=255)
                                                                                                                                                     agcGCCCGGCCGTACCTGGAGCAGGGTTACGACCACCTGCTACTGAGCTTCCATGGTTTGCCGGAACGGCACCTGAAGAAACTGTTTCCGAAGGGGGTCAAGCACGACCTGCGCGCTGCCGATTGCTGCCACGGTGCAAccgc         <  1:6903296/143‑1 (MQ=255)
                                                                                                                                                     agcGCCCGGCCGTACCTGGAGCAGGGTTACGACCACCTGCTACTGAGCTTCCATGGTTTGCCGGAACGGCACCTGAAGAAACTGTTTCCGAAGGGGGTCAAGCACGACCTGCGCGCTGCCGATGGCTGCCACGGTGCAACCGCCGaagga  >  2:5494968/1‑148 (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 < 12 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑