Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,212,475 G→A 5.4% intergenic (+39/+55) 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,4750GA 5.4% 287.7 / 4.2 92intergenic (+39/+55)PP_4593/PP_4594hypothetical protein/cystathionine gamma‑synthase
Reads supporting (aligned to +/- strand):  ref base G (43/44);  new base A (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 = 9.41e-01

GCCGCCGAAGGCAACGGCACCAACGACGAAGCCGTCAACCGCCTGCGCAAGGAAAGGTTGCTGATCAAGGACAAGATCGAACGGATCATTCATCCGCCCAAATCCTGACAGGCATCTATTGAATCCACCGGCCTCTTCGCGGGCAAGCCCGCTCCTACACCGTCGTAGGAGCGGGTTTACCCGCGAAAGGGCCTGCGGCCCTCACACAGCGACCAGCGCTTTCAACGCCCCATCCAGGTCCTCGACCAGGTCACGGTAATCCTCGATCCCTACCGACAACCGCAGCAGG  >  NC_002947/5212329‑5212617
                                                                                                                                                   |                                                                                                                                              
gccgccGAAGGCAACGGCACCAACGACGAAGCCGTCAACCGCCTGCGCAAGGAAAGGTTGCTGATCAAGGACAAGATCGAACGGATCATTCATCCGCCCAAATCCTGACAGGCATCTATTGAATCCACCGGCCTCTTCGCGGGCAAGc                                                                                                                                               <  1:167058/148‑1 (MQ=255)
   gccgAAGGCAACGGCACCAACGACGAAGCCGTCAACCGCCTGCGCAAGGAAAGGTTGCTGATCAAGGACAAGATCGAACGGATCATTCATCCGCCCAAATCCTGACAGGCATCTATTGAATCCACCGGCCTCTTCGCGGGCAAGCCCGc                                                                                                                                           >  5:91920/1‑149 (MQ=255)
   gccgAAGGCAACGGCACCAACGACGAAGCCGTCAACCGCCTGCGCAAGGAAAGGTTGCTGATCAAGGACAAGATCGAACGGATCATTCATCCGCCCAAATCCTGACAGGCATCTATTGAATCCACCGGCCTCTTCGCGGGCAAGCCCGc                                                                                                                                           <  6:91920/149‑1 (MQ=255)
    ccgAAGGCAACGGCACCAACGACGAAGCCGTCAACCGCCTGCGCAAGGAAAGGTTGCTGATCAAGGACAAGATCGAACGGATCATTCATCCGCCCAAATCCTGACAGGCATCTATTGAATCCACCGGCCTATTCGCGGGCAAGCCCGc                                                                                                                                           >  7:238885/1‑148 (MQ=255)
          gCAACGGCACCAACGACGAAGCCGTCAACCGCCTGCGCAAGGAAAGGTTGCTGATCAAGGACAAGATCGAACGGATCATTCATCCGCCCAAATCCTGACAGGCATCTATTGAATCCACCGGCCTCTTCGCGGGCAAGCCCGCTCCTaca                                                                                                                                    >  2:86665/1‑149 (MQ=255)
           cAACGGCACCAACGACGAAGCCGTCAACCGCCTGCGCAAGGAAAGGTTGCTGATCAAGGACAGATCGAACGGATCATTCATCCGCCCAAATCCTGACAGGCATCTATTGAATCCACCGGCCTCTTCGCGGGCAAGCCCGCTCCTACAcc                                                                                                                                  >  2:103570/1‑149 (MQ=255)
            aaCGGCACCAACGACGAAGCCGTCAACCGCCTGCGCAAGGAAAGGTTGCTGATCAAGGACAAGATCGAACGGATCATTCATCCGCCCAAATCCTGACAGGCATCTATTGAATCCACCGGCCTCTTCGCGGGCAAGCCCGCTCCTACAcc                                                                                                                                  >  2:283333/1‑149 (MQ=255)
               ggCACCAACGACGAAGCCGTCAACCGCCTGCGCAAGGAAAGGTTGCTGATCAAGGACAAGATCGAACGGATCATTCATCCGCCCAAATCCTGACAGGCATCTATTGAATCCACCGGCCTCTTCGCGGGCAAGCCCGCTCCTACACcgtc                                                                                                                               <  7:19125/149‑1 (MQ=255)
                 cACCAACGACGAAGCCGTCAACCGCCTGCGCAAGGAAAGGTTGCTGATCAAGGACAAGATCGAACGGATCATTCATCCGCCCAAATCCTGACAGGCATCTATTGAATCCACCGGCCTCTTCGCGGGCAAGCCCGCTCCTACACcgtcgt                                                                                                                             <  1:351008/149‑1 (MQ=255)
                   ccAACGACGAAGCCGTCAACCGCCTGCGCAAGGAAAGGTTGCTGATCATGGACAAGATCTAACGGATCATTCATCCGCCCAAATCCTTGACAGGCATCTATTGAATCCACCGGCCTCTTCGCGGGCAAGCCCGCTCCTACACCGTCGTa                                                                                                                            <  1:108921/149‑1 (MQ=255)
                   ccAACGACGAAGCCGTCAACCGCCTGCGCAAGGAAAGGTTGCTGATCAAGGACAAGATCGAACGGATCATTCATCCGCCCAAATCCTGACAGGCATCTATTGAATCCACCGGCCTCTTCGCGGGCAAGCCCGCTCCTACACCGTCGTAg                                                                                                                           <  1:86665/149‑1 (MQ=255)
                   ccAACGACGAAGCCGTCAACCGCCTGCGCAAGGAAAGGTTGCTGATCAAGGACAAGATCGAACGGATCATTCATCCGCCCAAATCCTGACAGGCATCTATTGAATCCACCGGCCTCTTCGCGGGCAAGCCCGCTCCTACACCGTCGTAg                                                                                                                           <  7:157335/149‑1 (MQ=255)
                   ccAACGACGAAGCCGTCAACCGCCTGCGCAAGGAAAGGTTGCTGATCAAGGACAAGATCGAACGGATCATTCATCCGCCCAAATCCTGACAGGCATCTATTGAATCCACCGGACTCTTCGCGGGCAAGCCCGCTCCTACACCGTCGTAg                                                                                                                           >  5:360740/1‑149 (MQ=255)
                     aaCGACGAAGCCGTCAACCGCCTGCGCAAGGAAAGGTTGCTGATCAAGGACAAGATCGAACGGATCATTCATCCGCCCAAATCCTGACAGGCATCTATTGAATCCACCGGCCTCTTCGCGGGCAAGCCCGCTCCTACACCGTCGTAg                                                                                                                           >  4:28048/1‑147 (MQ=255)
                       cgacgaAGCCGTCAACCGCCTGCGCAAGGAAAGGTTGCTGATCAAGGACAGATCGAACGGATCATTCATCCGCCCAAATCCTGACAGGCATCTATTGAATCCACCGGCCTCTTCGCGGGCAAGCCCGCTCCTACACCGTCGTAGGAGCg                                                                                                                      <  1:103570/149‑1 (MQ=255)
                         acgaAGCCGTCAACCGCCTGCGCAAGGAAAGGTTGCTGATCAAGGACAAGATCGAATGGATCATTCATCCGCCCAAATCCTGACAGGCATCTATTGAATCCACCGGCCTCTTCGCGGGCAAGCCCGCTCCTACACCGTCGTAGGAGCgg                                                                                                                     <  5:300553/149‑1 (MQ=255)
                         acgaAGCCGTCAACCGCCTGCGCAAGGAAAGGTTGCTGATCAAGGACAAGATCGAACGGATCATTCATCCGCCCAAATCCTGACAGGCATCTATTGAATCCACCGGCCTCTTCGCGGGCAAGCCCGCTCCTACACCGTCGTAGGAGCgg                                                                                                                     <  3:135831/149‑1 (MQ=255)
                          cgaAGCCGTCAACCGCCTGCGCAAGGAAAGGTTGCTGATCAAGGACAAGATCGAACGGATCATTCATCCGCCCAAATCCTGACAGGCATCTATTGAATCCACCGGCCTCTTCGCGGGCAAGCCCGCTCCTACACCGTCGTAGGAGCggg                                                                                                                    <  3:194897/149‑1 (MQ=255)
                           gaAGCCGTCAACCGCCTGCGCAAGGAAAGGTTGCTGATCAAGGACAAGATCGAACGGATCATTCATCCGCCCAAATCCTGACAGGCATCTATTGAATCCACCGGCCTCTTCGCGGGCAAGCCCGCTCCTACACCGTCGTAGGAGCGGGt                                                                                                                   >  3:327140/1‑149 (MQ=255)
                           gaAGCCGTCAACCGCCTGCGCAAGGAAAGGTTGCTGATCAAGGACAAGATCGAACGGATCATTCATCCGCCCAAATCCTGACAGGCATCTATTGAATCCACCGGCCTCTTCGCGGGCAAGCCCGCTCCTACACCGTCGTAGGAGCGGGt                                                                                                                   <  3:189720/149‑1 (MQ=255)
                            aaGCCGTCAACCGCCTGCGCAAGGAAAGGTTGCTGATCAAGGACAAGATCGAACGGATCATTCATCCGCCCAAATCCTGACAGGCATCTATTGAATCCACCGGCCTCTTCGCGGGCAAGCCCGCTCCTACACCGTCGTAGGAGCGGGtt                                                                                                                  <  4:47758/149‑1 (MQ=255)
                             aGCCGTCAACCGCCTGCGCAAGGAAAGGTTGCTGATCAAGGACAAGATCGAACGGATCATTCATCCGCCCAAATCCTGACAGGCATCTATTGAATCCACCGGCCTCTTCGCGGGCAAGCCCGCTCCTACACCGTCGTAGGAGCGGGtt                                                                                                                  >  6:225367/1‑148 (MQ=255)
                                 gTCAACCGCCTGCGCAAGGAAAGGTTGCTGATCAAGGACAAGATCGAACGGATCATTCATCCGCCCAAATCCTGACAGGCATCTATTGAATCCACCGGCCTCTTCGCGGGCAAGCCCGCTCCTACACCGTCGTAGGAGCGGGTTTAccc                                                                                                             <  4:327140/149‑1 (MQ=255)
                                 gTCAACCGCCTGCGCAAGGAAAGGTTGCTGATCAAGGACAAGATCGAACGGATCATTCATCCGCCCAAATCCTGACAGGCATCTATTGAATCCACCGGCCTCTTCGCGGGCAAGCCCGCTCCTACACCGTCGTAGGAGCGGGTTTAccc                                                                                                             <  8:336234/149‑1 (MQ=255)
                                  tCAACCGCCTGCGCAAGGAAAGGTTGCTGATCAAGGACAAGATCGAACGGATCATTCATCCGCCCAAATCCTGACAGGCATCTATTGAATCCACCGGCCTCTTCGCGGGCAAGCCCGCTCCTACACCGTCGTAGGAGCGGGTTTACCcg                                                                                                            <  2:174354/149‑1 (MQ=255)
                                  tCAACCGCCTGCGCAAGGAAAGGTTGCTGATCAAGGACAAGATCGAACGGATCATTCATCCGCCCAAATCCTGACAGGCATCTATTGAATCCACCGGCCTCTTCGCGGGCAAGCCCGCTCCTACACCGTCGTAGGAGCGGGTTTACAcg                                                                                                            >  5:45830/1‑149 (MQ=255)
                                                aaGGAAAGGTTGCTGATCAAGGACAAGATCGAACGGATCATTCATCCGCCCAAATCCTGACAGGCATCTATTGAATCCACCGGCCTCTTCGCGGGCAAGCCCGCTCCTACACCGTCGTAGGAGCGGGTGTACCCGCGAAAGGGCCTGCg                                                                                              >  8:349644/1‑149 (MQ=255)
                                                   gAAAGGTTGCTGATCAAGGACAAGATCGAACGGATCATTCATCCGCCCAAATCCTGACAGGCATCTATTGAATCCACCGGCCTCTTCGCGGGCAAGCCCGCTCCTACACCGTCGTAGGAGCGGGTTTACCCGCGAAAGGGCCTGCGGcc                                                                                           <  3:205369/149‑1 (MQ=255)
                                                   gAAAGGTTGCTGATCAAGGACAAGATCGAACGGATCATTCATCCGCCCAAATCCTGACAGGCATCTATTGAATCCACCGGCCTCTTCGCGGGCAAGCCCGCTCCTACACCGTCGTAGGAGCGGGTTTACCCGCGAAAGGGCCTGCGGcc                                                                                           <  7:349644/149‑1 (MQ=255)
                                                   gAAAGGTTGCTGATCAAGGACAAGATCGAACGGATCATTCATCCGCCCAAATCCTGACAGGCATCTATTGAATCCACCGGCCTCTTCGCGGGCAAGCCCGCTCCTACACCGTCGTAGGAGCGGGTTTACCCGCGAAAGGGCCTGCGGcc                                                                                           >  4:205369/1‑149 (MQ=255)
                                                         ttGCTGATCAAGGACAAGATCGAACCGATCATTCATCCGCCCAAATCCTGACAGGCATCTATTGAATCCACCGGCCTCTTCGCGGGCAAGCCCGCTCCTACACCGTCGTAGGAGCGGGTTTACCCGCGAAAGCGCCTGCGGCCCTcaca                                                                                     >  2:308504/1‑149 (MQ=255)
                                                           gCTGATCAAGGACAAGATCGAACGGATCATTCATCCGCCCAAATCCTGACAGGCATCTATTGAATCCACCGGCCTCTTCGCGGGCAAGCCCGCTCCTACACCGTCGTAGGAGCGGGTTTACCCGCGAAAGGGCCTGCGGCCCTcacaca                                                                                   >  7:199828/1‑149 (MQ=255)
                                                           gCTGATCAAGGACAAGATCGAACGGATCATTCATCCGCCCAAATCCTGACAGGCATCTATTGAATCCACCGGCCTCTTCGCGGGCAAGCCCGCTCCTACACCGTCGTAGGAGCGGGTTTACCCGCGAAAGGGCCTGCGGCCCTcaca                                                                                     >  8:101805/1‑147 (MQ=255)
                                                             tGATCAAGGACAAGATCGAACGGATCATTCATCCGCCCAAATCCTGACAGGCATCTATTGAATCCACCGGCCTCTTCGCGGGCAAGCCCGCTCCTACACCGTCGTAGGAGCGGGTTTACCCGCGAAAGGGCCTGc                                                                                               <  4:96440/135‑1 (MQ=255)
                                                             tGATCAAGGACAAGATCGAACGGATCATTCATCCGCCCAAATCCTGACAGGCATCTATTGAATCCACCGGCCTCTTCGCGGGCAAGCCCGCTCCTACACCGTCGTAGGAGCGGGTTTACCCGCGAAAGGGCCTGc                                                                                               >  3:96440/1‑135 (MQ=255)
                                                             tGATCAAGGACAAGATCGAACGGATCATTCATCCGCCCAAATCCTGACAGGCATCTATTGAATCCACCGGCCTCTTCGCGGGCAAGCCCGCTCCTACACCGTCGTAGGAGCGGGTTTACCCGCGAAAGGGCCGCAGGCCCtttcgc                                                                                    <  6:268857/146‑6 (MQ=255)
                                                             tGATCAAGGACAAGATCGAACGGATCATTCATCCGCCCAAATCCTGACAGGCATCTATTGAATCCACCGGCCTCTTCGCGGGCAAGCCCGCTCCTACACCGTCGTAGGAGCGGGTTTACCCGCGAAAGGGCCGCAGGCCCtttcgc                                                                                    >  5:268857/1‑141 (MQ=255)
                                                               aTCAAGGACAAGATCGAACGGATCATTCATCCGCCCAAATCCTGACAGGCATCTATTGAATCCACCGGCCTCTTCGCGGGTAAACCCGCTCCTACGACGGTGTAGGAGCGGGCTTGCCCGCGAAGAGGccggtgg                                                                                             >  1:88299/1‑130 (MQ=18)
                                                               aTCAAGGACAAGATCGAACGGATCATTCATCCGCCCAAATCCTGACAGGCATCTATTGAATCCACCGGCCTCTTCGCGGGTAAACCCGCTCCTACGACGGTGTAGGAGCGGGCTTGCCCGCGAAGAGGccggtgg                                                                                             <  2:88299/135‑6 (MQ=18)
                                                                 cAAGGACAAGATCGAACGGATCATTCATCCGCCCAAATCCTGACAGGCATCTATTGAATCCACCGGCCTCTTCGCGGGCAAGCCCGCTCCTACACCGTCGTAGGAGCGGGTTTACCCGCGAAAGGGCCTGCGGCCCTCACACAGCGAc                                                                              >  5:233313/1‑148 (MQ=255)
                                                                 cAAGGACAAGATCGAACGGATCATTCATCCGCCCAAATCCTGACAGGCATCTATTGAATCCACCGGCCTCTTCGCGGGCAAGCCCGCTCCTACACCGTCGTAGGAGCGGGTTTACCCGCGAAAGGGCCTGCGGCCCTCACACAGCGAc                                                                              >  3:256479/1‑148 (MQ=255)
                                                                  aaGGACAAGATCGAACGGATCATTCATCCGCCCAAATCCTGACAGGCATCTATTGAATCCACCGGCCTCTTCGCGGGCAAGCCCGCTCCTACACCGTCGTAGGAGCGGGTTTACCCGCGAAAGGGCCTGCGGCCCTCACACAGCGAc                                                                              >  4:88695/1‑147 (MQ=255)
                                                                    ggACAAGATCGAACGGATCATTCATCCGCCCAAATCCTGACAGGCATCTATTGAATCCACCGGCCTCTTCGCGGGCAAGCCCGCTCCTACACCGTCGTAGGAGCGGGTTTACCCGCGAAAGGGCCTGCGGCCCTCACACAGCGACCAgc                                                                          <  5:237260/149‑1 (MQ=255)
                                                                     gACAAGATCGAACGGATCATTCATCCGCCCAAATCCTGACAGGCATCTATTGAATCCACCGGCCTCTTCGCGGGCAAGCCCGCTCCTACACCGTCGTAGGAGCGGGTTTACCCGCGAAAGGGCCTGCGGCCCTCACACAGCGACCAgcg                                                                         >  8:346264/1‑149 (MQ=255)
                                                                       cAAGATCGAACGGATCATTCATCCGCCCAAATCCTGACAGGCATCTATTGAATCCACCGGCCTCTTCGCGGGCAAGCCCGCTCCTACACCGTCGTAGGAGCGGGTTTACCCGCGAAAGGGCCTGCGGCCCTCACACAGCGACCAGCGCt                                                                       >  1:9991/1‑149 (MQ=255)
                                                                       cAAGATCGAACGGATCATTCATCCGCCCAAATCCTGACAGGCATCTATTGAATCCACCGGCCTCTTCGCGGGCAAGCCCGCTCCTACACCGTCGTAGGAGCGGGTTTACCCGCGAAAGGGCCTGCGGCCCTCACACAGCGACCAGCGCt                                                                       >  7:206661/1‑149 (MQ=255)
                                                                          gATCGAACGGATCATTCATCCGCCCAAATCCTGACAGGCATCTATTGAATCCACCGGCCTCTTCGCGGGCAAGCCCGCTCCTACACCGTCGTAGGAGCGGGTTTACCCGCGAAAGGGCCTGCGGCCCTCACACAGCGACCAGCGCttt                                                                     >  6:358411/1‑148 (MQ=255)
                                                                          gATCGAACGGATCATTCATCCGCCCAAATCCTGACAGGCATCTATTGAATCCACCGGCCTCTTCGCGGGCAAGCCCGCTCCTACACCGTCGTAGGAGCGGGTTTACCCGCGAAAGGGCCTGCGGCCCTCACACAGCGACCAGCGCTTTc                                                                    <  8:206661/149‑1 (MQ=255)
                                                                               aaCGGATCATTCATCCGCCCAAATCCTGACAGGCATCTATTGAATCCACCGGCCTCTTCGCGGGCAAGCCCGCTCCTACACCGTCGTAGGAGCGGGTTTACCcgc                                                                                                           >  6:348465/1‑105 (MQ=255)
                                                                               aaCGGATCATTCATCCGCCCAAATCCTGACAGGCATCTATTGAATCCACCGGCCTCTTCGCGGGCAAGCCCGCTCCTACACCGTCGTAGGAGCGGGTTTACCcgc                                                                                                           <  5:348465/105‑1 (MQ=255)
                                                                                     tcattcatCCGCCCAAATCCTGACAGGCATCTATTGAATCCACCGGCCTCTTCGCGGGCAAGCCCGCTCCTACACCGTCGTAGGAGCGGGTTTACCCGCGAAAGGGCCTGCgg                                                                                             >  4:51681/1‑113 (MQ=255)
                                                                                     tcattcatCCGCCCAAATCCTGACAGGCATCTATTGAATCCACCGGCCTCTTCGCGGGCAAGCCCGCTCCTACACCGTCGTAGGAGCGGGTTTACCCGCGAAAGGGCCTGCgg                                                                                             >  4:303023/1‑113 (MQ=255)
                                                                                     tcattcatCCGCCCAAATCCTGACAGGCATCTATTGAATCCACCGGCCTCTTCGCGGGCAAGCCCGCTCCTACACCGTCGTAGGAGCGGGTTTACCCGCGAAAGGGCCTGCgg                                                                                             <  3:303023/113‑1 (MQ=255)
                                                                                     tcattcatCCGCCCAAATCCTGACAGGCATCTATTGAATCCACCGGCCTCTTCGCGGGCAAGCCCGCTCCTACACCGTCGTAGGAGCGGGTTTACCCGCGAAAGGGCCTGCgg                                                                                             <  3:51681/113‑1 (MQ=255)
                                                                                       attcatCCGCCCAAATCCTGACAGGCATCTATTGAATCCACCGGCCTCTTCGCGGGCAAGCCCGCTCCTACACCGTCGTAGGAGCGGGTTTACCCGCGAAAGGGCCTGCGGCCCTCACACAGCGACCAGCGCTTTCAACGCCCCATc                                                         >  2:114002/1‑147 (MQ=255)
                                                                                       attcatCCGCCCAAATCCTGACAGGCATCTATTGAATCCACCGGCCTCTTCGCGGGCAAGCCCGCTCCTACACCGTCGTAGGAGCGGGTTTACCCGCGAAAGGGCCTGCGGCCCTCACACAGCGACCAGCGCTTTCAACGCCCCATCCa                                                       <  5:10367/149‑1 (MQ=255)
                                                                                        ttcatCCGCCCAAATCCTGACAGGCATCTATTGAATCCACCGGCCTCTTCGCGGGCAAGCCCGCTCCTACACCGTCGTAGGAGCGGGTTTACCCGCGAAAGGGCCTGCGGCCCTCACACAGCGACCAGCGCTTTCAACGCCCCATCCAg                                                      <  6:183548/149‑1 (MQ=255)
                                                                                        ttcatCCGCCCAAATCCTGACAGGCATCTATTGAATCCACCGGCCTCTTCGCGGGCAAGCCCGCTCCTACACCGTCGTAGGAGCGGGTTTACCCGCGAAAGGGCCTGCGGCCCTCACACAGCGACCAGCGCTTTCAACGCCCCATCCAg                                                      >  5:113258/1‑149 (MQ=255)
                                                                                          catCCGCCCAAATCCTGACTGGCTTCTATTGAATCCACCGGCCTCTTCGCGGGCAAGGCCGCTCCTACACCGTCGTAGGTGCGGGTTTACCCGCGATAGGGCCTGCGGCCCTCACACAGCGACCAGCGCTTTCAACGCCCCATCCAgg                                                     <  1:308504/148‑1 (MQ=255)
                                                                                           atCCGCCCAAATCCTGACAGGCATCTATTGAATCCACCGGCCTCTTCGCGGGCAAGCCCGCTCCTACACCGTCGTAGGAGCGGGTTTACCCGCGAAAGGGCCTGCGGCCCTCACACAGCGACCAGCGCTTTCAACGCCCCATCCAGGTc                                                   >  8:63415/1‑149 (MQ=255)
                                                                                              ggCCCAAATCCTGACAGGCATCTATTGAATCCACCGGCCTCTTCGCGGGCAAGCCCGCTCCTACACCGTCGTAGGAGCGGGTTTACCCGCGAAAGGGCCTGCGGCCCTCACACAGCGACCAGCGCTTTCAACGCCCCATCCAGGTCCTc                                                <  3:88695/148‑1 (MQ=255)
                                                                                              cGCCCAAATCCTGACAGGCATCTATTGAATCCACCGGCCTCTTCGCGGGCAAGCCCGCTCCTACACCGTCGTAGGAGCGGGTTTACCCGCGAAAGGGCCTGCGGCCCTCACACAGCGACCAGCGCTTTCAACGCCCCATCCAGGTCCTc                                                <  1:283333/149‑1 (MQ=255)
                                                                                                cccAAATCCTGACAGGCATCTATTGAATCCACCGGCCTCTTCGCGGGCAAGCCCGCTCCTACACCGTCGTAGGAGCGGGTTTACCCGCGAAAGGGCCTGCGGCCCTCACACAGCGCCCAGCGCTTTCAACGCCCCCTCCAGGGCCTCGa                                              >  2:199009/1‑149 (MQ=255)
                                                                                                 ccAAATCCTGACAGGCATCTATTGAATCCACCGGCCTCTTCGCGGGCAAGCCCGCTCCTACACCGTCGTAGGAGCGGGTTTACCCGCGAAAGGGCCTGCGGCCCTCACACAGCGACCAGCGCTTTCAACGCCCCATCCAGGTCCTCGAc                                             >  5:96277/1‑149 (MQ=255)
                                                                                                  cAAATCCTGACAGGCATCTATTGAATCCACCGGCCTCTTCGCGGGCAAGCCCGCTCCTACACCGTCGTAGGAGCGGGTTTACCCGCGAAAGGGCCTGCGGCCCTCACACAGCGACCAGCGCTTTCAACGCCCCATCCAGGTCCTCGAc                                             >  3:49545/1‑148 (MQ=255)
                                                                                                   aaaTCCTGACAGGCATCTATTGAATCCACCGGCCTCTTCGCGGGCAAGCCCGCTCCTACACCGTCGTAGGAGCGGGTTTACCCGCGAAAGGGCCTGCGGCCCTCACACAGCGACCAGCGCTTTCAACGCCCCATCCAGGTCCTCGAc                                             >  6:84354/1‑147 (MQ=255)
                                                                                                        cTGACAGGCATCTATTGAATCCACCGGCCTCTTCGCGGGCTAGCCCGCTCCTACACCGTCGTAGGAGCGGGTTTACCCGCGAAAGGGCCTGCGGCCCTCACACAGCGACCAGCGCTTTCAACGCCCCATCCAGGTCCTCGACCAGGTCa                                      <  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                              <  7:101805/149‑1 (MQ=255)
                                                                                                                 cATCTATTGAATCCACCGGCCTCTTCGCGGGCAAGCCCGCTCCTACACCGTCGTAGGAGCGGGTTTACCCGCGAAAGGGCCTGCGGCCCTCACACAGCGACCAGCGCTTTCAACGCCCCATCCAGGTCCTCGACCAGGTCACGGTAATc                              <  2:9991/149‑1 (MQ=255)
                                                                                                                     tATTGAATCCACCGGCCTCTTCGCGGGCAAGCCCGCTCCTACACCGTCGTAGGAGCGGGTTTACCCGCCAAAGGGCCTGCGGCCCTCACACAGCGACCAGCGCATTCAACGCCCCAACCAGGTCCTCCACCAGCTCACGGTAATCCTc                           >  2:343347/1‑148 (MQ=38)
                                                                                                                       ttGAATCCACCGGCCTCTTCGCGGGCAAGCCCGCTCCTACACCGTCGTAGGAGCGGGTTTAccc                                                                                                             <  1:364129/64‑1 (MQ=255)
                                                                                                                       ttGAATCCACCGGCCTCTTCGCGGGCAAGCCCGCTCCTACACCGTCGTAGGAGCGGGTTTAccc                                                                                                             >  2:364129/1‑64 (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)
                                                                                                                                                   |                                                                                                                                              
GCCGCCGAAGGCAACGGCACCAACGACGAAGCCGTCAACCGCCTGCGCAAGGAAAGGTTGCTGATCAAGGACAAGATCGAACGGATCATTCATCCGCCCAAATCCTGACAGGCATCTATTGAATCCACCGGCCTCTTCGCGGGCAAGCCCGCTCCTACACCGTCGTAGGAGCGGGTTTACCCGCGAAAGGGCCTGCGGCCCTCACACAGCGACCAGCGCTTTCAACGCCCCATCCAGGTCCTCGACCAGGTCACGGTAATCCTCGATCCCTACCGACAACCGCAGCAGG  >  NC_002947/5212329‑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: