Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 87,959 G→C 8.8% R453R (CGC→CGG betC ← choline‑sulfatase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_00294787,9590GC 8.8% 173.8 / 6.4 57R453R (CGC→CGGbetCcholine‑sulfatase
Reads supporting (aligned to +/- strand):  ref base G (28/24);  new base C (2/3);  total (30/27)
Fisher's exact test for biased strand distribution p-value = 6.60e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.68e-01

GGCGGGCTGTGGATAACGTGCCTTGCGTTCGAGGTCATCGAGATCGATGTGGTTGCGCATGTACTGTTGGCTGGCGTCCACCAGTGGTTGGTGGTCCCAGCTTTTCAGCTTGCCGATGGCCAGCGCTTCGGCCACCAGCCGGCGGCGGCGCTGGCTGGCCAGTACCTGCTGGCGCAGGCTGGTGATGTCCCAGCGCTGTTGCGCTTCATCGACAAATGCCTGCAGCAGCACCTGGTGGTCCGGGCTGCCGGTGAGGTTCTCCCGCTCGTGCGGGTCGCGGCTCAGGTCATA  >  NC_002947/87812‑88102
                                                                                                                                                   |                                                                                                                                               
ggcgggcTGTGGATAACGTGCCTTGCGTTCGAGGTCATCGAGATCGATGTGGTTGCGCATGTACTGTTGGCTGGCGTCCACCAGTGGTTGGTGGTCCCAGCTTTTCAGCTTGCCGATGGCCAGCGCTTCGGCCACCAGCCggcggcggc                                                                                                                                                >  8:125015/1‑149 (MQ=255)
        gtgGATAACGTGCCTTGCGTTCGAGGTCATCGAGATCGATGTGGTTGCGCATGTACTGTTGGCTGGCGTCCACCAGTGGTTGGTGGTCCCAGCTTTTCAGCTTGCCGATGGCCAGCGCTTCGGCCACCAGCCGGCGGCGGCgctggctg                                                                                                                                        >  3:325023/1‑149 (MQ=255)
         tgGATAACGTGCCTTGCGTTCGAGGTCATCGAGATCGATGTGGTTGCGCATGTACTGTTGGCTGGCGTCCACCAGTGGTTGGTGGTCCCAGCTTTTCAGCTTGCCGATGGCCAGCGCTTCGGCCACCAGCCGGCGGCGGCgctggct                                                                                                                                         <  2:181645/147‑1 (MQ=255)
         tgGATAACGTGCCTTGCGTTCGAGGTCATCGAGATCGATGTGGTTGCGCATGTACTGTTGGCTGGCGTCCACCAGTGGTTGGTGGTCCCAGCTTTTCAGCTTGCCGATGGCCAGCGCTTCGGCCACCAGCCAGCGCCGCCgccggctg                                                                                                                                        <  5:98417/148‑1 (MQ=255)
           gATAACGTGCCTTGCGTTCGAGGTCATCGAGATCGATGTGGTTGCGCATGTACTGTTGGCTGGCGTCCACCAGTGGTTGGTGGTCCCAGCTTTTCAGCTTGCCGATGGCCAGCGCTTCGGCCACCAGCCGGCGGCGGCGCtggctggc                                                                                                                                      >  1:175842/1‑148 (MQ=255)
           gATAACGTGCCTTGCGTTCGAGGTCATCGAGATCGATGTGGTTGCGCATGTACTGTTGGCTGGCGTCCACCAGTGGTTGGTGGTCCCAGCTTTTCAGCTTGCCGATGGCCAGCGCTTCGGCCACCAGCCGGCGGCGGCGCtggctggc                                                                                                                                      <  2:175842/148‑1 (MQ=255)
             tAACGTGCCTTGCGTTCGAGGTCATCGAGATCGATGTGGTTGCGCATGTACTGTTGGCTGGCGTCCACCAGTGGTTGGTGGTCCCAGCTTTTCAGCTTGCCGATGGCCAGCGCTTCGGCCACCAGCCGGCGGCGGCGCTGGCTGGCCAg                                                                                                                                   >  1:183423/1‑149 (MQ=255)
               aCGTGCCTTGCGTTCGAGGTCATCGAGATCGATGTGGTTGCGCATGTACTGTTGGCTGGCGTCCACCAGTGGTTGGTGGTCCCAGCTTTTCAGCTTGCCGATGGCCAGCGCTCCGGCCACCAGCCGGCGGCGGCGCTGGCTGGCCAg                                                                                                                                   <  5:192509/147‑1 (MQ=255)
                  tGCCTGGCGTTCGAGTGCATCGAGATCGATGTGGTTGCGCATGTACTGTTGGCTGGCGTCCACCAGTGGGTGGTGGTCCCAGCTTTTCAGCTTGCCGATGGCCAGCGCTTCGGCCACCAGCCGGCGGCGGCGCTGGCTGGCCAGTACct                                                                                                                              <  6:349449/149‑1 (MQ=255)
                      ttGCGTTCGAGGTCATCGAGATCGATGTGGTTGCGCATGTACTGTTGGCTGGCGTCCACCAGTGGTTGGTGGTCTCAGCTTTTCAGCTTGCCGATGGCCAGCGCTTCGGCCACCAGCCGGCGGCGGCgctgg                                                                                                                                           <  6:12232/132‑1 (MQ=255)
                      ttGCGTTCGAGGTCATCGAGATCGATGTGGTTGCGCATGTACTGTTGGCTGGCGTCCACCAGTGGTTGGTGGTCTCAGCTTTTCAGCTTGCCGATGGCCAGCGCTTCGGCCACCAGCCGGCGGCGGCgctgg                                                                                                                                           >  5:12232/1‑132 (MQ=255)
                          gTTCGAGGTCATCGAGATCGATGTGGTTGCGCATGTACTGTTGGCTGGCGTCCACCAGTGGTTGGTGGTCCCAGCTTTTCAGCTTGCCGATGGCCAGCGCTTCGGCCACCAGCCGGCGGCGGCGCTGGCTGGCCAGTACCTGCTGgcgc                                                                                                                      <  2:183423/149‑1 (MQ=255)
                          gTTCGAGGTCATCGAGATCGATGTGGTTGCGCATGTACTGTTGGCTGGCGTCCACCAGTGGTTGGTGGTCCCAGCTTTTCAGCTTGCCGATGGCCAGCGCTTCGGCCACCAGCCGGCGGCGGCGCTGGCTGGCCAGTACCTGCTGgcg                                                                                                                       >  6:212385/1‑148 (MQ=255)
                                 gTCATCGAGATCGATGTGGTTGCGCATGTACTGTTGGCTGGCGTCCACCAGTGGTTGGTGGACCCAGCTTTTCAGCTTGCCGATGGCCAGCGCTTCGGCCACCAGCCGGCGGCGGCGCTGGCTGGCCAGTACCTGCTGGCGCAGGCTgg                                                                                                               >  7:69575/1‑149 (MQ=255)
                                  tCATCGAGATCGATGTGGTTGCGCATGTACTGTTGGCTGGCGTCCACCAGTGGTTGGTGGTCCCAGCTTTTCAGCTTGCCGATGGCCAGCGCTTCGGCCACCAGCCGGCGGCGGCGCTGGCTGGCCAGTACCTGCTGGCGCAGGCTGGt                                                                                                              <  4:325023/149‑1 (MQ=255)
                                  tCATCGAGATCGATGTGGTTGCGCATGTACTGTTGGCTGGCGTCCACCAGTGGTTGGTGGGCCCAGCTTTTCAGCTTGCCGATGGCCAGCGCTTCGGCCACCAGCCGGCGGCGGCGCTGGCTGGCCAGTACCTGCTGGCGCAGGCTGGt                                                                                                              <  8:69575/149‑1 (MQ=255)
                                          aTCGATGTGGTTGCGCATGTACTGTTGGCTGGCGTCCACCAGTGGTTGGTGGTCCCAGCTTTTCAGCTTGCCGATGGCCAGCGCTTCGGCCACCAGCCGGCGGCGGCGCTGGCTGGCCAGTACCTGCTGGCGCAGGCTGGTGATGTcc                                                                                                       >  8:64034/1‑148 (MQ=255)
                                          aTCGATGTGGTTGCGCATGTACTGTTGGCTGGCGTCCACCAGTGGTTGGTGGTCCCAGCTTTTCAGCTTGCCGATGGCCAGCGCTTCGGCCACCAGCCGGCGGCGGCGCTGGCTGGCCAGTACCTGCTGGCGCAGGCTGGTGATGTcc                                                                                                       >  6:224371/1‑148 (MQ=255)
                                             gATGTGGTTGCGCATGTACTGTTGGCTGGCGTCCACCAGTGGTTGGTGGTCCCAGCTTTTCAGCTTGCCGATGGCCAGCGCTTCGGCCACCAGCCGGCGGCGGCgctggct                                                                                                                                         >  8:256766/1‑111 (MQ=255)
                                             gATGTGGTTGCGCATGTACTGTTGGCTGGCGTCCACCAGTGGTTGGTGGTCCCAGCTTTTCAGCTTGCCGATGGCCAGCGCTTCGGCCACCAGCCGGCGGCGGCgctggct                                                                                                                                         <  7:256766/111‑1 (MQ=255)
                                                gtgGTTGCGCATGTACTGTTGGCTGGCGTCCACCAGTGGTTGGTGGTCCCAGCTTTTCAGCTTGCCGATGGCCAGCGCTTCGGCCACCAGCCGGCGGCGGCGCTGGCTGGCCAGTACCTGCTGGCGCAGGCTGGTGATGTCCCAGCGCt                                                                                                >  1:350532/1‑149 (MQ=255)
                                                gtgGTTGCGCATGTACTGTTGGCTGGCGTCCACCAGTGGTTGGTGGTCCCAGCTTTTCAGCTTGCCGATGGCCAGCGCTTCGGCCACCAGCCGGCGGCGGCGCTGGCAGGCCAGTACCTGCTGGCGCAGGCTGGTGATGTCCCAGCGCt                                                                                                >  2:153844/1‑149 (MQ=255)
                                                   gTTGCGCATGTACTGTTGGCTGGCGTCCACCAGTGGTTGGTGGTCCCAGCTTTTCAGCTTGCCGATGGCCAGCGCTTCGGCCACCAGCCGGCGGCGGCGCTGGCTGGCCAGTACCTGCTGGCGCAGGCTGGt                                                                                                              >  5:20423/1‑132 (MQ=255)
                                                   gTTGCGCATGTACTGTTGGCTGGCGTCCACCAGTGGTTGGTGGTCCCAGCTTTTCAGCTTGCCGATGGCCAGCGCTTCGGCCACCAGCCGGCGGCGGCGCTGGCTGGCCAGTACCTGCTGGCGCAGGCTGGt                                                                                                              <  6:20423/132‑1 (MQ=255)
                                                   gTTGCGCATGTACTGTTGGCTGGCGTCCACCAGTGGTTGGTGGTCCCAGCTTTTCAGCTTGCCGATGGCCAGCGCTTCGGCCACCAGCCGGCGGCGGCGCTGGCTGGCCAGTACCTGCTGGCGCAGGCTGGTGATGTCCCAGCGCTGtt                                                                                             <  7:125015/149‑1 (MQ=255)
                                                          aTGTACTGTTGGCTGGCGTCCACCAGTGGTTGGTGGTCCCAGCTTTTCAGCTTGCCGATGGCCAGCGCTTCGGCCACCAGCCGGCGGCGGCGCTGGCTGGCCAGTACCTGCTGGCGCAGGCTGGt                                                                                                              >  5:182532/1‑125 (MQ=255)
                                                          aTGTACTGTTGGCTGGCGTCCACCAGTGGTTGGTGGTCCCAGCTTTTCAGCTTGCCGATGGCCAGCGCTTCGGCCACCAGCCGGCGGCGGCGCTGGCTGGCCAGTACCTGCTGGCGCAGGCTGGt                                                                                                              <  6:182532/125‑1 (MQ=255)
                                                                  ttGGCTGGCGTCCACCAGTGGTTGGTGGTCCCAGCTTTTCAGCTTGCCGATGGCCAGCGCTTCGGCCACCAGCCGGCGGCGGCGCTGGCTGGCCAGTACCTGCTGGCGCAGGCTGGTGATGTCCCAGCGCTGTTGCGCTTCATCGACa                                                                               >  3:199174/1‑148 (MQ=255)
                                                                    ggctggcGTCCACCAGTGGTTGGTGGTCCCAGCTTTTCAGCTTGCCGATGGCCAGCGCTTCGGCCACCAGCCGGCGGCGGCGCTGGCTGGCCAGTACCTGCTGGCGCAGGCTGGTGATGTCCCAGCGCTGTTGCGCTTCATCGACaaa                                                                             <  2:350532/148‑1 (MQ=255)
                                                                     gctggcGTCCACCAGTGGTTGGTGGTCCCAGCTCTTCAGCTTGCCGATGGCCAGCGCTTCGGCCACCAGCCGGCGGCGGCGCTGGCTGGCCAGTACCTGCTGGCGCAGGCTGGTGATGTCCCAGCGCTGTTGCGCTTCATc                                                                                   >  6:249644/1‑141 (MQ=255)
                                                                     gctggcGTCCACCAGTGGTTGGTGGTCCCAGCTCTTCAGCTTGCCGATGGCCAGCGCTTCGGCCACCAGCCGGCGGCGGCGCTGGCTGGCCAGTACCTGCTGGCGCAGGCTGGTGATGTCCCAGCGCTGTTGCGCTTCATc                                                                                   <  5:249644/141‑1 (MQ=255)
                                                                      ctggcGTCCACCAGTGGTTGGTGGTCCCAGCTTTTCAGCTTGCCGATGGCCAGCGCTTCGGCCACCAGCCGGCGGCGGCGCTGGCTGGCCAGTACCTGCTGGCGCAGGCTGGTGATGTCCCAGCGCTGTTGCGCTTCATCGACAAATGc                                                                          <  4:199174/149‑1 (MQ=255)
                                                                         gcGTCCACCAGTGGTTGGTGGTCCCAGCTTTTCAGCTTGCCGATGGCCAGCGCTTCGGCCACCAGCCGGCGGCGGCGCTGGCTGGCCAGTACCTGCTGGCGCAGGCTGGTGATGTCCCAGCGCTGTTGCGCTTCATCGACAAATGCCTg                                                                       >  2:221163/1‑149 (MQ=255)
                                                                               accaGTGGTTGGTGGTCCCAGCTTTTCAGCTTGCCGATGGCCAGCGCTTCGGCCACCAGCCGGCGGCGGCGCTGGCTGGCCAGTACCTGCTGGCGCAGGCTGGTGATGTCCCAGCGCTGTTGCGCTTCATCGACAAATGCCTgcagcag                                                                 >  2:73962/1‑149 (MQ=255)
                                                                                      gttggtGGTCCCAGCTTTTCAGCTTGCCGATGGCCAGCGCTTCGGCCACCAGCCGGCGGCGGCGCTGGCTGGCCAGTACCTGCTGGCGCAGGCTGGTGATGTCCCAGCGCTGTTGCGCTTCATCGACAAATGCCTGCAGCAGCACCtgg                                                          <  5:224371/149‑1 (MQ=255)
                                                                                              tCCCAGCTTTTCAGCTTGCCGATGGCCAGCGCTTCGGCCACCAGCCGGCGGCGGCGCTGGCTGGCCAGTACCTGCTGGCGCAGGCTGGTGATGTCCCAGCGCTGTTGCGCTTCATCGACAAATGCCTGCAGCAGCACCTGGTGGTCCgg                                                  >  6:119499/1‑149 (MQ=255)
                                                                                              tCCCAGCTTTTCAGCTTGCCGATGGCCAGCGCTTCGGCCACCAGCCGGCGGCGGCGCTGGCTGGCCAGTACCTGCTGGCGCAGGCTGGTGATGTCCCAGCGCTGTTGCGCTTCATCGACAAATGCCTGCAGCAGCACCTGGTGGTCCgg                                                  <  5:212385/149‑1 (MQ=255)
                                                                                                ccAGCTTTTCAGCTTGCCGATGGCCAGCGCTTCGGCCACCAGCCGGCGGCGGTGCTGGCTGGCCAGTACCTGATGGCGCAGGCTGGTGATGTCCCAGCGCTGTTGCGCTTCATCGACAAATGCCTGCAGCAGCCCCTGGTGGTCCGGGc                                                >  1:156438/1‑149 (MQ=255)
                                                                                                ccAGCTTTTCAGCTTGCCGATGGCCAGCGCTTCGGCCACCAGCCGGCGGCGGCGCTGGCTGGCCAGTACCTGCTGGCGCAGGCTGGTGATGTCCCAGCGCTGTTGCGCTTCATCGACAAATGCCTGCAGCAGCACCTGGTGGTCCGGGc                                                >  7:235058/1‑149 (MQ=255)
                                                                                                     ttttCAGCTTGCCGATGGCCAGCGCTTCGGCCACCAGCCGGCGGCGGCGCTGGCTGGCCAGTACCTGCTGGCGCAGGCTGGTGATGTCCCAGCGCTGTTGCGCTTCATCGACAAATGCCTGCAGCAGCACCtg                                                           >  5:7810/1‑133 (MQ=255)
                                                                                                     ttttCAGCTTGCCGATGGCCAGCGCTTCGGCCACCAGCCGGCGGCGGCGCTGGCTGGCCAGTACCTGCTGGCGCAGGCTGGTGATGTCCCAGCGCTGTTGCGCTTCATCGACAAATGCCTGCAGCAGCACCtg                                                           <  6:7810/133‑1 (MQ=255)
                                                                                                     ttttCAGCTTGCCGATGGCCAGCGCTTCGGCCACCAGCCGGCGGCGGCGCTGGCTGGCCAGTACCTGCTGGCGCAGGCTGGTGATGTCCCAGCGCTGTTGCGCTTCATCGACAAATGCCTGCAGCAGCACCTGGTGGTCCGGGCTGCCg                                           >  5:338491/1‑149 (MQ=255)
                                                                                                          aGCTTGCCGATGGCCAGCGCTTCGGCCACCAGCCGGCGGCGGCGCTGGCTGGCCAGTACCTGCTGGCGCAGGCTGGTGATGTCCCAGCGCTGTTGCGCTTCATCGACAAATGCCTGCAGCAGCACCTGGTGGTCCGGGCTGCCGGTGAg                                      <  5:367540/149‑1 (MQ=255)
                                                                                                          aGCTTGCCGATGGCCAGCGCTTCGGCCACCAGCCGGCGGCGGCGCTGGCTGGCCAGTACCTGCTGGCGCAGGCTGGTGATGTCCCAGCGCTGTTGCGCTTCATCGACAAATGCCTGCAGCAGCACCTGGTGGTCCGGGCTGCCGGTGAg                                      >  7:228663/1‑149 (MQ=255)
                                                                                                                       ccAGCGCTTCGGCCACCAGCCGGCGGCGGCGCTGGCTGGCCAGTACCTGCTGGCGCAGGCTGGTGATGTCCCAGCGCTGTTGCGCCTCATCGACAAATGCCTGCAGCAGCACCTGGTGGTCCGGGCTGCCGGTGAGGTTCTCCCGCTCg                         >  5:240044/1‑149 (MQ=255)
                                                                                                                          gcgcTTGGGCCACGAGGCGGCGGCGGCGCTGGCTGGCCAGTACCTGCTGGCGCTGGCTGGTGATGTCCCAGCGCTGTTGCGCTTCATCGACAAATGCCTGCAGCAGCACCTGGTGGTCCTTGCTGCCGGTGAGGTTCTCCCGCTCGTGc                      <  5:367978/149‑1 (MQ=255)
                                                                                                                          gcgcTTCGGCCACCAGCCGGCGGCGGCGCTGGCTGGCCAGTACCTGCTGGCGCAGGCTGGTGATGTCCCAGCGCTGTTGCGCTTCATCGACAAATGCCTGCAGCAGCACCTGGTGGTCCGGGCTGCCGGTGAGGTTCTCCCGCTCGTGc                      <  8:228663/149‑1 (MQ=255)
                                                                                                                            gcTTCGGCCACCAGCCGGCGGCGGCGCTGGCTGGCCAGTACCTGCTGGCGCAGGCTGGTGATGTCCCAGCGCTGTTGCGCTTCATCGACAAATGCCTGCAGCAGCACCTGGTGGTCCGGGCTGCCGGTGAGGTTCTCCCGCTCGTGCgg                    >  7:123998/1‑149 (MQ=255)
                                                                                                                                  gCCACCAGCCGGCGGCGGCGCTGGCTGGCCAGTACCTGCTGGCGCAGGCTGGTGATGTCCCAGCGCTGTTGCGCTTCATCGACAAATGCCTGCAGCAGCACCTGGTGGTCCGGGCTGCCGGTGAGGTTCTCCCGCTCGTGCGGGTcgcg              >  2:146261/1‑149 (MQ=255)
                                                                                                                                  gCCACCAGCCGGCGGCGGCGCTGGCTGGCCAGTACCTGCTGGCGCAGGCTGGTGATGTCCCAGCGCTGTTGCGCTTCATCGACAAATGCCTGCAGCAGCACCTGGTGGTCCGGGCTGCCGGTGAGGTTCTCCCGCTCGTGCGGGTcgcg              >  4:217959/1‑149 (MQ=255)
                                                                                                                                     gccagccagCGCCGCCGCCGGCTGGCCAGTACCTGCTGGCGCAGGCTGGTGATGTCCCAGCGCTGTTGCGCTTCATCGACAAAtac                                                                          >  2:220596/2‑84 (MQ=38)
                                                                                                                                     gccagccagCGCCGCCGCCGGCTGGCCAGTACCTGCTGGCGCAGGCTGGTGATGTCCCAGCGCTGTTGCGCTTCATCGACAAAtac                                                                          <  1:220596/85‑3 (MQ=38)
                                                                                                                                       cagccagcGCCGCCGCCGGCTGGTCAGTACCTGCTGGCGCAGGCTGGTGATGTCCCAGCGCTGTTGCGCTTCATCGACAAATGCCTGCAGCAGCACCTGGTGGTCCGGGCTGCCGGTGa                                       >  6:336553/1‑119 (MQ=255)
                                                                                                                                       cagccagcGCCGCCGCCGGCTGGTCAGTACCTGCTGGCGCAGGCTGGTGATGTCCCAGCGCTGTTGCGCTTCATCGACAAATGCCTGCAGCAGCACCTGGTGGTCCGGGCTGCCGGTGa                                       <  5:336553/119‑1 (MQ=255)
                                                                                                                                       caGCCGGCGGCGGCGCTGGCTGGCCAGTACCTGCTGGCGCAGGCTGGTGATGTCCCAGCGCTGTTGCGCTTCATCGACAAATGCCTGCAGCAGCACCTGGTGGTCCGGGCTGCCGGTGAg                                      <  7:5948/120‑1 (MQ=255)
                                                                                                                                       caGCCGGCGGCGGCGCTGGCTGGCCAGTACCTGCTGGCGCAGGCTGGTGATGTCCCAGCGCTGTTGCGCTTCATCGACAAATGCCTGCAGCAGCACCTGGTGGTCCGGGCTGCCGGTGAg                                      >  8:5948/1‑120 (MQ=255)
                                                                                                                                       caGCCGGCGGCGGCGCTGGCTGGCCAGTACCTGCTGGCGCAGGCTGGTGATGTCCCAGCGCTGTTGCGCTTCATCGACAAATGCCTGCAGCAGCACCTGGTGGTCCGGGCTGCCGGTGAGGTTCTCCCGCTCGTGCGGGTCGCGGCTCa         <  8:235058/149‑1 (MQ=255)
                                                                                                                                          ccGGCGGCGGCGCTGGCTGGCCAGTACCTGCTGGCGCAGGCTGGTGATGTCCCAGCGCTGTTGCGCTTCATCGACAAATGCCTGCAGCAGCACCTGGTGGTCCGGGCTGCCGGTGAGGTTCTCCCGCTCGTGCGGGTCGCGGCTCAGGt      <  1:73962/149‑1 (MQ=255)
                                                                                                                                            ggcggcggcGCTGGCTGGCCAGTACCTGCTGGCGCAGGCTGGTGATGTCCCAGCGCTGTTGCGCTTCATCGACAAATGCCTGCAGCAGCACCTGGTGGTCCGGGCTGCCGGTGAGGTTCTCCCGCTCGTGCGGGTCGCGGCTCAGGTCa    <  1:221163/149‑1 (MQ=255)
                                                                                                                                             gcggcggcgCTGGCTGGCCAGTACCTGCTGGCGCAGGCTGGTGATGTCCCAGCGCTGTTGCGCTTCATCGACAAATGCCTGCAGCAGCACCTGGTGGTCCGGGCTGCCGGTGAGGTTCTCCCGCTCGTGCGGGTCGCGGCTCAGGTCa    >  2:107401/1‑148 (MQ=255)
                                                                                                                                              cggcggcgCTGGCAGGCCAGTACCTGCTGGCGCAGGCTGGTGATGTCCCAGCGCTGTTGCGCTTCATCGACAAATGCCTGCAGCAGCACCTGGTGGTCCGGGCTGCCGGTGAGGTTCTCCCGCTCGTGCGGGTCGCGGCTCAGGTCATa  <  1:153844/149‑1 (MQ=255)
                                                                                                                                                  ggcgCTGGCTGGCCAGTACCTGCTGGCGCAGGCTGGTGATGTCCCAGCGCTGTTGCGCTTCATc                                                                                   <  8:175058/64‑1 (MQ=255)
                                                                                                                                                  ggcgCTGGCTGGCCAGTACCTGCTGGCGCAGGCTGGTGATGTCCCAGCGCTGTTGCGCTTCATc                                                                                   >  7:175058/1‑64 (MQ=255)
                                                                                                                                                   |                                                                                                                                               
GGCGGGCTGTGGATAACGTGCCTTGCGTTCGAGGTCATCGAGATCGATGTGGTTGCGCATGTACTGTTGGCTGGCGTCCACCAGTGGTTGGTGGTCCCAGCTTTTCAGCTTGCCGATGGCCAGCGCTTCGGCCACCAGCCGGCGGCGGCGCTGGCTGGCCAGTACCTGCTGGCGCAGGCTGGTGATGTCCCAGCGCTGTTGCGCTTCATCGACAAATGCCTGCAGCAGCACCTGGTGGTCCGGGCTGCCGGTGAGGTTCTCCCGCTCGTGCGGGTCGCGGCTCAGGTCATA  >  NC_002947/87812‑88102

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