Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,901,761 A→T 9.8% D771V (GAC→GTC)  PP_1703 → assimilatory nitrate reductase/sulfite reductase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029471,901,7610AT 9.8% 228.7 / 13.0 89D771V (GAC→GTC) PP_1703assimilatory nitrate reductase/sulfite reductase
Reads supporting (aligned to +/- strand):  ref base A (45/35);  new base T (4/5);  total (49/40)
Fisher's exact test for biased strand distribution p-value = 7.26e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.32e-01

GACGCCCTTGATCTGAACACTGCCACTGTAGAGCCTGTCCGCATGCCGACCGCCACCCTGTCCCGTCTGCTTGGCCTGCACACCTTGCCGGTTCCGGTGCTGGCCGAAGA‑‑AGAGCGTCATTACCTGCAAGGTTTCCTGCTGGGCCTGGACCAGGCCCGCGCCGAGGGCGTACCCCGCCTGCCGGCCAGTGCGCCGCTGGCTGCCAGCCGTCGTTTGTTTGTCGACGGCTTGCTGGCCGGGCTGTTCGCCCAGCCTGCAACTGCACCCATCGCGGCACAGGCCGCGCCATGGCAT  >  NC_002947/1901613‑1901906
                                                                                                                                                        |                                                                                                                                                 
