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

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_002947863,3830AC14.3% 578.3 / 24.2 253S172R (AGC→CGC) hemHferrochelatase
Reads supporting (aligned to +/- strand):  ref base A (169/38);  new base C (4/39);  total (173/80)
Fisher's exact test for biased strand distribution p-value = 1.09e-19
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 277 total aligned reads displayed.

CACCGCTTTATCCACAGTTCGCCGACAGCACGGTGACCACTGTGATCGAGCAGGCAAAACAGACGGTCAGCGAGCACCAGTTACCGTTGCAGATGCGCGTGCTGCAGCCGTTCTATGAGCACCCGGCGTACATCGAGGCGCTGGCAGCCAGCGCCCGGCCGTACCTGGAGCAGGGTTACGACCACCTGCTACTGAGCTTCCATGGTTTGCCGGAACGGCACCTGAAGAAACTGTTTCCGAAGGGGGTCAAGCACGACCTGCGCGCTGCCGATTGCTGCCACGGTGCAACCGCCGAAGTA  >  NC_002947/863234‑863532
                                                                                                                                                     |                                                                                                                                                     
cACCGCTTTATCCACAGTTCGCCGACAGCACGGTGACCACTGTGATCGAGCAGGCAAAACAGACGGTCAGCGAGCACCAGTTACCGTTGCAGATGCGCGTGCTGCAGCCGTTCTATGAGCACCCGGCGTACATCGAGGCGCTGGcagcca                                                                                                                                                       >  2:1644794/1‑150 (MQ=255)
 aCCGCTTTATCCACAGTTCGCCGACAGCACGGTGACCACTGTGATCGAGCAGGCAAAACAGACGGTCAGCGAGCACCAGTTACCGTTGCAGATGCGCGTGCTGCAGCCGTTCTATGAGCACCCGGCGTACATCGAGGCGCTGGcagccag                                                                                                                                                      >  1:129229/1‑150 (MQ=255)
  ccGCTTTATCCACAGTTCGCCGACAGCACGGTGACCACTGTGATCGAGCAGGCAAAACAGACGGTCAGCGAGCACCAGTTACCGTTGCAGATGCGCGTGCTGCAGCCGTTCTATGAGCACCCGGCGTACATCGAGGCGCTGGcagccagc                                                                                                                                                     <  1:3769897/150‑1 (MQ=255)
   cGCTTTATCCACAGTTCGCCGACAGCACGGTGACCACTGTGATCGAGCAGGCAAAACAGACGGTCAGCGAGCACCAGTTACCGTTGCAGATGCGCGTGCTGCAGCCGTTCTATGAGCACCCGGCGTACATCGAGGCGCTGGCAGCCAgcg                                                                                                                                                    >  2:3937279/1‑150 (MQ=255)
      tttATCCACAGTTCGCCGACAGCACGGTGACCACTGTGATCGAGCAGGCAAAACAGACGGTCAGCGAGCACCAGTTACCGTTGCAGATGCGCGTGCTGCAGCCGTTCTATGAGCACCCGGCGTACATCGAGGCGCTGGCAGCCAGCGccc                                                                                                                                                 <  2:4169960/150‑1 (MQ=255)
         aTCCACCGTTCGCCGACAGCACGGTGACCACTGTGATCGAGCAGGCAAAACAGACGGTCAGCGAGCACCAGTTCCCGTGGCAGATGCGCGTGCTGCGGCCTTTGTTTGAGCACCCGGCGTACATCGCGGCGCTTGCAGCCTGCGCCCGGc                                                                                                                                              <  2:250536/150‑1 (MQ=255)
          tCCACAGTTCGCCGACAGCACGGTGACCACTGTGATCGAGCAGGCAAAACAGACGGTCAGCGAGCACCAGTTACCGTTGCAGATGCGCGTGCTGCAGCCGTTCTATGAGCACCCGGCGTACATCGAGGCGCTGGCAGCCAGCGCCCGGcc                                                                                                                                             >  2:6305216/1‑150 (MQ=255)
            cacaGTTCGCCGACAGCACGGTGACCACTGTGATCGAGCAGGCAAAACAGACGGTCAGCGAGCACCAGTTACCGTTGCATATGCGCGTGCTGCAGCCGTTCTCTGAGCACCCGGCGTACATCGAGGCGCTGGCAGCCAGCGCCCGGCCGt                                                                                                                                           >  2:5977202/1‑150 (MQ=255)
            cacaGTTCGCCGACAGCACGGTGACCACTGTGATCGAGCAGGCAAAACAGACGGTCAGCGAGCACCAGTTACCGTTGCAGATGCGCGTGCTGCAGCCGTTCTATGAGCACCCGGCGTACATCGAGGCGCTGGCAGCCAGCGCCCGGCCGt                                                                                                                                           >  2:92363/1‑150 (MQ=255)
            cacaGTTCGCCGACAGCACGGTGACCACTGTGATCGAGCAGGCAAAACAGACGGTCAGCGAGCACCAGTTACCGTTGCAGATGCGCGTGCTGCAGCCGTTCTATGAGCACCCGGCGTACATCGAGGCGCTGGCAGCCAGCGCCCGGCCGt                                                                                                                                           >  1:2602748/1‑150 (MQ=255)
            cacaGTTCGCCGACAGCACGGTGACCACTGACATCGAGCAGGCAAAACAGACGGTCAGCGAGCACCAGTTACCGTTGCAGATGCGCGTGCTGCAGCCGTTCTATGAGCACCCGGCGTACATCGAGGCGCTGGCAGCCCGCGCCCGGCCGt                                                                                                                                           >  2:3511362/1‑150 (MQ=255)
             acaGTTCGCCGACAGCACGGTGACCACTGTGATCGAGCAGGCAAAACAGACGGTCAGCGAGCACCAGTTACCGTTGCAGATGCGCGTGCTGCAGCCGTTCTATGAGCACCCGGCGTACATCGAGGCGCTGGCAGCCAGCGCCCGGCCGTa                                                                                                                                          >  2:1392991/1‑150 (MQ=255)
                gTTCGCCGACAGCTCGGTGACCTCTGTGTTCGAGCAGGAAAAACCGACGGTCAGCGAGCACGTGTTCCCGTTGCAGTTGCGCGTGCTGCAGCCGTTCTATGCGCACCGGGCGTACACCGAGGCGCTGGCCGCCAGCGCCCGGCCGTACCt                                                                                                                                       <  2:1490972/150‑1 (MQ=255)
                 ttCGCCGACAGCACGGTGACCACTGTGATCGAGCAGGCAAAACAGACGGTCAGCGAGCACCAGTTACCGTTGCAGATGCGCGTGCTGCAGCCGTTCTATGAGCACCCGGCGTACATCGAGGCGCTGGCAGCCAGCGCCCGGCCGTACCTg                                                                                                                                      >  2:1624278/1‑150 (MQ=255)
                    gCCGACAGCACGGTGACCACTGTGATCGAGCAGGCAAAACAGACGGTCAGCGAGCACCAGTTACCGTTGCAGATGCGCGTGCTGCAGCCGTTCTATGAGCACCCGGCGTACATCGAGGCGCTGGCAGCCAGCGCCCGGCCGTACCTGGAg                                                                                                                                   >  2:6072903/1‑150 (MQ=255)
                      cGACAGCACGGTGACCTCTGTGATCGAGCAGGCAAAACAGACGGTCAGCGAGCACCAGTTACCGTTGCAGATGCGCGTGCTGCAGCCGTTCTATGAGCCCCCGGCGTACATCGCGGCGCTGGCAGCCCGCGCCCGGCCGTACCTGGAGCa                                                                                                                                 <  2:5106742/150‑1 (MQ=255)
                      cGACAGCACGGTGACCACTGTGATCGAGCAGGCAAAACAGACGGTCAGCGAGCACCAGTTACCGTTGCAGATGCGCGTGCTGCAGCCGTTCTATGAGCACCCGGCGTACATCGAGGCGCTGGCAGCCAGCGCCCGGCCGTACCTGGAGCa                                                                                                                                 >  1:4593278/1‑150 (MQ=255)
                       gACAGCACGGTGACCACTGTGATCGAGCAGGCAAAACAGACGGTCAGCGAGCACCAGTTACCGTTGCAGATGCGCGTGCTGCAGCCGTTCTATGAGCACCCGGCGTACATCGAGGCGCTGGCAGCCAGCGCCCGGCCGTACCTGGAGCAg                                                                                                                                >  2:1548782/1‑150 (MQ=255)
                           gCACGGTGACCACTGTGATCGAGCAGGCAAAACAGACGGTCAGCGAGCACCAGTTACCGTTGCAGATGCGCGTGCTGCAGCCGTTCTATGAGCACCCGGCGTACATCGAGGCGCTGGCAGACAGCGCCCGGCCGTACCTGGAGCAGGGtt                                                                                                                            >  1:3087728/1‑150 (MQ=255)
                               ggTGACCACTGTGATCGAGCAGGCAAAACAGACGGTCAGCGAGCACCAGTTACCGTTGCAGATGCGCGTGCTGCAGCCGTTCTATGAGCACCCGGCGTACATCGAGGCGCTGGCAGCCAGCGCCCGGCCGTACCTGGAGCAGGGTTACGa                                                                                                                        >  1:233550/1‑150 (MQ=255)
                                   acccttgtgATCGAGCAGTCAAACCCGCCGGTCCGCGCGCACCTGTTCCGGTTGCTGTGGGGCGTGCTGCCGCCGTTCTTTGAGCCCCGGGCTTACAGCGCGGCGCTGGCCGCCAGCGCCCGGCCGTACCTGGAGCAGGGTTACGaccac                                                                                                                    <  2:670045/145‑1 (MQ=255)
                                     cACTGTGATCGAGCAGGCAAAACAGACGGTCAGCGAGCACCAGTTACCGTTGCAGATGCGCGTGCTGCAGCCGTTCTATGAGCACCCGGCGTACATCGAGGCGCTGGCAGCCCGCGCCCGGCCGTACCTGGAGCAGGGTTACGACCACCt                                                                                                                  >  2:214626/1‑150 (MQ=255)
                                      aCTGTGATCGAGCAGGCAAAACAGACGGTCAGCGAGCACCAGTTACCGTTGCAGATGCGCGTGCTGCAGCCGTTCTATGAGCACCCGGCGTACATCGAGGCGCTGGCAGCCAGCGCCCGGCCGTACCTGGAGCAGGGTTACGACCACCTg                                                                                                                 >  1:5330155/1‑150 (MQ=255)
                                        tgtgATCGAGCAGGCAAAACAGACGGTCAGCGAGCACCAGTTACCGTTGCAGATGCGCGTGCTGCAGCCGTTCTATGAGCACCCGGCGTACATCGAGGCGCTGGCAGCCAGCGCCCGGCCGTACCTGGAGCAGGGTTACGACCACCTGCt                                                                                                               >  2:4988799/1‑150 (MQ=255)
                                        tgtgATCGAGCAGGCAAAACAGACGGTCAGCGAGCACCAGTTACCGTTGCAGATGCGCGTGCTGCAGCCGTTCTATGAGCACCCGGCGTACATCGAGGCGCTGGCAGCCAGCGCCCGGCCGTACCTGGAGCAGGGGTACGACCACCTGCt                                                                                                               >  1:4056441/1‑150 (MQ=255)
                                          tgATCGAGCAGGCAAAACAGACGGTCAGCGAGCACCAGTTACCGTTGCAGATGCGCGTGCTGCAGCCGTTCTATGAGCACCCGGCGTACATCGAGGCGCTGGCAGCCAGCGCCCGGCCGTACCTGGAGCAGGGTTACGACCACCTGCTAc                                                                                                             >  1:3316319/1‑150 (MQ=255)
                                          agATCGAGCAGGCAAAACAGACGGTCAGCGAGCACCAGTTACCGTTGCAGATGCGCGTGCTGCAGCCGTTCTATGAGCACCCGGCGTACATCGAGGCGCTGGCAGCCAGCGCCCGGCCGTACCTGGAGCCGGGTTAAGACCACCAGCTAc                                                                                                             >  1:4939672/2‑150 (MQ=255)
                                             tCGAGCAGGCAAAACAGACGGTCAGCGAGCACCAGTTACCGTTGCAGATGCGCGTGCTGCAGCCGTTCTATGAGCACCCGGCGTACATCGAGGCGCTGGCAGCCAGCGCCCGGCCGTACCTGGAGCAGGGTTACGACCACCTGCTCCTGa                                                                                                          >  1:6744407/1‑150 (MQ=255)
                                             gcGAGCAGGCAACCCCGACGGTCCGCGCCCCCGTGTTCCGGTTGCCGCGGCGGGTGCGGCCGCCGTTCTCTGCGCCCCCGGCGTCCATCGCGGGGCTGGCAGCCAGCCCCCGGCCGTACCTGGAGCAGGGTTACGACCACCTGCTACTGa                                                                                                          <  2:281561/149‑1 (MQ=255)
                                                aGCAGGCAAAACAGACGGTCAGCGAGCACCAGTTACCGTTGCAGATGCGCGTGCTGCAGCCGTTCTATGAGCACCCGGCGTACATCGAGGCGCTGGCAGCCAGCGCCCGGCCGTACCTGGAGCAGGGTTACGACCACCTGCTACTGAGCt                                                                                                       >  2:1220729/1‑150 (MQ=255)
                                                    ggCAAACCCGACGGTCAGCGCGCACCCGCTCCGGTTGCTGCGGCGCGTGCTGCCGCCGTTCTATGCGCACCCGGCGTACCTCGCGGCGGTGGCAGCCAGCGCCCGGCCGTACCTGGAGCAGGGTTACGACCACCTGCTACTGAGCTTCCa                                                                                                   <  2:489967/150‑1 (MQ=255)
                                                     gCAAAACAGACGGTCAGCGAGCACCAGTTACCGTTGCAGATGCGCGTGCTGCAGCCGTTCTATGAGCACCCGGCGTACATCGAGGCGCTGGCAGCCAGCGCCCGGCCGTACCTGGAGCAGGGTTACGACCACCTGCTACTGAGCTTccct                                                                                                  >  1:3120905/1‑148 (MQ=255)
                                                      cAAAACAGACGGTCAGCGAGCACCAGTTACCGTTGCAGATGCGCGTGCTGCAGCCGTTCTATGAGCACCCGGCGTACATCGAGGCGCTGGCCGCCAGCGCCCGGCCGTACCTGGAGCAGGGTTACGACCACCTGCTACTGAGCTTCCATg                                                                                                 >  1:3722084/1‑150 (MQ=255)
                                                       aaaaCAGACGGTCAGCGAGCACCAGTTACCGTTGCAGATGCGCGTGCTGCAGCCGTTCTATGAGCACCCGGCGTACATCGAGGCGCTGGCAGCCAGCGCCCGGCCGTACCTGGAGCAGGGTTACGACCACCTGCTACTGAGCTTCCATgg                                                                                                >  1:5566512/1‑150 (MQ=255)
                                                       aaaaCAGACGGTCAGCGAGCACCAGTTACCGTTGCAGATGCGCGTGCTGCAGCCGTTCTATGAGCACCCGGCGTACATCGAGGCGCTGGCAGCCAGCGCCCGGCCGTACCTGGAGCAGGGTTACGACCACCTGCTACTGAGCTTCCATgg                                                                                                >  1:3580860/1‑150 (MQ=255)
                                                       aaaaCAGACGGTCAGCGAGCACCAGTTACCGTTGCAGATGCGCGTGCTGCAGCCGTTCTATGAGCACCCGGCGTACATCGAGGCGCTGGCAGCCAGCGCCCGGCCGTACCTGGAGCAGGGTTACGACCACCTGCTACTGAGCTTCCATgg                                                                                                >  2:5373512/1‑150 (MQ=255)
                                                       aaaaCAGACGGTCAGCGAGCACCAGTTACCGTTGCAGATGCGCGTGCTGCAGCCGTTCTATGAGCACCCGGCGTACATCGAGGCGCTGGCAGCCAGCGCCCGGCCGTACCTGGAGCAGGGTTACGACCACCTGCTACTGAGCTTCCATgg                                                                                                >  2:3168551/1‑150 (MQ=255)
                                                       aaaaCAGACGGTCAGCGAGCACCAGTTACCGTTGCAGATGCGCGTGCTGCAGCCGTTCTATGAGCACCCGGCGTACATCGAGGCGCTGGCAGCCAGCGCCCGGCCGTACCTGGAGCAGGGTTACGACCACCTGCTACTGAGCTTCCATgg                                                                                                >  2:5664514/1‑150 (MQ=255)
                                                       aaaaCAGACGGTCAGCGAGCACCAGTTACCGTTGCAGATGCGCGTGCTGCAGCCGTTCTATGAGCACCCGGCGTACATCGAGGCGCTGGCAGCCAGCGCCCGGCCGTACCTGGAGCAGGGTTACGACCACCTGCTACTGAGCTTCCATgg                                                                                                >  2:4665314/1‑150 (MQ=255)
                                                            aGACGGTCAGCGAGCACCAGTTACCGTTGCAGATGCGCGTGCTGCAGCCGTTCTATGAGCACCCGGCGTACATCGAGGCGCTGGCAGCCAGCGCCCGGCCGTACCTGGAGCAGGGTTACGACCACCAGCTACTGAGCTTCCATGGTTTGc                                                                                           >  1:2851551/1‑150 (MQ=255)
                                                             gcgggTCCGCGAGCACCTTTTCCGTTTGCTGTGGGGCGGGGTGCCGCCATTCTTTGCGCCCGCGGCTTACATGGCGGGGCTGGCCGCCCGCGCCCGGCCGTACCTGGAGCAGGGTTACGACCACCTGCTACTGAGCTTCCATGGTTTGcc                                                                                          <  2:1834288/147‑1 (MQ=255)
                                                                 gTCAGCGAGCACCTGTTACGGTTGCTGTTGGGCGTGCTGCCGCCTTTATTTGAGCACCCGGCTTACATCGAGGCGCTGGCAGCCAGCGCCCGGCCGTACCTGGAGCAGGGTTACGACCACCTGCTACTGAGCTTCCATGGTTTGCCGGaa                                                                                      <  2:2945534/150‑1 (MQ=255)
                                                                    aGCGAGCACCAGTTACCGTTGCAGATGCGCGTGCTGCAGCCGTTCTATGAGCACCCGGCGTACATCGAGGCGCTGGCAGCCAGCGCCCGGCCGTACCTGGAGCAGGGTTACGACCACCTGCTACTGAGCTTCCATGGTTTGCCGGAACgg                                                                                   >  1:1411836/1‑150 (MQ=255)
                                                                       gAGCACCAGTTACCGTTGCAGATGCGCGTGCTGCAGCCGTTCTATGAGCACCCGGCGTACATCGAGGCGCTGGCAGCCAGCGCCCGGCCGTACCTGGAGCAGGGTTACGACCACCTGCTACTGAGCTTCCATGGTTTGCCGGAACGGCAc                                                                                >  1:5179399/1‑150 (MQ=255)
                                                                       gAGCACCAGTTACCGTTGCAGATGCGCGTGCTGCAGCCGTTCTATGAGCACCCGGCGTACATCGAGGCGCTGGCAGCCAGCGCCCGGCCGTACCTGGAGCAGGGTTACGACCACCTGCTACTGAGCTTCCATGGTTTGCCGGAACGGCAc                                                                                >  2:6570485/1‑150 (MQ=255)
                                                                          ccccAGTTCCCTTTGCAGAGGCGCGGGCGGCAGCGGTTCTATGAGCCCCCGGCGTACATCGAGGCGCGGGCAGCCCGCGCCCGGCCGTACCTGGAGCAGGGTTACGACCACCTGCTACTGAGCTTCCATGGTTTGCCGGAACGGCACCTg                                                                             <  2:2656457/148‑1 (MQ=255)
                                                                          cACCAGTGCCCGTTGCCGCGGCGCGTGCTGCAGCCGGTCGCTCCGCCCCCGGCGCCCCGCGCGGCGCGGGCACCCCGCCCCCGGCCGTACCTGGAGCAGGGTTACGACCACCTGCTACTGAGCTTCCATGGTTTGCCGGAACGGCACCTg                                                                             <  1:5200883/150‑1 (MQ=255)
                                                                            cccgTTCCCGTTGCTGATGCGCGGGCTCCAGCCTTTTTATGCGCCCCCGGCGTACATCGCGGCGCGGGCCCCCAGCCCCCGGCCGTACCTGGAGCAGGGTTACGACCACCTGCTACTGAGCTTCCATGGTTTGCCGGAACGGCACCTgaa                                                                           <  2:2211227/147‑1 (MQ=255)
                                                                             cAGTTACCGTTGCAGATGCGCGTGCTGCAGCCGTTCTATGAGCACCCGGCGTACATCGAGGCGCTGGCAGCCAGCGCCCGGCCGTACCTGGAGCAGGGTTACGACCACCTGCTACTGAGCTTCCATGGTTTGCCGGAACGGCACCTgaag                                                                          >  1:5948125/1‑150 (MQ=255)
                                                                             cAGTTACCGTTGCAGATGCGCGTGCTGCAGCCGTTCTATGAGCACCCGGCGTACATCGAGGCGCTGGCAGCCAGCGCCCGGCCGTACCTGGAGCAGGGTTACGACCACCTGCTACTGAGCTTCCATGGTTTGCCGGAACGGCACCTgaag                                                                          >  2:2558051/1‑150 (MQ=255)
                                                                             cAGTTACCGTTGCAGATGCGCGTGCTGCAGCCGTTCTATGAGCACCCGGCGTACATCGAGGCGCTGGCAGCCAGCGCCCGGCCGTACCTGGAGCAGGGTTACGACCACCTGCTACTGAGCTTCCATGGTTTGCCGGAACGGCACCTgaag                                                                          >  2:3435033/1‑150 (MQ=255)
                                                                                gtcccGTTGCCGTGGGGCGTGCTGCAGCCGTTCTAGGCGCCCGCGGCGTCCCTCGCGGCGCGGGCAGCCCGCGCCCGGCCGTACCTGGAGCAGGGTTACGACCACCTGCTACTGAGCTTCCATGGTTTGCCGGAACGGCACCTGAAGaaa                                                                       <  1:1072055/147‑1 (MQ=255)
                                                                                  ccggTTGCTGCTGCGCGTGCTGCCGCTGTTCTATGCGCCCCCGGCGCACCTCGCGGCGCGGGCCCCCCGCCCCCGGCCGTACCTGGAGCAGGGTTACGACCACCTGCTACTGAGCTTCCATGGTTTGCCGGAACGGCACCTGAAGAAACt                                                                     <  1:92363/147‑1 (MQ=255)
                                                                                  aCCGTTGCCGATGCGCGTGCTGCAGCCGATCTAGGAGCACCCGGCGTACATCGCGGCGCTGGCAGCCCGCGCCCGGCCGTACCTGGAGCAGGGTTACGACCACCTGCTACTGAGCTTCCATGGTTTGCCGGAACGGCACCTGAAGAAACt                                                                     <  1:5709199/150‑1 (MQ=255)
                                                                                     gTTGCAGATGCGCGTGCTGCAGCCGTTCTATGAGCACCCGGCGTACATCGAGGCGCTGGCAGCCAGCGCCCGGCCGTACCTGGAGCAGGGTTACGACCACCTGCTACTGAGCTTCCATGGTTTGCCGGAACGGCACCTGAAGAAACTGtt                                                                  >  1:6080677/1‑150 (MQ=255)
                                                                                     gTTGCAGAGGCGGGTGCTGCAGCCGTTCTATGACCCCCCGGCGTACATCGTGGCGCTGGCAGCCAGCGCCCGGCCGTACCTGGAGCAGGGTTACGACCACCTGCTACTGAGCTTCCATGGTTTGCCGGAACGGCACCTGAAGAAACTGtt                                                                  <  2:5640941/150‑1 (MQ=255)
                                                                                      ttGCAGATGCGCGTGCTGCAGCCGTTCTATGAGCACCCGGCGTACATCGAGGCGCTGGCAGCCAGCGCCCGGCCGTACCTGGAGCAGGGTTACGACCACCTGCTACTGAGCTTCCATGGTTTGCCGGAACGGCACCAGAAGAAACTGttt                                                                 >  1:216658/1‑150 (MQ=255)
                                                                                        gCAGATGCGCGTGCTGCAGCCGTTCTATGAGCACCCGGCGTACATCGAGGCGCTGGCAGCCAGCGCCCGGCCGTACCTGGAGCAGGGTTACGACCACCTGCTACTGAGCTTCCATGGTTTGCCGGAACGGCACCTGAAGAAACTGTTTcc                                                               >  1:6404119/1‑150 (MQ=255)
                                                                                         cAGATGCGCGTGCTGCAGCCGTTCTATGAGCACCCGGCGTACATCGAGGCGCTGGCAGCCAGCGCCCGGCCGTACCTGGAGCAGGGTTACGACCACCTGCTACTGAGCTTCCATGGTTTGCCGGAACGGCACCTGAAGAAACTGTTTCCg                                                              >  1:231841/1‑150 (MQ=255)
                                                                                          aGATGCGCGTGCTGCAGCCGTTCTATGAGCACCCGGCGTACATCGAGGCGCTGGCAGCCAGCGCCCGGCCGTACCTGGAGCAGGGTTACGACCACCTGCTACTGAGCTTCCATGGTTTGCCGGAACGGCACCTGAAGAAACTGTTTCCGa                                                             >  2:6140806/1‑150 (MQ=255)
                                                                                           gctGCGCGTGCTCCCGCGGTGAAATGCGCCCCGGGCGTACAGCGCGGCGCGGCCCGCCCGCGCCCGGCCGTACCTGGAGCAGGGTTACGACCACCTGCTACTGAGCTTCCATGGTTTGCCGGAACGGCACCTGAAGAAACTGTTTCCGaa                                                            <  1:4259106/148‑1 (MQ=255)
                                                                                            gggggcgTGCTGCCGCCGTTGTTTGAGCCCCGGGCGTACTTCGAGGGGGGGGCAGCCCGCGCCCGGCCGTACCTGGAGCAGGGTTACGACCACCTGCTACTGAGCTTCCATGGTTTGCCGGAACGGCACCTGAAGAAACTGTTTCCGAAg                                                           <  2:1313909/146‑1 (MQ=255)
                                                                                             tgggggTGCTGCCGCGGTTCTCTGCGCCCCCGGCGTACCGCGCGGCGCCCCCCCCCCGCCCCCGGCCGTACCTGGAGCAGGGTTACGACCACCTGCTACTGAGCTTCCATGGTTTGCCGGAACGGCACCTGAAGAAACTGTTTCCGAAgg                                                          <  1:176549/145‑1 (MQ=255)
                                                                                              gcgcGTGCTGCAGCCGTTCTATGAGCACCCGGCGTACATCGAGGCGCTGGCAGCCAGCGCCCGGCCGTACCTGGAGCAGGGTTACGACCACCTGCTACTGAGCTTCCATGGTTTGCCGGAACGGCACCTGAAGAAACTGTTTCCGAAggg                                                         >  1:2502265/1‑150 (MQ=255)
                                                                                               ggcgTGCTGCAGCGTTTCTATGAGCCCCCGGCGTACATCGTGGCGGTGGCAGCCAGCGCCCGGCCGTACCTGGAGCAGGGTTACGACCACCTGCTACTGAGCTTCCATGGTTTGCCGGAACGGCACCTGAAGAAACTGTTTCCGAAgggg                                                        <  2:4636308/149‑1 (MQ=255)
                                                                                                gcgTGCTGCAGCCGTTCTATGAGCACCCGGCGTACATCGAGGCGCTGGCAGCCAGCGCCCGGCCGTACCTGGGGCAGGGTTACGACCACCTGCTACTGAGCTTCCATGGTTTGCCGGAACGGCTCCGGATGAAACTGTTTCCGAAggggg                                                       >  1:3003602/1‑150 (MQ=255)
                                                                                                gcgTGCTGCAGCCGTTCTATGAGCACCCGGCGTACATCGAGGCGCTGGCAGCCAGCGCCCGGCCGTACCTGGAGCAGGGTTACGACCACCTGCTACTGAGCTTCCATGGTTTGCCGGAACGGCACCTGAAGAAACTGTTTCCGAAggggg                                                       >  2:5434560/1‑150 (MQ=255)
                                                                                                    gctgctgctTTTCTATGAGCACCCGGCTTACCTGGTGGCGCGGGCAGCCAGCCCCCGGCCGTACCTGGAGCAGGGTTACGACCACCTGCTACTGAGCTTCCATGGTTTGCCGGAACGGCACCTGAAGAAACTGTTTCCGAAGGGGGTCaa                                                   <  1:4383441/150‑1 (MQ=255)
                                                                                                      tgcAGCCGTTCTATGAGCACCCGGCGTACATCGAGGCGCTGGCAGCCAGCGCCCGGCCGTACCTGGAGCAGGGTTACGACCACCTGCTACTGAGCTTCCATGGTTTGCCGGAACGGCACCTGAAGAAACTGTTTCCGAAGGGGGTCAAGc                                                 >  1:6489602/1‑150 (MQ=255)
                                                                                                      tgcAGCCGTTCTATGAGCACCCGGCGTACATCGAGGCGCTGGCAGCCAGCGCCCGGCCGTACCTGGAGCAGGGTTACGACCACCTGCTACTGAGCTTCCATGGTTTGCCGGAACGGCACCTGAAGAAACTGTTTCCGAAGGGGGTCAAGc                                                 >  2:2034768/1‑150 (MQ=255)
                                                                                                           ctgTTCTATGCGCCCCCGGCGTACCTCGCGGCGCTGGCAGCCAGCGCCCGGCCGTACCTGGAGCAGGGTTACGACCACCTGCTACTGAGCTTCCATGGTTTGCCGGAACGGCACCTGAAGAAACTGTTTCCGAAGGGGGTCAAGCACGAc                                            <  2:4056441/148‑1 (MQ=255)
                                                                                                           ccGTTCTATGCGCCCCCGGCGTCCCTCGCGGCGCGGGCAGCCCGCCCCCGGCCGTACCTGGAGCAGGGTTACGACCACCTGCTACTGAGCTTCCATGGTTTGCCGGAACGGCACCTGAAGAAACTGTTTCCGAAGGGGGTCAAGCACGAc                                            <  1:3966259/150‑1 (MQ=255)
                                                                                                               tCTTTGACCACCCGGCGAACATCGCGGCGCGGCCCCCCCGCGCCCGGCCGTACCTGGAGCAGGGTTACGACCACCTGCTACTGAGCTTCCATGGTTTGCCGGAACGGCACCTGAAGAAACTGTTTCCGAAGGGGGTCAAGCACGACCTgc                                        <  2:6501415/150‑1 (MQ=255)
                                                                                                                 tATGAGCACCCGGCGTACATCGAGGCGCTGGCAGCCAGCGCCCGGCCGTACCTGGAGCAGGGTTACGACCACCTGCTACTGAGCTTCCATGGTTTGCCGGAACGGCACCTGAAGAAACTGTTTCCGAAGGGGGTCAAGCACGACCTgcgc                                      >  2:5859670/1‑150 (MQ=255)
                                                                                                                 tATGAGCACCCGGCGTACATCGAGGCGCTGGCAGCCAGCGCCCGGCCGTACCTGGAGCAGGGTTACGACCACCTGCTACTGAGCTTCCATGGTTTGCCGGAACGGCACCTGAAGAAACTGTTTCCGAAGGGGGTCAAGCACGACCTgcgc                                      >  2:718485/1‑150 (MQ=255)
                                                                                                                        aCCCGGCGTACATCGAGGCGCTGGCAGCCAGCGCCCGGCCGTACCTGGAGCAGGGTTACGACAACCTGCTACTGAGCTTCCATGGTTTGCCGGAACGGCACCTGAAGAAACTGTTTCCGAAGGGGGTCAAGCACGACCTGCGCGCTGCCg                               >  1:5458576/1‑150 (MQ=255)
                                                                                                                          ccGGCTTACCTCGCGGCGCTGGCCGCCCGCGCCCGGCCGTACCTGGAGCAGGGTTACGACCACCTGCTACTGAGCTTCCATGGTTTGCCGGAACGGCACCTGAAGAAACTGTTTCCGAAGGGGGTCAAGCACGACCTGCGCGCTGCCGAt                             <  2:5218002/150‑1 (MQ=255)
                                                                                                                           cGGCGTACATCGAGGCGCTGGCAGCCAGCGCCCGGCCGTACCTGGAGCAGGGTTACGACCACCTGCTACTGAGCTTCCATGGTTTGCCGGAACGGCACCTGAAGAAACTGTTTCCGAAGGGGGTCAAGCACGACCTGCGCGCTGCCGAtt                            >  1:5071334/1‑150 (MQ=255)
                                                                                                                             gcttACATCGCGGCGGTGGCAGCCCGCGCCCGGCCGTACCTGGAGCAGGGTTACGACCACCTGCTACTGAGCTTCCATGGTTTGCCGGAACGGCACCTGAAGAAACTGTTTCCGAAGGGGGTCAAGCACGACCTGCGCGCTGCCGATtgc                          <  2:454167/147‑1 (MQ=255)
                                                                                                                             gCGTACATCGCGGCGCTGGCCGCCCGCGCCCGGCCGTACCTGGAGCAGGGTTACGACCACCTGCTACTGAGCTTCCATGGTTTGCCGGAACGGCACCTGAAGAAACTGTTTCCGAAGGGGGTCAAGCACGACCTGCGCGCTGCCGATtgc                          <  2:2882411/150‑1 (MQ=255)
                                                                                                                               gTACATCTAGGCGCTGGCAGCCAGCGCCCGGCCGTACCTGGAGCAGGGTTACGACCACCTGCTACTGAGCTTCCATGGTTTGCCGGAACGGCACCTGAAGAAACTGTTTCCGCGGGGGGTCAAGCACGACCTGCGCGCTGCCGCTtgctg                        >  1:5028251/1‑150 (MQ=255)
                                                                                                                               gTACATCTAGGCGCTGGCAGCCAGCGCCCGGCCGTACCTGGAGCAGGGTTACGACCACCTGCTACTGAGCTTCCATGGTTTGCCGGAACGGCACCTGAAGAAACTGTTTCCGAAGGGGGTCAAGCACGACCTGCGCGCTGCCGATtgctg                        >  1:2417931/1‑150 (MQ=255)
                                                                                                                               gTACATCGAGGCGCTGGCAGCCAGCGCCCGGCCGTACCTGGAGCAGGGTTACGACCACCTGCTACTGAGCTTCCATGGTTTGCCGGAACGGCACCTGAAGAAACTGTTTCCGAAGGGGGTCAAGCACGACCTGCGCGCTGCCGATcgctg                        >  1:653591/1‑150 (MQ=255)
                                                                                                                                tACATCGAGGCGCTGGCAGCCAGCGCCCGGCCGTACCTGGAGCAGGGTTACGACCACCTGCTACTGAGCTTCCATGGTTTGCCGGAACGGCACCTGAAGAAACTGTTTCCGAAGGGGGTCAAGCACGACCTGCGCGCTGCCGATtgctgc                       >  2:4245079/1‑150 (MQ=255)
                                                                                                                                 aCATCGAGGCGCTGGCAGCCAGCGCCCGGCCGTACCTGGAGCAGGGTTACGACCACCTGCTACTGAGCTTCCATGGTTTGCCGGAACGGCACCTGAAGAAACTGTTTCCGAAGGGGGTCAAGCACGACCTGCGCGCTGCCGATTGCTGcc                      >  2:2986498/1‑150 (MQ=255)
                                                                                                                                   ctCGTGGCGCTGGCAGCCAGCGCCCGGCCGTACCTGGAGCAGGGTTACGACCACCTGCTACTGAGCTTCCATGGTTTGCCGGAACGGCACCTGAAGAAACTGTTTCCGAAGGGGGTCAAGCACGACCTGCGCGCTGCCGATTGCTGCCAc                    <  1:6214528/149‑1 (MQ=255)
                                                                                                                                   agcgtggCGCTGGCCGCCCGCGCCCGGCCGTACCTGGAGCAGGGTTACGACCACCTGCTACTGAGCTTCCATGGTTTGCCGGAACGGCACCTGAAGAAACTGTTTCCGAAGGGGGTCAAGCACGACCTGCGCGCTGCCGATTGCTGCCAc                    <  2:3709604/145‑1 (MQ=255)
                                                                                                                                   aTCGAGGCGCTGGCAGCCAGCGCCCGGCCGTACCTGGAGCAGGGTTACGACCACCTGCTACTGAGCTTCCATGGTTTGCCGGAACGGCACCTGAAGAAACTGTTTCCGAAGGGGGTCAAGCACGACCTGCGCGCTGCCGATTGCTGCCAc                    >  1:1946982/1‑150 (MQ=255)
                                                                                                                                   aTCGAGGCGCTGGCAGCCAGCGCCCGGCCGTACCTGGAGCAGGGTTACGACCACCTGCTACTGAGCTTCCATGGTTTGCCGGAACGGCACCTGAAGAAACTGTTTCCGAAGGGGGTCAAGCACGACCTGCGCGCTGCCGATTGCTGCCAc                    >  1:3825499/1‑150 (MQ=255)
                                                                                                                                   aTCGAGGCGCTGGCAGCCAGCGCCCGGCCGTACCTGGAGCAGGGTTACGACCACCTGCTACTGAGCTTCCATGGTTTGCCGGAACGGCACCTGAAGAAACTGTTTCCGAAGGGGGTCAAGCACGACCTGCGCGCTGCCGATTGCTGCCAc                    >  2:207436/1‑150 (MQ=255)
                                                                                                                                    tCGAGGCGCTGGCAGCCAGCGCCCGGCCGAACCTGGAGCAGGGTTACGACCACCTGCTACTGAGCTTCCATGGTTTGCCGGAACGGCACCTGAAGAAACTGTTTCCGAAGGGGGTCAAGCACGACCTGCGCGCTGCCGATTGCTGCCACg                   <  2:6449334/150‑1 (MQ=255)
                                                                                                                                      gAGGCGCTGGCAGCCAGCGCCCGGCCGTACCTGGAGCAGGGTTACGACCACCTGCTACTGAGCTTCCATGGTTTGCCGGAACGGCACCTGAAGAAACTGTTTCCGAAGGGGGTCAAGCACGACCTGCGCGCTGCCGATTGCTGCCACGGt                 >  1:2772955/1‑150 (MQ=255)
                                                                                                                                         gcgcTGGCAGCCAGCGCCCGGCCGTACCTGGAGCAGGGTTACGACCACCTGCTACTGAGCTTCCATGGTTTGCCGGAACGGCACCTGAAGAAACTGTTTCCGAAGGGGGTCAAGCACGACCTGCGCGCTGCCGATTGCTGCCACGGTGCa              >  1:1647940/1‑150 (MQ=255)
                                                                                                                                           gcTGGCAGCCAGCGCCCGGCCGTACCTGGAGCAGGGTTACGACCACCTGCTACTGAGCTTCCATGGTTTGCCGGAACGGCACCTGAAGAAACTGTTTCCGAAGGGGGTCAAGCACGACCTGCGCGCTGCCGATTGCTGCCACGGTGCAAc            >  1:2120658/1‑150 (MQ=255)
                                                                                                                                             tGGCAGCCAGCGCCCGGCCGTACCTGGAGCAGGGTTACGACCACCTGCTACTGAGCTTCCATGGGTTGCCGGAACGGCACCTGAAGAAACTGTTTCCGAAGGGGGTCAAGCACGACCTGCGCGCTGCCGATTGCTGCCACGGTGCACccg          >  2:4432377/1‑150 (MQ=255)
                                                                                                                                              ccccgcccgcccCCGGCCGTACCTGGAGCAGGGTTACGACCACCTGCTACTGAGCTTCCATGGTTTGCCGGAACGGCACCTGAAGAAACTGTTTCCGAAGGGGGTCAAGCACGACCTGCGCGCTGCCGATTGCTGCCACGGTGCAAccgc         <  1:3900957/146‑1 (MQ=255)
                                                                                                                                               gCAGCCAGCGCCCGGCCGTACCTGGAGCAGGGTTACGACCACCTGCTACTGAGCTTCCATGGTTTGCCGGAACGGCACCTGAAGAAACTGTTTCCGAAGGGGGTCAAGCACGACCTGCGCGCTGCCGATTGCTGCCACGGTGCAAccgcc        >  1:1513629/1‑150 (MQ=255)
                                                                                                                                                 cgcccgcGCCCGGCCGTACCTGGAGCAGGGTTACGACCACCTGCTACTGAGCTTCCATGGTTTGCCGGAACGGCACCTGAAGAAACTGTTTCCGAAGGGGGTCAAGCACGACCTGCGCGCTGCCGATTGCTGCCACGGTGCAACCGCCGa      <  2:6089983/149‑1 (MQ=255)
                                                                                                                                                 agccagcGCCCGGCCGTACCTGGAGCAGGGTTACGACCACCTGCTACTGAGCTTCCATGGTTTGCCGGAACGGCACCTGAAGAAACTGTTTCCGAAGGGGGTCAAGCACGACCTGCGCGCTGCCGATTGCTGCCACGGTGCAACCGCCGa      >  2:1765582/1‑150 (MQ=255)
                                                                                                                                                    cagcGCCCGGCCGTACCTGGAGCAGGGTTACGACCACCTGCTACTGAGCTTCCATGGTTTGCCGGAACGGCACCTGAAGAAACTGTTTCCGAAGGGGGTCAAGCACGACCTGCGCGCTGCCGATTGCTGCCACGGTGCAACCGCCGaagt   <  2:2345771/150‑1 (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: ‑