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

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_002947863,3830AC12.3% 593.8 / 11.1 242S172R (AGC→CGC) hemHferrochelatase
Reads supporting (aligned to +/- strand):  ref base A (163/41);  new base C (3/35);  total (166/76)
Fisher's exact test for biased strand distribution p-value = 1.22e-17
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 266 total aligned reads displayed.

CACCGCTTTATCCACAGTTCGCCGACAGCACGGTGACCACTGTGATCGAGCAGGCAAAACAGACGGTCAGCGAGCACCAGTTACCGTTGCAGATGCGCGTGCTGCAGCCGTTCTATGAGCACCCGGCGTACATCGAGGCGCTGGCAGCCAGCGCCCGGCCGTACCTGGAGCAGGGTTACGACCACCTGCTACTGAGCTTCCATGGTTTGCCGGAACGGCACCTGAAGAAACTGTTTCCGAAGGGGGTCAAGCACGACCTGCGCGCTGCCGATTGCTGCCACGGTGCAACCGCCGAAGTA  >  NC_002947/863234‑863532
                                                                                                                                                     |                                                                                                                                                     
cACCGCTTTATCCACAGTTCGCCGACAGCACGGTGACCACTGTGATCGAGCAGGCAAAACAGACGGTCAGCGAGCACCAGTTACCGTTGCAGATGCGCGTGCTGCAGCCGTTCTATGAGCACCCGGCGTACATCGAGGCGCTGGcagcca                                                                                                                                                       >  1:5003966/1‑150 (MQ=255)
 aCCGCTTTATCCACAGTTCGCCGACAGCACGGTGACCACTGTGATCGAGCAGGCAAAACAGACGGTCAGCGAGCACCAGTTACCGTTGCAGATGCGCGTGCTGCAGCCGTTCTATGAGCACCCGGCGTACATCGAGGCGCTGGcagccag                                                                                                                                                      >  1:5593473/1‑150 (MQ=255)
 aCCGCTTTATCCACAGTTCGCCGACAGCACGGTGACCACTGTGATCGAGCAGGCAAAACAGACGGTCAGCGAGCACCAGTTACCGTTGCAGATGCGCGTGCTGCAGCCGTTCTATGAGCACCCGGCGTACATCGAGGCGCTGGcagccag                                                                                                                                                      >  2:2271931/1‑150 (MQ=255)
 aCCGCTTTATCCACAGTTCGCCGACAGCACGGTGACCACTGTGATCGAGCAGGCAAAACAGACGGTCAGCGAGCACCAGTTACCGTTGCAGATGCGCGTGCTGCAGCCGTTCTATGAGCACCCGGCGTACATCGAGGCGCTGGcagccag                                                                                                                                                      >  1:2882682/1‑150 (MQ=255)
 aCCGCTTTATCCACAGTTCGCCGACAGCACGGTGACCACTGTGATCGAGCAGGCAAAACAGACGGTCAGCGAGCACCAGTTACCGTTGCAGATGCGCGTGCTGCAGCCGTTCTATGAGCACCCGGCGTACATCGAGGCGCTGGcagccag                                                                                                                                                      >  2:1709503/1‑150 (MQ=255)
  ccGCTTTATCCACAGTTCGCCGACAGCACGGTGACCACTGTGATCGAGCAGGCAAAACAGACGGTCAGCGAGCACCAGTTACCGTTGCAGATGCGCGTGCTGCAGCCGTTCTATGAGCACCCGGCGTACATCGAGGCGCTGGcagccagc                                                                                                                                                     <  1:2987815/150‑1 (MQ=255)
  ccGCTTTATCCACAGTTCGCCGACAGCACGGTGACCACTGTGATCGAGCAGGCAAAACAGACGGTCAGCGAGCACCAGTTACCGTTGCAGATGCGCGTGCTGCAGCCGTTCTATGAGCACCCGGCGTACATCGAGGCGCTGGcagccagc                                                                                                                                                     <  1:4987017/150‑1 (MQ=255)
   cGCTTTATCCACAGTTCGCCGACAGCACGGTGACCACTGTGATCGAGCAGGCAAAACAGACGGTCAGCGAGCACCAGTTACCGTTGCAGATGCGCGTGCTGCAGCCGTTCTATGAGCACCCGGCGTACATCGAGGCGCTGGCAGCCAgcg                                                                                                                                                    >  1:1044350/1‑150 (MQ=255)
    gCTTTATCCACAGTTCGCCGACAGCACGGGGACCACTGTGATCGAGCAGGCAAAACAGACGGTCAGCGAGCACCAGTTACCGTTGCAGATGCGCGTGCTGCAGCCGTTCTATGAGCACCCGGCGTACATCGAGGCGCTGGCAGCCAgcgc                                                                                                                                                   <  1:1741440/150‑1 (MQ=255)
     cTTTATCCACAGTTCGCCGACAGCACGGTGACCACTGTGATCGAGCAGGCAAAACAGACGGTCAGCGAGCACCAGTTACCGTTGCAGATGCGCGTGCTGCAGCCGTTCTATGAGCACCCGGCGTACATCGAGGCGCTGGCAGCCAGCGcc                                                                                                                                                  >  2:1088315/1‑150 (MQ=255)
      tttATCCACAGTTCGCCGACAGCACGGTGACCACTGTGATCGAGCAGGCAAAACAGACGGTCAGCGAGCACCAGTTACCGTTGCAGATGCGCGTGCTGCAGCCGTTCTATGAGCACCCGGCGTACATCGAGGCGCTGGCAGCCAGCGccc                                                                                                                                                 <  1:895923/150‑1 (MQ=255)
      tttATCCACAGTTCGCCGACAGCACGGTGACCACTGTGATCGAGCAGGCAAAACAGACGGTCAGCGAGCACCAGTTACCGTTGCAGATGCGCGTGCTGCAGCCGTTCTATGAGCACCCGGCGTACATCGAGGCGCTGGCAGCCAGCGccc                                                                                                                                                 <  2:4500810/150‑1 (MQ=255)
       ttATCCACAGTTCGCCGACAGCACGGTGACCACTGTGATCGAGCAGGCAAAACAGACGGTCAGCGAGCACCAGTGACCGTTGCAGATGCGCGTGCTGCAGCCGTTCTATGAGCACCCGGCGTACATCGAGGCGCTGGCAGCCAGCGCCCg                                                                                                                                                <  2:4654117/150‑1 (MQ=255)
        tATCCACAGTTCGCCGACAGCACGGTGACCACTGTGATCGAGCAGGCAAAACAGACGGTCAGCGAGCACCAGTTACCGTTGCAGATGCGCGTGCTGCAGCCGTTCTATGAGCACCCGGCGTACATCGAGGCGCTGGCAGCCAGCGCCCgg                                                                                                                                               <  1:2684120/150‑1 (MQ=255)
            cacaGTTCGCCGACAGCACGGTGACCACTGTGATCGAGCAGGCAAAACAGACGGTCAGCGAGCACCAGTTACCGTTGCAGATGCGCGTGCTGCAGCCGTTCTATGAGCACCCGGCGTACATCGAGGCGCTGGCAGCCAGCGCCCGGCCGt                                                                                                                                           >  1:2207436/1‑150 (MQ=255)
            cacaGTTCGCCGACAGCACGGTGACCACTGTGATCGAGCAGGCAAAACAGACGGTCAGCGAGCACCAGTTACCGTTGCAGATGCGCGTGCTGCAGCCGTTCTATGAGCACCCGGCGTACATCGAGGCGCTGGCAGCCAGCGCCCGGCCGt                                                                                                                                           >  1:2648598/1‑150 (MQ=255)
             acaGTTCGCCGACAGCACGGTGACCACTGTGATCGAGCAGGCAAAACAGACGGTCAGCGAGCACCAGTTACCGTTGCAGATGCGCGTGCTGCAGCCGTTCTATGAGCACCCGGCGTACATCGAGGCGCTGGCAGCCAGCGCCCGGCCGTa                                                                                                                                          >  1:2383704/1‑150 (MQ=255)
             acaGTTCGCCGACAGCACGGTGACCACTGTGATCGAGCAGGCAAAACAGACGGTCAGCGAGCACCAGTTACCGTTGCAGATGCGCGTGCTGCAGCCGTTCTATGAGCACCCGGCGTACATCGAGGCGCTGGCAGCCAGCGCCCGGCCGTa                                                                                                                                          >  1:1137525/1‑150 (MQ=255)
             acaGTTCGCCGACAGCACGGTGACCACTGTGATCGAGCAGGCAAAACAGACGGTCAGCGAGCACCAGTTACCGTTGCAGATGCGCGTGCTGCAGCCGTTCTATGAGCACCCGGCGTACATCGAGGCGCTGGCAGCCAGCGCCCGGCCGTa                                                                                                                                          >  1:3494864/1‑150 (MQ=255)
             acaGTTCGCCGACAGCACGGTGACCACTGTGATCGAGCAGGCAAAACAGACGGTCAGCGAGCACCAGTTACCGTTGCAGATGCGCGTGCTGCAGCCGTTCTATGAGCACCCGGCGTACATCGAGGCGCTGGCAGCCAGCGCCCGGCCGTa                                                                                                                                          >  1:3279649/1‑150 (MQ=255)
             acaGTTCGCCGACAGCACGGTGACCACTGTGATCGAGCAGGCAAAACAGACGGTCAGCGAGCACCAGATACCGTTGCAGATGCGCGTGCTGCAGCCGTTCTATGAGCACCCGGCGTACATCGAGGCGCTGGCAGCCAGCGCCCGGCCGTa                                                                                                                                          >  2:4805578/1‑150 (MQ=255)
                 ttCGCCGACAGCACGGTGACCACTGTGATCGAGCAGGCAAAACAGACGGTCAGCGAGCACCAGTTACCGTTGCAGATGCGCGTGCTGCAGCCGTTCTATGAGCACCCGGCGTACATCGAGGCGCTGGCAGCCAGCGCCCGGCCGTACCTg                                                                                                                                      >  2:3707470/1‑150 (MQ=255)
                 tccGCCGACACCACGGTGACCAGTGTGTTCGAGCTGGCAAAACTGACGGTACGCGAGCACACGTTACCTTTGCAGATGCGCGTGCTGCAGCCGTTATATTAGCACCCGGCGTACATCGAGGCGCTGGCAGCCCGCGCCCGGCCGTACCTg                                                                                                                                      <  2:5469785/148‑1 (MQ=255)
                   cGCCGACAGCACGGTGACCACTGTGATCGATCTGGCAAAACAGTCGGTCTGCGAGCACCTTTTACCTTTGCAGTTGGGCGTGCTGCAGCCGTTATATGAGCACCCTGCGTACATCGCGGCGCTGGCAGCCAGCGCCCGGCCGTACCTGGa                                                                                                                                    <  2:1161363/150‑1 (MQ=255)
                    cccGACAGCCCGGCGACCTCTGTGAGCGAGCAGGCAAAACAGACGGTCAGCGAGCACCAGTTACGGTTGCAGATGGGCGTGCGGCCGCCGTTCTTGGCGCCCCCGGCTTACATCGCGGCGCTGGCACCCAGCGCCCGGCCGTACCTGGAg                                                                                                                                   <  2:4930229/149‑1 (MQ=255)
                     ccGACAGCACGGTGACCACTGTGATCGAGCAGGCAAAACAGACGGTCAGCGAGCACCAGTTACCGTTGCAGATGCGCGTGCTGCAGCCGTTCTATGAGCACCCGGCGTACATCGAGGCGCTGGCAGCCAGCGCCCGGCCGTACCTGGAGc                                                                                                                                  >  1:2947295/1‑150 (MQ=255)
                          aGCACGGTGACCACTGTGATCGAGCAGGCAAAACAGACGGTCAGCGAGCACCAGTTACCGTTGCAGATGCGCGTGCTGCAGCCGTTCTATGAGCACCCGGCGTACATCGAGGCGCTGGCAGCCAGCGCCCGGCCGTACCTGGAGCAGGGt                                                                                                                             >  1:5048952/1‑150 (MQ=255)
                          aGCACGGTGACCACTGTGATCGAGCAGGCAAAACAGACGGTCAGCGAGCACCAGTTACCGTTGCAGATGCGCGTGCTGCAGCCGTTCTATGAGCACCCGGCGTACATCGAGGCGCTGGCAGCCAGCGCCCGGCCGTACCGGGAGCAGGGt                                                                                                                             >  1:1971255/1‑150 (MQ=255)
                             agggTGACCCGTGTGATGGAGCAGTCAAAACAGACGTTATGCGCGCACCGGTTTCCGTTGCTGTGGGGGGTGCGGCAGCCGTTATATGAGCAGCCGGCTTACTTCGAGGCGCTGGCAGCCAGCGCCCGGCCGTACCTGGAGCAGGGTTAc                                                                                                                          <  2:4085353/148‑1 (MQ=255)
                              cGGTGACCACTGTGATCGAGCAGGCAAAACAGACGGTCAGCGAGCACCAGTTACCGTTGCAGATGCGCGTGCTGCAGCCGTTCTATGAGCACCCGGCGTACATCGAGGCGCTGGCAGCCAGCGCCCGGCCGTACCTGGAGCAGGGTTACg                                                                                                                         >  2:4056106/1‑150 (MQ=255)
                                gTGACCACTGTGATCGAGCAGGCAAAACAGACGGTCAGCGAGCACCAGTTACCGTTGCAGATGCGCGTGCTGCAGCCGTTCTATGAGCACCCGGCGTACATCGAGGCGCTGGCAGCCAGCGCCCGGCCGTACCTGGAGCAGGGTTACGac                                                                                                                       >  2:794788/1‑150 (MQ=255)
                                 tGACCACTGTGTTCGAGCAGCAAAACCAGCCGGCCAGCGAGCACCTGTGCCGGGTGCCGTTGCGGGTGCTGCCGCCGTTCTTTGAGCACCGGGGGTACATCGCGGCGCGGGCAGCCCGCGCCCGGCCGTACCTGGAGCAGGGTTACGacc                                                                                                                      <  2:1877534/150‑1 (MQ=255)
                                            aTCGAGCAGGCAAAACAGACGGTCAGCGAGCACCAGTTACCGTTGCAGATGCGCGTGCTGCAGCCGTTCTATGAGCACCCGGCGTACATCGAGGCGCTGGCAGCCAGCGCCCGGCCGTACCTGGAGCAGGGTTACGACCACCTGCTACTg                                                                                                           >  2:2383948/1‑150 (MQ=255)
                                               gAGCAGGCAAAACAGACGGTCAGCGAGCACCAGTTACCGTTGCAGATGCGCGTGCTGCAGCCGTTCTATGAGCACCCGGCGTACATCGAGGCGCTGGCAGCCAGCGCCCGGCCGTACCTGGAGCAGGGTTACGACCACCTGCTACTGAGc                                                                                                        >  2:5080427/1‑150 (MQ=255)
                                                aGCAGGCAACACCGCCGGCCGGCGCGCACCCGTGACGTTTGCCTCTGCGCGTGCTGCCGCGTTTATATGCCCACCCGGCGTACTCCGCGGCGCTGGCAGCCAGCGCCCGGCCGTACCTGGAGCAGGGTTACGACCACCTGCTACTGAGCt                                                                                                       <  2:5548170/150‑1 (MQ=255)
                                                aGCAGGCAAAACAGACGGTCAGCGAGCACCAGTTACCGTTGCAGATGCGCGTGCTGCAGCCGTTCTATGAGCACCCGGCGTACATCGAGGCGCTGGCAGCCAGCGCCCGGCCGTACCTGGAGCAGGGTTACGACCACCTGCTACTGAGCt                                                                                                       >  1:1832149/1‑150 (MQ=255)
                                                aGCAGGCAAAACAGACGGTCAGCGAGCACCAGTTACCGTTGCAGATGCGCGTGCTGCAGCCGTTCTATGAGCACCCGGCGTACATCGAGGCGCTGGCAGCCAGCGCCCGGCCGTACCTGGAGCAGGGTTACGACCACCTGCTACAGAGCt                                                                                                       >  1:2245659/1‑150 (MQ=255)
                                                     taaaaaCAGACGGTCCGCGAGCACCTTTCACGTTTGCTTTTGCGCTTGCTGCAGCCTTTTTTTTACCCCCCGGCTTACAGCGCGGGGCGGGCACCCCGCCCCCGGCCGTACCTGGAGCAGGGTTACGACCACCTGCTACTGAGCTTCCAt                                                                                                  <  2:1428658/148‑1 (MQ=255)
                                                      cAAAACTGCCGTTCAGCGCGCACCTTTTACCGGTGCAGAGGCGCGTGCGGCCGCCGTTCTCGGCCCCCCCGGCTTACTTCGTGGCCGTGGCCCCCAGCGCCCGGCCGTACCTGGAGCAGGGTTACGACCACCTGCTACTGAGCTTCCATg                                                                                                 <  2:1409862/150‑1 (MQ=255)
                                                       aaaaCAGACGGTCAGCGAGCACCAGTTACCGTTGCAGATGCGCGTGCTGCAGCCGTTCTATGAGCACCCGGCGTACATCGAGGCGCTGGCAGCCAGCGCCCGGCCGTACCTGGAGCAGGGTTACGACCACCTGCTACTGAGCTTCCATgg                                                                                                >  1:479485/1‑150 (MQ=255)
                                                       aaaaCAGACGGTCAGCGAGCACCAGTTACCGTTGCAGATGCGCGTGCTGCAGCCGTTCTATGAGCACCCGGCGTACATCGAGGCGCTGGCAGCCAGCGCCCGGCCGTACCTGGAGCAGGGGTACGACCACCTGCTACTGAGCTGCCATgg                                                                                                >  1:1064556/1‑150 (MQ=255)
                                                         aaCAGACGGTCAGCGAGCACCAGTTACCGTTGCAGATGCGCGTGCTGCAGCCGTTCTATGAGCACCCGGCGTACATCGAGGCGCTGGCAGCCAGCGCCCGGCCGTACCTGGAGCAGGGTTACGACCACCTGCTACTGAGCTTCCATGGtt                                                                                              >  1:853241/1‑150 (MQ=255)
                                                          cccgtcGGTCTGCGATCCCCAGGTCCCGTTGCAGATGGGGGTGGCGCCGCCGTAGACTGAGCCCCCGGCGTACCTCGAGGCGCTGGCAGCCCGCCCCCGGCCGTACCTGGAGCAGGGTTACGACCACCTGCTACTGAGCTTCCATGGttt                                                                                             <  2:4372424/145‑1 (MQ=255)
                                                          accgACGGCCAGCGAGCACCTTTTCCCGTTGCTGATGCGGGTGCTGCAGCCTTTCTATGAGCACCCGGCGTACATCGCGGCGCGGGCCGCCAGCGCCCGGCCGTACCTGGAGCAGGGTTACGACCACCTGCTACTGAGCTTCCATGGttt                                                                                             <  2:2134941/147‑1 (MQ=255)
                                                          aCAGACGGTCAGCGAGCACCAGTTACCGTTGCAGATGCGCGTGCTGCAGCCGTTCTATGAGCACCCGGCGTACATCGAGGCGCTGGCAGCCAGCGCCCGGCCGTACCTGGAGCAGGGTTACGACCACCTGCTACTGAGCTTCCATGGttt                                                                                             >  2:2676338/1‑150 (MQ=255)
                                                           cAGACGGTCAGCGAGCACCAGTTACCGTTGCAGATGCGCGTGCTGCAGCCGTTCTATGAGCACCCGGCGTACATCGAGGCGCTGGCAGCCAGCGCCCGGCCGTACCTGGAGCAGGGTTACGACCACCTGCTACTGAGCTTCCATGGTTTg                                                                                            >  1:455627/1‑150 (MQ=255)
                                                            aGACGGTCAGCGAGCACCAGTTACCGTTGCAGATGCGCGTGCTGCAGCCGTTCTATGAGCACCCGGCGTACATCGAGGCGCTGGCAGCCAGCGCCCGGCCGTACCTGGAGCAGGGTTACGACCACCTGCTACTGAGCTTCCATGGTTTGc                                                                                           >  1:1192320/1‑150 (MQ=255)
                                                             gtcGGTCAGCGAGCACCAGTTACCGTTGCAGATGCGCGTGCTGCAGCCGTTCTATGAGCACCCGGCTTACATCGCGGCGCGGGCCGCCCGCGCCCGGCCGTACCTGGAGCAGGGTTACGACCACCTGCTACTGAGCTTCCATGGTTTGcc                                                                                          <  1:4655566/148‑1 (MQ=255)
                                                                ggTCAGCGAGCACCAGTTACCGTTGCAGATGCGCGTGCTGCAGCCGTTCTATGAGCACCCGGCGTACATCGAGGCGCTGGCAGCCAGCGCCCGGCCGTACCTGGAGCAGGGTTACGACCACCTGCTACTGAGCTTCCATGGTTTGCCGGa                                                                                       >  1:2846664/1‑150 (MQ=255)
                                                                ggTCAGCGAGCACCAGTTACCGTTGCAGATGCGCGTGCTGCAGCCGTTCTATGAGCACCCGGCGTACATCGAGGCGCTGGCAGCCAGCGCCCGGCCGTACCTGGAGCAGGGTTACGACCACCTGCTACTGAGCTTCCATGGTTTGCCGGa                                                                                       >  1:4536999/1‑150 (MQ=255)
                                                                ggTCAGCGAGCACCAGTTACCGTTGCAGATGCGCGTGCTGCAGCCGTTCTATGAGCACCCGGCGTACATCGAGGCGCTGGCAGCCAGCGCCCGGCCGTACCTGGAGCAGGGTTACGACCACCTGCTACTGAGCTTCCATGGTTTGCCGGa                                                                                       >  2:4662444/1‑150 (MQ=255)
                                                                  tCAGCGAGCACCAGTTACCGTTGCAGATGCGCGTGCTGCAGCCGTTCTATGAGCACCCGGCGTACATCGAGGCGCTGGCAGCCAGCGCCCGGCCGTACCTGGAGCAGGGTTACGACCACCTGCTACTGAGCTTCCATGGTTTGCCGGAAc                                                                                     >  1:2141218/1‑150 (MQ=255)
                                                                  tCAGCGAGCACCAGTTACCGTTGCAGATGCGCGTGCTGCAGCCGTTCTATGAGCACCCGGCGTACATCGAGGCGCTGGCAGCCAGCGCCCGGCCGTACCTGGAGCAGGGTTACGACCACCTGCTACTGAGCTTCCATGGTTTGCCGGAAc                                                                                     >  2:1288757/1‑150 (MQ=255)
                                                                  tCAGCGACCACCCTTGACCGTTGCCGCGGCGCGGGCTGCCGCCTTTCTCTGCGCCCCCGGCGTACCGCGGGGCGCTGGCCGCCCGCCCCCGGCCGTACCTGGAGCAGGGTTACGACCACCTGCTACTGAGCTTCCATGGTTTGCCGGAAc                                                                                     <  2:1501399/150‑1 (MQ=255)
                                                                     gCGAGCACCAGTTACCGTTGCAGATGCGCGTGCTGCAGCCGTTCTATGAGCACCCGGCGTACATCGAGGCGCTGGCAGCCAGCGCCCGGCCGTACCTGGAGCAGGGTTACGACCACCTGCTACTGAGCTTCCATGGTTTGCCGGAACGGc                                                                                  >  2:3968861/1‑150 (MQ=255)
                                                                     gCGAGCACCAGTTACCGTTGCAGATGCGCGTGCTGCAGCCGTTCTATGAGCACCCGGCGTACATCGAGGCGCTGGCAGCCAGCGCCCGGCCGTACCTGGAGCAGGGTTACGACCACCTGCTACTGAGCTTCCATGGTTTGCCGGAACGGc                                                                                  >  2:2978615/1‑150 (MQ=255)
                                                                       gtccACCTGTGACCGTTGCAGTGGGGGGTGCTGCCGCCGTTCTTTGACCCCCCGTCTTACCGGGAGGCGCTGGCAGCCAGCGCCCGGCCGTACCTGGAGCAGGGTTACGACCACCTGCTACTGAGCTTCCATGGTTTGCCGGAACGGCAc                                                                                <  2:3225410/147‑1 (MQ=255)
                                                                            ccAGTTACCGTTGCAGATGCGCGTGCTGCAGCCGTTCTATGAGCACCCGGCGTACATCGAGGCGCTGGCAGCCAGCGCCCGGCCGTACCTGGAGCAGGGTTACGACCACCTGCTACTGAGCTTCCATGGTTTGCCGGAACGGCACCTgaa                                                                           >  1:249837/1‑150 (MQ=255)
                                                                             cAGTTACCGTTGCAGATGCGCGTGCTGCAGCCGTTCTATGAGCACCCGGCGTACATCGAGGCGCTGGCAGCCAGCGCCCGGCCGTACCTGGAGCAGGGTTACGACCACCTGCTACTGAGCTTCCATGGTTTGCCGGAACGGCACCTgaag                                                                          >  1:3233586/1‑150 (MQ=255)
                                                                                   ccGTTGCAGATGCGCGTGCTGCAGCCGTTCTATGAGCACCCGGCGTACATCGAGGCGCTGGCAGCCAGCGCCCGGCCGTACCTGGAGCAGGGTTACGACCACCTGCTACTGAGCTTCCATGGTTTGCCGGAACGGCACCTGAAGAAACTg                                                                    >  2:1584434/1‑150 (MQ=255)
                                                                                      ttGCAGATGCGCGTGCTGCAGCCGTTCTATGAGCACCCGGCGTACATCGAGGCGCTGGCAGCCAGCGCCCGGCCGTACCTGGAGCAGGGTTACGACCACCTGCTACTGAGCTTCCATGGTTTGCCGGAACGGCACCTGACGAAACTGttt                                                                 >  1:5442463/1‑150 (MQ=255)
                                                                                      ttGCAGATGCGCGTGCTGCAGCCGTTCTATGAGCACCCGGCGTACATCGAGGCGCTGGCAGCCAGCGCCCGGCCGTACCTGGAGCAGGGTTACGACCACCTGCTACTGAGCTTCCATGGTTTGCCGGAACGGCACCTGAAGAAACTGttt                                                                 >  2:559762/1‑150 (MQ=255)
                                                                                      ttGCAGATGCGCGTGCTGCAGCCGTTCTATGAGCACCCGGCGTACATCGAGGCGCTGGCAGCCAGCGCCCCGCCGTACCTGGAGCAGGGTGACGACCACCTGCTACTGAGCTTCCATGGTTTGCCGGAACGGCACCTGCAGAAACAGttt                                                                 >  1:1109626/1‑150 (MQ=255)
                                                                                      ttGCAGATGCGCGTGCGGCAGCCTTTCTTGGAGCACCCGGCTTACATCGAGGCGCTGGCAGCCCGCGCCCGGCCGTACCTGGAGCAGGGTTACGACCACCTGCTACTGAGCTTCCATGGTTTGCCGGAACGGCACCTGAAGAAACTGttt                                                                 <  2:1252905/150‑1 (MQ=255)
                                                                                       tgttgctGCGCGTGCGGCCGCCGTTCTAGGCGCACCCGGCGCACATGGCGGCGCTGGCCGCCCCCGCCCGGCCGTACCTGGAGCAGGGTTACGACCACCTGCTACTGAGCTTCCATGGTTTGCCGGAACGGCACCTGAAGAAACTGTTTc                                                                <  1:1299523/144‑1 (MQ=255)
                                                                                            ttgggcgTGCTGCCGGCGTTGCAGCAGCTCCGGGCGTCCTGCGTGGCGCTGGCCGCCCGCCCCCGGCCGTACCTGGAGCAGGGTTACGACCACCTGCTACTGAGCTTCAATGGTTTGCCGGAACGGCACCTGAAGAAACTGTTTCCGAAg                                                           <  1:5329632/146‑1 (MQ=255)
                                                                                            ttGCGCGGGCTGCCGCCGTTCTCCGCGCCCCCGGCGTACCCCGCGGCGCTGCCCGCCCCCCCCCGGCCGTACCTGGAGCAGGGTTACGACCACCTGCTACTGAGCTTCCATGGTTTGCCGGAACGGCACCTGAAGAAACTGTTTCCGAAg                                                           <  1:3450041/149‑1 (MQ=255)
                                                                                            tggcgcGTGCGTCCGCCTTTCTAGGAGCACCCGGCGTAACTCGCGGCGCTGGCCCCCACCCCCCGGCCGTACCTGGAGCAGGGTTACGACCACCTGCTACTGAGCTTCCATGGTTTGCCGGAACGGCACCTGAAGAAACTGTTTCCGAAg                                                           <  2:282304/148‑1 (MQ=255)
                                                                                                gcgTGCTGCAGCCGTTCTATGAGCACCCGGCGTACATCGAGGCGCTGGCAGCCAGCGCCCGGCCGTACCTGGAGCAGGGTTACGACCACCTGCTACTGAGCTTCCATGGTTTGCCGGAACGGCACCTGAAGAAACTGTTTCCGAAggggg                                                       >  1:1900902/1‑150 (MQ=255)
                                                                                                      tgcAGCCGTTCTATGAGCACCCGGCGTACATTGAGGCGCTGGCAGCCAGCGCCCGGCCGTACCTGGAGCAGGGTTACGACCACCTGCTACTGAGCTTCCATGGTTTGCCGGAACGGCACCTGAAGAAACTGTTTCCGAAGGGGGTCAAGc                                                 >  1:5359909/1‑150 (MQ=255)
                                                                                                      tgcAGCCGTTCTATGAGCACCCGGCGTACATCGAGGCGCTGGCAGCCAGCGCCCGGCCGTACCTGGAGCAGGGTTACGACCACCTGCTACTGAGCTTCCATGGTTTGCCGGAACGGCACCTGAAGAAACTGTTTCCGAAGGGGGTCAAGc                                                 >  2:943779/1‑150 (MQ=255)
                                                                                                       gccgCCGTTCTAGCCGCCCCCGGCGTACAGCGCGGCGCTGCCCGCCCCCCCCCGGCCGTACCTGGAGCAGGGTTACGACCACCTGCTACTGAGCTTCCATGGTTTGCCGGAACGGCACCTGAAGAAACTGTTTCCGAAGGGGGTCAAGCa                                                <  1:155153/147‑1 (MQ=255)
                                                                                                          gCCGTTCTATGAGCACCCGGCGTACATCGAGGCGCTGGCAGCCAGCGCCCGGCCGTACCTGGAGCAGGGTTACGACCACCTGCTACTGAGCTTCCATGGTTTGCCGGAACGGCACCTGAAGAAGCTGTTTCCGAAGGGGGTCAAGCACGa                                             >  1:3823953/1‑150 (MQ=255)
                                                                                                          gCCGTTCTATGAGCACCCGGCGTACAGCGCGGCGCTGGCCCCCCGCGCCCGGCCGTACCTGGAGCAGGGTTACGACCACCTGCTACTGAGCTTCCATGGTTTGCCGGAACGGCACCTGAAGAAACTGTTTCCGAAGGGGGTCAAGCACGa                                             <  2:697492/150‑1 (MQ=255)
                                                                                                             gTTCTATGAGCACCCGGCGTACATCGAGGCGCTGGCAGCCCGCGCCCGGCCGTACCTGGAGCAGGGTTACGACCACCTGCTACTGAGCTTCCATGGTTTGCCGGAACGGCACCTGAAGAAACTGTTTCCGAAGGGGGTCAAGCACGACCt                                          >  1:5489581/1‑150 (MQ=255)
                                                                                                             gTTCTATGAGCACCCGGCGTACATCGAGGCGCTGGCAGCCAGCGCCCGGCCGTACCTGGAGCAGGGTTACGACCACCTGCTACTGAGCTTCCATGGTTTGCCGGAACGGCACCTGAAGAAACTGTTTCCGAAGGGGGTCAAGCACGACCt                                          >  1:676126/1‑150 (MQ=255)
                                                                                                             gTTCTATGAGCACCCGGCGTACATCGAGGCGCTGGCAGCCAGCGCCCGGCCGTACCTGGAGCAGGGTTACGACCACCTGCTACTGAGCTTCCATGGTTTGCCGGAACGGCACCTGAAGAAACTGTTTCCGAAGGGGGTCAAGCACGACCt                                          >  1:510209/1‑150 (MQ=255)
                                                                                                                 tATGAGCACCCGGCGTACATCGAGGCGCTGGCAGCCAGCGCCCGGCCGTACCTGGAGCAGGGTTACGACCACCTGCTACTGAGCTTCCATGGTTTGCCGGAACGGCACCTGAAGAAACTGTTTCCGAAGGGGGTCAAGCACGACCTgcgc                                      >  2:425959/1‑150 (MQ=255)
                                                                                                                  aTGAGCACCCGGCGTACATCGAGGCGCTGGCAGCCAGCGCCCGGCCGTACCTGGAGCAGGGTTACGACCACCTGCTACTGAGCTTCCATGGTTTGCCGGAACGGCACCTGAAGAAACTGTTTCCGAAGGGGGTCAAGCACGACCTgcgcg                                     >  2:330548/1‑150 (MQ=255)
                                                                                                                   tGAGCACCCGGCGTACATCGAGGCGCTGGCAGCCAGCGCCCGGCCGTACCTGGAGCAGGGTTACGACCACCTGCTACTGAGCTTCCATGGTTTGCCGGAACGGCACCTGAAGAAACTGTTTCCGAAGGGGGTCAAGCACGACCTgcgcgc                                    >  2:4722758/1‑150 (MQ=255)
                                                                                                                    gAGCACCCGGCGTACATCGAGGCGCTGGCAGCCAGCGCCCGGCCGTACCTGGAGCAGGGTTACGACCACCTGCTACTGAGCTTCCATGGTTTGCCGGAACGGCACCTGAAGAAACTGTTTCCGAAGGGGGTCAAGCACGACCTGCGCGCt                                   >  1:4856878/1‑150 (MQ=255)
                                                                                                                     aGCACCCGGCGTACATCGAGGCGCTGGCAGCCAGCGCCCGGCCGTACCTGGAGCAGGGTTACGACCACCTGCTACTGAGCTTCCATGGTTTGCCGGAACGGCACCTGAAGAAACTGTTTCCGAAGGGGGTCAAGCACGACCTGCGCGCTg                                  >  1:784213/1‑150 (MQ=255)
                                                                                                                         cccGGCGTACATCGAGGCGCTGGCAGCCAGCGCCCGGCCGTACCTGGAGCAGGGTTACGACCACCTGCTACTGAGCTTCCATGGTTTGCCGGAACGGCACCTGAAGAAACTGTTTCCGAAGGGGGTCAAGCACGACCTGCGCGCTGCCGa                              >  1:2098618/1‑150 (MQ=255)
                                                                                                                         cccGGCGTACATCGAGGCGCTGGCAGCCAGCGCCCGGCCGTACCTGGAGCAGGGGTACGACCACCTGCTACTGAGCTTCCATGGTTTGCCGGAACGGCACCTGAAGAAACTGTTTCCGAAGGGGGTCAAGCACGACCTGCGCGCTGCCGa                              >  2:2509252/1‑150 (MQ=255)
                                                                                                                          ccGGCGTACATCGCGGCGCTGGCCGCCCGCGCCCGGCCGTACCTGGAGCAGGGTTACGACCACCTGCTACTGAGCTTCCATGGTTTGCCGGAACGGCACCTGAAGAAACTGTTTCCGAAGGGGGTCAAGCACGACCTGCGCGCTGCCGAt                             <  2:5037484/150‑1 (MQ=255)
                                                                                                                               gTACATCGAGGCGCTGGCAGCCAGCGCCCGGCCGTACCTGGAGCAGGGTTACGACCACCTGCTACTGAGCTTCCATGGTTTGCCGGAACGGCACCTGAAGAAACTGTTTCCGAAGGGGGTCAAGCACGACCTGCGCGCTGCCGATtgctg                        >  1:4286414/1‑150 (MQ=255)
                                                                                                                               gTACATCGAGGCGCTGGCAGCCAGCGCCCGGCCGTACCTGGAGCAGGGTTACGACCACCTGCTACTGAGCTTCCATGGTTTGCCGGAACGGCACCTGAAGAAACTGTTTCCGAAGGGGGTCAAGCACGACCTGCGCGCTGCCGATtgctg                        >  1:5203250/1‑150 (MQ=255)
                                                                                                                                   aTCGAGGCGCTGGCAGCCAGCGCCCGGCCGTACCTGGAGCAGGGTTACGACCACCTGCTACTGAGCTTCCATGGTTTGCCGGAACGGCACCTGAAGAAACTGTTTCCGAAGGGGGTCAAGCACGACCTGCGCGCTGCCGATTGCTGCCAc                    >  1:4133551/1‑150 (MQ=255)
                                                                                                                                    ccGAGGCGCCCCCAGCCACCCCCCGGCCGTACCTGGAGCAGGGTTACGACCACCTGCTACTGAGCTTCCATGGTTTGCCGGAACGGCACCTGAAGAAACTGTTTCCGAAGGGGGTCAAGCACGACCTGCGCGCTGCCGATTGCTGCCACg                   <  2:2538951/149‑1 (MQ=255)
                                                                                                                                      gAGGCGCGGGCCGCCCGCCCCCGGCCGTACCTGGAGCAGGGTTACGACCACCTGCTACTGAGCTTCCATGGTTTGCCGGAACGGCACCTGAAGAAACTGTTTCCGAAGGGGGTCAAGCACGACCTGCGCGCTGCCGATTGCTGCCACGGt                 <  2:4431001/150‑1 (MQ=255)
                                                                                                                                         gcgcTGGCTGCCAGCGCCTGGCCGTACCTGGAGCAGGGTTACGACCACCTGCTACTGAGCTTCCATGGTTTGCCGGAACGGCACCTGAAGAAACTGTTTCCGAAGGGGGTCAAGCACGACCTGCGCGCTGCCGATTGCTGCCACGGTGCa              >  1:3304501/1‑150 (MQ=255)
                                                                                                                                          cgcTGGCCGCCCGCGCCCGGCCGTACCTGGGGCAGGGTTACGACCACCTGCTACTGAGCTTCCATGGTTTGCCGGAACGGCACCTGAAGAAACTGTTTCCGAAGGGGGTCAAGCACGACCTGCGCGCTGCCGATTGCTGCCACGGTGCaa             <  2:2573544/150‑1 (MQ=255)
                                                                                                                                           gcTGGCAGCCAGCGCCCGGCCGTACCTGGAGCAGGGTTACGACCACCTGCTACTGAGCTTCCATGGTTTGCCGGAACGGCACCTGAAGAAACTGTTTCCGAAGGGGGTCAAGCACGACCTGCGCGCTGCCGATTGCTGCCACGGTGCAAc            >  1:5164805/1‑150 (MQ=255)
                                                                                                                                              ggCAGCCAGCGCCCGGCCGTACCTGGAGCAGGGTTACGACCACCTGCTACTGAGCTTCCATGGTTTGCCGGAACGGCACCTGAAGAAACTGTTTCCGAAGGGGGTCAAGCACGACCTGCGCGCTGCCGATTGCTGCCACGGTGCAAccgc         <  2:4536999/150‑1 (MQ=255)
                                                                                                                                              ggCAGCCAGCGCCCGGCCGTACCTGGAGCAGGGTTACGACCACCTGCTACTGAGCTTCCATGGTTTGCCGGAACGGCACCTGAAGAAACTGTTTCCGAAGGGGGTCAAGCACGACCTGCGCGCTGCCGATTGCTGCCACGGTGCAAccgc         >  2:5024006/1‑150 (MQ=255)
                                                                                                                                               gCAGCCAGCGCCCGGCCGTACCTGGAGCAGGGTTACGACCACCTGCTACTGAGCTTCCATGGTTTGCCGGAACGGCACCTGAAGAAACTGTTTCCGAAGGGGGTCAAGCACGACCTGCGCGCTGCCGATTGCTGCCACGGTGCAAccgcc        >  1:400253/1‑150 (MQ=255)
                                                                                                                                                 agccagcGCCCGGCCGTACCTGGAGCAGGGTTACGACCACCTGCTACTGAGCTTCCATGGTTTGCCGGAACGGCACCTGAAGAAACTGTTTCCGAAGGGGGTCAAGCACGACCTGCGCGCTGCCGATTGCTGCCACGGTGCAACCGCCGa      >  2:2902383/1‑150 (MQ=255)
                                                                                                                                                 agccaccCCCCGGCCGTACCTGGAGCAGGGTTACGACCACCTGCTACTGAGCTTCCATGGTTTGCCGGAACGGCACCTGAAGAAACTGTTTCCGAAGGGGGTCAAGCACGACCTGCGCGCTGCCGATTGCTGCCACGGTGCAACCGCCGa      <  2:3494864/150‑1 (MQ=255)
                                                                                                                                                  gccagcGCCTCGATGTACCTGGAGCAGGGTTACGACCACCTGCTACTGAGCTTCCATGGTTTGCCGGAACGGCACCTGAAGAAACTGTTTCCGAAGGGGGTCAAGCACGACCTGCGCGCTGCCGATTGCTGCCACGGTGCAACCGCCGaa     <  2:1814141/150‑1 (MQ=255)
                                                                                                                                                     agcGCCCGGCCGTACCTGGAGCAGGGTTACGACCACCTGCTACTGAGCTTCCATGGTTTGCCGGAACGGCACCTGAAGAAACTGTTTCCGAAGGGGGTCAAGCACGACCTGCGCGCTGCCGATTGCTGCCACGGTGCAACCGCCGaagta  >  2:79796/1‑150 (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: ‑