Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 86,060 A→G 5.8% W220R (TGG→CGG)  PP_0075 ← choline sulfate transporter

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_00294786,0600AG 5.8% 206.2 / 2.2 69W220R (TGG→CGG) PP_0075choline sulfate transporter
Reads supporting (aligned to +/- strand):  ref base A (34/30);  new base G (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 = 9.36e-01

CCCAGCATGCCTACCGCAACGGCGCTGGGCAGCAGACGCAACATGAGCTCCATGTCCAGCAAAGGCAGCGGGCTGAAAGGCGGGAGTTGCCCGGTGAACGCTGCCACCCGCGGTACATGGCCAAAGAAGCCCGGCAGCAGCCAGGCCAGCAAGCTGACCAGGATCAGGCTTATCAACAGGCTCGGCCAGCGTGGGCGCCAGAGTTTGAAGGCCGCACCGATCACGACGGTGGCCAGGCCCAGCAGCAAGGAAGGTAGGTCGACCTCGCCTGCGTGGCTGGCCAGGTCCTGC  >  NC_002947/85918‑86208
                                                                                                                                               |                                                                                                                                                     
cccAGCATGCCTACCGCAACGGCGCTGGGCAGCAGACGCAACATGAGCTCCATGTCCAGCAAAGGCAGCGGGCTGAAAGGCGGGAGTTGCCCGGTGAACGCTGCCACCCGCGGTACATGGCCAAAGAAGCCCGGCAGCAgccaggccag                                                                                                                                                 >  2:270783/1‑149 (MQ=255)
   aGCATGCCTACCCCACCGGCGCCGGGCACCACACGCAACATCACCTCCCTGTCCAGCAAAGGCACCGGGCTGAAAGGCGGGAGCTGCCCGGTGAACGCTGCCACCCCCGGTACATGGCCAAAGAAGCCCGGCAGCAGCCAGGCCAGCaa                                                                                                                                              <  8:241536/149‑1 (MQ=255)
      aTGCCTACCGCAACGGCGCTGGGCAGCAGACGCAACATGAGCTCCATGTCCAGCAAAGGCAGCGGGCTGAAAGGCGGGAGTTGCCCGGTGAACGCTGCCACTCGCGGTACATGGCCAAAGAAGCCCGGCAGCAGCCAGGCCAGCAAGCt                                                                                                                                           <  1:72327/149‑1 (MQ=255)
          cTACCGCAACGGCGCTGGGCAGCAGACGCAACATGAGCTCCATGTCCAGCAAAGGCAGCGGGCTGAAAGGCGGGAGTTGCCCGGTGAACGCTGCCACCCGCGGTACATGGCCAAAGAAGCCCGGCAGCAGCCAGGCCAGCAAGCTGAc                                                                                                                                        >  6:137566/1‑148 (MQ=255)
                 aaCGGCGCTGGGCAGCAGACGCAACATGAGCTCCATGTCCAGCAAAGGCAGCGGGCTGAAAGGCGGGAGTTGCCCGGTGAACGCTGCCACCCGCGGTACATGGCCAAAGAAGCCCGGCAGCAGCCAGGCCAGCAAGCTGACCAGGATCa                                                                                                                                <  7:86599/149‑1 (MQ=255)
                  aCGGCGCTGGGCAGCAGACGCAACATGAGCTCCATGTCCAGCAAAGGCAGCGGGCTGAAAGGCGGGAGTTGCCCGGTGAACGCTGCCACCCGCGGTACATGGCCAAAGAAGCCCGGCAGCAGCCAGGCCAGCAAGCTGACCAGGATCAg                                                                                                                               >  7:280479/1‑149 (MQ=255)
                   cGGCGCTGGGCAGCAGACGCAACATGAGCTCCATGTCCAGCAAAGGCAGCGGGCTGAAAGGCGGGAGTTGCCCGGTGAACGCTGCCACCCGCGGTACATGGCCAAAGAAGCCCGGCAGCAGCCAGGCCAGCAAGCTGACCAGGATCAgg                                                                                                                             >  5:134822/1‑149 (MQ=255)
                    ggCGCTGGGCAGCAGACGCAACATGAGCTCCATGTCCAGCAAAGGCAGCGGGCTGAAAGGCGGGAGTTGCCCGGTGAACGCTGCCACCCGCGGTACATGGCCAAAGAAGCCCGGCAGCAGCCAGGCCAGCAAGCTGACCAGGATCAGGc                                                                                                                            <  7:279456/149‑1 (MQ=255)
                      cgcTGGGCAGCAGACGCAACATGAGCTCCATGTCCAGCAAAGGCAGCGGGCTGAAAGGCGGGAGTTGCCCGGTGAACGCTGCCACCCGCGGTACATGGCCAAAGAAGCCCGGCAGCAGCCAGGCCAGCAAGCTGACCAGGATCAGGCtt                                                                                                                          <  6:134822/149‑1 (MQ=255)
                            gcagcaGACGCAACATGAGCTCCATGTCCAGCAAAGGCAGCGGGCTGAAAGGCGGGAGTTGCCCGGTGAACGCTGCCACCCGCGGTACATGGCCAAAGAAGCCCGGCAGCAGCCAGGCCAGCAAGCTGAcc                                                                                                                                       >  8:36028/1‑131 (MQ=255)
                            gcagcaGACGCAACATGAGCTCCATGTCCAGCAAAGGCAGCGGGCTGAAAGGCGGGAGTTGCCCGGTGAACGCTGCCACCCGCGGTACATGGCCAAAGAAGCCCGGCAGCAGCCAGGCCAGCAAGCTGAcc                                                                                                                                       <  7:36028/131‑1 (MQ=255)
                            gcagcaGACGCAACATGAGCTCCATGTCCAGCAAAGGCAGCGGGCTGAAAGGCGGGAGTTGCCCGGTGAACGCTGCCACCCGCGGTACATGGCCAAAGAAGCCCGGCAGCAGCCAGGCCAGCAAGCTGACCAGGATCAGGCTTATCaa                                                                                                                     >  1:114285/1‑148 (MQ=255)
                             cagcagACGCAACATGAGCTCCATGTCCAGCAAAGGCAGCGGGCTGAAAGGCGGGAGTTGCCCGGTGAACGCTGCCACCCGCGGTACATGGCCAAAGAAGCCCCGGCAGCAGCCAGGCCAGCAAGCTGACCAGGATCAGGCTTATCAAc                                                                                                                    <  1:75944/149‑1 (MQ=255)
                             cagcagACGCAACATGAGCTCCATGTCCAGCAAAGGCAGCGGGCTGAAAGGCGGGAGTTGCCCGGTGAACGCTGCCACCCGCGGTACATGGCCAAAGAAGCCCGGCAGCAGCCAGGCCAGCAAGCTGACCAGGATCAGGCTTATCAACa                                                                                                                   <  8:280479/149‑1 (MQ=255)
                              agcagACGCAACATGAGCTCCATGTCCAGCAAAGGCAGCGGGCTGAAAGGCGGGAGTTGCCCGGTGAACGCTGCCACCCGCGGTACATGGCCAAAGAAGCCCGGCAGCAGCCAGGCCAGCAAGCTGAc                                                                                                                                        >  4:178235/1‑128 (MQ=255)
                              agcagACGCAACATGAGCTCCATGTCCAGCAAAGGCAGCGGGCTGAAAGGCGGGAGTTGCCCGGTGAACGCTGCCACCCGCGGTACATGGCCAAAGAAGCCCGGCAGCAGCCAGGCCAGCAAGCTGAc                                                                                                                                        <  3:178235/128‑1 (MQ=255)
                                 agACGCAACATGAGCTCCATGTCCAGCAAAGGCAGCGGGCTGAAAGGCGGGAGTTGCCCGGTGAACGCTGCCACCCGCGGTACATGGCCAAAGAAGCCCGGCAGCAGCCAGGCCAGCAAGCTGACCAGGATCAGGCTTATCAACAGGCt                                                                                                               >  8:23258/1‑149 (MQ=255)
                                 agACGCAACATGAGCTCCATGTCCAGCAAAGGCAGCGGGCTGAAAGGCGGGAGTTGCCCGGTGAACGCTGCCACCCGCGGTACATGGCCAAAGAAGCCCGGCAGCAGCCAGGCCAGCAAGCTGACCAGGATCAGGCTTATCAACAGGCt                                                                                                               >  5:63552/1‑149 (MQ=255)
                                       aaCATGAGCTCCATGTCCAGCAAAGGCAGCGGGCTGAAAGGCGGGAGTTGCACGGTGAACGCTGCCACCCGCGGTACATGGCCAAAGAAGCCCGGCAGCAGCCAGGCCAGCAAGCTGACCAGGATCAGGCTTATCAACAGGCTCgga                                                                                                           >  8:132582/1‑146 (MQ=255)
                                        aCATGAGCTCCATGTCCAGCAAAGGCAGCGGGCTGAAAGGCGGGAGTTGCCCGGTGAACGCTGCCACCCGCGGTACATGGCCAAAGAAGCCCGGCAGCAGCCAGGCCAGCAAGCTGACCAGGATCAGGCTTATCAACAGGCTCGGCCAg                                                                                                        >  4:271062/1‑149 (MQ=255)
                                         cATGAGCTCCATGTCCAGCAAAGGCAGCGGGCTGAAAGGCGGGAGTTGCCCGGTGAACGCTGCCACCCGCGGTACATGGCCAAAGAAGCCCGGCAGCAGCCAGGCCAGCAAGCTGACCAGGATCAGGCTTATCAACAGGCTCGGCCAGc                                                                                                       <  2:114285/149‑1 (MQ=255)
                                           tGAGCTCCATGTCCAGCAAAGGCAGCGGGCTGAAAGGCGGGAGTTGCCCGGTGAACGCTGCCACCCGCGGTACATGGCCAAAGAAGCCCGGCAGCAGCCAGGCCAGCAAGCTGACCAGGATCAGGCTTATCAACAGGCTCGGCCAGCGt                                                                                                     >  7:4988/1‑149 (MQ=255)
                                              gCTCCATGTCCAGCAAAGGCAGCGGGCTGAAAGGCGGGAGTTGCCCGGTGAACGCTGCCACCCGCGGTACATGGCCAAAGAAGCCCGGCAGCAGCCAGGCCAGCAAGCTGACCAGGa                                                                                                                                   >  7:285580/1‑117 (MQ=255)
                                              gCTCCATGTCCAGCAAAGGCAGCGGGCTGAAAGGCGGGAGTTGCCCGGTGAACGCTGCCACCCGCGGTACATGGCCAAAGAAGCCCGGCAGCAGCCAGGCCAGCAAGCTGACCAGGa                                                                                                                                   <  8:285580/117‑1 (MQ=255)
                                              gCTCCATGTCCAGCAAAGGCAGCGGGCTGAAAGGCGGGAGTTGCCCGGTGAACGCTGCCACCCGCGGTACATGGCCAAAGAAGCCCGGCAGCAGCCAGGCCAGCAAGCTGACCAGGATCAGGCTTATCAAc                                                                                                                    <  2:186089/131‑1 (MQ=255)
                                              gCTCCATGTCCAGCAAAGGCAGCGGGCTGAAAGGCGGGAGTTGCCCGGTGAACGCTGCCACCCGCGGTACATGGCCAAAGAAGCCCGGCAGCAGCCAGGCCAGCAAGCTGACCAGGATCAGGCTTATCAAc                                                                                                                    >  1:186089/1‑131 (MQ=255)
                                                tCCATGTCCAGCAAAGGCAGCGGGCTGAAAGGCGGGAGTTGCCCGGTGAACGCTGCCACCCGCGGTACATGGCCAAAGAAGCCCGGCAGCAGCCAGGCCAGCAAGCTGACCAGGATCAGGCTTATCAACAGGCTCGGCCAGCGTGGgcg                                                                                                >  6:72446/1‑149 (MQ=255)
                                                 ccATGTCCAGCAAAGGCAGCGGGCTGAAAGGCGGGAGTTGCCCGGTGAACGCTGCCACCCGCGGTACATGGCCAAAGAAGCCCGGCAGCAGCCAGGCCAGCAAGCTGACCAGGATCAGGCTTATCAACAGGCTCGGCCAGCGTGGgtgc                                                                                               >  6:89242/1‑149 (MQ=255)
                                                  cATGTCCAGCAAAGGCAGCGGGCTGAAAGGCGGGAGTTGCCCGGTGAACGCTGCCACCCGCGGTACATGGCCAAAGAAGCCCGGCAGCAGCCAGGCCAGCAAGCTGACCAGGATCAGGCTTATCAACAGGCTCgg                                                                                                            <  6:67259/135‑1 (MQ=255)
                                                  cATGTCCAGCAAAGGCAGCGGGCTGAAAGGCGGGAGTTGCCCGGTGAACGCTGCCACCCGCGGTACATGGCCAAAGAAGCCCGGCAGCAGCCAGGCCAGCAAGCTGACCAGGATCAGGCTTATCAACAGGCTCgg                                                                                                            >  5:67259/1‑135 (MQ=255)
                                                      tCCAGCAAAGGCAGCGGGCTGAAAGGCGGGAGTTGCCCGGTGAACGCTGCCACCCGCGGTACATGGCCAAAGAAGCCCGGCAGCAGCCAGGCCAGCAAGCTGCCCAGGATCAGGCTTATCAACAGGCTCGGCCAGCGTGGGCGCCagag                                                                                          >  5:249891/1‑149 (MQ=255)
                                                        cAGCAAAGGCAGCGGGCTGAAAGGCGGGAGTTGCCCGGTGAACGCTGCCACCCGCGGTACATGGCCAAAGAAGCCCGGCAGCAGCCAGGCCAGCAAGCTGACCAGGATCAGGCTTATCAACAGGCTCGGCCAGCGTGGGTGCCAGAGt                                                                                         <  5:89242/148‑1 (MQ=255)
                                                        cAGCAAAGGCAGCGGGCTGAAAGGCGGGAGTTGCCCGGTGAACGCTGCCACCCGCGGTACATGGCCAAAGAAGCCCGGCAGCAGCCAGGCCAGCAAGCTGACCAGGATCAGGCTTATCAACAGGCTCGGCCAGCGTGGGCGCCAGAGt                                                                                         <  8:4988/148‑1 (MQ=255)
                                                          gCAAAGGCAGCGGGCTGAAAGGCGGGAGTTGCCCGGTGAACGCTGCCACCCGCGGTACATGGCCAAAGAAGCCCGGCAGCAGCCAGGCCAGCAAGCTGACCAGGATCAGGCTTATCAACAGGCTCGGCCAGCGTGGGCGCCAGAGTTTg                                                                                      <  7:261846/149‑1 (MQ=255)
                                                           cAAAGGCAGCGGGCTGAAAGGCGGGAGTTGCCCGGTGAACGCTGCCACCCGCGGTACATGGCCAAAGAAGCCCGGCAGCAGCCAGGCCAGCAAGCTGACCAGGATCAGGCTTATCAACAGGCTCGGCCAGCGTGGGCGCCAGAGTTTg                                                                                      >  7:193054/1‑148 (MQ=255)
                                                             aaGGCAGCGGGCTGAAAGGCGGGAGTTGCCCGGTGAACGCTGCCACCCGCGGTACATGGCCAAAGAAGCCCGGCAGCAGCCAGGCCAGCAAGCTGACCAGGATCAGGCTTATCAACAGGCTCGGCCAGCGTGGGCGCCAGAGTTTGAAg                                                                                   <  5:137566/149‑1 (MQ=255)
                                                                  aGCGGGCTGAAAGGCGGGAGTTGCCCGGTGAACGCTGCCACCCGCGGTACATGGCCAAAGAAGCCCGGCAGCAGCCAGGCCAGCAAGCTGACCAGGATCAGGCTTATCAACAGGCTCGGCCAGCGTGGGCGCCAGAGTTTGAAGGCCg                                                                               >  2:82444/1‑148 (MQ=255)
                                                                       gCTGAAAGGCGGGAGTTGCCCGGTGAACGCTGCCACCCGCGGTACATGGCCAAAGAAGCCCGGCAGCAGCCAGGCCAGCAAGCTGACCAGGATCAGGCTTATCAACAGGCTCGGCCAGCGTGGGCGCCAGAGTTTGAAGGCCGCACCGa                                                                         <  4:31732/149‑1 (MQ=255)
                                                                        cTGAAAGGCGGGAGTTGCCCGGTGAACGCTGCCACCCGCGGTACATGGCCAAAGAAGCCCGGCAGCAGCCAGGCCAGCAAGCTGACCAGGATCAGGCTTATCAACAGGCTCGGCCAGCGTGGGCGCCAGAGTTTGAAGGACGCACCGAt                                                                        >  7:299585/1‑149 (MQ=255)
                                                                           aaaGGCGGGAGTTGCCCGGGGAACGCTGCCACCCGCGGTACATGGCCAAAGAAGCCCGGCAGCAGCCAGGCCAGCAAGCTGACCAGGATCAGGCTTATCAACAGGCTCGGCCAGCGTGGGCGCCAGAGTTTGAAGGCCGCACCGATCac                                                                     <  8:252257/149‑1 (MQ=255)
                                                                                  ggAGTTGCCCGGTGAACGCTGCCACCCGCGGTACATGGCCAAAGAAGCCCGGCAGCAGCCAGGCCAGCAAGCTGACCAGGATCAGGCTTATCAACAGGCTCGGCCAGCGTGGGCGCCAGAGTTTGAAGGCCGCACCGATCACGACGGTg                                                              <  1:82444/149‑1 (MQ=255)
                                                                                         cccGGTGAACGCTGCCACCCGCGGTACATGGCCAAAGAAGCCCGGCAGCAGCCAGGCCAGCAAGCTGACCAGGATCAGGCTTATCAACAGGCTCGGCCAGCGTGGGCGCCAGAGTTTGAAGGCCGCACCGATCACGACGGTGGCCAGGc                                                       >  6:152958/1‑149 (MQ=255)
                                                                                          ccGGTTAAGGCTGCCACCCGCGGTACTTGGCAAACGACGCCCGGCAGCAGCCCGGCCTCCACGCTGACCTGCATCAGGCTTATCAACAGGCTCGGCCAGCGTGGGCGCCATAGTTTGATGGCCGCACCGATCACGCCGGTGGCCAGGcc                                                      <  8:193054/149‑1 (MQ=255)
                                                                                           cGGTGAACGCTGCCACCCGCGGTACATGGCCAAAGAAGCCCGGCAGCAGCCAGGCCAGCAAGCTGACCAGGATCAGGCTTATCAACAGGCTCGGCCAGCGTGGGCGCCAGAGTTTGAAGGCCGCACCGATCACGACGGTGGCCAGGccc                                                     >  2:114590/1‑149 (MQ=255)
                                                                                             gTGAACGCTGCCACCCGCGGTACATGGCCAAAGAAGCCCGGCAGCAGCCAGGCCAGCAAGCTGACCAGGATCAGGCTTATCAACAGGCTCGGCCAGCGTGGGCGCCAGAGTTTGAAGGCCGCACCGATCACGACGGTGGCCAGGCCcag                                                   >  3:180324/1‑149 (MQ=255)
                                                                                                   gCTGCCACCCGCGGTACATGGCCAAAGAAGCCCGGCAGCAGCCAGGCCAGCAAGCTGACCAGGATCAGGCTTATCAACAGGCTCGGCCAGCGTGGGCGCCAGAGTTTGAAGGCCGCACCGATCACGACGGTGGCCAGGCCCagcagc                                               >  2:170139/1‑147 (MQ=255)
                                                                                                   gCTGCCACCCGCGGTACATGGCCAAAGAAGCCCGGCAGCAGCCAGGCCAGCAAGCTGACCAGGATCAGGCTTATCAACAGGCTCGGCCAGCGTGGGCGCCAGAGTTTGAAGGCCGCACCGATCACGACGGTGGCCAGGCCCAGCAGCaa                                             <  1:170139/149‑1 (MQ=255)
                                                                                                    cTGCCACCCGCGGTACATGGCCAAAGAAGCCCGGCAGCAGCCAGGCCAGCAAGCTGACCAGGATCAGGCTTATCAACAGGCTCGGCCAGCGTGGGCGCCAGAGTTTGAAGGCCGCACCGATCACGACGGTGGCCAGGCCCAGCAGCaag                                            <  6:63552/149‑1 (MQ=255)
                                                                                                     tGCCACCCGCGGTACATGGCCAAAGAAGCCCGGCAGCAGCCAGGCCAGCAAGCTGACCAGGATCAGGCTTATCAACAGGCTCGGCCAGCGTGGGCGCCAGAGTTTGAAGGCCGCACCGATCACGACGGTGGCCAGGCCCAGCAGCaag                                            >  4:41075/1‑148 (MQ=255)
                                                                                                       ccACCCGCGGTACATGGCCAAAGAAGCCCGGCAGCAGCCAGGCCAGCAAGCTGACCAGGATCAGGCTTATCAACAGGCTCGGCCAGCGTGGGCGCCAGAGTTTGAAGGCCGCACCGATCACGACGGTGGCCAGGCCCAGCAGCaaggaa                                         <  5:72446/149‑1 (MQ=255)
                                                                                                       ccACCCGCGGTACATGGCCAAAGAAGCCCGGCAGCAGCCAGGCCAGCAAGCTGACCAGGATCAGGCTTATCAACAGGCTCGGCCAGCGTGGGCGCCAGAGTTTGAAGGCCGCACCGATCACGACGGTGGCCAGGCCCAGCAGCaagg                                           <  7:23258/147‑1 (MQ=255)
                                                                                                        cACCCGCGGTACATGGCCAAAGAAGCCCGGCAGCAGCCAGGCCAGCAAGCTGACCAGGATCAGGCTTATCAACAGGCTCGGCCAGCGTGGGCGCCAGAGTTTGAAGGCCGCACCGATCACGACGGTGGCCAGGCCCAGCAGCaaggaag                                        <  1:270783/149‑1 (MQ=255)
                                                                                                            cgcgGTACATGGCCAAAGAAGCCCGGCAGCAGCCAGGTCAGCAAGCTGACCAGGATCAGGCTTATCAACAGGCTCGGCCAGCGTGGGCGCCAGAGTTTGAAGGCCGCACCGATCACGACGGTGGCCAGGCCCAGCAGCAAGGAaggtag                                    <  6:279607/149‑1 (MQ=255)
                                                                                                            cgcgGTACATGGCCAAAGAAGCCCGGCAGCAGCCAGGCCAGCAAGCTGACCAGGATCAGGCTTATCAACAGGCTCGGCCAGCGTGGGCGCCAGAGTTTGAAGGCCGCACCGATCACGACGGTGGCCAGGCCCAGCAGCAAGGAaggtag                                    >  3:45028/1‑149 (MQ=255)
                                                                                                                   cATGGCCAAAGAAGCCCGGCAGCAGCCAGGCCAGCAAGCTGACCAGGATCAGGCTTATCAACAGGCTCGGCCAGCGTGGGCGCCAGAGTTTGAAGGCCGCACCGATCACGACGGTGGCCAGGCCCAGCAGCAAGGAAGGTAGGTCGAc                              <  3:271062/148‑1 (MQ=255)
                                                                                                                       gCCAAAGAAGCCCGGCAGCAGCCAGGCCAGCAAGCTGACCAGGATCAGGCTTATCAACAGGCTCGGCCAGCGTGGGCGCCAGAGTTTGAAGGCCGCACCGATCACGACGGTGGCCAGGCCCAGCAGCAAGGAAGGTAGGTCGACCTCGc                         >  3:287184/1‑149 (MQ=255)
                                                                                                                        ccAAAGAAGCCCGGCAGCAGCCAGGCCAGCAAGCTGACCAGGATCAGGCTTATCAACAGGCTCGGCCAGCGTGGGCGCCAGAGTTTGAAGGCCGCACCGATCACGACGGTGGCCAGGCCCAGCAGCAAGGAAGGTAGGTCGACCTCGcc                        <  2:317001/149‑1 (MQ=255)
                                                                                                                        caaaaGAAGCCCGGCAGCAGCCATGCCAGCAAGCTCAACAGGATCAGCCTTATCAACAGCCTCAGCCAGCGTGGGCGCCAGAGTTTGAACGCCGAACCCCTCACGCAGGTCGCCAGGCCCAGCAGCAAGCAACGTAGGACCACCTCGc                         >  7:126750/3‑148 (MQ=255)
                                                                                                                         cAAAGAAGCCCGGCAGCAGCCAGGCCAGCAAGCTGACCAGGATCAGTGCTTATCAACAGGCTCGGCCAGCGTGGGCGCCAGAGTTTGAAGGCCGCACCGATCACGACGGTGGCCAGGCCCAGCAGCAAGGAAGGTAGGTCGACCTCGcc                        >  1:70854/1‑149 (MQ=255)
                                                                                                                         cAAAGAAGCCCGGCAGCAGCCAGGCCAGCAAGCTGACCAGGATCAGGCTTATCAACAGGCTCGGCCAGCGTGGGCGCCAGAGTTTGAAGGCCGCACCGATCACGCCGGTTGCCAGGCCCAGCAGCAAGGAagggaggt                                  >  4:30386/1‑138 (MQ=255)
                                                                                                                         cAAAGAAGCCCGGCAGCAGCCAGGCCAGCAAGCTGACCAGGATCAGGCTTATCAACAGGCTCGGCCAGCGTGGGCGCCAGAGTTTGAAGGCCGCACCGATCACGACGGTGGCCAGGCCCAGCAGCaag                                            >  6:55038/1‑128 (MQ=255)
                                                                                                                         cAAAGAAGCCCGGCAGCAGCCAGGCCAGCAAGCTGACCAGGATCAGGCTTATCAACAGGCTCGGCCAGCGTGGGCGCCAGAGTTTGAAGGCCGCACCGATCACGACGGTGGCCAGGCCCAGCAGCaag                                            <  5:55038/128‑1 (MQ=255)
                                                                                                                         cAAAGAAGCCCGGCAGCAGCCAGGCCAGCAAGCTGACCAGGATCAGGCTTATCAACAGGCTCGGCCAGCGTGGGCGCCAGAGTTTGAAGGCCGCACCGATCACGACGGTGGCCAGGCCCAGCAGCAAGGAaggtaggt                                  <  3:30386/138‑1 (MQ=255)
                                                                                                                          ttgctggCCTGGCTGCTGCCGGGCCAGCAAGCTGACCAGGATCAGGCTTATCAACAGGCTCGGCCAGCGTGGGCGCCAGAGTTTGAAGGCCGCACCGATCACGACGGTGGCCAGGCCCAGCAGCaaggaag                                        >  5:85673/7‑131 (MQ=255)
                                                                                                                          ttgctggCCTGGCTGCTGCCGGGCCAGCAAGCTGACCAGGATCAGGCTTATCAACAGGCTCGGCCAGCGTGGGCGCCAGAGTTTGAAGGCCGCACCGATCACGACGGTGGCCAGGCCCAGCAGCaaggaag                                        <  6:85673/125‑1 (MQ=255)
                                                                                                                          aaaGAAGCCCGGCAGCAGCCAGGCCAGCAAGCTGACCAGGATCAGGCTTATCAACAGGCTCGGCCAGCGTGGGCGCCAGAGTTTGAAGGCCGCACCGATCACGACGGTGGCCAGGCCCAGCAGCAAGGAAGGTAGGTCGACCTCGCCTg                      >  3:89283/1‑149 (MQ=255)
                                                                                                                              tggCCTGGCTGCTGCCGGGCCAGCAAGCTGACCAGGATCAGGCTTATCAACAGGCTCGGCCAGCGTGGGCGCCAGAGTTTGAAGGCCGCACCGATCACGACGGTGGCCAGGCCCAGCAGCAAGGAaggt                                      >  5:248755/3‑129 (MQ=255)
                                                                                                                              tggCCTGGCTGCTGCCGGGCCAGCAAGCTGACCAGGATCAGGCTTATCAACAGGCTCGGCCAGCGTGGGCGCCAGAGTTTGAAGGCCGCACCGATCACGACGGTGGCCAGGCCCAGCAGCAAGGAaggt                                      <  6:248755/127‑1 (MQ=255)
                                                                                                                                  cccGGCAGCAGCCAGGCCAGCAAGCTGACCAGGATCAGGCTTATCAACAGGCTCGGCCAGCGTGGGCGCCAGAGTTTGAAGGCCGCACCGATCACGACGGTGGCCAGGCCCAGCAGCAAGGAAGGTAGGTCGACCTCGCCTGCGtggct               >  3:7300/1‑149 (MQ=255)
                                                                                                                                          cagcCAGGCCAGCAAGCTGACCAGGATCAGGCTTATCAACAGGCTCGCCCAGCGTGGGCGCCAGACTTTGAAGGCCGCACCGATCACGACGGTGGCCAGGACCAGCAGCAAGGAAGGTAGGTCGACCTCGCCTGCGTCGCTGGCCAGGt       >  8:30842/1‑149 (MQ=255)
                                                                                                                                            gccaggccagCAAGCTGACCAGGATCAGGCTTATCAACAGGCTCGGCCAGCGTGGGCGCCAGAGTTTGAAGGCCGCACCGATCACGACGGTGGCCAGGCCCAGCAGCAAGGAAGGTAGGTCGACCTCGCCTGCGTGGCTGGCCAGGTcc     <  4:180324/149‑1 (MQ=255)
                                                                                                                                             ccaggccagCAAGCTGACCAGGATCAGGCTTATCAACAGGCTCGGCCAGCGTGGGCGCCAGAGTTTGAAGGCCGCACCGATCACGACGGTGGCCAGGCCCAGCAGCAAGGAAGGTAGGTCGACCTCGCCTGCGTGGCTGGCCAGGTCCt    <  5:152958/149‑1 (MQ=255)
                                                                                                                                             ccaggccagCAAGCTGACCAGGATCAGGCTTATCAACAGGCTCGGCCAGCGTGGGCGCCAGAGTTTGAAGGCCGCACCGATCACGACGGTGGCCAGGCCCAGCAGCAAGGAAGGTAGGTCGACCTCGCCTGCGTGGCTGGCCAGGTCCt    >  6:318248/1‑149 (MQ=255)
                                                                                                                                               aggccagCAAGCTGACCAGGATCAGGCTTATCAACAGGCTCGGCCAGCGTGGGCGCCAGAGTTTGAAGGCCGCACCGATCACGACGGTGGCCAGGCCCAGCAGCAAGGAAGGTAGGTCGACCTCGCCTGCGTGGCTGGCCAGGTCCTGc  >  7:134850/1‑149 (MQ=255)
                                                                                                                                               |                                                                                                                                                     
CCCAGCATGCCTACCGCAACGGCGCTGGGCAGCAGACGCAACATGAGCTCCATGTCCAGCAAAGGCAGCGGGCTGAAAGGCGGGAGTTGCCCGGTGAACGCTGCCACCCGCGGTACATGGCCAAAGAAGCCCGGCAGCAGCCAGGCCAGCAAGCTGACCAGGATCAGGCTTATCAACAGGCTCGGCCAGCGTGGGCGCCAGAGTTTGAAGGCCGCACCGATCACGACGGTGGCCAGGCCCAGCAGCAAGGAAGGTAGGTCGACCTCGCCTGCGTGGCTGGCCAGGTCCTGC  >  NC_002947/85918‑86208

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: