Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 2,750,020 A→T 5.4% F157Y (TTC→TAC)  PP_2406 ← quinate/shikimate dehydrogenase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029472,750,0200AT 5.4% 232.2 / 2.5 74F157Y (TTC→TAC) PP_2406quinate/shikimate dehydrogenase
Reads supporting (aligned to +/- strand):  ref base A (39/30);  new base T (2/2);  total (41/33)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.61e-01

GGCCATGCCCACCGGCGTGGTGTTGACCAGCCCGTCCGCTTCGGCCAAGGCTGTGGCCAGGTCGGTGCCAAGCACGGCGCGCTCCGCGCCGAAATGGGTGTTCAGGTTGTCCACCAACGCCTGCGCACGGGTCGCATCCACTTCGAACAGCACCAGCCGTTCTACCCCTTCACCCAGCAGGGCATGGGCCACGGCCGAGCCAGCGCCACCGGCACCCATCTGCACCACCTGGCGCCGTGCAACGTCGGGCAGGCCACGGCGCAAGCCTTCGGCAAAGCCCAGGCAATCGGTGT  >  NC_002947/2749875‑2750167
                                                                                                                                                 |                                                                                                                                                   
ggCCATGCCCACCGGCGTGGTGTTGACCAGCCCGTCCGCTTCGGCCAAGGCTGTGGCCAGGTCGGTGCCAAGCACGGCGCGCTCCGCGCCGAAATGGGTGTTCAGGTTGTCCACCAACGCCTGCGCACGGGTCGCATCCACTTCGaa                                                                                                                                                    <  2:197843/147‑1 (MQ=255)
   cATGCCCACCGGCGTGGTGTTGACCAGCCCGTCCGCTTCGGCCAAGGCTGTGGCCAGGTCGGTGCCAAGCACGGCGCGCTCCGCGCCGAAATGGGTGTTCAGGTTGTCCACCAACGCCTGCGCACGGGTCGCATCCACTTCGAACAGCa                                                                                                                                               >  7:362106/1‑149 (MQ=255)
      gCCCACCGGCGTGGTGTTGACCAGCCCGTCCGCTTCGGCCAAGGCTGTGGCCAGGTCGGTGCCAAGCACGGCGCGCTCCGCGCCGAAATGGGTGTTCAGGTTGTCCACCAACGCCTGCGCACGGGTCGCATCCACTTCGAACAGCAc                                                                                                                                              >  4:191442/1‑147 (MQ=255)
         cACCGGCGTGGTGTTGACCAGCCCGTCCGCTTCGGCCAAGGCTGTGGCCAGGTCGGTGCCAAGCACGGCGCCCTCCGCGCCGAAATGGGTGTTCAGGTTGTCCACCAACGCCTGCGCACGGGTCGCATCCACTTCGAACAGCACCAGcc                                                                                                                                         <  1:128383/149‑1 (MQ=255)
          aCCGGCGTGGTGTTGACCAGCCCGTCCGCTTCGGCCAAGGCTGTGGCCAGGTCGGTGCCAAGCACGGCGCGCTCCGCGCCGAAATGGGTGTTCAGGTTGT‑CACCAACGCCTGCGCACGGGTCGCATCCACTTCGAACAGCACCAGCCGt                                                                                                                                       >  7:75169/1‑149 (MQ=255)
                gtggtgTTGACCAGCCCGTCCGCTTCGGCCAAGGCTGTGGCCAGGTCGGTGCCAAGCACGGCGCGCTCCGCGCCGAAATGGGTGTTCAGGTTGTCCACCAACGCCTGCGCACGGGTCGCATCCACTTCGAACAGCACCAGCCGTTCTAc                                                                                                                                  <  8:49058/149‑1 (MQ=255)
                     gtTGACCAGCCCGTCCGCTTCGGCCAAGGCTGTGGCCAGGTCGGTGCCAAGCACGGCGCGCTCCGCGCCGAAATGGGTGTTCAGGTTGTCCACCAACGCCTGCGCACGGGTCGCATCCACTTCGAACAGCACCAGCCGTTCTACCCCtt                                                                                                                             >  4:303999/1‑149 (MQ=255)
                        gACCAGCCCGTCCGCTTCGGCCAAGGCTGTGGCCAGGTCGGTGCCAAGCACGGCGCGCTCCGCGCCGAAATGGGTGTTCAGGTTGTCCACCAACGCCTGCGCACGGGTCGCATCCACTTCGAACAGCACCAGCCGTTCTACCCCTTCa                                                                                                                           <  4:115987/148‑1 (MQ=255)
                        gACCAGCCCGTCCGCTTCGGCCAAGGCTGTGGCCAGGTCGGTGCCAAGCACGGCGCGCTCCGCGCCGAAATGGGTGTTCAGGTTGTCCACCAACGCCTGCGCACGGGTCGCATCCACTTCGAACAGCACCAGCCGTTCTACCCCTTCAc                                                                                                                          >  2:411640/1‑149 (MQ=255)
                        gACCAGCCCGTCCGCTTCGGCCAAGGCTGTGGCCAGGTCGGTGCCAAGCACGGCGCGCTCCGCGCCGAAATGGGTGTTCAGGTTGTCCACCAACGCCTGCGCACGGGTCGCATCCACTTCGAACAGCACCAGCCGTTCTACCCCTTCAc                                                                                                                          >  1:368819/1‑149 (MQ=255)
                         aCCAGCCCGTCCGCTTCGGCCAAGGCTGTGGCCAGGTCGGTGCCAAGCACGGCGCGCTCCGCGCCGAAATGGGTGTTCAGGTTGTCCACCAACGCCTGCGCACGGGTCGCATCCACTTCGTAc                                                                                                                                                   >  4:326247/1‑123 (MQ=255)
                         aCCAGCCCGTCCGCTTCGGCCAAGGCTGTGGCCAGGTCGGTGCCAAGCACGGCGCGCTCCGCGCCGAAATGGGTGTTCAGGTTGTCCACCAACGCCTGCGCACGGGTCGCATCCACTTCGTAc                                                                                                                                                   >  4:320631/1‑123 (MQ=255)
                         aCCAGCCCGTCCGCTTCGGCCAAGGCTGTGGCCAGGTCGGTGCCAAGCACGGCGCGCTCCGCGCCGAAATGGGTGTTCAGGTTGTCCACCAACGCCTGCGCACGGGTCGCATCCACTTCGTAc                                                                                                                                                   <  3:326247/123‑1 (MQ=255)
                         aCCAGCCCGTCCGCTTCGGCCAAGGCTGTGGCCAGGTCGGTGCCAAGCACGGCGCGCTCCGCGCCGAAATGGGTGTTCAGGTTGTCCACCAACGCCTGCGCACGGGTCGCATCCACTTCGTAc                                                                                                                                                   <  3:320631/123‑1 (MQ=255)
                         aCCAGCCCGTCCGCTTCGGCCAAGGCTGTGGCCAGGTCGGTGCCAAGCACGGCGCGCTCCGCGCCGAAATGGGTGTTCAGGTTGTCCACCAACGCCTGCGCACGGGTCGCATCCACTTCGAACAGCACCAGCCGTTCTACCCCTTCAc                                                                                                                          <  8:362106/148‑1 (MQ=255)
                          ccAGCCCGTCCGCTTCGGCCAAGGCTGTGGCCAGGTCGGTGCCAAGCACGGCGCGCTCCGCGCCGAAATGGGTGTTCAGGTTGTCCACCAACGCCTGCGCACGGGTCGCATCCCCTTCGAACAGCACCAGCCGTTCTACCCCTTCAccc                                                                                                                        >  4:224956/1‑149 (MQ=255)
                                 gTCCGCTTCGGCCAAGGCTGTGGCCAGGTCGGTGCCAAGCACGGCGCGCTCCGCGCCGAAATGGGTGTTCAGGTTGTCCACCAACGCCTGCGCACGGGTCGCATCCACTTCGAACAGCACCAGCCGTTCTACCCCTTCACCCAGCAgg                                                                                                                  >  5:383794/1‑148 (MQ=255)
                                 gTCCGCTTCGGCCAAGGCTGTGGCCAGGTCGGTGCCAAGCACGGCGCGCTCCGCGCCGAAATGGGTGTTCAGGTTGTCCACCAACGCCTGCGCACGGGTCGCATCCACTTCGAACAGCACCAGCCGTTCTACCCCTTCACCCAGCAgg                                                                                                                  <  6:383794/148‑1 (MQ=255)
                                     gCTTCGGCCAAGGCTGTGGCCAGGTCGGTGCCAAGCACGGCGCGCTCCGCGCCGAAATGGGTGTTCAGGTTGTCCACCAACGCCTGCGCACGGGTCGCATCCACTTCGAACAGCACCAGCCGTTCTACCCCTTCACCCAGCAGGGCATg                                                                                                             <  3:191442/149‑1 (MQ=255)
                                     gCTTCGGCCAAGGCTGTGGCCAGGTCGGTGCCAAGCACGGCGCGCTCCGCGCCGAAATGGGTGTTCAGGTTGT‑CACCAACGCCTGCGCACGGGTCGCATCCACTTCGAACAGCACCAGCCGTTCTACCCCTTCACCCAGCAGGGCATgg                                                                                                            <  8:75169/149‑1 (MQ=255)
                                          ggCCAAGGCTGTGGCCAGGTCGGTGCCAAGCACGGCGCGCTCCGCGCCGAAATGGGTGTTCAGGTTGTCCACCAACGCCTGCGCACGGGTCGCATCCACTTCGAACAGCACCAGCCGTTCTACCCCTTCACCCAGCAGGGCATGGGCCa                                                                                                        >  5:311880/1‑149 (MQ=255)
                                           gCCAAGGCTGTGGCCAGGTCGGTGCCAAGCACGGCGCGCTCCGCGCCGAAATGGGTGTTCAGGTTGTCCACCAACGCCTGCGCACGGGTCGCATCCCCTTCGAACAGCACCAGCCGTTCTACCCCTTCACCCAGCAGGGCATGGGCCAc                                                                                                       >  7:275191/1‑149 (MQ=255)
                                           gCCAAGGCTGTGGCCAGGTCGGTGCCAAGCACGGCGCGCTCCGCGCCGAAATGGGTGTTCAGGTTGTCCACCAACGCCTGCGCACGGGTCGCATCCACTTCGAACAGCACCAGCCGTTCTAc                                                                                                                                  <  8:223171/122‑1 (MQ=255)
                                           gCCAAGGCTGTGGCCAGGTCGGTGCCAAGCACGGCGCGCTCCGCGCCGAAATGGGTGTTCAGGTTGTCCACCAACGCCTGCGCACGGGTCGCATCCACTTCGAACAGCACCAGCCGTTCTAc                                                                                                                                  >  7:223171/1‑122 (MQ=255)
                                           gCCAAGGCTGTGGCCAGGTCGGTGCCAAGCACGGCGCGCTCCGCGCCGAAATGGGTGTTCAGGTTGTCCACCAACGCCTGCGCACGGGTCGCATCCACTTCGAACAGCACCAGCCGTTCTACCCCTTCACCCAGCAGGGCATGGGCCAc                                                                                                       >  2:189868/1‑149 (MQ=255)
                                            ccAAGGCTGTGGCCAGGTCGGTGCCAAGCACGGCGCGCTCCGCGCCGAAATGGGTGTTCAGGTTGTCCACCAACGCCTGCGCACGGGTCGCATCCACTTCGAACAGCACCAGCCGTTCTACCCCTTCAccc                                                                                                                        >  2:3129/1‑131 (MQ=255)
                                            ccAAGGCTGTGGCCAGGTCGGTGCCAAGCACGGCGCGCTCCGCGCCGAAATGGGTGTTCAGGTTGTCCACCAACGCCTGCGCACGGGTCGCATCCACTTCGAACAGCACCAGCCGTTCTACCCCTTCAccc                                                                                                                        <  1:3129/131‑1 (MQ=255)
                                                ggCTGTGGCCAGGTCGGTGCCAAGCACGGCGCGCTCCGCGCCGAAATGGGTGTTCAGGTTGTCCACCAACGCCTGCGCACGGGTCGCATCCACTTCGAACAGCACCAGCCGTTCTACCCCTTCACCCAGCAGGGCATGGGCCACGGCCg                                                                                                  >  1:198166/1‑149 (MQ=255)
                                                 gCTGTGGCCAGGTCGGTGCCAAGCACGGCGCGCTCCGCGCCGAAATGGGTGTTCAGGTTGTCCACCAACGCCTGCGCACGGGTCGCATCCACTTCGAACAGCACCAGCCGTTCTACCCCTTCACCCAGCAGGGCATGGGCCACGGCCGa                                                                                                 <  3:117366/149‑1 (MQ=255)
                                                         catgTCGGTGGAAAGCACGGAGCTCTGCGCGGCGAAATGGGTGTTCAGGTTGTCCACCAAAGCTTGCGCTCGTGTCGCATCAACTTCCAACAGCACCAGGCGTTCTACCCCTTCACCCAGCAGGGCATGGGCCACGGCCGAGCCAgcgc                                                                                         <  6:311880/146‑1 (MQ=255)
                                                          aGGTCGGTGCCAAGCACGGCGCGCTCCGCGCCGAAATGGGTGTTCAGGTTGTCCACCAACGCCTGCGCACGGGTCGCATCCACTTCGAACAGCACCAGCCGTTCTACCCCTTCACCCAGCAGGGCATGGGCCACGGCCGAGCCAgcgc                                                                                         >  2:289555/1‑148 (MQ=255)
                                                             tCGGTGCCAAGCACGGCGCGCTCCGCGCCGAAATGGGTGTTCAGGTTGTCCACCAACGCCTGCGCACGGGTCGCATCCACTTCGAACAGCACCAGCCGTTCTACCCCTTCACCCAGCAGGGCATGGGCCACGGCCGAGCCAGCGCCAcc                                                                                     <  2:368819/149‑1 (MQ=255)
                                                                 tGCCAAGCACGGCGCGCTCCGCGCCGAAATGGGTGTTCAGGTTGTCCACCAACGCCTGCGCACGGGTCGCATCCACTTCGAACAGCACCAGCCGTTCTACCCCTTCACCCAGCAGGGCATGGGCCACGGCCGAGCCAGCGCCACCGGCa                                                                                 >  3:177907/1‑149 (MQ=255)
                                                                   ccAAGCACGGCGCGCTCCGCGCCGAAATGGGTGTTCAGGTTGTCCACCAACGCCTGCGCACGGGTCGCATCCACTTCGAACAGCACCAGCCGTTCTACCCCTTCACCCAGCAGGGCATGGGCCACGGCCGAGCCAGCGCCACCGGCAc                                                                                <  3:303999/148‑1 (MQ=255)
                                                                    cAAGCACGGCGCGCTCCGCGCCGAAATGGGTGTTCAGGTTGTCCACCAACGCCTGCGCACGGGTCGCATCCACTTCGAACAGCACCAGCCGTTCTACCCCTTCACCCAGCAGGGCATGGGCCACGGCCGAGCCAGCGCCACCGGCAcc                                                                               >  2:58967/1‑148 (MQ=255)
                                                                      aGCACGGCGCGCTCCGCGCCGAATTGGGTGTTCAGGTTGTCCACCAACGCCTGCGCACGGGTCGCATCCACTTCGCACAGCACCAGCCGTTCTACCCCTTCACCCAGCAGGGCATGGGCCACGGCCGAGCCAGCGCCACCGGCACCCAt                                                                            <  7:164088/149‑1 (MQ=255)
                                                                            gcgcgcTCCGCGCCGAAATGGGTGTTCAGGTTGTCCACCAACGCCTGCGCACGGGTCGCATCCACTTCGAACAGCACCAGCCGTTCTACCCCTTCACCCAGCAGGGCATGGGCCACGGCCGAGCCAGCGCCACCGGCACCCATCTGcac                                                                      <  2:198166/149‑1 (MQ=255)
                                                                                  tCCGCGCCGAAATGGGTGTTCAGGTTGTCCACCAACGCCTGCGCACGGGTCGCATCCACTTCGAACAGCACCAGCCGTTCTACCCCTTCACCCAGCAGGGCATGGGCCACGGCCGAGCCAGCGCCACCGGCACCCATCTGCACCACCTg                                                                >  7:336065/1‑149 (MQ=255)
                                                                                  tCCGCGCCGAAATGGGTGTTCAGGTTGTCCACCAACGCCTGCGCACGGGTCGCATCCACTTCGAACAGCACCAGCCGTTCTACCCCTTCACCCAGCAGGGCATGGGCCACGGCCGAGCCAGCGCCACCGGCACCCATCTGCACCACCTg                                                                >  1:102803/1‑149 (MQ=255)
                                                                                    cgcgCCGAAATGGGTGTTCAGGTTGTCCACCAACGCCTGCGCACGGGTCGCATCCACTTCGAACAGCACCAGCCGTTCTACCCCTTCACCCAGCAGGGCATGGGCCACGGCCGAGCCAGCGCCACCGGCACCCATCTGCACCACCTGgc                                                              >  3:224367/1‑149 (MQ=255)
                                                                                        ccGAAATGGGTGTTCAGGTTGTCCACCAACGCCTGCGCACGGGTCGCATCCACTTCGAACAGCACCAGCCGTTCTACCCCTTCACCCAGCAGGGCATGGGCCACGGCCGAGCCAGCGCCACCGGCACCCATCTGCACCACCTGGCGCCg                                                          <  4:224367/149‑1 (MQ=255)
                                                                                         cGAAATGGGTGTTCAGGTTGTCCACCAACGCCTGCGCACGGGTCGCATCCACTTCGAACAGCACCAGCCGTTCTACCCCTTCACCCAGCAGGGCATGGGCCACGGCCGAGCCAGCGCCACCGGCACCCATCTGCACCACCTGGCGCCGt                                                         <  4:186832/149‑1 (MQ=255)
                                                                                               gggTGTTCAGGTTGTCCACCAACGCCTGCGCACGGGTCGCATCCACTTCGAACAGCACCAGCCGTTCTACCCCTTCACCCAGCAGGGCATGGGCCACGGCCGAGCCAGCGCCACCGGCACCCATCTGCACCACCTGGCGCCGTGCAACg                                                   >  4:402693/1‑149 (MQ=255)
                                                                                                    ttCAGGTTGTCCACCAACGCCTGCGCACGGGTCGCATCCACTTCGAACAGCACCAGCCGTTCTACCCCTTCACCCAGCAGGGCATGGGCCACGGCCGAGCCAGCGCCACCGGCACCCATCTGCACCACCTGGCGCCGTGCAAc                                                    >  4:359839/1‑143 (MQ=255)
                                                                                                    ttCAGGTTGTCCACCAACGCCTGCGCACGGGTCGCATCCACTTCGAACAGCACCAGCCGTTCTACCCCTTCACCCAGCAGGGCATGGGCCACGGCCGAGCCAGCGCCACCGGCACCCATCTGCACCACCTGGCGCCGTGCAACGTCggg                                              >  6:17301/1‑149 (MQ=255)
                                                                                                    ttCAGGGTGTCCACCAACGCCTGCGCACGGGTCGCATCCACTTCGAACAGCACCAGCCGTTCTACCCCTTCACCCAGCAGGGCATGGGCCACGGCCGAGCCAGCGCCACCGGCACCCATCTGCACCACCTGGCGCCGTGCAAc                                                    <  3:359839/143‑1 (MQ=255)
                                                                                                      cAGGTTGTCCACCAACGCCTGCGCACGGGTCGCATCCACTTCGAACAGCACCAGCCGTTCTACCCCTTCACCCAGCAGGGCATGGGCCACGGCCGAGCCAGCGCCACCGGCACCCATCTGCACCACCTGGCGCCGTGCAACGTCggg                                              >  5:36515/1‑147 (MQ=255)
                                                                                                               caccaACGCCTGCGCACGGGTCGCATCCACTTCGAACAGCACCAGCCGTTCTACCCCTTCACCCAGCAGGGCATGGGCCACGGCCGAGCCAGCGCCACCGGCACCCATCTGCACCACCTGGCGCCGTGCAACGTCGGGCAGGCCACgg                                    <  2:102803/148‑1 (MQ=255)
                                                                                                               caccaACGCCTGCGCACGGGTCGCATCCACTTCGAACAGCACCAGCCGTTCTACCCCTTCACCCAGCAGGGCATGGGCCACGGCCGAGCCAGCGCCACCGGCACCCATCTGCACCACCTGGCGCCGTGCAACGTCGGGCAGGCCACGgc                                   >  6:211004/1‑149 (MQ=255)
                                                                                                                accaACGCCTGCGCACGGGTCGCATCCACTTCGAACAGCACCAGCCGTTCTACCCCTTCACCCAGCAGGGCATGGGCCACGGCCGAGCCAGCGCCACCGGCACCCATCTGCACCACCTGGCGCCGTGCAACGTCGGGCAGGCCACGgcg                                  >  2:254249/1‑149 (MQ=255)
                                                                                                                accaACGCCTGCGCACGGGTCGCATCCACTTCGAACAGCACCAGCCGTTCTACCCCTTCACCCAGCAGGGCATGGGCCACGGCCGAGCCAGCGCCACCGGCACCCATCTGCACCACCTGGCGCCGTGCAACGTCGGGCAGGCCACGgcg                                  >  5:235023/1‑149 (MQ=255)
                                                                                                                 ccaACGCCTGCGCACGGGTCGCATCCACTTCGAACAGCACCAGCCGTTCTACCCCTTCACCCAGCAGGGCATGGGCCACGGCCGAGCCAGCGCCACCGGCACCCATCTGCACCACCTGGCGCCGTGCAACGTCGGGCAGGCCACGgcg                                  >  8:30591/1‑148 (MQ=255)
                                                                                                                   aaCGCCTGCGCACGGGTCGCATCCACTTCGAACAGCACCAGCCGTTCTACCCCTTCACCCAGCAGGGCATGGGCCACGGCCGAGCCAGCGCCACCGGCACCCATCTGCACCACCTGGCGCCGTGCAACGTCGGGCAGGCCACGGCGCaa                               <  3:402693/149‑1 (MQ=255)
                                                                                                                   aaCGCCTGCGCACGGGTCGCATCCACTTCGAACAGCACCAGCCGTTCTACCCCTTCACCCAGCAGGGCATGGGCCACGGCCGAGCCAGCGCCACCGGCACCCATCTGCACCACCTGGCGCCGTGCAACGTCGGGCAGGCCACGGCGCaa                               <  1:189868/149‑1 (MQ=255)
                                                                                                                      gCCTGCGCACGGGTCGCATCCACTTCGAACAGCACCAGCCGTTCTACCCCTTCACCCAGCAGGGCATGGGCCACGGCCGAGCCAGCGCCACCGGCACCCATCTGCACCACCTGGCGCCGTGCAACGTCGGGCAGGCCACGGCGCAAGcc                            >  4:148744/1‑149 (MQ=255)
                                                                                                                        cTGCGCACGGGTCGCATCCACTTCGAACAGCACCAGCCGTTCTACCCCTTCACCCAGCAGGGCATGGGCCACGGCCGAGCCAGCGCCACCGGCACCCATCTGc                                                                        >  6:207687/1‑103 (MQ=255)
                                                                                                                        cTGCGCACGGGTCGCATCCACTTCCAACAGCACCAGCCGTTCTACCCCTTCACCCAGCAGGGCATGGGCCACGGCCGAGCCAGCGCCACCGGCACCCATCTGc                                                                        <  5:207687/103‑1 (MQ=255)
                                                                                                                           cgcACGGGTCGCATCCACTTCGAACAGCACCAGCCGTTCTACCCCTTCACCCAGCAGGGCATGGGCCACGGCCGAGCCAGCGCCACCGGCACCCATCTGCACCACCTGGCGCCGTGCAACGTCGGGCAGGCCACGGCGCAAGCCTTCgg                       <  5:17301/149‑1 (MQ=255)
                                                                                                                            gcACGGGTCGCATCCACTTCGAACAGCACCAGCCGTTCTACCCCTTCACCCAGCAGGGCATGGGCCACGGCCGAGCCa                                                                                             >  6:400778/1‑78 (MQ=255)
                                                                                                                              aCGGGTCGCATCCACTTCGAACAGCACCAGCCGTTCTACCCCTTCACCCAGCAGGGCATGGGCCACGGCCGAGCCAGCGCCACCGGCACCCATCTGCACCACCTGGCGCCGTGCAACGTCGGGCAGGCCACGGCGCAAGCCTTCGGCa                     >  4:329880/1‑148 (MQ=255)
                                                                                                                                   tCGCATCCACTTCGAACAGCACCAGCCGTTCTACCCCTTCACCCAGCAGGGCATGGGCCACGGCCGAGCCAGCGCCACCGGCACCCATCTGCACCACCTGGCGCCGTGCAACGTCGGGCAGGCCACGGCGCAAGCCTTCGGCAAAGcc                >  2:100744/1‑148 (MQ=255)
                                                                                                                                     gCATCCACTTCGAACAGCACCAGCCGTTCTACCCCTTCACCCAGCAGGGCATGGGCCACGGCCGAGCCAGCGCCACCGGCACCCATCTGCACCACCTGGCGCCGTGCAACGTCGGGCAGGCCACGGCGCAAGCCTTCGGCAAAGCCCAg             <  1:58967/149‑1 (MQ=255)
                                                                                                                                      cATCCACTTCGAACAGCACCAGCCGTTCTACCCCTTCACCCAGCAGGGCATGGGCCACGGCCGAGCCAGCGCCACCGGCACCCATCTGCACCACCTGGCGCCGTGCAACGTCGGGCAGGCCACGGCGCAAGCCTTCGGCAAAGCCCAgg            >  2:399792/1‑149 (MQ=255)
                                                                                                                                       aTCCACTTCGAACAGCACCAGCCGTTCTACCCCTTCACCCAGCAGGGCATGGGCCACGGCCGAGCCAGCGCCACCGGCACCCATCTGCACCACCTGGCGCCGTGCAACGTCGGGCAGGCCACGGCGCAAGCCTTCGGCAAAGCCCAgg            >  4:261193/1‑148 (MQ=255)
                                                                                                                                       aTCCACTTCGAACAGCACCAGCCGTTCTACCCCTTCACCCAGCAGGGCATGGGCCACGGCCGAGCCAGCGCCACCGGCACCCATCTGCACCACCTGGCGCCGTGCAACGTCGGGCAGGCCACGGCGCAAGCCTTCGGCAAAGCCCAGGc           <  1:399792/149‑1 (MQ=255)
                                                                                                                                          cACTTCGAACAGCACCAGCCGTTCTACCCCTTCACCCAGCAGGGCATGGGCCACGGCCGAGCCAGCGCCACCGGCACCCATCTGCACCACCTGGCGCCGTGCAACGTCGGGCAGGCCACGGCGCAAGCCTTCGGCAAAGCCCAg             >  6:200112/1‑144 (MQ=255)
                                                                                                                                          cACTTCGAACAGCACCAGCCGTTCTACCCCTTCACCCAGCAGGGCATGGGCCACGGCCGAGCCAGCGCCACCGGCACCCATCTGCACCACCTGGCGCCGTGCAACGTCGGGCAGGCCACGGCGCAAGCCTTCGGCAAAGCCCAg             <  5:200112/144‑1 (MQ=255)
                                                                                                                                             ttCGAACAGCACCAGCCGTTCTACCCCTTCACCCAGCAGGGCATGGGCCACGGCCGAGCCa                                                                                             <  4:12486/61‑1 (MQ=255)
                                                                                                                                             ttCGAACAGCACCAGCCGTTCTACCCCTTCACCCAGCAGGGCATGGGCCACGGCCGAGCCa                                                                                             >  3:12486/1‑61 (MQ=255)
                                                                                                                                               cGAACAGCACCAGCCGTTCTACCCCTTCACCCAGCAGGGCATGGGCCACGGCCGAGCCAGCGCCACCGGCACCCATCTGCACCACCTGGCGCCGTGCAACGTCGGGCAGGCCACGGCGCAAGCCTTCGGCAAAGCCCAg             <  4:131086/139‑1 (MQ=255)
                                                                                                                                               cGAACAGCACCAGCCGTTCTACCCCTTCACCCAGCAGGGCATGGGCCACGGCCGAGCCAGCGCCACCGGCACCCATCTGCACCACCTGGCGCCGTGCAACGTCGGGCAGGCCACGGCGCAAGCCTTCGGCAAAGCCCAg             >  3:131086/1‑139 (MQ=255)
                                                                                                                                               cGAACAGCACCAGCCGTTCTACCCCTTCACCCAGCAGGGCATGGGCCACGGCCGAGCCAGCGCCACCGGCACCCATCTGCACCACCTGGCGCCGTGCAACGTCGGGCAGGCCACGGCGCAAGCCTTCGGCAAAGCCCAGGCAATc       >  1:186360/1‑145 (MQ=255)
                                                                                                                                               cGAACAGCACCAGCCGTTCTACCCCTTCACCCAGCAGGGCATGGGCCACGGCCGAGCCAGCGCCACCGGCACCCATCTGCACCACCTGGCGCCGTGCAACGTCGGGCAGGCCACGGCGCAAGCCTTCGGCAAAGCCCAGGCAATc       <  2:186360/145‑1 (MQ=255)
                                                                                                                                               cGAACAGCACCAGCCGTTCTACCCCTTCACCCAGCAGGGCATGGGCCACGGCCGAGCCAGCGCCACCGGCACCCATCTGCACCACCTGGCGCCGTGCAACGTCGGGCAGGCCACGGCGCAAGCCTTCGGCAAAGCCCAGGCAATCGgt    <  6:235023/148‑1 (MQ=255)
                                                                                                                                                gAACAGCACCAGCCGTTCTACCCCTTCACCCAGCAGGGCATGGGCCACGGCCGAGCCAGCGCCACCGGCACCCATCTGCACCACCTGGCGCCGTGCAACGTCGCGCAGGCCACGGCGCAAGCCTTCGGCAAAGCCCAGGCAATCGgtgt  <  1:289555/149‑1 (MQ=255)
                                                                                                                                                 |                                                                                                                                                   
GGCCATGCCCACCGGCGTGGTGTTGACCAGCCCGTCCGCTTCGGCCAAGGCTGTGGCCAGGTCGGTGCCAAGCACGGCGCGCTCCGCGCCGAAATGGGTGTTCAGGTTGTCCACCAACGCCTGCGCACGGGTCGCATCCACTTCGAACAGCACCAGCCGTTCTACCCCTTCACCCAGCAGGGCATGGGCCACGGCCGAGCCAGCGCCACCGGCACCCATCTGCACCACCTGGCGCCGTGCAACGTCGGGCAGGCCACGGCGCAAGCCTTCGGCAAAGCCCAGGCAATCGGTGT  >  NC_002947/2749875‑2750167

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