Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 86,053 A→T 5.8% L222Q (CTG→CAG)  PP_0075 ← choline sulfate transporter

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_00294786,0530AT 5.8% 212.0 / 2.8 69L222Q (CTG→CAG) PP_0075choline sulfate transporter
Reads supporting (aligned to +/- strand):  ref base A (34/31);  new base T (2/2);  total (36/33)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.73e-01

CCCAGCATGCCTACCGCAACGGCGCTGGGCAGCAGACGCAACATGAGCTCCATGTCCAGCAAAGGCAGCGGGCTGAAAGGCGGGAGTTGCCCGGTGAACGCTGCCACCCGCGGTACATGGCCAAAGAAG‑CCCGGCAGCAGCCAGGCCAGCAAGCTGACCAGGATCAG‑GCTTATCAACAGGCTCGGCCAGCGTGGGCGCCAGAGTTTGAAGGCCGCACCGATCACGACGGTGGCCAGGCCCAGCAGCAAGGAAGGTAGGTCGACCTCGCCTGCGTGGCT  >  NC_002947/85918‑86195
                                                                                                                                        |                                                                                                                                               
cccAGCATGCCTACCGCAACGGCGCTGGGCAGCAGACGCAACATGAGCTCCATGTCCAGCAAAGGCAGCGGGCTGAAAGGCGGGAGTTGCCCGGTGAACGCTGCCACCCGCGGTACATGGCCAAAGAAG‑CCCGGCAGCAgccaggccag                                                                                                                                    >  2:270783/1‑149 (MQ=255)
   aGCATGCCTACCCCACCGGCGCCGGGCACCACACGCAACATCACCTCCCTGTCCAGCAAAGGCACCGGGCTGAAAGGCGGGAGCTGCCCGGTGAACGCTGCCACCCCCGGTACATGGCCAAAGAAG‑CCCGGCAGCAGCCAGGCCAGCaa                                                                                                                                 <  8:241536/149‑1 (MQ=255)
      aTGCCTACCGCAACGGCGCTGGGCAGCAGACGCAACATGAGCTCCATGTCCAGCAAAGGCAGCGGGCTGAAAGGCGGGAGTTGCCCGGTGAACGCTGCCACTCGCGGTACATGGCCAAAGAAG‑CCCGGCAGCAGCCAGGCCAGCAAGCt                                                                                                                              <  1:72327/149‑1 (MQ=255)
          cTACCGCAACGGCGCTGGGCAGCAGACGCAACATGAGCTCCATGTCCAGCAAAGGCAGCGGGCTGAAAGGCGGGAGTTGCCCGGTGAACGCTGCCACCCGCGGTACATGGCCAAAGAAG‑CCCGGCAGCAGCCAGGCCAGCAAGCTGAc                                                                                                                           >  6:137566/1‑148 (MQ=255)
                 aaCGGCGCTGGGCAGCAGACGCAACATGAGCTCCATGTCCAGCAAAGGCAGCGGGCTGAAAGGCGGGAGTTGCCCGGTGAACGCTGCCACCCGCGGTACATGGCCAAAGAAG‑CCCGGCAGCAGCCAGGCCAGCAAGCTGACCAGGATCa                                                                                                                   <  7:86599/149‑1 (MQ=255)
                  aCGGCGCTGGGCAGCAGACGCAACATGAGCTCCATGTCCAGCAAAGGCAGCGGGCTGAAAGGCGGGAGTTGCCCGGTGAACGCTGCCACCCGCGGTACATGGCCAAAGAAG‑CCCGGCAGCAGCCAGGCCAGCAAGCTGACCAGGATCAg                                                                                                                  >  7:280479/1‑149 (MQ=255)
                   cGGCGCTGGGCAGCAGACGCAACATGAGCTCCATGTCCAGCAAAGGCAGCGGGCTGAAAGGCGGGAGTTGCCCGGTGAACGCTGCCACCCGCGGTACATGGCCAAAGAAG‑CCCGGCAGCAGCCAGGCCAGCAAGCTGACCAGGATCAg‑g                                                                                                                >  5:134822/1‑149 (MQ=255)
                    ggCGCTGGGCAGCAGACGCAACATGAGCTCCATGTCCAGCAAAGGCAGCGGGCTGAAAGGCGGGAGTTGCCCGGTGAACGCTGCCACCCGCGGTACATGGCCAAAGAAG‑CCCGGCAGCAGCCAGGCCAGCAAGCTGACCAGGATCAG‑Gc                                                                                                               <  7:279456/149‑1 (MQ=255)
                      cgcTGGGCAGCAGACGCAACATGAGCTCCATGTCCAGCAAAGGCAGCGGGCTGAAAGGCGGGAGTTGCCCGGTGAACGCTGCCACCCGCGGTACATGGCCAAAGAAG‑CCCGGCAGCAGCCAGGCCAGCAAGCTGACCAGGATCAG‑GCtt                                                                                                             <  6:134822/149‑1 (MQ=255)
                            gcagcaGACGCAACATGAGCTCCATGTCCAGCAAAGGCAGCGGGCTGAAAGGCGGGAGTTGCCCGGTGAACGCTGCCACCCGCGGTACATGGCCAAAGAAG‑CCCGGCAGCAGCCAGGCCAGCAAGCTGAcc                                                                                                                          <  7:36028/131‑1 (MQ=255)
                            gcagcaGACGCAACATGAGCTCCATGTCCAGCAAAGGCAGCGGGCTGAAAGGCGGGAGTTGCCCGGTGAACGCTGCCACCCGCGGTACATGGCCAAAGAAG‑CCCGGCAGCAGCCAGGCCAGCAAGCTGAcc                                                                                                                          >  8:36028/1‑131 (MQ=255)
                            gcagcaGACGCAACATGAGCTCCATGTCCAGCAAAGGCAGCGGGCTGAAAGGCGGGAGTTGCCCGGTGAACGCTGCCACCCGCGGTACATGGCCAAAGAAG‑CCCGGCAGCAGCCAGGCCAGCAAGCTGACCAGGATCAG‑GCTTATCaa                                                                                                        >  1:114285/1‑148 (MQ=255)
                             cagcagACGCAACATGAGCTCCATGTCCAGCAAAGGCAGCGGGCTGAAAGGCGGGAGTTGCCCGGTGAACGCTGCCACCCGCGGTACATGGCCAAAGAAGCCCCGGCAGCAGCCAGGCCAGCAAGCTGACCAGGATCAG‑GCTTATCAAc                                                                                                       <  1:75944/149‑1 (MQ=255)
                             cagcagACGCAACATGAGCTCCATGTCCAGCAAAGGCAGCGGGCTGAAAGGCGGGAGTTGCCCGGTGAACGCTGCCACCCGCGGTACATGGCCAAAGAAG‑CCCGGCAGCAGCCAGGCCAGCAAGCTGACCAGGATCAG‑GCTTATCAACa                                                                                                      <  8:280479/149‑1 (MQ=255)
                              agcagACGCAACATGAGCTCCATGTCCAGCAAAGGCAGCGGGCTGAAAGGCGGGAGTTGCCCGGTGAACGCTGCCACCCGCGGTACATGGCCAAAGAAG‑CCCGGCAGCAGCCAGGCCAGCAAGCTGAc                                                                                                                           >  4:178235/1‑128 (MQ=255)
                              agcagACGCAACATGAGCTCCATGTCCAGCAAAGGCAGCGGGCTGAAAGGCGGGAGTTGCCCGGTGAACGCTGCCACCCGCGGTACATGGCCAAAGAAG‑CCCGGCAGCAGCCAGGCCAGCAAGCTGAc                                                                                                                           <  3:178235/128‑1 (MQ=255)
                                 agACGCAACATGAGCTCCATGTCCAGCAAAGGCAGCGGGCTGAAAGGCGGGAGTTGCCCGGTGAACGCTGCCACCCGCGGTACATGGCCAAAGAAG‑CCCGGCAGCAGCCAGGCCAGCAAGCTGACCAGGATCAG‑GCTTATCAACAGGCt                                                                                                  >  5:63552/1‑149 (MQ=255)
                                 agACGCAACATGAGCTCCATGTCCAGCAAAGGCAGCGGGCTGAAAGGCGGGAGTTGCCCGGTGAACGCTGCCACCCGCGGTACATGGCCAAAGAAG‑CCCGGCAGCAGCCAGGCCAGCAAGCTGACCAGGATCAG‑GCTTATCAACAGGCt                                                                                                  >  8:23258/1‑149 (MQ=255)
                                       aaCATGAGCTCCATGTCCAGCAAAGGCAGCGGGCTGAAAGGCGGGAGTTGCACGGTGAACGCTGCCACCCGCGGTACATGGCCAAAGAAG‑CCCGGCAGCAGCCAGGCCAGCAAGCTGACCAGGATCAG‑GCTTATCAACAGGCTCgga                                                                                              >  8:132582/1‑146 (MQ=255)
                                        aCATGAGCTCCATGTCCAGCAAAGGCAGCGGGCTGAAAGGCGGGAGTTGCCCGGTGAACGCTGCCACCCGCGGTACATGGCCAAAGAAG‑CCCGGCAGCAGCCAGGCCAGCAAGCTGACCAGGATCAG‑GCTTATCAACAGGCTCGGCCAg                                                                                           >  4:271062/1‑149 (MQ=255)
                                         cATGAGCTCCATGTCCAGCAAAGGCAGCGGGCTGAAAGGCGGGAGTTGCCCGGTGAACGCTGCCACCCGCGGTACATGGCCAAAGAAG‑CCCGGCAGCAGCCAGGCCAGCAAGCTGACCAGGATCAG‑GCTTATCAACAGGCTCGGCCAGc                                                                                          <  2:114285/149‑1 (MQ=255)
                                           tGAGCTCCATGTCCAGCAAAGGCAGCGGGCTGAAAGGCGGGAGTTGCCCGGTGAACGCTGCCACCCGCGGTACATGGCCAAAGAAG‑CCCGGCAGCAGCCAGGCCAGCAAGCTGACCAGGATCAG‑GCTTATCAACAGGCTCGGCCAGCGt                                                                                        >  7:4988/1‑149 (MQ=255)
                                              gCTCCATGTCCAGCAAAGGCAGCGGGCTGAAAGGCGGGAGTTGCCCGGTGAACGCTGCCACCCGCGGTACATGGCCAAAGAAG‑CCCGGCAGCAGCCAGGCCAGCAAGCTGACCAGGa                                                                                                                      >  7:285580/1‑117 (MQ=255)
                                              gCTCCATGTCCAGCAAAGGCAGCGGGCTGAAAGGCGGGAGTTGCCCGGTGAACGCTGCCACCCGCGGTACATGGCCAAAGAAG‑CCCGGCAGCAGCCAGGCCAGCAAGCTGACCAGGa                                                                                                                      <  8:285580/117‑1 (MQ=255)
                                              gCTCCATGTCCAGCAAAGGCAGCGGGCTGAAAGGCGGGAGTTGCCCGGTGAACGCTGCCACCCGCGGTACATGGCCAAAGAAG‑CCCGGCAGCAGCCAGGCCAGCAAGCTGACCAGGATCAG‑GCTTATCAAc                                                                                                       <  2:186089/131‑1 (MQ=255)
                                              gCTCCATGTCCAGCAAAGGCAGCGGGCTGAAAGGCGGGAGTTGCCCGGTGAACGCTGCCACCCGCGGTACATGGCCAAAGAAG‑CCCGGCAGCAGCCAGGCCAGCAAGCTGACCAGGATCAG‑GCTTATCAAc                                                                                                       >  1:186089/1‑131 (MQ=255)
                                                tCCATGTCCAGCAAAGGCAGCGGGCTGAAAGGCGGGAGTTGCCCGGTGAACGCTGCCACCCGCGGTACATGGCCAAAGAAG‑CCCGGCAGCAGCCAGGCCAGCAAGCTGACCAGGATCAG‑GCTTATCAACAGGCTCGGCCAGCGTGGgcg                                                                                   >  6:72446/1‑149 (MQ=255)
                                                 ccATGTCCAGCAAAGGCAGCGGGCTGAAAGGCGGGAGTTGCCCGGTGAACGCTGCCACCCGCGGTACATGGCCAAAGAAG‑CCCGGCAGCAGCCAGGCCAGCAAGCTGACCAGGATCAG‑GCTTATCAACAGGCTCGGCCAGCGTGGgtgc                                                                                  >  6:89242/1‑149 (MQ=255)
                                                  cATGTCCAGCAAAGGCAGCGGGCTGAAAGGCGGGAGTTGCCCGGTGAACGCTGCCACCCGCGGTACATGGCCAAAGAAG‑CCCGGCAGCAGCCAGGCCAGCAAGCTGACCAGGATCAG‑GCTTATCAACAGGCTCgg                                                                                               <  6:67259/135‑1 (MQ=255)
                                                  cATGTCCAGCAAAGGCAGCGGGCTGAAAGGCGGGAGTTGCCCGGTGAACGCTGCCACCCGCGGTACATGGCCAAAGAAG‑CCCGGCAGCAGCCAGGCCAGCAAGCTGACCAGGATCAG‑GCTTATCAACAGGCTCgg                                                                                               >  5:67259/1‑135 (MQ=255)
                                                      tCCAGCAAAGGCAGCGGGCTGAAAGGCGGGAGTTGCCCGGTGAACGCTGCCACCCGCGGTACATGGCCAAAGAAG‑CCCGGCAGCAGCCAGGCCAGCAAGCTGCCCAGGATCAG‑GCTTATCAACAGGCTCGGCCAGCGTGGGCGCCagag                                                                             >  5:249891/1‑149 (MQ=255)
                                                        cAGCAAAGGCAGCGGGCTGAAAGGCGGGAGTTGCCCGGTGAACGCTGCCACCCGCGGTACATGGCCAAAGAAG‑CCCGGCAGCAGCCAGGCCAGCAAGCTGACCAGGATCAG‑GCTTATCAACAGGCTCGGCCAGCGTGGGTGCCAGAGt                                                                            <  5:89242/148‑1 (MQ=255)
                                                        cAGCAAAGGCAGCGGGCTGAAAGGCGGGAGTTGCCCGGTGAACGCTGCCACCCGCGGTACATGGCCAAAGAAG‑CCCGGCAGCAGCCAGGCCAGCAAGCTGACCAGGATCAG‑GCTTATCAACAGGCTCGGCCAGCGTGGGCGCCAGAGt                                                                            <  8:4988/148‑1 (MQ=255)
                                                          gCAAAGGCAGCGGGCTGAAAGGCGGGAGTTGCCCGGTGAACGCTGCCACCCGCGGTACATGGCCAAAGAAG‑CCCGGCAGCAGCCAGGCCAGCAAGCTGACCAGGATCAG‑GCTTATCAACAGGCTCGGCCAGCGTGGGCGCCAGAGTTTg                                                                         <  7:261846/149‑1 (MQ=255)
                                                           cAAAGGCAGCGGGCTGAAAGGCGGGAGTTGCCCGGTGAACGCTGCCACCCGCGGTACATGGCCAAAGAAG‑CCCGGCAGCAGCCAGGCCAGCAAGCTGACCAGGATCAG‑GCTTATCAACAGGCTCGGCCAGCGTGGGCGCCAGAGTTTg                                                                         >  7:193054/1‑148 (MQ=255)
                                                             aaGGCAGCGGGCTGAAAGGCGGGAGTTGCCCGGTGAACGCTGCCACCCGCGGTACATGGCCAAAGAAG‑CCCGGCAGCAGCCAGGCCAGCAAGCTGACCAGGATCAG‑GCTTATCAACAGGCTCGGCCAGCGTGGGCGCCAGAGTTTGAAg                                                                      <  5:137566/149‑1 (MQ=255)
                                                                  aGCGGGCTGAAAGGCGGGAGTTGCCCGGTGAACGCTGCCACCCGCGGTACATGGCCAAAGAAG‑CCCGGCAGCAGCCAGGCCAGCAAGCTGACCAGGATCAG‑GCTTATCAACAGGCTCGGCCAGCGTGGGCGCCAGAGTTTGAAGGCCg                                                                  >  2:82444/1‑148 (MQ=255)
                                                                       gCTGAAAGGCGGGAGTTGCCCGGTGAACGCTGCCACCCGCGGTACATGGCCAAAGAAG‑CCCGGCAGCAGCCAGGCCAGCAAGCTGACCAGGATCAG‑GCTTATCAACAGGCTCGGCCAGCGTGGGCGCCAGAGTTTGAAGGCCGCACCGa                                                            <  4:31732/149‑1 (MQ=255)
                                                                        cTGAAAGGCGGGAGTTGCCCGGTGAACGCTGCCACCCGCGGTACATGGCCAAAGAAG‑CCCGGCAGCAGCCAGGCCAGCAAGCTGACCAGGATCAG‑GCTTATCAACAGGCTCGGCCAGCGTGGGCGCCAGAGTTTGAAGGACGCACCGAt                                                           >  7:299585/1‑149 (MQ=255)
                                                                           aaaGGCGGGAGTTGCCCGGGGAACGCTGCCACCCGCGGTACATGGCCAAAGAAG‑CCCGGCAGCAGCCAGGCCAGCAAGCTGACCAGGATCAG‑GCTTATCAACAGGCTCGGCCAGCGTGGGCGCCAGAGTTTGAAGGCCGCACCGATCac                                                        <  8:252257/149‑1 (MQ=255)
                                                                                  ggAGTTGCCCGGTGAACGCTGCCACCCGCGGTACATGGCCAAAGAAG‑CCCGGCAGCAGCCAGGCCAGCAAGCTGACCAGGATCAG‑GCTTATCAACAGGCTCGGCCAGCGTGGGCGCCAGAGTTTGAAGGCCGCACCGATCACGACGGTg                                                 <  1:82444/149‑1 (MQ=255)
                                                                                         cccGGTGAACGCTGCCACCCGCGGTACATGGCCAAAGAAG‑CCCGGCAGCAGCCAGGCCAGCAAGCTGACCAGGATCAG‑GCTTATCAACAGGCTCGGCCAGCGTGGGCGCCAGAGTTTGAAGGCCGCACCGATCACGACGGTGGCCAGGc                                          >  6:152958/1‑149 (MQ=255)
                                                                                          ccGGTTAAGGCTGCCACCCGCGGTACTTGGCAAACGACG‑CCCGGCAGCAGCCCGGCCTCCACGCTGACCTGCATCAG‑GCTTATCAACAGGCTCGGCCAGCGTGGGCGCCATAGTTTGATGGCCGCACCGATCACGCCGGTGGCCAGGcc                                         <  8:193054/149‑1 (MQ=255)
                                                                                           cGGTGAACGCTGCCACCCGCGGTACATGGCCAAAGAAG‑CCCGGCAGCAGCCAGGCCAGCAAGCTGACCAGGATCAG‑GCTTATCAACAGGCTCGGCCAGCGTGGGCGCCAGAGTTTGAAGGCCGCACCGATCACGACGGTGGCCAGGccc                                        >  2:114590/1‑149 (MQ=255)
                                                                                             gTGAACGCTGCCACCCGCGGTACATGGCCAAAGAAG‑CCCGGCAGCAGCCAGGCCAGCAAGCTGACCAGGATCAG‑GCTTATCAACAGGCTCGGCCAGCGTGGGCGCCAGAGTTTGAAGGCCGCACCGATCACGACGGTGGCCAGGCCcag                                      >  3:180324/1‑149 (MQ=255)
                                                                                                   gCTGCCACCCGCGGTACATGGCCAAAGAAG‑CCCGGCAGCAGCCAGGCCAGCAAGCTGACCAGGATCAG‑GCTTATCAACAGGCTCGGCCAGCGTGGGCGCCAGAGTTTGAAGGCCGCACCGATCACGACGGTGGCCAGGCCCagcagc                                  >  2:170139/1‑147 (MQ=255)
                                                                                                   gCTGCCACCCGCGGTACATGGCCAAAGAAG‑CCCGGCAGCAGCCAGGCCAGCAAGCTGACCAGGATCAG‑GCTTATCAACAGGCTCGGCCAGCGTGGGCGCCAGAGTTTGAAGGCCGCACCGATCACGACGGTGGCCAGGCCCAGCAGCaa                                <  1:170139/149‑1 (MQ=255)
                                                                                                    cTGCCACCCGCGGTACATGGCCAAAGAAG‑CCCGGCAGCAGCCAGGCCAGCAAGCTGACCAGGATCAG‑GCTTATCAACAGGCTCGGCCAGCGTGGGCGCCAGAGTTTGAAGGCCGCACCGATCACGACGGTGGCCAGGCCCAGCAGCaag                               <  6:63552/149‑1 (MQ=255)
                                                                                                     tGCCACCCGCGGTACATGGCCAAAGAAG‑CCCGGCAGCAGCCAGGCCAGCAAGCTGACCAGGATCAG‑GCTTATCAACAGGCTCGGCCAGCGTGGGCGCCAGAGTTTGAAGGCCGCACCGATCACGACGGTGGCCAGGCCCAGCAGCaag                               >  4:41075/1‑148 (MQ=255)
                                                                                                       ccACCCGCGGTACATGGCCAAAGAAG‑CCCGGCAGCAGCCAGGCCAGCAAGCTGACCAGGATCAG‑GCTTATCAACAGGCTCGGCCAGCGTGGGCGCCAGAGTTTGAAGGCCGCACCGATCACGACGGTGGCCAGGCCCAGCAGCaaggaa                            <  5:72446/149‑1 (MQ=255)
                                                                                                       ccACCCGCGGTACATGGCCAAAGAAG‑CCCGGCAGCAGCCAGGCCAGCAAGCTGACCAGGATCAG‑GCTTATCAACAGGCTCGGCCAGCGTGGGCGCCAGAGTTTGAAGGCCGCACCGATCACGACGGTGGCCAGGCCCAGCAGCaagg                              <  7:23258/147‑1 (MQ=255)
                                                                                                        cACCCGCGGTACATGGCCAAAGAAG‑CCCGGCAGCAGCCAGGCCAGCAAGCTGACCAGGATCAG‑GCTTATCAACAGGCTCGGCCAGCGTGGGCGCCAGAGTTTGAAGGCCGCACCGATCACGACGGTGGCCAGGCCCAGCAGCaaggaag                           <  1:270783/149‑1 (MQ=255)
                                                                                                            cgcgGTACATGGCCAAAGAAG‑CCCGGCAGCAGCCAGGTCAGCAAGCTGACCAGGATCAG‑GCTTATCAACAGGCTCGGCCAGCGTGGGCGCCAGAGTTTGAAGGCCGCACCGATCACGACGGTGGCCAGGCCCAGCAGCAAGGAaggtag                       <  6:279607/149‑1 (MQ=255)
                                                                                                            cgcgGTACATGGCCAAAGAAG‑CCCGGCAGCAGCCAGGCCAGCAAGCTGACCAGGATCAG‑GCTTATCAACAGGCTCGGCCAGCGTGGGCGCCAGAGTTTGAAGGCCGCACCGATCACGACGGTGGCCAGGCCCAGCAGCAAGGAaggtag                       >  3:45028/1‑149 (MQ=255)
                                                                                                                   cATGGCCAAAGAAG‑CCCGGCAGCAGCCAGGCCAGCAAGCTGACCAGGATCAG‑GCTTATCAACAGGCTCGGCCAGCGTGGGCGCCAGAGTTTGAAGGCCGCACCGATCACGACGGTGGCCAGGCCCAGCAGCAAGGAAGGTAGGTCGAc                 <  3:271062/148‑1 (MQ=255)
                                                                                                                       gCCAAAGAAG‑CCCGGCAGCAGCCAGGCCAGCAAGCTGACCAGGATCAG‑GCTTATCAACAGGCTCGGCCAGCGTGGGCGCCAGAGTTTGAAGGCCGCACCGATCACGACGGTGGCCAGGCCCAGCAGCAAGGAAGGTAGGTCGACCTCGc            >  3:287184/1‑149 (MQ=255)
                                                                                                                        ccAAAGAAG‑CCCGGCAGCAGCCAGGCCAGCAAGCTGACCAGGATCAG‑GCTTATCAACAGGCTCGGCCAGCGTGGGCGCCAGAGTTTGAAGGCCGCACCGATCACGACGGTGGCCAGGCCCAGCAGCAAGGAAGGTAGGTCGACCTCGcc           <  2:317001/149‑1 (MQ=255)
                                                                                                                        caaaaGAAG‑CCCGGCAGCAGCCATGCCAGCAAGCTCAACAGGATCAG‑CCTTATCAACAGCCTCAGCCAGCGTGGGCGCCAGAGTTTGAACGCCGAACCCCTCACGCAGGTCGCCAGGCCCAGCAGCAAGCAACGTAGGACCACCTCGc            >  7:126750/3‑148 (MQ=255)
                                                                                                                         cAAAGAAG‑CCCGGCAGCAGCCAGGCCAGCAAGCTGACCAGGATCAGTGCTTATCAACAGGCTCGGCCAGCGTGGGCGCCAGAGTTTGAAGGCCGCACCGATCACGACGGTGGCCAGGCCCAGCAGCAAGGAAGGTAGGTCGACCTCGcc           >  1:70854/1‑149 (MQ=255)
                                                                                                                         cAAAGAAG‑CCCGGCAGCAGCCAGGCCAGCAAGCTGACCAGGATCAG‑GCTTATCAACAGGCTCGGCCAGCGTGGGCGCCAGAGTTTGAAGGCCGCACCGATCACGCCGGTTGCCAGGCCCAGCAGCAAGGAagggaggt                     >  4:30386/1‑138 (MQ=255)
                                                                                                                         cAAAGAAG‑CCCGGCAGCAGCCAGGCCAGCAAGCTGACCAGGATCAG‑GCTTATCAACAGGCTCGGCCAGCGTGGGCGCCAGAGTTTGAAGGCCGCACCGATCACGACGGTGGCCAGGCCCAGCAGCaag                               >  6:55038/1‑128 (MQ=255)
                                                                                                                         cAAAGAAG‑CCCGGCAGCAGCCAGGCCAGCAAGCTGACCAGGATCAG‑GCTTATCAACAGGCTCGGCCAGCGTGGGCGCCAGAGTTTGAAGGCCGCACCGATCACGACGGTGGCCAGGCCCAGCAGCaag                               <  5:55038/128‑1 (MQ=255)
                                                                                                                         cAAAGAAG‑CCCGGCAGCAGCCAGGCCAGCAAGCTGACCAGGATCAG‑GCTTATCAACAGGCTCGGCCAGCGTGGGCGCCAGAGTTTGAAGGCCGCACCGATCACGACGGTGGCCAGGCCCAGCAGCAAGGAaggtaggt                     <  3:30386/138‑1 (MQ=255)
                                                                                                                          ttgctgg‑CCTGGCTGCTGCCGGGCCAGCAAGCTGACCAGGATCAG‑GCTTATCAACAGGCTCGGCCAGCGTGGGCGCCAGAGTTTGAAGGCCGCACCGATCACGACGGTGGCCAGGCCCAGCAGCaaggaag                           <  6:85673/125‑1 (MQ=255)
                                                                                                                          ttgctgg‑CCTGGCTGCTGCCGGGCCAGCAAGCTGACCAGGATCAG‑GCTTATCAACAGGCTCGGCCAGCGTGGGCGCCAGAGTTTGAAGGCCGCACCGATCACGACGGTGGCCAGGCCCAGCAGCaaggaag                           >  5:85673/7‑131 (MQ=255)
                                                                                                                          aaaGAAG‑CCCGGCAGCAGCCAGGCCAGCAAGCTGACCAGGATCAG‑GCTTATCAACAGGCTCGGCCAGCGTGGGCGCCAGAGTTTGAAGGCCGCACCGATCACGACGGTGGCCAGGCCCAGCAGCAAGGAAGGTAGGTCGACCTCGCCTg         >  3:89283/1‑149 (MQ=255)
                                                                                                                              tgg‑CCTGGCTGCTGCCGGGCCAGCAAGCTGACCAGGATCAG‑GCTTATCAACAGGCTCGGCCAGCGTGGGCGCCAGAGTTTGAAGGCCGCACCGATCACGACGGTGGCCAGGCCCAGCAGCAAGGAaggt                         <  6:248755/127‑1 (MQ=255)
                                                                                                                              tgg‑CCTGGCTGCTGCCGGGCCAGCAAGCTGACCAGGATCAG‑GCTTATCAACAGGCTCGGCCAGCGTGGGCGCCAGAGTTTGAAGGCCGCACCGATCACGACGGTGGCCAGGCCCAGCAGCAAGGAaggt                         >  5:248755/3‑129 (MQ=255)
                                                                                                                                  cccGGCAGCAGCCAGGCCAGCAAGCTGACCAGGATCAG‑GCTTATCAACAGGCTCGGCCAGCGTGGGCGCCAGAGTTTGAAGGCCGCACCGATCACGACGGTGGCCAGGCCCAGCAGCAAGGAAGGTAGGTCGACCTCGCCTGCGtggct  >  3:7300/1‑149 (MQ=255)
                                                                                                                                        |                                                                                                                                               
CCCAGCATGCCTACCGCAACGGCGCTGGGCAGCAGACGCAACATGAGCTCCATGTCCAGCAAAGGCAGCGGGCTGAAAGGCGGGAGTTGCCCGGTGAACGCTGCCACCCGCGGTACATGGCCAAAGAAG‑CCCGGCAGCAGCCAGGCCAGCAAGCTGACCAGGATCAG‑GCTTATCAACAGGCTCGGCCAGCGTGGGCGCCAGAGTTTGAAGGCCGCACCGATCACGACGGTGGCCAGGCCCAGCAGCAAGGAAGGTAGGTCGACCTCGCCTGCGTGGCT  >  NC_002947/85918‑86195

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