Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,212,472 C→T 5.4% intergenic (+36/+58) PP_4593 → / ← PP_4594 hypothetical protein/cystathionine gamma‑synthase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,212,4720CT 5.4% 285.4 / 3.7 92intergenic (+36/+58)PP_4593/PP_4594hypothetical protein/cystathionine gamma‑synthase
Reads supporting (aligned to +/- strand):  ref base C (43/44);  new base T (3/2);  total (46/46)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.35e-01

TGGTAGCCGCCGAAGGCAACGGCACCAACGACGAAGCCGTCAACCGCCTGCGCAAGGAAAGGTTGCTGATCAAGGACAAGATCGAACGGATCATTCATCCGCCCAAATCC‑TGACAGGCATCTATTGAATCCACCGGCCTCTTCGCGGGCAAGCCCGCTCCTACACCGTCGTAGGAGCGGGTTTACCCGCGAAAGGGCCTGCGGCCCTCACACAGCGACCAGCGCTTTCAACGCCCCATCCAGGTCCTCGACCAGGTCACGGTAATCCTCGATCCCTACCGACAACCGCAGCAGG  >  NC_002947/5212324‑5212617
                                                                                                                                                     |                                                                                                                                                 
tggtAGCCGCCGAAGGCAACGGCACCAACGACGAAGCCGTCAACCGCCTGCGCGAGGAAAGGTTGCTGATCAAGGACAAGATCGAACGGATCATTCATCCGCCCAAATCC‑TGACAGGCATCTATTGAATCCACCGGCCTCTTCGCGGGc                                                                                                                                                   <  7:347193/149‑1 (MQ=255)
 ggtAGCCGCCGAAGGCAACGGCACCAACGACGAAGCCGTCAACCGCCTGCGCAAGGAAAGGTTGCTGATCAAGGACAAGATCGAACGGATCATTCATCCGCCCAAATCC‑TGACAGGCATCTATTGAATCCACCGGCCTCTTCGCGGGc                                                                                                                                                   >  2:108921/1‑148 (MQ=255)
  gtAGCCGCCGAAGGCAACGGCACCAACGACGAAGCCGTCAACCGCCTGCGCAAGGAAAGGTTGCTGATCAAGGACAAGATCGAACGGATCATTCATCCGCCCAAATCC‑TGACAGGCATCTATTGAATCCACCGGCCTCTTCGCGGGc                                                                                                                                                   >  6:10367/1‑147 (MQ=255)
     gccgccGAAGGCAACGGCACCAACGACGAAGCCGTCAACCGCCTGCGCAAGGAAAGGTTGCTGATCAAGGACAAGATCGAACGGATCATTCATCCGCCCAAATCC‑TGACAGGCATCTATTGAATCCACCGGCCTCTTCGCGGGCAAGc                                                                                                                                               <  1:167058/148‑1 (MQ=255)
        gccgAAGGCAACGGCACCAACGACGAAGCCGTCAACCGCCTGCGCAAGGAAAGGTTGCTGATCAAGGACAAGATCGAACGGATCATTCATCCGCCCAAATCC‑TGACAGGCATCTATTGAATCCACCGGCCTCTTCGCGGGCAAGCCCGc                                                                                                                                           >  5:91920/1‑149 (MQ=255)
        gccgAAGGCAACGGCACCAACGACGAAGCCGTCAACCGCCTGCGCAAGGAAAGGTTGCTGATCAAGGACAAGATCGAACGGATCATTCATCCGCCCAAATCC‑TGACAGGCATCTATTGAATCCACCGGCCTCTTCGCGGGCAAGCCCGc                                                                                                                                           <  6:91920/149‑1 (MQ=255)
         ccgAAGGCAACGGCACCAACGACGAAGCCGTCAACCGCCTGCGCAAGGAAAGGTTGCTGATCAAGGACAAGATCGAACGGATCATTCATCCGCCCAAATCC‑TGACAGGCATCTATTGAATCCACCGGCCTATTCGCGGGCAAGCCCGc                                                                                                                                           >  7:238885/1‑148 (MQ=255)
               gCAACGGCACCAACGACGAAGCCGTCAACCGCCTGCGCAAGGAAAGGTTGCTGATCAAGGACAAGATCGAACGGATCATTCATCCGCCCAAATCC‑TGACAGGCATCTATTGAATCCACCGGCCTCTTCGCGGGCAAGCCCGCTCCTaca                                                                                                                                    >  2:86665/1‑149 (MQ=255)
                cAACGGCACCAACGACGAAGCCGTCAACCGCCTGCGCAAGGAAAGGTTGCTGATCAAGGAC‑AGATCGAACGGATCATTCATCCGCCCAAATCC‑TGACAGGCATCTATTGAATCCACCGGCCTCTTCGCGGGCAAGCCCGCTCCTACAcc                                                                                                                                  >  2:103570/1‑149 (MQ=255)
                 aaCGGCACCAACGACGAAGCCGTCAACCGCCTGCGCAAGGAAAGGTTGCTGATCAAGGACAAGATCGAACGGATCATTCATCCGCCCAAATCC‑TGACAGGCATCTATTGAATCCACCGGCCTCTTCGCGGGCAAGCCCGCTCCTACAcc                                                                                                                                  >  2:283333/1‑149 (MQ=255)
                    ggCACCAACGACGAAGCCGTCAACCGCCTGCGCAAGGAAAGGTTGCTGATCAAGGACAAGATCGAACGGATCATTCATCCGCCCAAATCC‑TGACAGGCATCTATTGAATCCACCGGCCTCTTCGCGGGCAAGCCCGCTCCTACACcgtc                                                                                                                               <  7:19125/149‑1 (MQ=255)
                      cACCAACGACGAAGCCGTCAACCGCCTGCGCAAGGAAAGGTTGCTGATCAAGGACAAGATCGAACGGATCATTCATCCGCCCAAATCC‑TGACAGGCATCTATTGAATCCACCGGCCTCTTCGCGGGCAAGCCCGCTCCTACACcgtcgt                                                                                                                             <  1:351008/149‑1 (MQ=255)
                        ccAACGACGAAGCCGTCAACCGCCTGCGCAAGGAAAGGTTGCTGATCATGGACAAGATCTAACGGATCATTCATCCGCCCAAATCCTTGACAGGCATCTATTGAATCCACCGGCCTCTTCGCGGGCAAGCCCGCTCCTACACCGTCGTa                                                                                                                            <  1:108921/149‑1 (MQ=255)
                        ccAACGACGAAGCCGTCAACCGCCTGCGCAAGGAAAGGTTGCTGATCAAGGACAAGATCGAACGGATCATTCATCCGCCCAAATCC‑TGACAGGCATCTATTGAATCCACCGGCCTCTTCGCGGGCAAGCCCGCTCCTACACCGTCGTAg                                                                                                                           <  1:86665/149‑1 (MQ=255)
                        ccAACGACGAAGCCGTCAACCGCCTGCGCAAGGAAAGGTTGCTGATCAAGGACAAGATCGAACGGATCATTCATCCGCCCAAATCC‑TGACAGGCATCTATTGAATCCACCGGCCTCTTCGCGGGCAAGCCCGCTCCTACACCGTCGTAg                                                                                                                           <  7:157335/149‑1 (MQ=255)
                        ccAACGACGAAGCCGTCAACCGCCTGCGCAAGGAAAGGTTGCTGATCAAGGACAAGATCGAACGGATCATTCATCCGCCCAAATCC‑TGACAGGCATCTATTGAATCCACCGGACTCTTCGCGGGCAAGCCCGCTCCTACACCGTCGTAg                                                                                                                           >  5:360740/1‑149 (MQ=255)
                          aaCGACGAAGCCGTCAACCGCCTGCGCAAGGAAAGGTTGCTGATCAAGGACAAGATCGAACGGATCATTCATCCGCCCAAATCC‑TGACAGGCATCTATTGAATCCACCGGCCTCTTCGCGGGCAAGCCCGCTCCTACACCGTCGTAg                                                                                                                           >  4:28048/1‑147 (MQ=255)
                            cgacgaAGCCGTCAACCGCCTGCGCAAGGAAAGGTTGCTGATCAAGGAC‑AGATCGAACGGATCATTCATCCGCCCAAATCC‑TGACAGGCATCTATTGAATCCACCGGCCTCTTCGCGGGCAAGCCCGCTCCTACACCGTCGTAGGAGCg                                                                                                                      <  1:103570/149‑1 (MQ=255)
                              acgaAGCCGTCAACCGCCTGCGCAAGGAAAGGTTGCTGATCAAGGACAAGATCGAATGGATCATTCATCCGCCCAAATCC‑TGACAGGCATCTATTGAATCCACCGGCCTCTTCGCGGGCAAGCCCGCTCCTACACCGTCGTAGGAGCgg                                                                                                                     <  5:300553/149‑1 (MQ=255)
                              acgaAGCCGTCAACCGCCTGCGCAAGGAAAGGTTGCTGATCAAGGACAAGATCGAACGGATCATTCATCCGCCCAAATCC‑TGACAGGCATCTATTGAATCCACCGGCCTCTTCGCGGGCAAGCCCGCTCCTACACCGTCGTAGGAGCgg                                                                                                                     <  3:135831/149‑1 (MQ=255)
                               cgaAGCCGTCAACCGCCTGCGCAAGGAAAGGTTGCTGATCAAGGACAAGATCGAACGGATCATTCATCCGCCCAAATCC‑TGACAGGCATCTATTGAATCCACCGGCCTCTTCGCGGGCAAGCCCGCTCCTACACCGTCGTAGGAGCggg                                                                                                                    <  3:194897/149‑1 (MQ=255)
                                gaAGCCGTCAACCGCCTGCGCAAGGAAAGGTTGCTGATCAAGGACAAGATCGAACGGATCATTCATCCGCCCAAATCC‑TGACAGGCATCTATTGAATCCACCGGCCTCTTCGCGGGCAAGCCCGCTCCTACACCGTCGTAGGAGCGGGt                                                                                                                   >  3:327140/1‑149 (MQ=255)
                                gaAGCCGTCAACCGCCTGCGCAAGGAAAGGTTGCTGATCAAGGACAAGATCGAACGGATCATTCATCCGCCCAAATCC‑TGACAGGCATCTATTGAATCCACCGGCCTCTTCGCGGGCAAGCCCGCTCCTACACCGTCGTAGGAGCGGGt                                                                                                                   <  3:189720/149‑1 (MQ=255)
                                 aaGCCGTCAACCGCCTGCGCAAGGAAAGGTTGCTGATCAAGGACAAGATCGAACGGATCATTCATCCGCCCAAATCC‑TGACAGGCATCTATTGAATCCACCGGCCTCTTCGCGGGCAAGCCCGCTCCTACACCGTCGTAGGAGCGGGtt                                                                                                                  <  4:47758/149‑1 (MQ=255)
                                  aGCCGTCAACCGCCTGCGCAAGGAAAGGTTGCTGATCAAGGACAAGATCGAACGGATCATTCATCCGCCCAAATCC‑TGACAGGCATCTATTGAATCCACCGGCCTCTTCGCGGGCAAGCCCGCTCCTACACCGTCGTAGGAGCGGGtt                                                                                                                  >  6:225367/1‑148 (MQ=255)
                                      gTCAACCGCCTGCGCAAGGAAAGGTTGCTGATCAAGGACAAGATCGAACGGATCATTCATCCGCCCAAATCC‑TGACAGGCATCTATTGAATCCACCGGCCTCTTCGCGGGCAAGCCCGCTCCTACACCGTCGTAGGAGCGGGTTTAccc                                                                                                             <  8:336234/149‑1 (MQ=255)
                                      gTCAACCGCCTGCGCAAGGAAAGGTTGCTGATCAAGGACAAGATCGAACGGATCATTCATCCGCCCAAATCC‑TGACAGGCATCTATTGAATCCACCGGCCTCTTCGCGGGCAAGCCCGCTCCTACACCGTCGTAGGAGCGGGTTTAccc                                                                                                             <  4:327140/149‑1 (MQ=255)
                                       tCAACCGCCTGCGCAAGGAAAGGTTGCTGATCAAGGACAAGATCGAACGGATCATTCATCCGCCCAAATCC‑TGACAGGCATCTATTGAATCCACCGGCCTCTTCGCGGGCAAGCCCGCTCCTACACCGTCGTAGGAGCGGGTTTACCcg                                                                                                            <  2:174354/149‑1 (MQ=255)
                                       tCAACCGCCTGCGCAAGGAAAGGTTGCTGATCAAGGACAAGATCGAACGGATCATTCATCCGCCCAAATCC‑TGACAGGCATCTATTGAATCCACCGGCCTCTTCGCGGGCAAGCCCGCTCCTACACCGTCGTAGGAGCGGGTTTACAcg                                                                                                            >  5:45830/1‑149 (MQ=255)
                                                     aaGGAAAGGTTGCTGATCAAGGACAAGATCGAACGGATCATTCATCCGCCCAAATCC‑TGACAGGCATCTATTGAATCCACCGGCCTCTTCGCGGGCAAGCCCGCTCCTACACCGTCGTAGGAGCGGGTGTACCCGCGAAAGGGCCTGCg                                                                                              >  8:349644/1‑149 (MQ=255)
                                                        gAAAGGTTGCTGATCAAGGACAAGATCGAACGGATCATTCATCCGCCCAAATCC‑TGACAGGCATCTATTGAATCCACCGGCCTCTTCGCGGGCAAGCCCGCTCCTACACCGTCGTAGGAGCGGGTTTACCCGCGAAAGGGCCTGCGGcc                                                                                           <  3:205369/149‑1 (MQ=255)
                                                        gAAAGGTTGCTGATCAAGGACAAGATCGAACGGATCATTCATCCGCCCAAATCC‑TGACAGGCATCTATTGAATCCACCGGCCTCTTCGCGGGCAAGCCCGCTCCTACACCGTCGTAGGAGCGGGTTTACCCGCGAAAGGGCCTGCGGcc                                                                                           >  4:205369/1‑149 (MQ=255)
                                                        gAAAGGTTGCTGATCAAGGACAAGATCGAACGGATCATTCATCCGCCCAAATCC‑TGACAGGCATCTATTGAATCCACCGGCCTCTTCGCGGGCAAGCCCGCTCCTACACCGTCGTAGGAGCGGGTTTACCCGCGAAAGGGCCTGCGGcc                                                                                           <  7:349644/149‑1 (MQ=255)
                                                              ttGCTGATCAAGGACAAGATCGAACCGATCATTCATCCGCCCAAATCC‑TGACAGGCATCTATTGAATCCACCGGCCTCTTCGCGGGCAAGCCCGCTCCTACACCGTCGTAGGAGCGGGTTTACCCGCGAAAGCGCCTGCGGCCCTcaca                                                                                     >  2:308504/1‑149 (MQ=255)
                                                                gCTGATCAAGGACAAGATCGAACGGATCATTCATCCGCCCAAATCC‑TGACAGGCATCTATTGAATCCACCGGCCTCTTCGCGGGCAAGCCCGCTCCTACACCGTCGTAGGAGCGGGTTTACCCGCGAAAGGGCCTGCGGCCCTcacaca                                                                                   >  7:199828/1‑149 (MQ=255)
                                                                gCTGATCAAGGACAAGATCGAACGGATCATTCATCCGCCCAAATCC‑TGACAGGCATCTATTGAATCCACCGGCCTCTTCGCGGGCAAGCCCGCTCCTACACCGTCGTAGGAGCGGGTTTACCCGCGAAAGGGCCTGCGGCCCTcaca                                                                                     >  8:101805/1‑147 (MQ=255)
                                                                  tGATCAAGGACAAGATCGAACGGATCATTCATCCGCCCAAATCC‑TGACAGGCATCTATTGAATCCACCGGCCTCTTCGCGGGCAAGCCCGCTCCTACACCGTCGTAGGAGCGGGTTTACCCGCGAAAGGGCCTGc                                                                                               <  4:96440/135‑1 (MQ=255)
                                                                  tGATCAAGGACAAGATCGAACGGATCATTCATCCGCCCAAATCC‑TGACAGGCATCTATTGAATCCACCGGCCTCTTCGCGGGCAAGCCCGCTCCTACACCGTCGTAGGAGCGGGTTTACCCGCGAAAGGGCCTGc                                                                                               >  3:96440/1‑135 (MQ=255)
                                                                  tGATCAAGGACAAGATCGAACGGATCATTCATCCGCCCAAATCC‑TGACAGGCATCTATTGAATCCACCGGCCTCTTCGCGGGCAAGCCCGCTCCTACACCGTCGTAGGAGCGGGTTTACCCGCGAAAGGGCCGCAGGCCCtttcgc                                                                                    <  6:268857/146‑6 (MQ=255)
                                                                  tGATCAAGGACAAGATCGAACGGATCATTCATCCGCCCAAATCC‑TGACAGGCATCTATTGAATCCACCGGCCTCTTCGCGGGCAAGCCCGCTCCTACACCGTCGTAGGAGCGGGTTTACCCGCGAAAGGGCCGCAGGCCCtttcgc                                                                                    >  5:268857/1‑141 (MQ=255)
                                                                    aTCAAGGACAAGATCGAACGGATCATTCATCCGCCCAAATCC‑TGACAGGCATCTATTGAATCCACCGGCCTCTTCGCGGGTAAACCCGCTCCTACGACGGTGTAGGAGCGGGCTTGCCCGCGAAGAGGccggtgg                                                                                             <  2:88299/135‑6 (MQ=18)
                                                                    aTCAAGGACAAGATCGAACGGATCATTCATCCGCCCAAATCC‑TGACAGGCATCTATTGAATCCACCGGCCTCTTCGCGGGTAAACCCGCTCCTACGACGGTGTAGGAGCGGGCTTGCCCGCGAAGAGGccggtgg                                                                                             >  1:88299/1‑130 (MQ=18)
                                                                      cAAGGACAAGATCGAACGGATCATTCATCCGCCCAAATCC‑TGACAGGCATCTATTGAATCCACCGGCCTCTTCGCGGGCAAGCCCGCTCCTACACCGTCGTAGGAGCGGGTTTACCCGCGAAAGGGCCTGCGGCCCTCACACAGCGAc                                                                              >  3:256479/1‑148 (MQ=255)
                                                                      cAAGGACAAGATCGAACGGATCATTCATCCGCCCAAATCC‑TGACAGGCATCTATTGAATCCACCGGCCTCTTCGCGGGCAAGCCCGCTCCTACACCGTCGTAGGAGCGGGTTTACCCGCGAAAGGGCCTGCGGCCCTCACACAGCGAc                                                                              >  5:233313/1‑148 (MQ=255)
                                                                       aaGGACAAGATCGAACGGATCATTCATCCGCCCAAATCC‑TGACAGGCATCTATTGAATCCACCGGCCTCTTCGCGGGCAAGCCCGCTCCTACACCGTCGTAGGAGCGGGTTTACCCGCGAAAGGGCCTGCGGCCCTCACACAGCGAc                                                                              >  4:88695/1‑147 (MQ=255)
                                                                         ggACAAGATCGAACGGATCATTCATCCGCCCAAATCC‑TGACAGGCATCTATTGAATCCACCGGCCTCTTCGCGGGCAAGCCCGCTCCTACACCGTCGTAGGAGCGGGTTTACCCGCGAAAGGGCCTGCGGCCCTCACACAGCGACCAgc                                                                          <  5:237260/149‑1 (MQ=255)
                                                                          gACAAGATCGAACGGATCATTCATCCGCCCAAATCC‑TGACAGGCATCTATTGAATCCACCGGCCTCTTCGCGGGCAAGCCCGCTCCTACACCGTCGTAGGAGCGGGTTTACCCGCGAAAGGGCCTGCGGCCCTCACACAGCGACCAgcg                                                                         >  8:346264/1‑149 (MQ=255)
                                                                            cAAGATCGAACGGATCATTCATCCGCCCAAATCC‑TGACAGGCATCTATTGAATCCACCGGCCTCTTCGCGGGCAAGCCCGCTCCTACACCGTCGTAGGAGCGGGTTTACCCGCGAAAGGGCCTGCGGCCCTCACACAGCGACCAGCGCt                                                                       >  7:206661/1‑149 (MQ=255)
                                                                            cAAGATCGAACGGATCATTCATCCGCCCAAATCC‑TGACAGGCATCTATTGAATCCACCGGCCTCTTCGCGGGCAAGCCCGCTCCTACACCGTCGTAGGAGCGGGTTTACCCGCGAAAGGGCCTGCGGCCCTCACACAGCGACCAGCGCt                                                                       >  1:9991/1‑149 (MQ=255)
                                                                               gATCGAACGGATCATTCATCCGCCCAAATCC‑TGACAGGCATCTATTGAATCCACCGGCCTCTTCGCGGGCAAGCCCGCTCCTACACCGTCGTAGGAGCGGGTTTACCCGCGAAAGGGCCTGCGGCCCTCACACAGCGACCAGCGCttt                                                                     >  6:358411/1‑148 (MQ=255)
                                                                               gATCGAACGGATCATTCATCCGCCCAAATCC‑TGACAGGCATCTATTGAATCCACCGGCCTCTTCGCGGGCAAGCCCGCTCCTACACCGTCGTAGGAGCGGGTTTACCCGCGAAAGGGCCTGCGGCCCTCACACAGCGACCAGCGCTTTc                                                                    <  8:206661/149‑1 (MQ=255)
                                                                                    aaCGGATCATTCATCCGCCCAAATCC‑TGACAGGCATCTATTGAATCCACCGGCCTCTTCGCGGGCAAGCCCGCTCCTACACCGTCGTAGGAGCGGGTTTACCcgc                                                                                                           <  5:348465/105‑1 (MQ=255)
                                                                                    aaCGGATCATTCATCCGCCCAAATCC‑TGACAGGCATCTATTGAATCCACCGGCCTCTTCGCGGGCAAGCCCGCTCCTACACCGTCGTAGGAGCGGGTTTACCcgc                                                                                                           >  6:348465/1‑105 (MQ=255)
                                                                                          tcattcatCCGCCCAAATCC‑TGACAGGCATCTATTGAATCCACCGGCCTCTTCGCGGGCAAGCCCGCTCCTACACCGTCGTAGGAGCGGGTTTACCCGCGAAAGGGCCTGCgg                                                                                             >  4:51681/1‑113 (MQ=255)
                                                                                          tcattcatCCGCCCAAATCC‑TGACAGGCATCTATTGAATCCACCGGCCTCTTCGCGGGCAAGCCCGCTCCTACACCGTCGTAGGAGCGGGTTTACCCGCGAAAGGGCCTGCgg                                                                                             <  3:303023/113‑1 (MQ=255)
                                                                                          tcattcatCCGCCCAAATCC‑TGACAGGCATCTATTGAATCCACCGGCCTCTTCGCGGGCAAGCCCGCTCCTACACCGTCGTAGGAGCGGGTTTACCCGCGAAAGGGCCTGCgg                                                                                             >  4:303023/1‑113 (MQ=255)
                                                                                          tcattcatCCGCCCAAATCC‑TGACAGGCATCTATTGAATCCACCGGCCTCTTCGCGGGCAAGCCCGCTCCTACACCGTCGTAGGAGCGGGTTTACCCGCGAAAGGGCCTGCgg                                                                                             <  3:51681/113‑1 (MQ=255)
                                                                                            attcatCCGCCCAAATCC‑TGACAGGCATCTATTGAATCCACCGGCCTCTTCGCGGGCAAGCCCGCTCCTACACCGTCGTAGGAGCGGGTTTACCCGCGAAAGGGCCTGCGGCCCTCACACAGCGACCAGCGCTTTCAACGCCCCATc                                                         >  2:114002/1‑147 (MQ=255)
                                                                                            attcatCCGCCCAAATCC‑TGACAGGCATCTATTGAATCCACCGGCCTCTTCGCGGGCAAGCCCGCTCCTACACCGTCGTAGGAGCGGGTTTACCCGCGAAAGGGCCTGCGGCCCTCACACAGCGACCAGCGCTTTCAACGCCCCATCCa                                                       <  5:10367/149‑1 (MQ=255)
                                                                                             ttcatCCGCCCAAATCC‑TGACAGGCATCTATTGAATCCACCGGCCTCTTCGCGGGCAAGCCCGCTCCTACACCGTCGTAGGAGCGGGTTTACCCGCGAAAGGGCCTGCGGCCCTCACACAGCGACCAGCGCTTTCAACGCCCCATCCAg                                                      >  5:113258/1‑149 (MQ=255)
                                                                                             ttcatCCGCCCAAATCC‑TGACAGGCATCTATTGAATCCACCGGCCTCTTCGCGGGCAAGCCCGCTCCTACACCGTCGTAGGAGCGGGTTTACCCGCGAAAGGGCCTGCGGCCCTCACACAGCGACCAGCGCTTTCAACGCCCCATCCAg                                                      <  6:183548/149‑1 (MQ=255)
                                                                                               catCCGCCCAAATCC‑TGACTGGCTTCTATTGAATCCACCGGCCTCTTCGCGGGCAAGGCCGCTCCTACACCGTCGTAGGTGCGGGTTTACCCGCGATAGGGCCTGCGGCCCTCACACAGCGACCAGCGCTTTCAACGCCCCATCCAgg                                                     <  1:308504/148‑1 (MQ=255)
                                                                                                atCCGCCCAAATCC‑TGACAGGCATCTATTGAATCCACCGGCCTCTTCGCGGGCAAGCCCGCTCCTACACCGTCGTAGGAGCGGGTTTACCCGCGAAAGGGCCTGCGGCCCTCACACAGCGACCAGCGCTTTCAACGCCCCATCCAGGTc                                                   >  8:63415/1‑149 (MQ=255)
                                                                                                   ggCCCAAATCC‑TGACAGGCATCTATTGAATCCACCGGCCTCTTCGCGGGCAAGCCCGCTCCTACACCGTCGTAGGAGCGGGTTTACCCGCGAAAGGGCCTGCGGCCCTCACACAGCGACCAGCGCTTTCAACGCCCCATCCAGGTCCTc                                                <  3:88695/148‑1 (MQ=255)
                                                                                                   cGCCCAAATCC‑TGACAGGCATCTATTGAATCCACCGGCCTCTTCGCGGGCAAGCCCGCTCCTACACCGTCGTAGGAGCGGGTTTACCCGCGAAAGGGCCTGCGGCCCTCACACAGCGACCAGCGCTTTCAACGCCCCATCCAGGTCCTc                                                <  1:283333/149‑1 (MQ=255)
                                                                                                     cccAAATCC‑TGACAGGCATCTATTGAATCCACCGGCCTCTTCGCGGGCAAGCCCGCTCCTACACCGTCGTAGGAGCGGGTTTACCCGCGAAAGGGCCTGCGGCCCTCACACAGCGCCCAGCGCTTTCAACGCCCCCTCCAGGGCCTCGa                                              >  2:199009/1‑149 (MQ=255)
                                                                                                      ccAAATCC‑TGACAGGCATCTATTGAATCCACCGGCCTCTTCGCGGGCAAGCCCGCTCCTACACCGTCGTAGGAGCGGGTTTACCCGCGAAAGGGCCTGCGGCCCTCACACAGCGACCAGCGCTTTCAACGCCCCATCCAGGTCCTCGAc                                             >  5:96277/1‑149 (MQ=255)
                                                                                                       cAAATCC‑TGACAGGCATCTATTGAATCCACCGGCCTCTTCGCGGGCAAGCCCGCTCCTACACCGTCGTAGGAGCGGGTTTACCCGCGAAAGGGCCTGCGGCCCTCACACAGCGACCAGCGCTTTCAACGCCCCATCCAGGTCCTCGAc                                             >  3:49545/1‑148 (MQ=255)
                                                                                                        aaaTCC‑TGACAGGCATCTATTGAATCCACCGGCCTCTTCGCGGGCAAGCCCGCTCCTACACCGTCGTAGGAGCGGGTTTACCCGCGAAAGGGCCTGCGGCCCTCACACAGCGACCAGCGCTTTCAACGCCCCATCCAGGTCCTCGAc                                             >  6:84354/1‑147 (MQ=255)
                                                                                                             c‑TGACAGGCATCTATTGAATCCACCGGCCTCTTCGCGGGCTAGCCCGCTCCTACACCGTCGTAGGAGCGGGTTTACCCGCGAAAGGGCCTGCGGCCCTCACACAGCGACCAGCGCTTTCAACGCCCCATCCAGGTCCTCGACCAGGTCa                                      <  8:342403/149‑1 (MQ=255)
                                                                                                                 aCAGGCATCTATTGAATCCACCGGCCTCTTCGCGGGCAAGCCCGCTCCTACACCGTCGTAGGAGCGGGTTTACCCGCGAAAGGGCCTGCGGCCCTCACACAGCGACCAGCGCTTTCAACGCCCCATCCAGGTCCTCGACCAGGTCACgg                                   <  6:113258/149‑1 (MQ=255)
                                                                                                                 aCAGGCATCTATTGAATCCACCGGCCTCTTCGCGGGCAAGCCCGCTCCTACACCGTCGTAGGAGCGGGTTTACCCGCGAAAGGGCCTGCGGCCCTCACACAGCGACCAGCGCTTTCAACGCCCCATCCAGGTCCTCGACCAGGTCACgg                                   <  7:346264/149‑1 (MQ=255)
                                                                                                                 aCAGGCATCTATTGAATCCACCGGCCTCTTCGCGGGCAAGCCCGCTCCTACACCGTCGTAGGAGCGGGTTTACCCGCGAAAGGGCCTGCGGCCCTCACACAGCGACCAGCGCTTTCAACGCCCCATCCAGGTCCTCGACCAGGTCACg                                    <  6:233313/148‑1 (MQ=255)
                                                                                                                   aGGCATCTATTGAATCCACCGGCCTCTTCGCGGGCAAGCCCGCTCCTACACCGTCGTAGGAGCGGGTTTACCCGCGAAAGGGCCTGCGGCCCTCACACAGCGACCAGCGCTTTCAACGCCCCATCCAGGTCCTCGACCAGGTCACGGt                                  >  5:387353/1‑148 (MQ=255)
                                                                                                                     gCATCTATTGAATCCACCGGCCTCTTCGCGGGCAAGCCCGCTCCTACACCGTCGTAGGAGCGGGTTTACCCGCGAAAGGGCCTGCGGCCCTCACACAGCGACCAGCGCTTTCAACGCCCCATCCAGGTCCTCGACCAGGTCACGGTAAt                               >  3:284057/1‑149 (MQ=255)
                                                                                                                      cATCTATTGAATCCACCGGCCTCTTCGCGGGCAAGCCCGCTCCTACACCGTCGTAGGAGCGGGTTTACCCGCGAAAGGGCCTGCGGCCCTCACACAGCGACCAGCGCTTTCAACGCCCCATCCAGGTCCTCGACCAGGTCACGGTAATc                              <  2:9991/149‑1 (MQ=255)
                                                                                                                      cATCTATTGAATCCACCGGCCTCTTCGCGGGCAAGCCCGCTCCTACACCGTCGTAGGAGCGGGTTTACCCGCGAAAGGGCCTGCGGCCCTCACACAGCGACCAGCGCTTTCAACGCCCCATCCAGGTCCTCGACCAGGTCACGGTAATc                              <  7:101805/149‑1 (MQ=255)
                                                                                                                          tATTGAATCCACCGGCCTCTTCGCGGGCAAGCCCGCTCCTACACCGTCGTAGGAGCGGGTTTACCCGCCAAAGGGCCTGCGGCCCTCACACAGCGACCAGCGCATTCAACGCCCCAACCAGGTCCTCCACCAGCTCACGGTAATCCTc                           >  2:343347/1‑148 (MQ=38)
                                                                                                                            ttGAATCCACCGGCCTCTTCGCGGGCAAGCCCGCTCCTACACCGTCGTAGGAGCGGGTTTAccc                                                                                                             >  2:364129/1‑64 (MQ=255)
                                                                                                                            ttGAATCCACCGGCCTCTTCGCGGGCAAGCCCGCTCCTACACCGTCGTAGGAGCGGGTTTAccc                                                                                                             <  1:364129/64‑1 (MQ=255)
                                                                                                                            ttGAATCCACCGGCCTCTTCGCGGGCAAGCCCGCTCCTACACCGTCGTAGGAGCGGGTTTACCCGCGAAAGGGCCTGCGGCCCTCACACAGCGACCAGCGCTTTCAACGCCCCATCCAGGTCCTCGACCAGGTCACGGTAATc                              <  7:250606/143‑1 (MQ=255)
                                                                                                                            ttGAATCCACCGGCCTCTTCGCGGGCAAGCCCGCTCCTACACCGTCGTAGGAGCGGGTTTACCCGCGAAAGGGCCTGCGGCCCTCACACAGCGACCAGCGCTTTCAACGCCCCATCCAGGTCCTCGACCAGGTCACGGTAATc                              >  8:250606/1‑143 (MQ=255)
                                                                                                                            ttGAATCCACCGGCCTCTTCGCGGGCAAGCCCGCTCCTACACCGTCGTAGGAGCGGGTTTACCCGCGAAAGGGCCTGCGGCCCTCACACAGCGACCAGCGCTTTCAACGCCCCATCCAGGTCCTCGACCAGGTCACGGTAATCCTCGa                         <  3:28048/148‑1 (MQ=255)
                                                                                                                             tGAATCCACCGGCCTCTTCGCGGGCAAGCCCGCTCCTACACCGTCGTAGGAGCGGGTTTACCCGCGAAAGGGCCTGCGGCCCTCACACAGCGACCAGCGCTTTCAACGCCCCATCCAGGTCCTCGACCAGGTCACGGTAATCCTCGAt                        <  1:199009/148‑1 (MQ=255)
                                                                                                                                   cACCGGCCTCTTCGCGGGCAAGCCCGCTCCTACACCGTCGTAGGAGCGGGTTTAACCGCGAAAGGGCCTGCGGCCCTCACACCGCGACCAGCGCTATCAACGCCCCATCCAGGTCCTCGACCAGGACACGGTAATCCTCGATCCCTAcc                 >  6:77481/1‑149 (MQ=255)
                                                                                                                                    aCCGGCCTCTTCGCGGGCAAGCCCGCTCCTACACCGTCGTAGAAGCGGGTTTACCCGCGAAAGGGCCTGCGGCCCTCACACAGCGACCAGCGCTTTCAACGCCCCATCCAGGTCCTCGACCAGGTCACGGTAATCCTCGATCCCTAcc                 >  4:54870/1‑148 (MQ=255)
                                                                                                                                          ctctTCGCGGGCAAGCCCGCTCCTACACCGTCGTAGGAGCGGGTTTACCCGCGAAAGGGCCTGCGGCCCTCACACAGCGACCAGCGCTTTCAACGCCCCATCCAGGGCCTCGACCAGGTCACGGTAATCCTCGATCCCTACCGACAccc          >  6:391006/1‑149 (MQ=255)
                                                                                                                                           tctTCGCGGGCAAGCCCGCTCCTACACCGTCGTAGGAGCGGGTTTACCCGCGAAAGGGCCTGCGGCCCTCACACAGCGACCAGCGCTTTCAACGCCCCATCCAGGTCCTCGACCAGGTCACGGTAATCCTCGATCCCTACCGACAACCg         >  6:250568/1‑149 (MQ=255)
                                                                                                                                            tttCGCGGGTAAACCCGCTCCTACGACGGTGTAGGAGCGGGTTTACCCGCGAAAGGGCCTGCGGCCCTCACACAGCGACCAGCGCTTTCAACGCCCCATCCAGGTCCTCGACCAGGTCACGGTAATCCTCGATCCCTACCGACAACCg         >  8:206823/2‑148 (MQ=38)
                                                                                                                                             ttCGCGGGCAAGCCCGCTCCTACACCGTCGTAGGAGCGGGTTTACCCGCGAAAGGGCCTGCGGCCCTCACACAGCGACCAGCGCTTTCAACGCCCCATCCAGGTCCTCGACCAGGTCACGGTAATCCTCGATCCCTACCGACAACCg         <  5:225367/147‑1 (MQ=255)
                                                                                                                                               cgcgGGCAAGCCCGCTCCTACACCGTCGTAGGAGCGGGTTTACCCGCGAAAGGGCCTGCGGCCCTCACACAGCGACCAGCGCTTTCAACGCCCCATCCAGGTCCTCGACCAGGTCACGGTAATCCTCGATCCCTACCGACAACCgca       <  5:358411/147‑1 (MQ=255)
                                                                                                                                                 cgGGTAAACCCGCTCCTACGACGGTGTAGGAGCGGGCTTGCCCGCGAAAGGGCCTGCGGCCCTCACACAGCGACCAGCGCTTTCAACGCCCCATCCAGGTCCTCGACCAGGTCACGGTAATCCTCGATCCCTACCGACAACCGcagcag   >  3:160673/1‑149 (MQ=37)
                                                                                                                                                 cgGGCAAGCCCGCTCCTACACCGTCGTAGGAGCGGGTTTACCCGCGAAAGGGCCTGCGGCCCTCACACAGCGACCAGCGCTTTCAACGCCCCATCCAGGTCCTCGACCAGGTCACGGTAATCCTCGATCCCTACCGACAACCGcagcag   <  1:343347/149‑1 (MQ=255)
                                                                                                                                                  gggTAAACCCGCTCCTACGACGGTGTAGGAGCGGGCTTGCCCGCGAAAGGGCCTGCGGCCCTCACACAGCGACCAGCGCTTTCAACGCCCCATCCAGGTCCTCGACCAGGTCACGGTAATCCTCGACCCCTACCGACAACCGCAGCAgg  <  2:93104/149‑1 (MQ=37)
                                                                                                                                                  gggCAAGCCCGCTCCTACACCGTCGTAGGAGCGGGTTTACCCGCGAAAGGGCCTGCGGCCCTCACACAGCGACCAGCGCTTTCAACGCCCCATCCAGGTCCTCGACCAGGTCACGGTAATCCTCGATCCCTACCGACAACCGCAGCAgg  <  7:63415/149‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                 
TGGTAGCCGCCGAAGGCAACGGCACCAACGACGAAGCCGTCAACCGCCTGCGCAAGGAAAGGTTGCTGATCAAGGACAAGATCGAACGGATCATTCATCCGCCCAAATCC‑TGACAGGCATCTATTGAATCCACCGGCCTCTTCGCGGGCAAGCCCGCTCCTACACCGTCGTAGGAGCGGGTTTACCCGCGAAAGGGCCTGCGGCCCTCACACAGCGACCAGCGCTTTCAACGCCCCATCCAGGTCCTCGACCAGGTCACGGTAATCCTCGATCCCTACCGACAACCGCAGCAGG  >  NC_002947/5212324‑5212617

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