gACGCCCTTGATCTGAACACTGCCACTGTAGAGCCTGTCCGCATGCCGACCGCCACCCTGTCCCGTCTGCTTGGCCTGCACACCTTGCCGGTTCCGGTGCTGGCCGAAGA‑‑AGAGCGTCATTACCTGCAAGGTTTCCTGCTGGGCCTGGa                                                                                                                                                   <  4:204430/149‑1 (MQ=255)
 aCGCCCTTGATCTGAACACTGCCACTGTAGAGCCTGTCCGCATGCCGACCGCCACCCTGTCCCGTCTGCTTGGCCTGCACACCTTGCCGGTTCCGGTGCTGGCCGAAGA‑‑AGAGCGTCATTACCTGCAAGGTTTCCTGCTGGGCCTGGAc                                                                                                                                                  <  8:98372/149‑1 (MQ=255)
  cGCCCTTGATCTGAACACTGCCACTGTAGAGCCTGTCCGCATGCCGACCGCCACCCTGTCCCGTCTGCTTGGCCTGCACACCTTGCCGGTTCCGGTGCTGGCCGAAGA‑‑AGAGCGTCATTACCTGCAAGGTTTCCTGCTGGGCCTGGAcc                                                                                                                                                 <  5:133869/149‑1 (MQ=255)
  cGCCCTTGATCTGAACACTGCCACTGTAGAGCCTGTCCGCATGCCGACCGCCACCCTGTCCCGTCTGCTTGGCCTGCACACCTTGCCGGTTCCGGTGCTGGCCGAAGA‑‑AGAGCGTCATTACCTGCAAGGTTTCCTGCTGGGCCTGGAcc                                                                                                                                                 <  7:273818/149‑1 (MQ=255)
   gCCCTTGATCTGAACACTGCCACTGTAGAGCCTGTCCGCATGCCGACCGCCACCCTGTCCCGTCTGCTTGGCCTGCACACCTTGCCGGTTCCGGTGCTGGCCGAAGA‑‑AGAGCGTCATTACCTGCAAGGTTTCCTGCTGGGCCTGGTCCa                                                                                                                                                <  6:160494/149‑1 (MQ=255)
   gCCCTTGATCTGAACACTGCCACTGTAGAGCCTGTCCGCATGCCGACCGCCACCCTGTCCCGTCTGCTTGGCCTGCACACCTTGCCGGTTCCGGCGCTGGCCGAAGA‑‑AGAGCGTCATTACCTGCAAGGTTTCCTGCTGGGCCTGGAc                                                                                                                                                  >  2:157787/1‑147 (MQ=255)
    cccTTGATCTGAACACTGCCACTGTAGAGCCTGTCCGCATGCCGACCGCCACCCTGTCCCGTCTGCTTGGCCTGCACACCTTGCCGGTTCCGGTGCTGGCCGAAGA‑‑AGAGCGTCATTACCTGCAAGGTTTCCTGCTGGGCCTGGACCAg                                                                                                                                               >  2:198243/1‑149 (MQ=255)
    cccTTGATCTGAACACTGCCACTGTAGAGCCTGTCCGCATGCCGACCGCCACCCTGTCCCGTCTGCTTGGCCTGCACACCTTGCCGGTTCCGGTGCTGGCCGAAGA‑‑AGAGCGTCATTACCTGCAAGGTTTCCTGCTGGGCCTGGACCAg                                                                                                                                               <  1:46773/149‑1 (MQ=255)
             tGAACACTGCCACTGTAGAGCCTGTCCGCATGCCGACCGCCACCCTGTCCCGTCTGCTTGGCCTGCACACCTTGCCGGTTCCGGTGCTGGCCGAAGA‑‑AGAGCGTCATTACCTGCAAGGTTTCCTGCTGGGCCTGGACCAGGCCCGCGcc                                                                                                                                      >  6:201152/1‑149 (MQ=255)
               aaCCCTGCCACTGTAGAGCCTGTCCGCATGCCGACCGCCACCCTGTCCCGTCTGCTTGGCCTGCACACCTTGCCGGTTCCGGTGCTGGCCGAAGA‑‑AGAGCGTCATTACCTGCAAGGTTTCCTGCTGGGCCTGGTCCAGGcccagca                                                                                                                                       <  6:291287/146‑5 (MQ=255)
                       cACTGTAGAGCCTGTCCGCATGCCGACCGCCACCCTGTCCCGTCTGCTTGGCCTGCACACCTTGCCGGTTCCGGTGCTGGCCGAAGA‑‑AGAGCGTCATTACCTGCAAGGTTTCCTGCTGGGCCTGGACCAGGCCCGCGCCGAGGGCGTAc                                                                                                                            <  1:32598/149‑1 (MQ=255)
                       cACTGTAGAGCCTGTCCGCATGCCGACCGCCACCCTGTCCCGTCTGCTTGGCCTGCACACCTTGCCGGTTCCGGTGCTGGCCGAAGA‑‑AGAGCGTCATTACCTGCAAGGTTTCCTGCTGAGCCTGGACCAGGCCCGCGCCGAGGGCGTAc                                                                                                                            >  6:110177/1‑149 (MQ=255)
                        aCTGTAGAGCCTGTCCGCATGCCGACCGCCACCCTGTCCCGTCTGCTTGGCCTGCACACCTTGCCGGTTCCGGTGCTGGCCGAAGA‑‑AGAGCGTCATTACCTGCAAGGTTTCCTGCTGGGCCTGGTCCAGGcccagca                                                                                                                                       <  4:257646/137‑5 (MQ=255)
                        aCTGTAGAGCCTGTCCGCATGCCGACCGCCACCCTGTCCCGTCTGCTTGGCCTGCACACCTTGCCGGTTCCGGTGCTGGCCGAAGA‑‑AGAGCGTCATTACCTGCAAGGTTTCCTGCTGGGCCTGGTCCAGGcccagca                                                                                                                                       >  3:257646/1‑133 (MQ=255)
                               agCCTGTCCGCATGCCGACCGCCACCCTGTCCCGTCTGCTTGGCCTGCACACCTTGCCGGTTCCGGTGCTGGCCGAAGA‑‑AGAGCGTCATTACCTGCAAGGTTTCCTGCTGGGCCTGGACCAGGCCCGCGCCGAGGGCGTACCCCGCCTg                                                                                                                    >  8:198011/1‑149 (MQ=255)
                                 ccTGTCCGCATGCCGACCGCCACCCTGTCCCGTCTGCTTGGCCTGCACACCTTGCCGGTTCCGGTGCTGGCCGAAGA‑‑AGAGCGTCATTACCTGCAAGGTTTCCTGCTGGGCCTGGACCAGGCCCGCGCCGAGGGCGTACCCCgtacga                                                                                                                   >  8:242976/1‑143 (MQ=255)
                                  cTGTCCGCATGCCGACCGCCACCCTGTCCCGTCTGCTTGGCCTGCACACCTTGCCGGTTCCGGTGCTGGCCGAAGA‑‑AGAGCGTCATTACCTGCAAGGTTTCCTGCTGGGCCTGGACCAGGCCCGCGCCGAGGGCGTACACCGCCTGCCg                                                                                                                 >  7:196088/1‑149 (MQ=255)
                                   tGTCCGCATGCCGACCGCCACCCTGTCCCGTCTGCTTGGCCTGCACACCTTGCCGGTTCCGGTGCTGGCCGAAGA‑‑AGAGCGTCATTACCTGCAAGGTTTCCTGCTGGGCCTGGACCAGGCCCGCGCCGAGGGCGTACCCCGCCTGCCgg                                                                                                                >  8:10202/1‑149 (MQ=255)
                                    gTCCGCAGCCCGACCGCCACCCGGGCCCGTCTGGTTGGGCTGCACACGTTGGCGGTTCCGGTGCTGGCCTAATA‑‑AGTGGGTCATTACCTGCAAGGTTTCCTGCTGGGCCTGGACCCGGCCCGCGCCGAGGGCGTACCCCGCCTGCCGGc                                                                                                               <  2:369245/149‑1 (MQ=255)
                                            gCCGACCGCCACCCTGTCCCGTCTGCTTGGCCTGCACACCTTGCCGGTTCCGGTGCTGGCCGAAGAGTGGAGCGTAATTACCGGCAAGGTTTCCTGCGGGGCCGGGACCCGGCCCCCCCCGCGGGCGTACCCCCCCTGCCGGCCCGTgc                                                                                                         >  1:303428/1‑149 (MQ=255)
                                              cGACCGCCACCCTGTCCCGTCTGCTTGGCCTGCACACCTTGCCGGTTCCGGTGCTGGCCGAAGA‑‑AGAGCGTCATTACCTGCAAGGTTTCCTGCTGGGCCTGGACCAGGCCCGCGCCGAGGGCGTACCCCGCCTGCCGGCCAGTGcgccg                                                                                                     >  7:97702/1‑149 (MQ=255)
                                                aCCGCCACCCTGTCCCGTCTGCTTGGCCTGCACACCTTGCCGGTTCCGGTGCTGGCCGAAGA‑‑AGAGCGTCATTACCTGCAAGGTTTCCTGCTGGGCCTGGTCCAGGccc                                                                                                                                           >  5:53591/1‑109 (MQ=255)
                                                aCCGCCACCCTGTCCCGTCTGCTTGGCCTGCACACCTTGCCGGTTCCGGTGCTGGCCGAAGA‑‑AGAGCGTCATTACCTGCAAGGTTTCCTGCTGGGCCTGGTCCAGGccc                                                                                                                                           <  6:53591/109‑1 (MQ=255)
                                                 ccGCCACCCTGTCCCGTCTGCTTGGCCTGCACACCTTGCCGGTTCCGGTGCTGGCCGAAGA‑‑AGAGCGTCATTACCTGCAAGGTTTCCTGCTGGGCCTGGACCAGGCCCGCGCCGAGGGCGTACCCCGCCTGCCGGCCAGTGCGCCgctg                                                                                                  <  8:97702/149‑1 (MQ=255)
                                                   gCCACCCTGTCCCGTCTGCTTGGCCTGCACACCTTGCCGGTTCCGGTGCTGGCCGAAGA‑‑AGAGCGTACATTACATCCAAGGTTTCATCATGGCCATGGACCAGGCCCCCCCCAAGGGGGTCCCCCGCATGCCGGCCATTGCGCCgctg                                                                                                  >  1:303759/1‑148 (MQ=255)
                                                    ccACCCTGTCCCGTCTGCTTGGCCTGCACACCTTGCCGGTTCCGGTGCTGGCCGAAGA‑‑AGAGCGTCATTACCTGCAAGGTTTCCTGCTGGGCCTGGACCAGGCCCGCGCCGAGGGCGTACCCCGCCTGCCGGCCAGTGCGCCgctggct                                                                                               >  8:45964/1‑149 (MQ=255)
                                                        ccTGTCCCGTCTGCTTGGCCTGCACACCTTGCCGGTTCCGGTGCTGGCCGAAGA‑‑AGAGCGTCATTACCTGCAAGGTTTCCTGCTGGGCCTGGACCAGGCCCGCGCCGAGGGCGTACCCCGCCTGCCGGCCAGTGCGCCGCTGGCTGc                                                                                             >  1:14548/1‑147 (MQ=255)
                                                                   tGCTTGGCCTGCACACCTTGCCGGTTCCGGTGCTGGCCGAAGA‑‑AGAGCGTCATTACCTGCAAGGTTTCCTGCTGGGCCTGGACCAGGCCCGCGCCGAGGGCGTACCCCGCCTGCCGGCCAGTGCGCCGCTGGCTGCCAGCCGTCgtttg                                                                                >  7:200310/1‑149 (MQ=255)
                                                                    gCTTGGCCTGCACACCTTGCCGGTTCCGGTGCTGGCCGAAGA‑‑AGAGCGTCATTACCTGCAAGGTTTCCTGCTGGGCCTGGACCAGGCCCGCGCCGAGGGCGTACCCCGCCTGCCGGCCAGTGCGCCGCTGGCTGCCAGCCGTCgtttgt                                                                               >  5:212553/1‑149 (MQ=255)
                                                                    gCTTGGCCTGCACACCTTGCCGGTGCCGGTGCTGGCCGAAGA‑‑AGAGCGTCATTACCTGCAAGGTTTCCTGCTGGGCCTGGACCAGGCCCGCGCCGAGGGCGTACCCCGCCTGCCGGCCAGTGCGCCGCTGGCTGCCAGCCGTCgtttgt                                                                               <  2:303428/149‑1 (MQ=255)
                                                                     cTTGGCCTGCACACCTTGCCGGTTCCGGTGCTGGCCGAAGA‑‑AGAGCGTCATTACCTGCAAGGTTTCGTGCTGGGCCTGGACCAGGCCCGCGCCGAGGGCGTACCCCGCCTGCCGGCCAGTGCGCCGCTGGCTGCCAGCCGTCgttt                                                                                 <  3:229188/146‑1 (MQ=255)
                                                                     cTTGGCCTGCACACCTTGCCGGTTCCGGTGCTGGCCGAAGA‑‑AGAGCGTCATTACCTGCAAGGTTTCCTGCTGGGCCTGGAcc                                                                                                                                                 <  5:164055/82‑1 (MQ=255)
                                                                     cTTGGCCTGCACACCTTGCCGGTTCCGGTGCTGGCCGAAGA‑‑AGAGCGTCATTACCTGCAAGGTTTCCTGCTGGGCCTGGAcc                                                                                                                                                 >  6:164055/1‑82 (MQ=255)
                                                                     cTTGGCCTGCACACCTTGCCGGTTCCGGTGCTGGCCGAAGA‑‑AGAGCGTCATTACCTGCAAGGTTTCCTGCTGGGCCTGGACCAGGCCCGCGCCGAGGGCGTACCCCGCCTGCCGGCCAGTGCGCCGCTGGCTGCCAGCCGTCgtttgtt                                                                              <  5:325272/149‑1 (MQ=255)
                                                                     cTTGGCCTGCACACCTTGCCGGTTCCGGTGCTGGCCGAAGA‑‑AGAGCGTCATTACCTGCAAGGTTTCCTGCTGGGCCTGGACCAGGCCCGCGCCGAGGGCGTACCCCGCCTGCCGGCCAGTGCGCCGCTGGCTGCCAGCCGTCgtttgt                                                                               <  7:15927/148‑1 (MQ=255)
                                                                      ttGGCCTGCACACCTTGCCGTTTCGGGTGCTGGCCGATGA‑‑ATAGCGTGATTACCTGCATGGTTTCGTGTTGGGCGTGTTCCAGGCCCGCGCCGAGTGCGTACCCCGCCTGCCGGCCAGTGCGCCGCTGGCTGCCAGCCGTCgtttgttt                                                                             <  2:14548/149‑1 (MQ=255)
                                                                      ttGGCCTGCACACCTTGCCGGTTCCGGTGCTGGCCGAAGA‑‑AGAGCGTCATTACCTGCAAGGTTTCCTGCTGGGCCTGGACCAGGCCCGCGCCGAgg                                                                                                                                  <  2:31158/96‑1 (MQ=255)
                                                                      ttGGCCTGCACACCTTGCCGGTTCCGGTGCTGGCCGAAGA‑‑AGAGCGTCATTACCTGCAAGGTTTCCTGCTGGGCCTGGACCAGGCCCGCGCCGAgg                                                                                                                                  >  1:31158/1‑96 (MQ=255)
                                                                      ttGGCCTGCACACCTTGCCGGTTCCGGTGCTGGCCGAAGA‑‑AGAGCGTCATTACCTGCAAGGTTTCCTGCTGGGCCTGGACCAGGCCCGCGCCGAGGGCGTACCCCGCCTGCCGGCCAGTGCGCCGCTGGCTGCCAGCCGTCgtttgttt                                                                             >  6:324246/1‑149 (MQ=255)
                                                                        ggCCTGCACACCTTGCCGGTTCCGGTGCTGGCCGAAGA‑‑AGAGCGTCATTACCTGCAAGGTTTCCTGCTGGGCCTGGACCAGGCCCGCGCCGAGGGCGTACCCCGCCTGCCGGCCAGTGCGCCGCTGGCTGCCAGCCGTCgtttgttt                                                                             <  5:324246/147‑1 (MQ=255)
                                                                        ggCCTGCACACCTTGCCGGTTCCGGTGCTGGCCGAAGA‑‑AGAGCGTCATTACCTGCAAGGTTTCCTGCTGGGCCTGGACCAGGCCCGCGCCGAGGGCGTACCCCGCCTGCCGGCCAGTGCGCCGCTGGCTGCCAGCCGTCgtttgttt                                                                             <  1:198243/147‑1 (MQ=255)
                                                                         gCCTGCACACCTTGCCGGTTCCGGCGCTGGCCGAAGA‑‑TGAGCGTCATTACCTGCAAGGTTTCCTGCTGGGCGTGGACCAGGCCCGCGCCGAGGGCGTACCCCGCCTGCCGGCCAGTGCGCCGCTGGCTGCCAGCCGTCGTttgtttgt                                                                           <  1:157787/148‑1 (MQ=255)
                                                                           cTGCACACCTTGCCGGTTCCGGTGCTGGCCGAAGA‑‑AGAGCGTCATTACCTGCAAGGTTTCCTGCTGGGCCTGGACCAGGCCCGCGCCGAGGGCGTACCCCGCCTGCCGGCCAGTGCGCCGCTGGCTGCCAGCCGTCGTTTGTTTGTCGa                                                                        <  6:212553/149‑1 (MQ=255)
                                                                            tGCACACCTTGCCGGTTCCGGTGCTGGCCGAAGA‑‑AGAGCGTCATTACCTGCAAGGTTTCCTGCTGGGCCTGGACCAGGCCCGCGCCGAGGGCGTACCCCGCCTGCCGGCCAGTGCGCCGCTGGCTGCCAGCCGTCGTTTGTTTGTCGAc                                                                       >  6:135812/1‑149 (MQ=255)
                                                                            tGCACACCTTGCCGGTTCCGGTGCTGGCCGAAGA‑‑AGAGCGTCATTACCTGCAAGGTTTCCTGCTGGGCCTGGACCAGGCCCGCGCCGAGGGCGTACCCCGCCTGCCGGCCAGTGCGCCGCTGGCTGCCAGCCGTCGATTGTTTGTCGAc                                                                       >  5:159323/1‑149 (MQ=255)
                                                                              cacaCCTTGCCGGTTCCGGTGCTGGCCGAAGA‑‑AGAGCGTCATTACCCTGCAAGGTTTCCTGCTGGGCCTGGACCAGGCCCGCGCCGAg                                                                                                                                   <  1:263666/88‑1 (MQ=255)
                                                                              cacaCCTTGCCGGTTCCGGTGCTGGCCGAAGA‑‑AGAGCGTCATTACCTGCAAGGTTTCCTGCTGGGCCTGGACCAGGCCCGCGCCGAg                                                                                                                                   >  2:263666/1‑87 (MQ=255)
                                                                              cacaCCTTGCCGGTTCCGGTGCTGGCCGAAGA‑‑AGAGCGTCATTACCTGCAAGGTTTCCTGCTGGGCCTGGACCAGGCCCGCGCCGAGGGCGTACCCCGCCTGCCGGCCAGTGCGCCGCTGGCTGCCAGCCGTCGTTTGTTTGTCGACgg                                                                     >  2:104534/1‑149 (MQ=255)
                                                                                 acCTTGCCGGTTCCGGTGCTGGCCGAAGA‑‑AGAGCGTCATTACCTGCAAGGTTTCCTGCTGGGCCTGGACCAGGCCCGCGCCGAGGGCGTCCCCCGCCTGCCGGCCAGTGCGCCGCTGGCTGCCAGCCGTCGTTTGTTTGTCGACGGCtt                                                                  >  2:117744/1‑149 (MQ=255)
                                                                                 acCTTGCCGGTTCCGGTGCTGGCCGAAGA‑‑AGAGCGTCATTACCTGCAAGGTTTCCTGCTGGGCCTGGACCAGGCCAGCGCCGAGGGCGTACCCCGCCTGCCGGCCAGTGCGCCGCTGGCTGCCAGCCGTCGTTTGTTTGTCGACGGCtt                                                                  >  2:7946/1‑149 (MQ=255)
                                                                                   cTTGCCGGTTCCGGTGCTGGCCGAAGA‑‑AGAGCGTCATTACCTGCAAGGTTTCCTGCTGGGCCTGGACCAGGCCCGCGCCGAGGGCGTACCCCGCCTGCCGGCCAGTGCGCCGCTGGCTGCCAGCCGTCGTTTGTTTGTCGACGGCTTGc                                                                >  6:34413/1‑149 (MQ=255)
                                                                                   cTTGCCGGTTCCGGTGCTGGCCGAAGA‑‑AGAGCGTCATTACCTGCAAGGTTTCCTGCTGGGCCTGGACCAGGCCCGCGCCGAGGGCGTACCCCGCCTGCCGGCCAGTGCGCCGCTGGCTGCCAGCCGTCGTTTGTTTGTCGACGGCTTGc                                                                >  6:159419/1‑149 (MQ=255)
                                                                                        cGGTTCCGGTGCTGGCCGAAGA‑‑AGAGCGTCATTACCTGCAAGGTTTCCTGCTGGGCCTGGACCAGGCCCGCGCCGAGGGCGTACCCCGCCTGCCGGCCAGTGCGCAGCTGGCTGCCAGCCGTCGTTTGTTTGTCGACGGCTTGCTGGcc                                                           >  1:344853/1‑149 (MQ=255)
                                                                                              ctgTGCTGGCCGAAGA‑‑ATAGCTTCATTACGTGCAAGGTTTCATGCTGGGCCTGTACCAGGCCCGCGCCGAGGGCGTACCCCGCCTGCCGGCCAGTGCGCCGCTGGCTGCCAGCCGTCGTTTGTTTGTCGACGGCTTGCTGGCCGGGCt                                                      <  2:344853/146‑1 (MQ=255)
                                                                                               ggTGCTGGCCGAAGA‑‑AGAGCGTCATTACCTGCAAGGTTTCCTGCTGGGCCTGGACCAGGCCCGCGCCGAGGGCGTACCCCGCCTGCCGGCCAGTGCGCCGCTGGCTGCCAGCCGTCGTTTGTTTGTCGACGGCTTGCTGGCCGGTCTgc                                                    >  3:154190/1‑148 (MQ=255)
                                                                                                 tGCTGGCCGAAGA‑‑AGAGCGTCATTACCTGCAAGGTTTCCTGCTGGGCCTGGACCAGGCCCGCGCCGAGGGCGTACCCCGCCTGCCGGCCAGTGCGCCGCTGGCTGCCAGCCGTCGTTTGTTTGTCGACGGCTTGCTGGCCGGGCTGTTc                                                  <  4:160054/149‑1 (MQ=255)
                                                                                                    tGGCCGAAGA‑‑AGAGCGTCATTACCTGCAAGGTTTCCTGCTGGGCCTGGACCAGGCCCGCGCCGAGGGCGTACCCCGCCTGCCGGCCAGTGCGCCGCTGGCTGCCAGCCGTCGTTTGTTTGTCGACGGCTTGCTGGCCGGGCTGTTCGcc                                               >  8:283140/1‑149 (MQ=255)
                                                                                                     ggCCGAAGA‑‑AGAGCGTCATTACCTGCAAGGTTTCCTGCTGGGCCTGGACCAGGCCCGCGCCGAGGGCGTACCCCGCCTGCCGGCCAGTGCGCCGCTGGCTGCCAGCCGTCGTTTGTTTGTCGACGGCTTGCTGGCCGGGCTGTTCGccc                                              <  5:201152/149‑1 (MQ=255)
                                                                                                      gCCGAAGA‑‑AGAGCGTCATTACCTGCAAGGTTTCCTGCTGGGCCTGGACCAGGCCCGCGCCGAGGGCGTACCCCGCCTGCCGGCCAGTGCGCCGCTGGCTGCCAGCCGTCGTTTGTTTGTCGACGGCTTGCTGGCCGGGCTGTTCGccc                                              <  8:200310/148‑1 (MQ=255)
                                                                                                      gCCGAAGA‑‑AGAGCGTCATTACCTGCAAGGTTTCCTGCTGGGCCTGGACCAGGCCCGCGCCGAGGGCGTACCCCGCCTGCCGGCCAGTGCGCCGCTGGCTGCCAGCCGTCGTTTGTTTGTCGACGGCTTGCTGGCCGGGCTGTTCGCCCa                                             <  7:45964/149‑1 (MQ=255)
                                                                                                      gCCGAAGA‑‑AGAGCGTCATTACCTGCAAGGTTTCCTGCTGGGCCTGGACCAGGCCCGCGCCGAGGGCGTACCCCGCCTGCCGGCCAGTGCGCCGCTGGCTGCCAGCCGTCGTTTGTTTGTCGACGGCTTGCTGGCCGGGCTGTTCGCCCa                                             <  7:198011/149‑1 (MQ=255)
                                                                                                        cGAAGA‑‑AGAGCGTCATTACCTGCAAGGTTTCCTGCTGGGCCTGGACCAGGCCCGCGCCGAGGGCGTACCCCGCCTGCCGGCCAGTGCGCCGCTGGCTGCCAGCCGTCGTTTGTTTGTCGACGGCTTGCTGGCCGGGCTGTTCGCCCAGc                                           >  2:7583/1‑149 (MQ=255)
                                                                                                         gaaga‑‑aGAGCGTCATTACCTGCAAGGTTTCCTGCTGGGCCTGGACCAGGCCCGCGCCGAGGGCGTACCCCGCCTGCCGGCCAGTGCGCCGCTGGCTGCCAGCCGTCGTTTGTTTGTCGACGGCTTGCTGGCCGGGCTGTTCGCCCAGcc                                          >  2:344189/1‑149 (MQ=255)
                                                                                                          aaga‑‑agAGCGTCATTACCTGCAAGGTTTCCTGCTGGGCGTGGACCAGGCCCGCGCCGAGGGCGTACCCCGCCTGCCGGCCAGTGCGCCGCTGGCTGCCAGCCGTCGTTTGTTTGTCGACGGCTTGCTGGCCGGGCTGTTCGCCCAGCCt                                         <  1:7946/149‑1 (MQ=255)
                                                                                                           aga‑‑agaGCGTCATTACCTGCAAGGTTTCCTGCTGGGCCTGGACCAGGCCCGCGCCGAGGGCGTACCCCGCCTGCCGGCCAGTGCGCCGCTGGCTGCCAGCCGTCGTTTGTTTGTCGACGGCTTGCTGGCCGGGCTGTTCGCCCAGCCTg                                        <  8:196088/149‑1 (MQ=255)
                                                                                                           aga‑‑agaGCGTCATTACCTGCAAGGTTTCCTGCTGGGCCTGGACCAGGCCCGCGCCGAGGGCGTACCCCGCCTGCCGGCCAGTGCGCCGCTGGCTGCCAGCCGTCGTTTGTTTGTCGACGGCTTGCTGGCCGGGCTGTTCGCCCAGCCTg                                        <  5:159419/149‑1 (MQ=255)
                                                                                                                          ttACCTGCAAGGTTTCCTGCTGGGCCTGGACCAGGCCCGCGCCGAGGGCGTACCCCGCCTGCCGGCCAGTGCGCCGCTGGCTGCCAGCCGTCGTTTGTTTGTCGACGGCTTGCTGGCCGGGCTGTTCGCCCAGCCTGCAACTGCACCCa                            >  3:46594/1‑149 (MQ=255)
                                                                                                                          ttACCTGCAAGGTTTCCTGCTGAGCCTGGACCAGGCCCGCGCCGAGGGCGTACCCCGCCTGCCGGCCAGTGCGCCGCTGGCTGCCAGCCGTCGTTTGTTTGTCGACGGCTTGCTGGCCGGGCTGTTCGCCCAGCCTGCAACTGCACCCa                            <  5:110177/149‑1 (MQ=255)
                                                                                                                              ccTGCAAGGTTTCCTGCTGGGCCTGGACCAGGCCCGCGCCGAGGGCGTACCCCGCCTGCCGGCCAGTGCGCCGCTGGCTGCCAGCCGTCGTTTGTTTGTCGACGGCTTGCTGGCCGGGCTGTTCGCCCAGCCTGCAACTGCACCCATcg                         <  5:34413/149‑1 (MQ=255)
                                                                                                                               cTGCAAGGTTTCCTGCTGGGCCTGGACCAGGCCCGCGCCGAGGGCGTACCCCGCCTGCCGGCCAGTGCGCCGCTGGCTGCCAGCCGTCGTTTGTTTGTCGACGGCTTGCTGGCCGGGCTGTTCGCCCAGCCTGCAACTGCACCCATcgc                        >  7:235798/1‑149 (MQ=255)
                                                                                                                                tGCAAGGTTTCCTGCTGGGCCTGGACCAGGCCCGCGCCGAGGGCGTACCCCGCCTGCCGGCCAGTGCGCCGCTGGCTGCCAGCCGTCGTTTGTTTGTCGACGGCTTGCTGGCCGGGCTGTTCGCCCAGCCTGCAACTGCACCCATcgcg                       >  3:11134/1‑149 (MQ=255)
                                                                                                                                tGCAAGGTTTCCTGCTGGGCCTGGACCAGGCCCGCGCCGAGGGCGTACCCCGCCTGCCGGCCAGTGCGCCGCTGGCTGCCAGCCGTCGTTTGTTTGTCGACGGCTTGCTGGCCGGGCTGTTCGCCCAGCCTGCAACTGCACCCATcgcg                       >  2:152327/1‑149 (MQ=255)
                                                                                                                                tGCAAGGTTTCCTGCTGGGCCTGGACCAGGCCCGCGCCGAGGGCGTAACCCGCCTGCCGCCCAGTGAGCCGCTGGCTGCCAGCCGTCGTTTGTTTGTCGACGGCTTGCTGGCCGGGCTGTTCGCCCAGCCTGCAACTGCACCCATcgcg                       >  6:313774/1‑149 (MQ=255)
                                                                                                                                 gCAAGGTTTCCTGCTGGGCCTGGACCAGGCCCGCGCCGAGGGCGTACCCCGCCTGCCGGCCAGTGCGCCGCTGGCTGCCAGCCGTCGTTTGTTTGTCGACGGCTTGCTGGCCGGGCTGTTCGCCCAGCCTGCAACTGCCCCCATCGCgg                      >  8:43191/1‑149 (MQ=255)
                                                                                                                                  cAAGGTTTCCTGCTGGGCCTGGACCAGGCCCGCGCCGAGGGCGTACCCCGCCTGCCGGCCAGTGCGCCGCTGGCTGCCAGCCGTCGTTTGGTTGTCGACGGCTTGCTGGCCGGGCTGTTCGCCCAGCCTGCAACTGCACCCATCGCgg                      <  7:10202/148‑1 (MQ=255)
                                                                                                                                       tttCCTGCTGGGCCTGGACCAGGCCCGCGCCGAGGGCGTACCCCGCCTGCCGGCCAGTGCGCCGCTGGCTGCCAGCCGTCGTTTGTTTGTCGACGGCTTGCTGGCCGGGCTGTTCGCCCAGCCTGCAACTGCACCCATCGCGGCACAgg                <  1:7583/149‑1 (MQ=255)
                                                                                                                                            gcgcgggCCTGGTCCAGGCCCGCGCCGAGGGCGTACCCCGCCTGCCGGCCAGTGCGCCGCTGGCTGCCAGCCGTCGTTTGTTTGTCGACGGCTTGCTGGCCGGGCTGTTCGCCCAGCCTGCAACTGCACCCATCGCGGCACAGGCcgcg           >  5:184397/5‑149 (MQ=255)
                                                                                                                                             gctgGGCCTGGACCAGGCCCGCGCCGAGGGCGTACCCCGCCTGCCGGCCAGTGCGCCGCTGGCTGCCAGCCGTCGTTTGTTTGTCGACGGCTTGCTGGCCggg                                                        <  7:267838/103‑1 (MQ=255)
                                                                                                                                             gctgGGCCTGGACCAGGCCCGCGCCGAGGGCGTACCCCGCCTGCCGGCCAGTGCGCCGCTGGCTGCCAGCCGTCGTTTGTTTGTCGACGGCTTGCTGGCCggg                                                        >  8:267838/1‑103 (MQ=255)
                                                                                                                                              gcgggCCTGGTCCAGGCCCGCGCCGAGGGCGTACCCCGCCTGCCGGCCAGTGCGCCGCTGGCTGCCAGCCGTCGTTTGTTTGTCGACGGCTTGCTGGCCGGGCTGTTCGCCCAGCCTGCAACTGCACCCATCGCGGCACAGGCCgcgc          >  8:194847/3‑148 (MQ=255)
                                                                                                                                              ctgGGCCTGGACCAGGCCCGCGCCGAGGGCGTACCCCGCCTGCCGGCCAGTGCGCCGCTGGCTGCCAGCCGTCGTTTGTTTGTCGACGGCTTGCTGGCCGGGCTGTTCGCCCAGCCTGCAACTGCACCCATCGCGGCACAGGCcgcg           <  5:135812/147‑1 (MQ=255)
                                                                                                                                                gggCCTGGACCAGGCCCGCGCCGAGGGCGTACCCCGCCTGCCGGCCAGTGCGCCGCTGGCTGCCAGCCGTCGTTTGTTGGTCGACGGCTTGCTGGCCGGGCTGTTCGCCCAGCCTGCAACTGCACCCATCGCGGCTCAGGCCGCGCCAt       >  5:373342/1‑149 (MQ=255)
                                                                                                                                                 ggCCTGGACCAGGCCCGCGCCGAGGGCGTACCCCGCCTGCCGGCCAGTGCGCCGCTGGCTGCCAGCCGTCGTTTGTTTGTCGACGGCTTGCTGGCCGGGCTGTTCGCCCAGCCTGCAACTGCACCCATCGCGGCACAGGCCGCGCCATg      >  6:375381/1‑149 (MQ=255)
                                                                                                                                                   ccTGGACCAGGCCCGCGCCGAGGGCGTACCCCGCCTGCCGGCCAGTGCGCCGCTGGCTGCCAGCCGTCGTTTGTTTGTCGACGGCTTGCTGGCCGGGCTGTTCGCCCAGCCTGCAACTGCACCCATCGCgg                      >  8:220579/1‑131 (MQ=255)
                                                                                                                                                   ccTGGACCAGGCCCGCGCCGAGGGCGTACCCCGCCTGCCGGCCAGTGCGCCGCTGGCTGCCAGCCGTCGTTTGTTTGTCGACGGCTTGCTGGCCGGGCTGTTCGCCCAGCCTGCAACTGCACCCATCGCgg                      <  7:220579/131‑1 (MQ=255)
                                                                                                                                                   ccTGGACCAGGCCCGCGCCGAGGGCGTACCCCGCCTGCCGGCCAGTGCGCCGCTGGCTGCCAGCCGTCGTTTGTTTGTCGACGGCTTGCTGGCCGGGCTGTTCGCCCAGCCTGCAACTGCACCCATCGCGGCACAGGCCGCGCCATgg     >  3:268757/1‑148 (MQ=255)
                                                                                                                                                   ccTGGACCAGGCCCGCGCCGAGGGCGTACCCCGCCTGCCGGCCAGTGCGCCGCTGGCTGCCAGCCGTCGTTTGTTTGTCGACGGCTTGCTGGCCGGGCTGTTCGCCCAGCCTGCAACTGCACCCATCGCGGCACAGGCCGCGCCATgg     >  2:43592/1‑148 (MQ=255)
                                                                                                                                                    cTGGACCAGGCCCGCGCCGAGGGCGTACCCCGCCTGCCGGCCAGTGCGCCGCTGGCTGCCAGCcgtc                                                                                     <  4:216932/67‑1 (MQ=255)
                                                                                                                                                    cTGGACCAGGCCCGCGCCGAGGGCGTACCCCGCCTGCCGGCCAGTGCGCCGCTGGCTGCCAGCcgtc                                                                                     >  3:216932/1‑67 (MQ=255)
                                                                                                                                                     tGGACCAGGCCCGCGCCGAGGGCGTACCCCGCCTGCCGGCCAGTGCGCCGCTGGCTGCCAGCCGTCGTTTGTTTGTCGACGGCTTGCTGGCCGGGCTGTTCGCCAAGCCTGCAACTGCACCCATCGCGGCACAGGCCGCGCCATGGCAt  >  2:17079/1‑149 (MQ=255)
                                                                                                                                                        |                                                                                                                                                 
GACGCCCTTGATCTGAACACTGCCACTGTAGAGCCTGTCCGCATGCCGACCGCCACCCTGTCCCGTCTGCTTGGCCTGCACACCTTGCCGGTTCCGGTGCTGGCCGAAGA‑‑AGAGCGTCATTACCTGCAAGGTTTCCTGCTGGGCCTGGACCAGGCCCGCGCCGAGGGCGTACCCCGCCTGCCGGCCAGTGCGCCGCTGGCTGCCAGCCGTCGTTTGTTTGTCGACGGCTTGCTGGCCGGGCTGTTCGCCCAGCCTGCAACTGCACCCATCGCGGCACAGGCCGCGCCATGGCAT  >  NC_002947/1901613‑1901906

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: