Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 88,653 T→A 5.1% Q222L (CAG→CTG)  betC ← choline‑sulfatase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_00294788,6530TA 5.1% 278.5 / 4.6 92Q222L (CAG→CTG) betCcholine‑sulfatase
Reads supporting (aligned to +/- strand):  ref base T (38/49);  new base A (2/3);  total (40/52)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.45e-01

TCTGCAGCAATTGGCCGATGTTGTCATCGATATAGCTGCACGCGCCGAAGTAGGCGCGGCGGGCATCGCGAACCTTATCCACAGGCAGTGGCTTGTTCCAAAGGTCGTAGACCTTCAGCAGGCGCTGTGAATGCGGGTCGAGTTCTGCCTGGCTGAACTCGGCACGGGGCATGGGGATATCCACACCTTCGTAGCGGTCCCAGTAGCGTTTGGCAATGGTGTAGGGGTCATGCGGGTGGGTCATCGAAACGGTCAGGCAAAATGGCCGGCCATCGTTTTCGCGCACATGGTCGTACAGG  >  NC_002947/88504‑88802
                                                                                                                                                     |                                                                                                                                                     
tCTGCAGCAATTGGCCGATGTTGTCATCGATATAGCTGCACGCGCCGAAGTAGGCGCGGCGGGCATCGCGAACCTTATCCACAGGCAGTGGCTTGTTCCAAAGGTCGTAGACCTTCAGCAGGCGCTGTGAATGCGGGTCGAGTTCTGCCt                                                                                                                                                       >  1:707043/1‑150 (MQ=255)
tCTGCAGCAATTGGCCGATGTTGTCATCGATATAGCTGCACGCGCCGAAGTAGGCGCGGCGGGCATCGCGAACCTTATCCACAGGCAGTGGCTTGTTCCAAAGGTCGTAGACCTTCAGCAGGCGCTGTGAATGCGGGTCGAGTTCTGCCt                                                                                                                                                       >  1:1833821/1‑150 (MQ=255)
    cagcaATTGGCCGATGTTGTCATCGATATAGCTGCACGCGCCGAAGTAGGCGCGGCGGGCATCGCGAACCTTATCCACAGGCAGTGGCTTGTTCCAAAGGTCGTATACCTTCAGCAGGCGCTGTGAATGCGGGTCGAGTTCTGCCTGGCt                                                                                                                                                   >  2:576731/1‑150 (MQ=255)
         aTTGGCCGATGTTGTCATCGATATAGCTGCACGCGCCGAAGTAGGCGCGGCGGGCATCGCGAACCTTATCCACAGGCAGTGGCTTGTTCCAAAGGTCGTAGACCTTCAGCAGGCGCTGTGAATGCGGGTCGAGTTCAGCCAGGCAGAACt                                                                                                                                              <  2:328105/150‑1 (MQ=255)
          ttGGCCGATGTTGTCATCGATATAGCTGCACGCGCCGAAGTAGGCGCGGCGGGCATCGCGAACCTTATCCACAGGCAGTGGCTTGTTCCAAAGGTCGTAGACCTTCAGCAGGCGCTGTGAATGCGGGTCGAGTTCTGCCTGGCTGAACTc                                                                                                                                             >  1:1848767/1‑150 (MQ=255)
               cGATGTTGTCATCGATATAGCTGCACGCGCCGAAGTAGGCGCGGCGGGCATCGCGAACCTTATCCACAGGCAGTGGCTTGTTCCAAAGGTCGTAGACCTTCAGCAGGCGCTGTGAATGCGGGTCGAGTTCTGCCTGGCTGAACTCGGCAc                                                                                                                                        >  1:662611/1‑150 (MQ=255)
                 aTGTTGTCATCGATATAGGTGGACGCGCCTAAGTAGGCGCGGCGGGCATCGCGAACCTTATCCACAGGCAGTGGCTTGTTCCAAAGGTCGTAGACCTTCAGCAGGCGCTGTGAATGCGGGTCGAGTTCTGCCTGGCTGAACTCGGCACgg                                                                                                                                      <  2:1885075/150‑1 (MQ=255)
                  tgttgtCATCGATATAGCTGCACGCGCCGAAGTAGGCGCGGCGGGCATCGCGAACCTTATCCACAGGCAGTGGCTTGTTCCAAAGGTCGTAGACCTTCAGCAGGCGCTGTGAATGCGGGTCGAGTTCTGCCTGGCCGAACTCGGCACggg                                                                                                                                     >  1:1038748/1‑150 (MQ=255)
                       tCATCGATATAGCTGCACGCGCCGAAGTAGGCGCGGCGGGCATCGCGAACCTTATCCACAGGCAGTGGCTTGTTCCAAAGGTCGTAGACCTTCAGCAGGCGCTGTGAATGCGGGTCGAGTTCTGCCTGGCTGAACTCGGCACGGGGCATg                                                                                                                                >  2:826451/1‑150 (MQ=255)
                       tCATCGATATAGCTGCACGCGCCGAAGTAGGCGCGGCGGGCATCGCGAACCTTATCCACAGGCAGTGGCTTGTTCAAAAGGTCGTAGACCTTCAGCAGGCGCTGTGAATGCGGGTCGAGTTCTGCCTGGCTGAACTCGGCACGGGGCATg                                                                                                                                >  2:1207531/1‑150 (MQ=255)
                         aTCGATATAGCTGCACGCGCCGAAGTAGGCGCGGCGGGCATCGCGAACCTTATCCACAGGCAGTGGCTTGTTCCAAAGGTCGTAGACCTTCAGCAGGCGCTGTGAATGCGGGTCGAGTTCTGCCTGGCTGAACTCTGCACGGGGCATggg                                                                                                                              >  1:650822/1‑150 (MQ=255)
                         aTCGATATAGCTGCACGCGCCGAAGTAGGCGCGGCGGGCATCGCGAACCTTATCCACAGGCAGTGGCTTGTTCCAAAGGTCGTAGACCTTCAGCAGGCGCTGTGAATGCGGGTCGAGTTCTGCCTGGCTGAACTCGGCACGGGGCATggg                                                                                                                              >  1:36015/1‑150 (MQ=255)
                         aTCGATATAGCTGCACGCGCCGAAGTAGGCGCGGCGGGCATCGCGAACCTTATCCACAGGCAGTGGCTTGTTCCAAAGGTCGTAGACCTTCAGCAGGCGCTGTGAATGCGGGTCGAGTTCTGCCTGGCTGAACTCGGCACGGGGCATggg                                                                                                                              >  2:108365/1‑150 (MQ=255)
                          tCGATATAGCTGCACGCGCCGAAGTAGGCGCGGCGGGCATCGCGAACCTTATCCACAGGCAGTGGCTTGTTCCAAAGGTCGTAGACCTTCAGCAGGCGCTGTGAATGCGGGTCGAGTTCTGCCTGGCTGAACTCGGCACGGGGCATgggg                                                                                                                             <  1:893837/150‑1 (MQ=255)
                              tataGCTGCACGCGCCGAAGTAGGCGCGGCGGGCATCGCGAACCTTATCCACAGGCAGTGGCTTGTTCCAAAGGTCGTAGACCTTCAGCAGGCGCTGTGAATGCGGGTCGAGTTCTGCCTGGCTGAACTCGGCACGGGGCATGGGGatct                                                                                                                         <  2:976792/150‑3 (MQ=255)
                              tataGCTGCACGCGCCGAAGTAGGCGCGGCGGGCATCGCGAACCTTATCCACAGGCAGTGGCTTGTTCCAAAGGTCGTAGACCTTCAGCAGGCGCTGTGAATGCGGGTCGAGTTCTGCCTGGCTGAACTCGGCACGGGGCATGGGGatat                                                                                                                         <  2:1833821/150‑1 (MQ=255)
                               ataGCTGCACGCGCCGAAGTAGGCGCGGCGGGCATCGCGAACCTTATCCACCGGCAGTGGCTTGTTCCAAAGGTCGTAGACCTTCAGCAGGCGCTGTGAATCCGGGTCGAGTTCTGCCTGGCTGAACTCGGCACGGGGCATGGGGATATc                                                                                                                        <  2:1797391/150‑1 (MQ=255)
                               ataGCTGCACGCGCCGAAGTAGGCGCGGCGGGCATCGCGAACCTTATCCACAGGCAGTGGCTTGTTCCAAAGGTCGTAGACCTTCAGCAGGCGCTGTGAATGCGGGTCGAGTTCTGCCTGGCTGAACTCGGCACGGGGCATGGGGATATc                                                                                                                        >  1:1656659/1‑150 (MQ=255)
                                 aGCTGCACGCGCCGAAGTAGGCGCGGCGGGCATCGCGAACCTTATCCACAGGCAGTGGCTTGTTCCAAAGGTCGTAGACCTTCAGCAGGCGCTGTGAATGCGGGTCGAGTTCTGCCTGGCTGAACTCGGCACGGGGCATGGGGATATCca                                                                                                                      <  1:1302573/150‑1 (MQ=255)
                                 aGCTGCACGCGCCGAAGTAGGCGCGGCGGGCATCGCGAACCTTATCCACAGGCAGTGGCTTGTTCCAAAGGTCGTAGACCTTCAGCAGGCGCTGTGAATGCGGGTCGAGTTCTGCCTGGCTGAACTCGGCACGGGGCATGGGGATATCca                                                                                                                      >  1:515731/1‑150 (MQ=255)
                                  gCTGCACGCGCCGAAGTAGGCGCGGCGGGCATCGCGAACCTTATCCACAGGCAGTGGCTTGTTCCAAAGGTCGTAGACCTTCAGCAGGCGCTGTGAATGCGGGTCGAGTTCTGCCTGGc                                                                                                                                                    <  2:1425526/119‑1 (MQ=255)
                                  gCTGCACGCGCCGAAGTAGGCGCGGCGGGCATCGCGAACCTTATCCACAGGCAGTGGCTTGTTCCAAAGGTCGTAGACCTTCAGCAGGCGCTGTGAATGCGGGTCGAGTTCTGCCTGGc                                                                                                                                                    >  1:1425526/1‑119 (MQ=255)
                                   cTGCACGCGCCGAAGTAGGCGCGGCGGGCATCGCGAACCTTATCCACAGGCAGTGGCTTGTTCCAAAGGTCGTAGACCTTCAGCAGGCGCTGTGAATGCGGGTCGAGTTCTGCCTGGCTGAACTCg                                                                                                                                            <  1:949874/126‑1 (MQ=255)
                                   cTGCACGCGCCGAAGTAGGCGCGGCGGGCATCGCGAACCTTATCCACAGGCAGTGGCTTGTTCCAAAGGTCGTAGACCTTCAGCAGGCGCTGTGAATGCGGGTCGAGTTCTGCCTGGCTGAACTCg                                                                                                                                            >  2:949874/1‑126 (MQ=255)
                                    tgtaCGCGCCGAAGTAGGCGCGGCGGGCATCGCGAACCTTATCCACAGGCAGTGGCTTGTTCCAAAGGTCGTAGACCTTCAGCAGGCGATGTGAATGCGGGTCGAGTTCTGCCTGGCTGAACTCGGCACGGGGCATGGGGATATCCacac                                                                                                                   <  2:624764/147‑1 (MQ=255)
                                    tGCACGCGCCGAAGTAGGCGCGGCGGGCATCGCGAACCTTATCCACAGGCAGTGGCTTGTTCCAAAGGTCGTAGACCTTCAGCAGGCGCTGTGAATGCGGGTCGAGTTCTGCCTGGCTGAACTCGGCACGGGGCATGGGGATATCCacac                                                                                                                   <  1:108365/150‑1 (MQ=255)
                                    tGCACGCGCCGAAGTAGGCGCGGCGGGCATCGCGAACCTTATCCACAGGCAGTGGCTTGTTCCAAAGGTCGTAGACCTTCAGCAGGCGCTGTGAATGCGGGTCGAGTTCTGCCTGGCTGAACTCGGCACGGGGCATGGGGATATCCacac                                                                                                                   <  1:1325485/150‑1 (MQ=255)
                                      cACGCGCCGAAGTAGGCGCGGCGGGCATCGCGAACCTTATCCACAGGCAGTGGCTTGTTTCAAAGGTCGTAGACCTTCAGCAGGCGCTGTGAATGCGGGTCGAGTTCTGCCTGGCTGAACTCGGCACGGGGCATGGGGATATCCACACCt                                                                                                                 >  1:876556/1‑150 (MQ=255)
                                      cACGCGCCGAAGTAGGCGCGGCGGGCATCGCGAACCTTATCCACAGGCAGTGGCTTGTTCCAAAGGTCGTAGACCTTCAGCAGGCGCTGTGAATGCGGGTCGAGTTCTGCCTGGCTGAACTCGGCACGGGGCATGGGGATATCCACACCt                                                                                                                 >  2:694780/1‑150 (MQ=255)
                                       aCGCGCCGAAGTAGGCGCGGCGGGCATCGCGAACCTTATCCACAGGCAGTGGCTTGTTCCAAAGGTCGTAGACCTTCAGCAGGCGCTGTGAATGCGGGTCGAGTTCTGCCTGGCTGAACTc                                                                                                                                             >  2:2112748/1‑121 (MQ=255)
                                       aCGCGCCGAAGTAGGCGCGGCGGGCATCGCGAACCTTATCCACAGGCAGTGGCTTGTTCCAAAGGTCGTAGACCTTCAGCAGGCGCTGTGAATGCGGGTCGAGTTCTGCCTGGCTGAACTc                                                                                                                                             <  1:2112748/121‑1 (MQ=255)
                                       aCGCGCCGAAGTAGGCGCGGCGGGCATCGCGAACCTTATCCACAGGCAGTGGCTTGTTCCAAAGGTCGTAGACCTTCAGCAGGCGCTGTGAATGCGGGTCGAGTTCTGCCTGGCTGAACTCGGCa                                                                                                                                         <  1:232444/125‑1 (MQ=255)
                                       aCGCGCCGAAGTAGGCGCGGCGGGCATCGCGAACCTTATCCACAGGCAGTGGCTTGTTCCAAAGGTCGTAGACCTTCAGCAGGCGCTGTGAATGCGGGTCGAGTTCTGCCTGGCTGAACTCGGCa                                                                                                                                         >  2:232444/1‑125 (MQ=255)
                                         gcgcCGAAGTAGGCGCGGCGGGCATCGCGAACCTTATCCACAGGCAGTGGCTTGTTCCAAAGGTCGTAGACCTTCAGCAGGCGCTGTGAATGCGGGTCGAGTTCTGCCTGGCTGAACTCGGCACGGGGCATGGGGATATCCACACCTTCg                                                                                                              >  2:537216/1‑150 (MQ=255)
                                           gcCGAAGTAGGCGCGGCGGGCATCGCGAACCTTATCCACAGGCAGTGGCTTGTTCCAAAGGTCGTAGACCTTCAGCAGGCGCTGTGAATGCGGGTCGAGTTCTGCCTGGCTGAACTCGGCACGGGGCATGGGGATATCCACACCTTCGTa                                                                                                            >  2:902550/1‑150 (MQ=255)
                                           gcCGAAGTAGGCGCGGCGGGCATCGCGAACCTTATCCACAGGCAGTGGCTTGTTCCAAAGGTCGTAGACCTTCAGCAGGCGCTGTGAATGCGGGTCGAGTTCTGCCTGGCTGAACTCGGCACGGGGCATGGGGATATCCAAACCTTCGTa                                                                                                            <  1:1693063/150‑1 (MQ=255)
                                            ccGAAGTAGGCGCGGCGGGCATCGCGAACCTTATCCACAGGCAGTGGCTTGTTCCAAAGGTAGTAGACCTTCAGCCGGCGCTGTGAATGCGGGTCGAGTTCTGCCTGGCTGAACTCGGCACGGGGCATGGGGATATCCACACCTTCGTAg                                                                                                           <  1:1449827/150‑1 (MQ=255)
                                            ccGAAGTAGGCGCGGCGGGCATCGCGAACCTTATCCACAGGCAGTGGCTTGTTCAAAAGGTCGTAGACCTTCAGCAGGCGCTGTGAATGCGGGTCGAGTTCTGCCTGGCTGAACTCGGCACGGGGCATGGGGATATCCACACCTTCGTAg                                                                                                           <  1:1207531/150‑1 (MQ=255)
                                            cagAAGTAGGCGCGGCGGGCATCGCGAACCTTATCCACAGGCAGTGGCTTGTTCCAAAGGTCGTAGACCTTCAGCAGGCGCTGTGAATGCGGGTCGAGTTCTGCCTGGCTGAACTCGGCACGGGGCATGGGGATATCCACACCTTCGTAg                                                                                                           <  2:1344202/148‑1 (MQ=255)
                                             cGAAGTAGGCGCGGCGGGCATCGCGAACCTTATCCACAGGCAGTGGCTTGTTCCAAAGGTCGTAGACCTTCAGCAGGCGCTGTGAATGCGGGTCGAGTTCTGCCTGGCTGAACTCGGCACGGGGCATGGGGATATCCACACCTTCGTAGc                                                                                                          <  2:515731/150‑1 (MQ=255)
                                                 gTAGGCGCGGCGGGCATCGCGAACCTTATCCACAGGCAGTGGCTTGTTCCAAAGGTCGTAGACCTTCAGCAGGCGCTGTGAATGCGGGTCGAGTTCTGCCTGGCTGAACTCGGCACGGGGCATGGGGATATCCACACCTTCGTAGCGGTc                                                                                                      <  1:565898/150‑1 (MQ=255)
                                                  tAGGCGCGGCGGGCATCGCGAACCTTATCCACAGGCAGTGGCTTGTTCCAAAGGTCGTAGACCTTCAGCAGGCGCTGTGAATGCGGGTCGAGTTCTGCCTGGCTGAACTCGGCACGGGGCATGGGGATATCCACACCTTCGTAGCGGTcc                                                                                                     <  1:537216/150‑1 (MQ=255)
                                                    ggCGCGGCGGGCATCGCGAACCTTATCCACAGGCAGTGGCTTGTTCCAAAGGTCGTAGACCTTCAGCAGGCGCTGTGAATGCTGGTCGAGTTCTGCCTGGCTGAACTCGGCACGGGGCATGGGGATATCCACACCTTCGTAGCGGTCCCa                                                                                                   >  1:1086663/1‑150 (MQ=255)
                                                    ggCGCGGCGGGCATCGCGAACCTTATCCACAGGCAGTGGCTTGTTCCAAAGGTCGTAGACCTTCAGCAGGCGCTGTGAATGCGGGTCGAGTTCTGCCTGGCTGAACTCGGCACGGGGCATGGGGATATCCACACCTTCGTAGCGGTCCCa                                                                                                   <  1:979909/150‑1 (MQ=255)
                                                    ggCGCGGCGGGCATCGCGAACCTTATCCACAGGCAGTGGCTTGTTCCAAAGGTCGTAGACCTTCAGCAGGCGCTGTGAATGCGGGTCGAGTTCTGCCTGGCTGAACTCGGCACGGGGCATGGGGATATCCACACCTTCGTAGCGGTCCCa                                                                                                   >  2:278076/1‑150 (MQ=255)
                                                       gcggcgGGCATCGCGAACCTTATCCACAGGCAGTGGCTTGTTCCAAAGGTCGTAGACCTTCAGCAGGCGCTGTGAATGCGGGTCGAGTTCTGCCTGGCTGAACTCGGCACGGGGCATGGGGATATCCACACCTTCGTAGCGGTCCCAGTa                                                                                                <  1:2175642/150‑1 (MQ=255)
                                                       gcggcgGGCATCGCGAACCTTATCCACAGGCAGTGGCTTGTTCCAAAGGTCGTAGACCTTCAGCAGGCGCTGTGAATGCGGGTCGAGTTCTGCCTGGCTGAACTCGGCACGGGGCATGGGGATATCCACACCTTCGTAGCGGTCCCAGTa                                                                                                >  2:558950/1‑150 (MQ=255)
                                                              gCATCGCGCAACTTATCCACAGGCCGTGGCTTTTTCCAAAGGTCTTAGACCTTCAGCAGGCGCTGTGAATGCTGGTCGAGTTCTGCCTGGCTGAACTCGGCACGGGGCATGGGGATATCCACACCTTCGTAGCGGTCCCAGTAGCGTTTg                                                                                         <  2:1086663/150‑1 (MQ=255)
                                                              gCATCGCGAACCTTATCCACAGGCAGTGGCTTGTTCCAAAGGTCGTAGACCTTCAGCAGGCGCTGTGAATGCGGGTCGAGTTCTGCCTGGCTGAACTCGGCACGGGGCATGGGGATATCCACACCTTCGTAGCGGTCCCAGTAGCGTTTg                                                                                         <  1:1771471/150‑1 (MQ=255)
                                                                    cgAACCTTATCCACAGGCAGTGGCTTGTTCCAAAGGTCGTAGACCTTCAGCAGGCGCTGTGAATGCGGGTCGAGTTCTGCCTGGCTGAACTCGGCACGGGGCATGGGGATATCCACACCTTCGTAGCGGTCCCAGTAGCGTTTGGCAATg                                                                                   <  1:1878754/150‑1 (MQ=255)
                                                                     gAACCTTATCCACAGGCAGTGGCTTGTTCCAAAGGTCGTAGACCTTCAGCAGGCGCTGTGCATGCGGTTCGAGTTCTGCCTGGCTGAACTCGGCACGGGGCATGGGGATATCCACACCTTCGTAGCGGTCCCAGTAGCGTTTGGCAATgg                                                                                  <  2:662611/150‑1 (MQ=255)
                                                                      aaCCTTATCCACAGGCAGTGGCTTGTTCCAAAGGTCGTAGACCTTCAGCAGGCGCTGTGAATGCGGGTCGAGTTCTGCCTGGCTGAACTCGGCACGGGGCATGGGGATATCCACACCTTCGTAGCGGTCCCAGTAGCGTTTGGCAATGgt                                                                                 <  2:1038748/150‑1 (MQ=255)
                                                                       aCCTTATCCACAGGCAGTGGCTTGTTCCAAAGGTCGTAGACCTTCAGCAGGCGCTGTGAATGCGGGTCGAGTTCTGCCTGGCTGAACTCGGCACGGGGCATGGGGATATCCACACCTTCGTAGCGGTCCCAGTAGCGTTTGGCAATGgtg                                                                                <  1:963324/150‑1 (MQ=255)
                                                                          ttATCCACAGGCAGTGGCTTGTTCCAAAGGTCGTAGACCTTCAGCAGGCGCTGTGAATGCGGGTCGAGTTCTGCCTGGCTGAACTCGGCACGGGGCATGGGGATATCCACACCTTCGTAGCGGTCCCAGTAGCGTTTGGCAATGGTGTAg                                                                             <  2:1769740/150‑1 (MQ=255)
                                                                           tATCCACAGGCAGTGGCTTGTTCCAAAGGTCGTAGACCTTCAGCAGGCGCTGTGAATGCGGGTCGAGTTCTGCCTGGCTGAACTCGGCACGGGGCATGGGGATATCCACACCTTCGTAGCGGTCCCAGTAGCGTTTGGCAATGGTGTAgg                                                                            <  1:782524/150‑1 (MQ=255)
                                                                             tCCACAGGCAGTGGCTTGTTCCAAAGGTCGTAGACCTTCAGCAGGCGCTGTGAATGCGGGTCGAGTTCTGCCTGGCTGAACTCGGCACGGGGCATGGGGATATCCACACCTTCGTAGCGGTCCCAGTAGCGTTTGGCAATGGTGTAgggg                                                                          <  1:1666448/150‑1 (MQ=255)
                                                                                     cAGTGGCTTGTTCCAAAGGTCGTAGACCTTCAGCAGGCGCTGTGAATGCGGGTCGAGTTCTGCCTGGCTGAACTCGGCACGGGGCATGGGGATATCCACACCTTCGTAGCGGTCCCAGTAGCGTTTGGCAATGGTGTAGGGGTCATGCgg                                                                  <  2:707043/150‑1 (MQ=255)
                                                                                        tGGCTTGTTCCAAAGGTCGTAGACCTTCAGCAGGCGCTGTGAATGCGGGTCGAGTTCTGCCTGGCTGAACTCGGCACGGGGCATGGGGATATCCACACCTTCGTAGCGGTCCCAGTAGCGTTTGGCAATGGTGTAGGGGTCATGCgggtg                                                               >  1:514982/1‑150 (MQ=255)
                                                                                         ggCTTGTTCCAAAGGTCGTAGACCTTCAGCAGGCGCTGTGAATGCGGGTCGAGTTCTGCCTGGCTGAACTCGGCACGGGGCATGGGGATATCCACACCTTCGTAGCGGTCCCAGTAGCGTTTGGCAATGGTGTAGGGGTCATGCgggtgg                                                              <  1:897255/150‑1 (MQ=255)
                                                                                            ttGTTCCAAAGGTCGTAGACCTTCAGCAGGCGCTGTGAATGCGGGTCGAGTTCTGCCTg                                                                                                                                                      >  1:2135403/1‑59 (MQ=255)
                                                                                            ttGTTCCAAAGGTCGTAGACCTTCAGCAGGCGCTGTGAATGCGGGTCGAGTTCTGCCTg                                                                                                                                                      <  2:2135403/59‑1 (MQ=255)
                                                                                            ttGTTCCAAAGGTCGTAGACCTTCAGCAGGCGCTGTGAATGCGGGTCGAGTTCTGCCTGGCTGAACTCGGCACGGGGCATGGGGATATCCACACCTTCGTAGCGGTCCCAGTAGCGTTTGGCAATGGTGTAGGGGTCATGCGGGTGGGTc                                                           >  1:959660/1‑150 (MQ=255)
                                                                                              gTTCCAAAGGTCGTAGACCTTCAGCAGGCGCTGTGAATGCGGGTCGAGTTCTGCCTGGCTGAACTCGGCACGGGGCATGGGGATATCCACACCTTCGTAGCGGTCCCAGTAGCGTTTGGCAATGGTGTAGGGGTCATGCGGGTGGGTCAt                                                         <  2:1656659/150‑1 (MQ=255)
                                                                                               tttcAAAGGTCGTAGACCTTCAGTAGGCTCTGTGAAGGCGGGTCGAGTTCTTCCTGGCTGAACTCGGCACGGGGCATGGGGATATCCACACCTTCGTAGCGGTCCCAGTAGCGTTTGGCAATGGTGTAGGGGTCATGCGGGTGGGTCATc                                                        <  2:876556/147‑1 (MQ=255)
                                                                                               ttCCAAAGGTCGTAGACCTTCAGCAGGCGCTGTGAATGCGGGTCGAGTTCTGCCTGGCTGAACTCGGCACGGGGCATGGGGATATCCACACCTTCGTAGCGGTCCCAGTAGCGTTTGGCAATGGTGTAGGGGTCATGCGGGTGGGTCATc                                                        >  2:1096404/1‑150 (MQ=255)
                                                                                                  cAAAGATCGTAGACCTTCAGCAGGCGCTGTGAATGCGGGTCGAGTTCTGCCTGGCTGAACTCGGCACGGGGCATGGGGATATCCACACCTTCGTAGCGGTCCCAGTa                                                                                                <  1:1150132/107‑1 (MQ=255)
                                                                                                  cAAAGATCGTAGACCTTCAGCAGGCGCTGTGAATGCGGGTCGAGTTCTGCCTGGCTGAACTCGGCACGGGGCATGGGGATATCCACACCTTCGTAGCGGTCCCAGTa                                                                                                >  2:1150132/1‑107 (MQ=255)
                                                                                                   aacggTCGTAGACCTTCAGCAGGCGCTGTGAATGCGGGTCGAGTTCTGCCTGGCTGAACTCGGCACGGGGCATGGGGATATCCACACCTTCGTAGCGGTCCCAGTAGCGTTTGGCAATGGTGTAGGGGTCATGCGGGTGGGTCATCGaaa                                                    >  2:1842690/4‑150 (MQ=255)
                                                                                                   aaaGGTCGTAGACCTTCAGCAGGCGCTGTGAATGCGGGTCGAGTTCTGCCTGGCTGAACTCGGCACGGGGCATGGGGATATCCACACCTTCGTAGCGGTCCCAGTAGCGTTTGGCAATGGTGTAGGGGTCATGCGGGTGGGTCATCGaaa                                                    <  1:576731/150‑1 (MQ=255)
                                                                                                   aaaGGTCGTAGACCTTCAGCAGGCGCTGTGAATGCGGGTCGAGTTCTGCCTGGCTGAACTCGGCACGGGGCATGGGGATATCCACACCTTCGTAGCGGTCCCAGTAGCGTTTGGCAATGGTGTAGGGGTCATGCGGGTGGGTCATCGaaa                                                    <  1:1842690/150‑1 (MQ=255)
                                                                                                        tCGTAGACCTTCAGCAGGCGCTGTGAATGCGGGTCGAGTTCTGCCTGGCTGAACTCGGCACGGGGCATGGGGATATCCACACCTTCGTAGCGGTCCCAGTAGCGTTTGGCAATGGTGTAGGGGTCATGCGGGTGGGTCATCGAAACGGTc                                               <  1:535420/150‑1 (MQ=255)
                                                                                                         cGTAGACCTTCAGCAGGCGCTGTGAATGCGGGTCGAGTTCTGCCTGGCTGAACTCGGCACGGGGCATGGGGATATCCACACCTTCGTAGCGGTCCCAGTAGCGTTTGGCAATGGTGTAGGGGTCATGCGGGTGGGTCCTCGAAACGGTCa                                              <  1:902550/150‑1 (MQ=255)
                                                                                                            aGACCTTCAGCAGGCGCTGTGAATGCGGGTCGAGTTCTGCCTGGCTGAATTCGTCACGGGGCATGGGGATATCCACACCTTCGTAGCGGTCCCAGTAGCGTTTGGCAATTGTGTAGGGGTCATGCGGGTGGGTCATCGAAACggagatcg                                           >  1:1227242/1‑144 (MQ=255)
                                                                                                                cTTCAGCAGGCGCTGTGAATGCGGGTCGAGTTCTGCCTGGCTGAACTCGGCACGGGGCATGGGGATATCCACACCTTCGTAGCGGTCCCAGTAGCGtt                                                                                           <  2:1862159/98‑1 (MQ=255)
                                                                                                                cTTCAGCAGGCGCTGTGAATGCGGGTCGAGTTCTGCCTGGCTGAACTCGGCACGGGGCATGGGGATATCCACACCTTCGTAGCGGTCCCAGTAGCGtt                                                                                           >  1:1862159/1‑98 (MQ=255)
                                                                                                                   cagcagGCGCTGTGAATGCGGGTCGAGTTCTGCCTGGCTGAACTCGGCACGGGGCATGGGGATATCCACACCTTCGTAGCGGTCCCAGTAGCGTTTGGCAATGGTGTAGGGGTCATGCGGGTGGGTCATCGAAACGGTCAGGCAAAATgg                                    >  2:881113/1‑150 (MQ=255)
                                                                                                                     gcagGCGCTGTGAATGCGGGTCGAGTTCTGCCTGGCTGAACTCGGCCAGAGTCAGAGGGATATCCACACCTTCGTAGCGGTCCCAGTAGCGTTTGGCAATGGTGTAGGGGTCATGCGGGTGGGTCATCGAAACGGGCAGGCAAAATgtgc                                  >  2:327143/1‑147 (MQ=255)
                                                                                                                            cTGTGAATGCGGGTCGAGTTCTGCCTGGCTGAACTCGGCACGGGGCATGGGGATATCCACACCTTCGTAGCGGTCCCAGTAGCGTTTGGCAATGGTGTAGGGGTCATGCGGGTGGGTCATCGAAACGGTCAGGCAAAATGGCCGGCCATc                           <  2:36015/150‑1 (MQ=255)
                                                                                                                                  aTGCGGGTCGAGTTCTGCCTGGCTGAACTCGGCACGGGGCATGGGGATATCCACACCTTCGTAGCGGTCCCAGTAGCGTTTGGCAATGGTGTAGGGGTCATGCGGGTGGGTCATCGAAACGGTCAGGCAAAATGGCCGGCCa                             >  1:1396581/1‑142 (MQ=255)
                                                                                                                                  aTGCGGGTCGAGTTCTGCCTGGCTGAACTCGGCACGGGGCATGGGGATATCCACACCTTCGTAGCGGTCCCAGTAGCGTTTGGCAATGGTGTAGGGGTCATGCGGGTGGGTCATCGAAACGGTCAGGCAAAATGGCCGGCCa                             <  2:1396581/142‑1 (MQ=255)
                                                                                                                                       ggTCGAGTTCTGCCTGGCTGAACTCGGCACGGGGCATGGGGATATCCACACCTTCGTAGCGGTCCCAGTAGCGTTTGGCAATGGTGTAGGGGTCATGCGGGTGGGTCATCGAAACGGTCAgg                                            >  2:441491/1‑122 (MQ=255)
                                                                                                                                       ggTCGAGTTCTGCCTGGCTGAACTCGGCACGGGGCATGGGGATATCCACACCTTCGTAGCGGTCCCAGTAGCGTTTGGCAATGGTGTAGGGGTCATGCGGGTGGGTCATCGAAACGGTCAgg                                            <  1:441491/122‑1 (MQ=255)
                                                                                                                                        gTCGAGTTCTGCCTGGCTGAACTCGGCACGGGGCATGGGGATATCCACACCTTCGTAGCGGTCCCAGTAGCGttt                                                                                          >  2:1401491/1‑75 (MQ=255)
                                                                                                                                        gTCGAGTTCTGCCTGGCTGAACTCGGCACGGGGCATGGGGATATCCACACCTTCGTAGCGGTCCCAGTAGCGttt                                                                                          <  1:1401491/75‑1 (MQ=255)
                                                                                                                                        gTCGAGTTCTGCCTGGCTGAACTCGGCACGGGGCATGGGGATATCCACACCTTCGTAGCGGTCCCAGTAGCGTTTGGCAATg                                                                                   >  1:829857/1‑82 (MQ=255)
                                                                                                                                        gTCGAGTTCTGCCTGGCTGAACTCGGCACGGGGCATGGGGATATCCACACCTTCGTAGCGGTCCCAGTAGCGTTTGGCAATg                                                                                   <  2:829857/82‑1 (MQ=255)
                                                                                                                                          cGAGTTCAGCCAGGCAGAACTCGGCACGGGGCATGGGGATATCCACACCTTCGTAGCGGTCCCAGTAGCGTTTGGCAATGGTGTAgggg                                                                          <  2:1824891/89‑1 (MQ=39)
                                                                                                                                          cGAGTTCAGCCAGGCAGAACTCGGCACGGGGCATGGGGATATCCACACCTTCGTAGCGGTCCCAGTAGCGTTTGGCAATGGTGTAgggg                                                                          >  1:1824891/1‑89 (MQ=39)
                                                                                                                                          cGAGTTCAGCCAGGCAGAACTCGGCACGGGGCATGGGGATATCCACACCTTCGTAGCGGTCCCAGTAGCGTTTGGCAATGGTGTAGGGGTCATGCGGGTGGGTCATCGAAAc                                                   >  2:1958991/1‑112 (MQ=255)
                                                                                                                                          cGAGTTCAGCCAGGCAGAACTCGGCACGGGGCATGGGGATATCCACACCTTCGTAGCGGTCCCAGTAGCGTTTGGCAATGGTGTAGGGGTCATGCGGGTGGGTCATCGAAAc                                                   <  1:1958991/112‑1 (MQ=255)
                                                                                                                                              ttCTGCCTGGCTTAACTCGGCACGGGGCATGGGGATATCCACACCTTCTTAGCGGTCCCAGTAGCGTTTGGCAATGGTGTAGGGGTCATGCGGGTGGGTCATCGAAACGGTCAGGCAAAATGGCCGGCCATCGTTTTCGCGCACATGGTc         <  1:1274376/150‑1 (MQ=255)
                                                                                                                                              ttCTGCCTGGCTGAACTCGGCACGGGGCATGGGGATATCCACACCTTCGTAGCGGTCCCAGTAGCGTTTGGCAATGGTGTAGGGGTCATGCGGGTGGGTCATCGAAACGGTCAGGCAAAATGGCCGGCCATCGTTTTCGCGCACATGGTc         >  2:1053321/1‑150 (MQ=255)
                                                                                                                                                 tGCCTGGCTGAACTCGGCACGGGGCATGGGGATATCCACACCTTCGTAGCGGTCCCAGTAGCGTTTGGCAATGGTGTAGGGGTCATGCGGGTGGGTCATCGAAACGGTCAGGCAAAATGGCCGGCCATCGTTTTCGCGCACATGGTCGTa      <  2:514982/150‑1 (MQ=255)
                                                                                                                                                 tGCCTGGCTGAACTCGGCACGGGGCATGGGGATATCCACACCTTCGTAGCGGTCCCAGTAGCGTTTGGCAATGGTGTAGGGGTCATGCGGGTGGGTCATCGAAACGGTCAGGCAAAATGGCCGGCCATCGTTTTCGCGCACATGGTCGTa      >  2:1763757/1‑150 (MQ=255)
                                                                                                                                                     tGGCTGAACTCGGCACGGGGCATGGGGATATCCACACCTTCGTAGCGGTCCCAGTAGCGTTTTTCAATGGTGTAGGGGTCATGCGGGTGGGTCATCGAAACGGTCAGGCAAAATGGCCGGCCATCGTTTTCGCGCACATGGTCGTACAgg  <  2:1824129/150‑1 (MQ=255)
                                                                                                                                                     tGGCTGAACTCGGCACGGGGCATGGGGATATCCACACCTTCGTAGCGGTCCCAGTAGCGTTTGGCAATGGTGTAGGGGTCATGCGGGTGGGTCATCGAAACGGTCAGGCAAAATGGCCGGCCATCGTTTTCGCGCACATGGTCGTACAgg  <  1:382022/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
TCTGCAGCAATTGGCCGATGTTGTCATCGATATAGCTGCACGCGCCGAAGTAGGCGCGGCGGGCATCGCGAACCTTATCCACAGGCAGTGGCTTGTTCCAAAGGTCGTAGACCTTCAGCAGGCGCTGTGAATGCGGGTCGAGTTCTGCCTGGCTGAACTCGGCACGGGGCATGGGGATATCCACACCTTCGTAGCGGTCCCAGTAGCGTTTGGCAATGGTGTAGGGGTCATGCGGGTGGGTCATCGAAACGGTCAGGCAAAATGGCCGGCCATCGTTTTCGCGCACATGGTCGTACAGG  >  NC_002947/88504‑88802

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