Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,791,405 C→G 8.0% G283G (GGC→GGG PP_3353 → sulfatase‑modifying protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029473,791,4050CG 8.0% 233.8 / 7.5 75G283G (GGC→GGGPP_3353sulfatase‑modifying protein
Reads supporting (aligned to +/- strand):  ref base C (31/38);  new base G (2/4);  total (33/42)
Fisher's exact test for biased strand distribution p-value = 6.89e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.64e-01

TGGTTTTACCGGTACATCGCCCGTGGGTTGCTTCCCGGCCAACGGCTTTGGCCTGTTCGATGCCGGCGGCAATGTCTGGGAGCTGACTCGCACGGGCTATCGGCCAGGCCATGACGCACAGCGCGACGCCAAGCTCGACCCCTCAGGCCCGGCCCTGAGTGACAGCTTCGACCCGGCAGACCCCGGCGTGCCGGTGGCGGTAATCAAAGGCGGCTCGCACCTGTGTTCGGCGGACCGCTGCATGCGCTACCGCCCCTCGGCACGCCAGCCGCAGCCGGTGTTCATGACGACCTCGC  >  NC_002947/3791258‑3791553
                                                                                                                                                   |                                                                                                                                                    
tGGTTTTACCGGTCCATCGCCCGTGGGTTGCTTCCCGGCCAACGGCTTTGGCCTGTTCGATGCCGGCGGCAATGTCTGGGAGCTGACTCGCACGGGCTATCGGCCAGGCCATGACGCACAGCGCGACGCCAAGCTCGACCCCTCAggc                                                                                                                                                      <  2:213144/148‑1 (MQ=255)
           gTACATCGCCCGTGGGTTGCTTCCCGGCCAACGGCTTTGGCCTGTTCGATGCCGGCGGCAATGTCTGGGAGCTGACTCGCACGGGCTATCGGCCAGGCCATGACGCACAGCGCGACGCCAAGCTCGACCCCTCAGGCCCGGCCCTGAGt                                                                                                                                          >  3:141833/1‑149 (MQ=255)
              cATCGCCCGTGGGTTGCTTCCCGGCCAACGGCTTTGGCCTGTTCGATGCCGGCGGCAATGTCTGGGAGCTGACTCGCACGGGCTATCGGTCAGGCCATGACGCACAGCGCGACGCCAAGCTCGACCCCTCAGGCCCGGCCCTGAGTGa                                                                                                                                        <  5:358148/148‑1 (MQ=255)
              cATCGCCCGTGGGTTGCTTCCCGGCCAACGGCTTTGGCCTGTTCGATGCCGGCGGCAATGTCTGGGAGCTGACTCGCACGGGCTATCGGCCAGGCCATGACGCACAGCGCGACGCCAAGCTCGACCCCTCAGGGCCGGGCCTGAgggg                                                                                                                                        <  1:11335/148‑4 (MQ=255)
              cATCGCCCGTGGGTTGCTTCCCGGCCAACGGCTTTGGCCTGTTCGATGCCGGCGGCAATGTCTGGGAGCTGACTCGCACGGGCTATCGGCCAGGCCATGACGCACAGCGCGACGCCAAGCTCGACCCCTCAGGCGCGGCGCTGAGGGAc                                                                                                                                       >  8:95545/1‑149 (MQ=255)
               aTCGCCCGTGGGTTGCTTCCCGGCCAACGGCTTTGGCCTGTTCGATGCCGGCGGCAATGTCTGGGAGCTGACTCGCACGGGCTATCGGCCAGGCCATGACGCACAGCGCGACGCCAAGCTCGACCCCTCAGGCCCGGCCCTGAGTGa                                                                                                                                        >  6:215321/1‑147 (MQ=255)
               aTCGCCCGTGGGTTGCTTCCCGGCCAACGGCTTTGGCCTGTTCGATGCCGGCGGCAATGTCTGGGAGCTGACTCGCACGGGCTATCGGCCAGGCCATGACGCACAGCGCGACGCCAAGCTCGACCCCTCAGGCCCGGCCCTGAGTGACa                                                                                                                                      <  3:268644/149‑1 (MQ=255)
                 cGCCCGTGGGTTGCTTCCCGGCCAACGGCTTTGGCCTGTTCGATGCCGGCGGCAATGTCTGGGAGCTGACTCGCACGGGCTATCGGCCAGGCCATGACGCACAGCGCGACGCCAAGCTCGACCCCTCAGGGCCGGGCCTGAggg                                                                                                                                         >  2:23907/1‑142 (MQ=255)
                 cGCCCGTGGGTTGCTTCCCGGCCAACGGCTTTGGCCTGTTCGATGCCGGCGGCAATGTCTGGGAGCTGACTCGCACGGGCTATCGGCCAGGCCATGACGCACAGCGCGACGCCAAGCTCGACCCCTCAGGGCCGGGCCTGAggg                                                                                                                                         <  1:23907/144‑3 (MQ=255)
                   cccGTGGGTTGCTTCCCGGCCAACGGCTTTGGCCTGTTCGATGCCGGCGGCAATGTCTGGGAGCTGACTCGCACGGGCTATCGGCCAGGCCATGACGCACAGCGCGACGCCAAGCTCGACCCCTCAGGCCCGGCCCTGAGTGACAGCtt                                                                                                                                  <  5:246964/149‑1 (MQ=255)
                              cTTCCCGGCCAACGGCTTTGGCCTGTTCGATGCCGGCGGCAATGTCTGGGAGCTGACTCGCACGGGCTATCGGCCAGGCCATGACGCACAGCGCGACGCCAAGCTCGACCCCTCAGGCCCGGCCCTGAGTGACAGCTTCGACCCgg                                                                                                                          >  7:344682/1‑146 (MQ=255)
                              cTTCCCGGCCAACGGCTTTGGCCTGTTCGATGCCGGCGGCAATGTCTGGGAGCTGACTCGCACGGGCTATCGGCCAGGCCATGACGCACAGCGCGACGCCAAGCTCGACCCCTCAGGCCCGGCCCTGAGTGACAGCTTCGACCCgg                                                                                                                          >  2:303614/1‑146 (MQ=255)
                                tCCCGGCCAACGGCTTTGGCCTGTTCGATGCCGGCGGCAATGTCTGGGAGCTGACTCGCACGGGCTATCGGCCAGGCCATGACGCACAGCGCGACGCCAAGCTCGACCCCTCAGGCCCGGCCCTGAGt                                                                                                                                          <  8:247905/128‑1 (MQ=255)
                                tCCCGGCCAACGGCTTTGGCCTGTTCGATGCCGGCGGCAATGTCTGGGAGCTGACTCGCACGGGCTATCGGCCAGGCCATGACGCACAGCGCGACGCCAAGCTCGACCCCTCAGGCCCGGCCCTGAGt                                                                                                                                          >  7:247905/1‑128 (MQ=255)
                                 cccGGCCAACGGCTTTGGCCTGTTCGATGCCGGCGGCAATGTCTGGGAGCTGACTCGCACGGGCTATCGGCCAGGCCATGACGCACAGCGCGACGCCAAGCTCGACCCCTCAGGCCCGGCCCTGAGTGACAGCTTCGACCCGGCAGAcc                                                                                                                    >  6:132930/1‑149 (MQ=255)
                                  ccGGCCAACGGCTTTGGCCTGTTCGATGCCGGCGGCAATGTCTGGGAGCTGACTCGCACGGGCTATCGGCCAGGCCATGACGCACAGCGCGACGCCAAGCTCGACCCCTCAGGCCCGGCCCTGAGTGACAGCTTCGACCCGGCAGAccc                                                                                                                   >  5:35044/1‑149 (MQ=255)
                                  ccGGCCAACGGCTTTGGCCTGTTCGATGCCGGCGGCAATGTCTGGGAGCTGACTCGCACGGGCTATCGGCCAGGCCATGACGCACAGCGCGACGCCAAGCTCGACCCCTCAGGCCCGGCCCTGAGTGACAGCTTCGACCCGGCAGAccc                                                                                                                   >  5:283982/1‑149 (MQ=255)
                                  ccGGCCAACGGCTTTGGCCTGTTCGATGCCGGCGGCAATGTCTGGGAGCTGACTCGCACGGGCTATCGGCCAGGCCATGACGCACAGCGCGACGCCAAGCTCGACCCCTCAGGCCCGGCCCTGAGTGACAGCTTCGACCCGGCAGAccc                                                                                                                   >  3:234487/1‑149 (MQ=255)
                                   cGGCCAACGGCTTTGGCCTGTTCGATGCCGGCGGCAATGTCTGGGAGCTGACTCGCACGGGCTATCGGCCAGGCCATGACGCACAGCGCGACGCCAAGCTCGACCCCTCAGGCCCGGCCCTGAGTGACAGCTTCGACCCGGCAGAcccc                                                                                                                  >  2:110735/1‑149 (MQ=255)
                                      ccAACGGCTTTGGCCTGTTCGATGCCGGCGGCAATGTCTGGGAGCTGACTCGCACGGGCTATCGGCCAGGCCATGACGCACAGCGCGACGCCAAGCTCGACCCCTCAGGCCCGGCCCTGAGTGACAGCTTCGACCCGGCAGACCCCGGc                                                                                                               <  5:132930/149‑1 (MQ=255)
                                        aaCGGCTTTGGCCTGTTCGATGCCGGCGGCAATGTCTGGGAGCTGACTCGCACGGGCTATCGGCCAGGCCATGACGCACAGCGCGACGCCAAGCTCGACCCCTCAGGCCCGGCCCTGAGTGACAGCTTCGAccc                                                                                                                            >  8:188249/1‑134 (MQ=255)
                                        aaCGGCTTTGGCCTGTTCGATGCCGGCGGCAATGTCTGGGAGCTGACTCGCACGGGCTATCGGCCAGGCCATGACGCACAGCGCGACGCCAAGCTCGACCCCTCAGGCCCGGCCCTGAGTGACAGCTTCGAccc                                                                                                                            <  7:188249/134‑1 (MQ=255)
                                        aaCGGCTTTGGCCTGTTCGATGCCGGCGGCAATGTCTGGGAGCTGACTCGCACGGGCTATCGGCCAGGCCATGACGCACAGCGCGACGCCAAGCTCGACCCCTCAGGCCCGGCCCTGAGTGACAGCTTCGACCCGGCAGACCCCgggcg                                                                                                             <  8:269122/149‑4 (MQ=255)
                                        aaCGGCTTTGGCCTGTTCGATGCCGGCGGCAATGTCTGGGAGCTGACTCGCACGGGCTATCGGCCAGGCCATGACGCACAGCGCGACGCCAAGCTCGACCCCTCAGGCCCGGCCCTGAGTGACAGCTTCGACCCGGCAGACCCCGGCGt                                                                                                             <  4:141833/149‑1 (MQ=255)
                                        aaCGGCTTTGGCCTGTTCGATGCCGGCGGCAATGTCTGGGAGCTGACTCGCACGGGCTATCGGCCAGGCCATGACGCACAGCGCGACGCCAAGCTCGACCCCTCAGGCCCGGCCCTGAGTGACAGCTTCGACCCGGCAGACCCCGGCGt                                                                                                             >  2:7280/1‑149 (MQ=255)
                                           ggCTTTGGCCTGTTCGATGCCGGCGGCAATGTCTGGGAGCTGACTCGCACGGGCTATCGGCCAGGCCATGACGCACAGCGCGACGCCAAGCTCGACCCCTCAGGCCCGGCCCTGAGTGACAGCTTCGACCCGGCAGACCCCGGCGTGcc                                                                                                          <  2:173334/149‑1 (MQ=255)
                                                  gCCTGTTCGATGCCGGCGGCAATGTCTGGGAGCTGACTCGCACGGGCTATCGGCCAGGCCATGACGCACAGCGCGACGCCAAGCTCGACCCCTCAGGCCCGGCCCTGAGTGACAGCTTCGACCCGGCAGACCCCGGCGTGCCGGTGGc                                                                                                    <  1:110735/148‑1 (MQ=255)
                                                  gCCTGTTCGATGCCGGCGGCAATGTCTGGGAGCTGACTCGCACGGGCTATCGGCCAGGCCATGACGCACAGCGCGACGCCAAGCTCGACCCCTCAGGCCCGGCCCTGAGTGACAGCTTCGACCCGGCAGACCCCGGCGTGCCGGTGGCg                                                                                                   <  8:344682/149‑1 (MQ=255)
                                                  gCCTGTTCGATGCCGGCGGCAATGTCTGGGAGCTGACTCGCACGGGCTATCGGCCAGGCCATGACGCACAGCGCGACGCCAAGCTCGACCCCTCAGGCCCGGCCCTGAGTGACAGCTTCGACCCGGCAGACCCCGGCGTGCCGGTGGCg                                                                                                   >  6:53968/1‑149 (MQ=255)
                                                     tGTTCGATGCCGGCGGCAATGTCTGGGAGCTGACTCGCACGGGCTATCGGCCAGGCCATGACGCACAGCGCGACGCCAAGCTCGACCCCTCAGGCCCGGCCCTGAGTGACAGCTTCGACCCGGCAGACCCCGGCGTGCCGGTGGCGGt                                                                                                 >  5:236299/1‑148 (MQ=255)
                                                       ttCGATGCCGGCGGCAATGTCTGGGAGCTGACTCGCACGGGCTATCGGCCAGGCCATGACGCACAGCGCGACGCCAAGCTCGACCCCTCAGGCCCGGCCCTGAGTGACAGCTTCGACCCGGCAGACCCCGGCGTGCCGGTGGCGGTaa                                                                                               <  6:313288/148‑1 (MQ=255)
                                                         cGATGCCGGCGGCAATGTCTGGGAGCTGACTCGCACGGGCTATCGGCCAGGCCATGACGCACAGCGCGACGCCAAGCTCGACCCCTCAGGCCCGGCCCTGAGTGACAGCTTCGACCCGGCAGACCCCGGCGTGCCGGTGGCGGTAATc                                                                                             <  4:368907/148‑1 (MQ=255)
                                                            tGCCGGCGGCAATGTCTGGGAGCTGACTCGCACGGGCTATCGGCCAGGCCATGACGCACAGCGCGACGCCAAGCTCGACCCCTCAGGCCCGGCCCTGAGTGACAGCTTCGACCCGGCAGACCCCGGCGTGCCGGTgg                                                                                                     <  3:331648/137‑1 (MQ=255)
                                                            tGCCGGCGGCAATGTCTGGGAGCTGACTCGCACGGGCTATCGGCCAGGCCATGACGCACAGCGCGACGCCAAGCTCGACCCCTCAGGCCCGGCCCTGAGTGACAGCTTCGACCCGGCAGACCCCGGCGTGCCGGTgg                                                                                                     >  4:331648/1‑137 (MQ=255)
                                                            tGCCGGCGGCAATGTCTGGGAGCTGACTCGCACGGGCTATCGGCCAGGCCATGACGCACAGCGCGACGCCAAGCTCGACCCCTCAGGCCCGGCCCTGAGTGACAGCTTCGACCCGGCAGACCCCGGCGTGCCGGTGGCGGTAATCAAAg                                                                                         >  4:267667/1‑149 (MQ=255)
                                                              ccGGCGGCAATGTCTGGGAGCTGACTCGCACGGGCTATCGGCCAGGCCATGACGCACAGCGCGACGCCAAGCTCGACCCCTCAGGCCCGGCCCTGAGTGACAGCTTCGACCCGGCAGACCCCGGCGTGCCGGTGGCGGTAATCAAAggc                                                                                       <  5:301860/149‑1 (MQ=255)
                                                                ggcggcAATGTCTGGGAGCTGACTCGCACGGGCTATCGGCCAGGCCATGACGCACAGCGCGACGCCAAGCTCGACCCCTCAGGCCCGGCCCTGAGTGACAGCTTCGACCCGGCAGACCCCGGCGTGCCGGTgg                                                                                                     >  5:19367/1‑133 (MQ=255)
                                                                ggcggcAATGTCTGGGAGCTGACTCGCACGGGCTATCGGCCAGGCCATGACGCACAGCGCGACGCCAAGCTCGACCCCTCAGGCCCGGCCCTGAGTGACAGCTTCGACCCGGCAGACCCCGGCGTGCCGGTgg                                                                                                     <  6:19367/133‑1 (MQ=255)
                                                                ggcggcAATGTCTGGGAGATGACTCGCACGGGCTATCGGCCAGGCCATGACGCACAGCGCGACGCCAAGCTCGACCCCTCAGGCCCGGCCCTGAGTGACAGCTTCGACCCGGCAGACCCCGGCGTGCCGGTGGCGGTAATCAAAggcgg                                                                                     <  2:68135/149‑1 (MQ=255)
                                                                 gcggcAATGTCTGGGAGCTGACTCGCACGGGCTATCGGCCAGGCCATGACGCACAGCGCGACGCCAAGCTCGACCCCTCAGGCCCGGCCCTGAGTGACAGCTTCGACCCGGCAGACCCCGGCGTGCCGGTGGCGGTAATCAAAggcggc                                                                                    <  4:242788/149‑1 (MQ=255)
                                                                   ggcAATGTCTGGGAGCTGACTCGCACGGGCTATCGGCCAGGCCATGACGCACAGCGCGACGCCAAGCTCGACCCCTCAGGCCCGGCCCTGAGTGACAGCTTCGACCCGGCAGACCCCGGCGTGCCGGTGGCGGTAATCAAAGGCGGCTc                                                                                  <  1:38375/149‑1 (MQ=255)
                                                                   ggcAATGTCTGGGAGCTGACTCGCACGGGCTATCGGCCAGGCCATGACGCACAGCGCGACGCCAAGCTCGACCCCTCAGGCCCGGCCCTGAGTGACAGCTTCGACCCGGCAGACCCCGGCGTGCCGGTGGCGGTAATCAAAGGCGGCTc                                                                                  <  1:148623/149‑1 (MQ=255)
                                                                   ggcAATGTCTGGGAGCTGACTCGCACGGGCTATCGGCCAGGCCATGACGCACAGCGCGACGCCAAGCTCGACCCCTCAGGCCCGGCCCTGAGTGACAGCTTCGACCCGGCAGACCCCGGCGTGCCGGTGGCGGTAATCAAAGGCGGCTc                                                                                  <  6:275744/149‑1 (MQ=255)
                                                                    gcAATGTCTGGGAGCTGACTCGCACGGGCTATCGGCCAGGCCATGACGCACAGCGCGACGCCAAGCTCGACCCCTCAGGCCCGGCCCTGAGTGACAGCTTCGACCCGGCAGACCCCGGCGTGCCGGTGGCGGTAATCAAAGGCGGCTCg                                                                                 >  6:46366/1‑149 (MQ=255)
                                                                      aaTGTCTGGGAGCTGACTCGCACGGGCTATCGGCCAGGCCATGACGCACAGCGCGACGCCAAGCTCGACCCCTCAGGCCCGGCCCTGAGTGACAGCTTCGACCCGGCAGACCCCGGCGTGCCGGTGGCGGTAATCAAAGGCGGCTCGCa                                                                               >  3:52176/1‑149 (MQ=255)
                                                                           cTGGGAGCTGACTCGCACGGGCTATCGGCCAGGCCATGACGCACAGCGCGACGCCAAGCTCGACCCCTCAGGCCCGGCCCTGAGTGACAGCTTCGACCCGGCAGACCCCGGCGTGCCGGTGGCGGTAATCAAAGGCGGCTCGCACCtgt                                                                          <  5:215321/149‑1 (MQ=255)
                                                                              ggAGCTGACTCGCACGGGCTATCGGCCAGGCCATGACGCACAGCGCGACGCCAAGCTCGACCCCTCAGGCCCGGCCCTGAGTGACAGCTTCGACCCGGCAGACCCCGGCGTGCCGGTGGCGGTAATCAAAGGCGGCTCGCACCTgtgt                                                                        <  6:236299/148‑1 (MQ=255)
                                                                              ggAGCTGACTCGCACGGGCTATCGGCCAGGCCATGACGCACAGCGCGACGCCAAGCTCGACCCCTCAGGCCCGGCCCTGAGTGACAGCTTCGACCCGGCAGACCCCGGCGTGCCGGTGGCGGTAATCAAAGGCGGCTCGCACCTGTGtt                                                                       <  8:181436/149‑1 (MQ=255)
                                                                               gAGCTGACTCGCACGGGCTATCGGCCAGGCCATGACGCACAGCGCGACGCCAAGCTCGACCCCTCAGGCCCGGCCCTGAGTGACAGCTTCGACCCGGCAGACCCCGGCGTGCCGGTGGCGGTAATCAAAGGCGGCTCGCACCTGTGTTc                                                                      <  3:267667/149‑1 (MQ=255)
                                                                                  cTGACTCGCTCGGGCTATCGGCCAGGCCATGACGCACAGCGCGACGCCACGCTCGACCCCTCAGGCCCGGCCCTGAGTGACAGCTTCGACCCGGCAGACCCCGGCGTGCCGGTGGCGGTAATCAAAGGCGGCTCGCACCTGTGTTcggc                                                                   <  2:360108/149‑1 (MQ=255)
                                                                                     aCTCGCACGGGCTATCGGCCAGGCCATGACGCACAGCGCGACGCCAAGCTCGACCCCTCAGGCCCGGCCCTGAGTGACAGCTTCGACCCGGCAGACCCCGGCGTGCCGGTGGCGGTAATCAAAGGCGGCTCGCACCTGTGTTCGGCGGa                                                                >  7:339175/1‑149 (MQ=255)
                                                                                       tCGCACGGGCTATCGGCCAGGCCATGACGCACAGCGCGACGCCAAGCTCGACCCCTCAGGCCCGGCCCTGAGTGACAGCTTCGACCCGGCAGACCCCGGCGTGCCGGTGGCGGTAATCAAAg                                                                                         <  5:292889/122‑1 (MQ=255)
                                                                                       tCGCACGGGCTATCGGCCAGGCCATGACGCACAGCGCGACGCCAAGCTCGACCCCTCAGGCCCGGCCCTGAGTGACAGCTTCGACCCGGCAGACCCCGGCGTGCCGGTGGCGGTAATCAAAg                                                                                         >  6:292889/1‑122 (MQ=255)
                                                                                       tCGCACGGGCTATCGGCCAGGCCATGACGCACAGCGCGACGCCAAGCTCGACCCCTCAGGCCCGGCCCTGAGTGACAGCTTCGACCCGGCAGACCCCGGCGTGCCGGTGGCGGTAATCAAAGGCGGCTCGCACCTGTGTTCGGCGGAcc                                                              >  5:316308/1‑149 (MQ=255)
                                                                                            cGGGCTATCGGCCAGGCCATGACGCACAGCGCGACGCCAAGCTCGACCCCTCAGGCCCGGCCCTGAGTGACAGCTTCGACCCGGCAGACCCCGGCGTGCCGGTGGCGGTAATCAAAGGCGGCTCGCACCTGTGTTCGGCGGACCGCTGc                                                         <  3:33379/149‑1 (MQ=255)
                                                                                                 tATCGGCCAGGCCATGACGCACAGCGCGACGCCAAGCTCGACCCCTCAGGCCCGGCCCTGAGTGACAGCTTCGACCCCGCAGACCCCGGCGTGCCGGTGGCGGTAATCAAAGGCGGCACGCACCTGTGTTCGGCCGACCGCTGCATgc                                                     >  3:236562/1‑148 (MQ=255)
                                                                                                   tgggccaggcgaTGACGCACAGCGCGACGCGAAGCTCGACCCCTCAGGGCCGGCCCTGAGTGAGAGCTTCGACCCGGCAGACCGCGGCGTGCCGGTGGCGGTAATCAAAGGCGGCTCGCACCTGTGTTCGGCGGACCGCTGCATgc                                                     <  4:236562/144‑1 (MQ=255)
                                                                                                    cGGCCAGGCCATGACGCACAGCGCGACGCCAAGCTCGACCCCTCAGGCCCGGCCCTGAGTGACAGCTTCGACCCGGCAGACCCCGGCGTGCCGGTGGCGGTAATCAAAGGCGGCTCGCACCTGTGTTCGGCGGACCGCTGCAt                                                       >  7:148209/1‑143 (MQ=255)
                                                                                                    cGGCCAGGCCATGACGCACAGCGCGACGCCAAGCTCGACCCCTCAGGCCCGGCCCTGAGTGACAGCTTCGACCCGGCAGACCCCGGCGTGCCGGTGGCGGTAATCAAAGGCGGCTCGCACCTGTGTTCGGCGGACCGCTGCAt                                                       <  8:148209/143‑1 (MQ=255)
                                                                                                       ccaggccaTGACGCACAGCGCGACGCCAAGCTCGACCCCTCAGGCCCGGCCCTGAGTGACAGCTTCGACCCGGCAGACCCCGGCGTGCCGGTGGCGGTAATCAAAGGCGGCTCGCACCTGTGTTCGGCGGACCGCTGCATGCGCTACCg                                              <  5:46366/149‑1 (MQ=255)
                                                                                                               tGACGCACAGCGCGACGCCAAGCTCGACCCCTCAGGCCCGGCCCTGAGTGACAGCTTCGACCCGGCAGACCCCGGCGTGCCGGTGGCGGTAATCAAAGGCGGCTCGCACCTGTGTTCGGCGGACCGCTGCATGCGCTACCGCCCCTCgg                                      <  1:7280/149‑1 (MQ=255)
                                                                                                               tGACGCACAGCGCGACGCCAAGCTCGACCCCTCAGGCCCGGCCCTGAGTGACAGCTTCGACCCGGCAGACCCCGGCGTGCCGGTGGCGGTAATCAAAGGCGGCTCGCACCTGTGTTCGGCGGACCGCTGCATGCGCTACCGCCCCTCgg                                      >  8:185833/1‑149 (MQ=255)
                                                                                                                 aCGCACAGCGCGACGCCAAGCTCGACCCCTCAGGCCCGGCCCTGAGTGACAGCTTCGACCCGGCAGACCCCGGCGTGCCGGTGGCGGTAATCAAAGGCGGCTCGCACCTGTGTTCGGCGGACCGCTGCATGCGCTACCGCCCCTCGGc                                     <  7:185833/148‑1 (MQ=255)
                                                                                                                         cgcgACGCCAAGCTCGACCCCTCAGGCCCGGCCCTGAGTGACAGCTTCGACCCGGCAGACCCCGGCGTGCCGGTGGCGGTAATCAAAGGCGGCTCGCACCTGTGTTCGGCGGACCGCTGCATGCGCTACCGCCCCTCGGCACGCcagcc                            >  8:152517/1‑149 (MQ=255)
                                                                                                                                ccAAGCTCGACCCCTCAGGCCCGGCCCTGAGTGACAGCTTCGACCCGGCAGACCCCGGCGTGCCGGTGGCGGTAATCAAAGGCGGCTCGCACCTGTGTTCGGCGGACCGCTGCATGCGCTACCGCCCCTCGGCACGCcagccgcagccg                     <  6:283982/149‑1 (MQ=255)
                                                                                                                                ccAAGCTCGACCCCTCAGGCCCGGCCCTGAGTGACAGCTTCGACCCGGCAGACCCCGGCGTGCCGGTGGCGGTAATCAAAGGCGGCTCGCACCTGTGTTCGGCGGACCGCTGCATGCGCTACCGCCCCTCGGCACGCcagccgcagccg                     <  8:339175/149‑1 (MQ=255)
                                                                                                                                ccAAGCTCGACCCCTCAGGCCCGGCCCTGAGTGACAGCTTCGACCCGGCAGACCCCGGCGTGCCGGTGGCGGTAATCAAAGGCGGCTCGCACCTGTGTTCGGCGGACCGCTGCATGCGCTACCGCCCCTCGGCACGCcagccgcagccg                     <  4:234487/149‑1 (MQ=255)
                                                                                                                                 cAAGCTCGACCCCTCAGGCCCGGCCCTGAGTGACAGCTTCGACCCGGCAGACCCCGGCGTGCCGGTGGCGGTAATCAAAGGCGGCTCGCACCTGTGTTCGGCGGACCGCTGCATGCGCTACCGCCCCTCGGCACGCCAGCCGCAGCCgg                    >  6:142573/1‑149 (MQ=255)
                                                                                                                                 cAAGCTCGACCCCTCAGGCCCGGCCCTGAGTGACAGCTTCGACCCGGCAGACCCCGGCGTGCCGGTGGCGGTAATCAAAGGCGGCTCGCACCTGTGTTCGGCGGACCGCTGCATGCGCTACCGCCCCTCGGCACGCCAGCCGCAGCCgg                    >  3:19794/1‑149 (MQ=255)
                                                                                                                                  aaGCTCAACCCCTCAGGCCCGGCCCTGAGTGACAGCTTCGACCCGGCAGACCCCGGCGTGCCGGTGGCGGTAATCAAAGGCGGCTCGCACCTGTGTTCGGCGGACCGCTGCATGCGCTACCGCCCCTCGGCACGCCAGCCGCAGCCGgt                   >  6:203680/1‑149 (MQ=255)
                                                                                                                                   aGCTCGACCCCTCAGGCCCGGCCCTGAGTGACAGCTTCGACCCGGCAGACCCCGGCGTGCCGGTGGCGGTAATCAAAGGCGGCTCGCACCTGTGTTCGGCGGACCGCTGCATGCGCTACCGCCCCTCGGCACGCCAGCCGCAGCCGgtg                  <  6:35044/149‑1 (MQ=255)
                                                                                                                                        tgtcacTCAGGGCCGGGCCTGAGTGACAGCTTCGACCCGGCAGACCCCGGCGTGCCGGTGGCGGTAATCAAAGGCGGCTCGCACCTGTGTTCGGCGGACCGCTGCATGCGCTACCGCCCCTCGGCACGCCAGCCGCAGCCGGTGtt                >  6:66336/6‑146 (MQ=255)
                                                                                                                                        tgtcacTCAGGGCCGGGCCTGAGTGACAGCTTCGACCCGGCAGACCCCGGCGTGCCGGTGGCGGTAATCAAAGGCGGCTCGCACCTGTGTTCGGCGGACCGCTGCATGCGCTACCGCCCCTCGGCACGCCAGCCGCAGCCGGTGtt                <  5:66336/141‑1 (MQ=255)
                                                                                                                                        gACCCCTCAGGCCCGGCCCTGAGTGACAGCTTCGACCCGGCAGACCCCGGCGTGCCGGTGGCGGTAATCAAAGGCGGCTCGCACCTGTGTTCGGCGGACCGCTGCATGCGCTACCGCCCCTCGGCACGCCAGCCGCAGCCGGTGTTc               <  4:19794/147‑1 (MQ=255)
                                                                                                                                          ccccTCAGGCCCGGCCCTGAGTGACAGCTTCGACCCGGCAGACCCCGGCGTGCCGGTGGCGGTAATCAAAGGCGGCTCGCACCTGTGTTCGGCGGACCGCTGCATGCGCTACCGCCCCTCGGCACGCCAGCCTCAGCCGGTGTTCATga           <  5:50021/149‑1 (MQ=255)
                                                                                                                                             cTCAGGCCCGGCCCTGAGTGACAGCTTCGACCCGGCAGACCCCGGCGTGCCGGTGGCGGTAATCAAAGGCGGCTCGCACCTGTGTTCGGCGGACCGCTGCATGCGCTACCGCCCCTCGGCACGCCAGCCGCAGCCGGTGTTCATgacga        >  4:53761/1‑149 (MQ=255)
                                                                                                                                                   cccggcccTGAGTGACAGCTTCGACCCGGCAGACCCCGGCGTGCCGGTGGCGGTAATCAAAGGCGGCTCGCACCTGTGTTCGGCGGACCGCTGCATGCGCTACCGCCCCTCGGCACGCCAGCCGCAGCCGGTGTTCATGACGACCTCGc  <  7:152517/149‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                    
TGGTTTTACCGGTACATCGCCCGTGGGTTGCTTCCCGGCCAACGGCTTTGGCCTGTTCGATGCCGGCGGCAATGTCTGGGAGCTGACTCGCACGGGCTATCGGCCAGGCCATGACGCACAGCGCGACGCCAAGCTCGACCCCTCAGGCCCGGCCCTGAGTGACAGCTTCGACCCGGCAGACCCCGGCGTGCCGGTGGCGGTAATCAAAGGCGGCTCGCACCTGTGTTCGGCGGACCGCTGCATGCGCTACCGCCCCTCGGCACGCCAGCCGCAGCCGGTGTTCATGACGACCTCGC  >  NC_002947/3791258‑3791553

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: