Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,763,850 G→T 6.2% G409V (GGC→GTC)  PP_5057 → M23/M37 family peptidase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,763,8500GT 6.2% 263.8 / 5.8 81G409V (GGC→GTC) PP_5057M23/M37 family peptidase
Reads supporting (aligned to +/- strand):  ref base G (45/31);  new base T (3/2);  total (48/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 = 8.77e-01

GTGGTATTCGCGGACTGGTTGCGCGGTGCCGGGCTTCTGGTCATTCTCGACCATGGCAACGGTTACCTGAGCCTGTACGGCCACAACCAGAGCCTGCTCAAGAATGCCGGTGATATCGTCAAGGCCGGTGAAGCCATTTCCACCGTTGGCGACAGCGGTGGTCAGGACGCGGCAGGCTTGTATTTCGCTATCCGCCAGCAGGGCCGGCCTACCGACCCTTCGCAGTGGTGCCGGGGCTAGTCAAATTTCTACGGCGCAGC‑CTGAGTGCTGCGCCGATGCTCCCCAGGGAGCGCTTAC  >  NC_002947/5763702‑5763998
                                                                                                                                                    |                                                                                                                                                     
gtgGTATTCGCGGACTGGTTGCGCGGTGCCGGGCTTCTGGTCATTCTCGACCATGGCAACGGTTACCTGAGCCTGTACGGCCACAACCAGAGCCTGCTCAAGAATGCCGGTGATATCGTCAAGGCCGGTGAAGCCATTTCCACCGTTgg                                                                                                                                                       <  8:3993/149‑1 (MQ=255)
 tgGTATTCGCGGACTGGTTGCGCGGTGCCGGGCTTCTGGTCATTCTCGACCATGGCAACGGTTACCTGAGCCTGTACGGCCACAACCAGAGCCTGCTCAAGAATGCCGGTGATATCGTCAAGGCCGGTGAAGCCATTTCCACCGTTGGc                                                                                                                                                      >  6:232075/1‑149 (MQ=255)
 tgGTATTCGCGGACTGGTTGCGCGGTGCCGGGCTTCTGGTCATTCTCGACCATGGCAACGGTTACCTGAGCCTGTACGGCCACAACCAGAGCCTGCTCAAGAATGCCGGTGATATCGTCAAGGCCGGTGAAGCCATTTCCACCGTTGGc                                                                                                                                                      >  2:88357/1‑149 (MQ=255)
 tgGTATTCGCGGACTGGTTGCGCGGTGCCGGGCTTCTGGTCATTCTCGACCATGGCAACGGTTACCTGAGCCTGTACGGCCACAACCAGAGCCTGCTCAAGAATGCCGGTGATATCGTCAAGGCCGGTGAAGCCATTTCCACCGTTGGc                                                                                                                                                      <  3:160062/149‑1 (MQ=255)
 tgGTATTCGCGGACTGGGTGCGCGGTGCCGGGCTTCTGGTCATTCTCGACCATGGCAACGGTTACCTGAGCCTGTACGGCCACAACCAGAGCCTGCTCAAGAATGCCGGTGATATCGTCAAGGCCGGTGAAGCCATTTCCACCGTTGGc                                                                                                                                                      <  4:261750/149‑1 (MQ=255)
      ttCGCGGACTGGTTGCGCGGTGCCGGGCTTCTGGTCATTCTCGACCATGGCAACGGTTACCTGAGCCTGTACGGCCACAACCAGAGCCTGCTCAAGAATGCCGGTGATATCGTCAAGGCCGGTGAAGCCATTTCCACCGTTGGCGACAg                                                                                                                                                 <  5:298275/149‑1 (MQ=255)
       tCGCGGACTGGTTGGGCGGTGCCGGGCTTCTGGTCATTCTCGACCATGGCAACGGTTACCTGAGCCTGTACGGCCACAACCAGAGCCTGCTCAAGAATGCCGGTGATATCGTCAAGGCCGGTGAAGCCATTTCCACCGTTGGCGACAGc                                                                                                                                                >  7:347949/1‑149 (MQ=255)
       tCGCGGACTGGTTGCGCGGTGCCGGGCTTCTGGTCATTCTCGACCATGGCAACGGTTACCTGAGCCTGTACGGCCACAACCAGAGCCTGCTCAAGAATGCCGGTGATATCGTCAAGGCCGGTGAAGCCATTTCCACCGTTGGCGACAGc                                                                                                                                                >  2:261406/1‑149 (MQ=255)
           ggACTGGTTGCGCGGTGCCGGGCTTCTGGTCATTCTCGACCATGGCAACGGTTACCTGAGCCTGTACGGCCACAAGCAGAGCCTGCTCAAGAATGCCGGTGATATCGTCAAGGCCGGTGAAGCCATTTCCACCGTTGGCGACAGCggtg                                                                                                                                            <  5:232075/149‑1 (MQ=255)
           ggACTGGTTGCGCGGTGCCGGGCTTCTGGTCATTCTCGACCATGGCAACGGTTACCTGAGCCTGTACGGCCACAACCAGAGCCTGCTCAAGAATGCCGGTGATATCGTCAAGGCCGGTGAAGCCATTTCCACCGTTGGCGACAGCggtg                                                                                                                                            >  3:228276/1‑149 (MQ=255)
             aCTGGTTGCGCGGTGCCGGGCTTCTGGTCATTCTCGACCATGGCAACGGTTACCTGAGCCTGTACGGCCACAACCAGAGCCTGCTCAAGAATGCCGGTGATATCGTCAAGGCCGGTGAAGCCATTTCCACCGTTGGCGACAGCggtg                                                                                                                                            <  8:78665/147‑1 (MQ=255)
                      gcgGTGCCGGGCTTCTGGTCATTCTCGACCATGGCAACGGTTACCTGAGCCTGTACGGCCACAACCAGAGCCTGCTCAAGAATGCCGGTGATATCGTCAAGGCCGGTGAAGCCATTTCCACCGTTGGCGACAGCGGTGGTCAGGAcgcg                                                                                                                                 <  3:269975/149‑1 (MQ=255)
                              gggCTTCTGGTCATTCTCGACCATGGCAACGGTTACCTGAGCCTGTACGGCCACAACCAGAGCCTGCTCAAGAATGCCGGTGATATCGTCAAGGCCGGTGAAGCCATTTCCACCGTTGGCGACAGCGGTGGTCAGGACGCGGCAGGCtt                                                                                                                         >  6:107994/1‑149 (MQ=255)
                                gCTTCTGGTCATTCTCGACCATGGCAACGGTTACCTGAGCCTGTACGGCCACAACCAGAGCCTGCTCAAGAATGCCGGTGATATCGTCAAGGCCGGTGAAGCCATTTCCACCGTTGGCGACAGCGGTGGTCAGGACGCGGCAGGCTTGt                                                                                                                       <  5:107994/149‑1 (MQ=255)
                                gCTTCTGGTCATTCTCGACCATGGCAACGGTTACCTGAGCCTGTACGGCCACAACCAGAGCCTGCTCAAGAATGCCGGTGATATCGTCAAGGCCGGTGAAGCCATTTCCACCGTTGGCGACAGCGGTGGTCAGGACGCGGCAGGCTTGt                                                                                                                       <  3:13399/149‑1 (MQ=255)
                                gCTTCGGGGGATTCTCGACCATGGCAACGGTTACCTGAGCCTGTACGGCCACAACCAGAGCCTGCTCAAGAATGCCTGTGATATCCTCAAGGCCGGTGAAGCCATTTCCACCGTTGGCGACAGCGGTGGTCAGGACGCGGCAGGCTTGt                                                                                                                       <  2:107654/149‑1 (MQ=255)
                                      ggTCATTCTCGACCATGGCAACGGTTACCTGAGCCTGTACGGCCACAACCAGAGCCTGCTCAAGAATGCCGGTGATATCGTCAAGGCCGGTGAAGCCATTTCCACCGTTGGCGACAGCGGTGGTCAGGACGCGGCAGGCTTGTATTTCg                                                                                                                 >  2:57169/1‑149 (MQ=255)
                                       gTCATTCTCGACCATGGCAACGGTTACCTGAGCCTGTACGGCTACAACCAGAGCCTGCTCAAGAATGCCGGTGATATCGTCAAGGCCGGTGAAGCCATTTCCACCGTTGGCGACAGCGGTGGTCAGGACGCGGCAGGCTTGTATTTCg                                                                                                                 <  7:170234/148‑1 (MQ=255)
                                       gTCATTCTCGACCATGGCAACGGTTACCTGAGCCTGTACGGCCACAACCAGAGCCTGCTCAAGAATGCCGGTGATATCGTCAAGGCCGGTGAAGCCATTTCCACCGTTGGCGACAGCGGTGGTCAGGACGCGGCAGGCTTGTa                                                                                                                      <  4:332881/143‑1 (MQ=255)
                                       gTCATTCTCGACCATGGCAACGGTTACCTGAGCCTGTACGGCCACAACCAGAGCCTGCTCAAGAATGCCGGTGATATCGTCAAGGCCGGTGAAGCCATTTCCACCGTTGGCGACAGCGGTGGTCAGGACGCGGCAGGCTTGTa                                                                                                                      >  3:332881/1‑143 (MQ=255)
                                       gTCATTCTCGACCATGGCAACGGTTACCTGAGCCTGTACGGCCACAACCAGAGCCTGCTCAAGAATGCCGGTGATATCGTCAAGGCCGGTGAAGCCATTTCCACCGTTGGCGACAGCGGTGGTCAGGACGCGGCAGGCTTGTATTTCg                                                                                                                 <  8:347949/148‑1 (MQ=255)
                                       gTCATTCTCGACCATGGCAACGGTTACCTGAGCCTGTACGGCCACAACCAGAGCCTGCTCAAGAATGCCGGTGATATCGTCAAGGCCGGTGAAGCCATTTCCACCGTTGGCGACAGCGGTGGTCAGGACGCGGCAGGCTTGTATTTCGc                                                                                                                >  5:180279/1‑149 (MQ=255)
                                       gTCATTCTCGACCATGGCAACGGTTACCTGAGCCTGTACGGCCACAACCAGAGCCTGCTCAAGAATGCCGGTGATATCGTCAAGGCCGGTGAAGCCATTTCCACCGTTGGCGACAGCGGTGGTCAGGACGCGGCAGGCTTGTATTTCGc                                                                                                                >  7:168021/1‑149 (MQ=255)
                                       gTCATTCTCGACCATGGCAACGGTTACCTGAGCCTGTACGGCCACAACCAGAGCCTGCTCAAGAATGCCGGTGATATCGTCAAGGCCGGTGAAGCCATTTCCACCGTTGGCGACAGCGGTGGTCAGGACGCGGCAGGCTTGTATTTCGc                                                                                                                >  8:73634/1‑149 (MQ=255)
                                         cATTCTCGACCATGGCAACGGTTACCTGAGCCTGTACGGCCACAACCAGAGCCTGCTCAAGAATGCCGGTGATATCGTCAAGGCCGGTGAAGCCATTTCCACCGCTGTCGCCAACggtg                                                                                                                                            <  3:261626/119‑1 (MQ=255)
                                         cATTCTCGACCATGGCAACGGTTACCTGAGCCTGTACGGCCACAACCAGAGCCTGCTCAAGAATGCCGGTGATATCGTCAAGGCCGGTGAAGCCATTTCCACCGCTGTCGCCAACggtg                                                                                                                                            >  4:261626/1‑119 (MQ=255)
                                          aTTCTCGACCATGGCAACGGTTACCTGAGCCTGTACGGCCACAACCAGAGCCTGCTCAAGAATGCCGGTGATATCGTCAAGGCCGGTGAAGCCATTTCCACCGTTGGCGACAGCGGTGGTCAGGACGCGGCAGGCTTGTATTTCGCt                                                                                                               >  8:346996/1‑147 (MQ=255)
                                          aTTCTCGACCATGGCAACGGTTACCTGAGCCTGTACGGCCACAACCAGAGCCTGCTCAAGAATGCCGGTGATATCGTCAAGGCCGGTGAAGCCATTTCCACCGTTGGCGACAGCGGTGGTCAGGACGCGGCAGGCTTGTATTTCGCt                                                                                                               <  7:346996/147‑1 (MQ=255)
                                            tctcGACCATGGCAACGGTTACCTGAGCCTGTACGGCCACAACCAGAGCCTGCTCAAGAATGCCGGTGATATCGTCAAGGCCGGTGAAGCCATTTCCACCGTTGGCGACAGCGGTGGTCAGGACGCGGCAGGCTTGTATTTCGCTATcc                                                                                                           >  2:192574/1‑149 (MQ=255)
                                                gACCATGGCAACGGTTACCTGAGCCTGTACGGCCACAACCAGAGCCTGCTCAAGAATGCCGGTGATATCGTCAAGGCCGGTGAAGCCATTTCCACCGTTGGCGACAGCGGTGGTCAGGACGCGGCAGGCt                                                                                                                          >  1:289541/1‑130 (MQ=255)
                                                gACCATGGCAACGGTTACCTGAGCCTGTACGGCCACAACCAGAGCCTGCTCAAGAATGCCGGTGATATCGTCAAGGCCGGTGAAGCCATTTCCACCGTTGGCGACAGCGGTGGTCAGGACGCGGCAGGCt                                                                                                                          <  2:289541/130‑1 (MQ=255)
                                                    aTGGCAACGGTTACCTGAGCCTGTACGGCCACAACCAGAGCCTGCTCAAGAATGCCGGTGATATCGTCAAGGCCGGTGAAGCCATTTCCACCGTTGGCGACAGCGGTGGTCAGGACGCGGCAGGCTTGTATTTCGCTATCCGCcagcag                                                                                                   >  5:201216/1‑149 (MQ=255)
                                                      ggCAACGGTTACCTGAGCCTGTACGGCCACAACCAGAGCCTGCTCAAGAATGCCGGTGATATCGTCAAGGCCGGTGAAGCCATTTCCACCGTTGGCGACAGCGGTGGTCAGGACGCGGCAGGCTTGTATTTCGCTATCCGCCAGCAggg                                                                                                 >  4:61097/1‑149 (MQ=255)
                                                       gCAACGGTTACCTGAGCCTGTACGGCCACAACCAGAGCCTGCTCAAGAATGCCGGTGATATCGTCAAGGCCGGTGAAGCCATTTCCACCGTTGGCGACAGCGGTGGTCAGGACGCGGCAGGCTTGTATTTCGCTATCCGCCAGCAggg                                                                                                 <  6:201216/148‑1 (MQ=255)
                                                        cAACGGTTACCTGAGCCTGTACGGCCACAACCAGAGCCTGCTCAAGAATGCCGGTGATATCGTCAAGGCCGGTGAAGCCATTTCCACCGTTGGCGACAGCGGTGGTCAGGACGCGGCAGGCTTGTATTTCGCTATCCGCCAGCAGggcc                                                                                               <  7:73634/149‑1 (MQ=255)
                                                            ggTTACCTGAGCCTGTACGGCCACAACCAGAGCCTGCTCAAGAATGCCGGTGATATCGTCAAGGCCGGTGAAGCCATTTTCACCGTTGGCGACAGCGGTGGTCa                                                                                                                                        <  6:2951/104‑1 (MQ=255)
                                                            ggTTACCTGAGCCTGTACGGCCACAACCAGAGCCTGCTCAAGAATGCCGGTGATATCGTCAAGGCCGGTGAAGCCATTTACACCGTTGGCGACAGCGGTGGTCa                                                                                                                                        >  5:2951/1‑104 (MQ=255)
                                                                    gAGCCTGTACGGCCACAACCAGAGCCTGCTCAAGAATGCCGGTGATATCGTCAAGGCCGGTGAAGCCATTTCCACCGTTGGCGACAGCGGTGGTCAGGACGCGGCAGGCTTGTATTTCGCTATCCGCCAGCAGGGCCGGCCTACCGAcc                                                                                   >  4:259677/1‑149 (MQ=255)
                                                                    gAGCCTGTACGGCCACAACCAGAGCCTGCTCAAGAATGCCGGTGATATCGTCAAGGCCGGTGAAGCCATTTCCACCGTTGGCGACAGCGGTGGTCAGGACGCGGCAGGCTTGTATTTCGCTATCCGCCAGCAGGGCCGGCCTACCGAcc                                                                                   >  3:48063/1‑149 (MQ=255)
                                                                          gTACGGCCACAACCAGAGCCTGCTCAAGAATGCCGGTGATATCGTCAAGGCCGGTGAAGCCATTTCCACCGTTGGCGACAGCGGTGGTCAGGACGCGGCAGGCTTGTATTTCGCTATCCGCCAGCAGGGCCGGCCTACCGACCCTTCg                                                                              >  5:77723/1‑148 (MQ=255)
                                                                           tACGGCCACAACCAGAGCCTGCTCAAGAATGCCGGTGATATCGTCAAGGCCGGTGAAGCCATTTCCACCGTTGGCGACAGCGGTGGTCAGGACGCGGCAGGCTTGTATTTCGCTATCCGCCAGCAGGGCCGGCCTACCGACCCTTCGc                                                                             <  1:88357/148‑1 (MQ=255)
                                                                            aCGGCCACAACCAGAGCCTGCTCAAGAATGCCGGTGATATCGTCAAGGCCGGTGAAGCCATTTCCACCGTTGGCGACAGCGGTGGTCAGGACGCGGCAGGCtt                                                                                                                         >  1:202048/1‑103 (MQ=255)
                                                                            aCGGCCACAACCAGAGCCTGCTCAAGAATGCCGGTGATATCGTCAAGGCCGGTGAAGCCATTTCCACCGTTGGCGACAGCGGTGGTCAGGACGCGGCAGGCtt                                                                                                                         <  2:202048/103‑1 (MQ=255)
                                                                              ggCCACAACCAGAGCCTGCTCAAGAATGCCGGTGATATCGTCAAGGCCGGTGAAGCCATTTCCACCGTTGGCGACAGCGGTGGTCAGGACGCGGCAGGCTTGTATTTCGCTATCCGCCAGCAGGGCCGGCCTACCGACCCTTCGCAgtg                                                                         >  4:176633/1‑149 (MQ=255)
                                                                              ggCCACAACCAGAGCCTGCTCAAGAATGCCGGTGATATCGTCAAGGCCGGTGAAGCCATTTCCACCGTTGGCGACAGCGGTGGTCAGGACGCGGCAGGCTTGTATTTCGCTATCCGCCAGCAGGGCCGGCCTACCGACCCTTCGCAgtg                                                                         <  8:168021/149‑1 (MQ=255)
                                                                                caacaACCAGAGCCAGCTCAAGAATGCCGGTGATATCGTCAACGCCAGTTAAGCCATTTACACAGTTGGCGACAGCGCTCCTCAGGACGCGTCAGGCTAGTATTTCGCTATCCGCCACCAGCGCCGTCCAACCGCCCCTTCGCAgtgct                                                                       >  8:294120/3‑147 (MQ=255)
                                                                                 cacaACCAGAGCCTGCTCAAGAATGCCGGTGATATCGTCAAGGCCGGTGAAGCCATTTCCACCGTTGGCGACAGCGGTGGTCAGGACGCGGCAGGCTTGTATTTCGCTATCCGCCAGAAGGGCCGGCCTACCGCCCCTACGCActggtg                                                                      >  3:75792/1‑149 (MQ=255)
                                                                                  acaACCAGAGCCTGCTCAAGAATGCCGGTGATATCGTCAAGGCCGGTGAAGCCATTTCCACCGTTGGCGACAGCGGTGGTCAGGACGCGGCAGGCTTGTATTTCGCTATCCGCCAGCAGGGCCGGCCTACCGAACCTTCGCAGTGGTGc                                                                     >  4:212994/1‑149 (MQ=255)
                                                                                  acaACCAGAGCCTGCTCAAGAATGCCGGTGATATCGTCAAGGCCGGTGAAGCCATTTCCACCGTTGGCGACAGCGGTGGTCAGGACGCGGCAGGCTTGTATTTCGCTATCAGCCAGCAGGGCCGGCCTACCGACCCTTCGCAgtggttc                                                                     >  5:247896/1‑147 (MQ=255)
                                                                                   caACCAGAGCCTGCTCAAGAATGCCGGTGATATCGTCAAGGCCGGTGAAGCCATTTCCACCGTTGGCGACAGCGGTGGTCAGGACGCGGCAGGCTTGTATTTCg                                                                                                                 <  7:142518/104‑1 (MQ=255)
                                                                                   caACCAGAGCCTGCTCAAGAATGCCGGTGATATCGTCAAGGCCGGTGAAGCCATTTCCACCGTTGGCGACAGCGGTGGTCAGGACGCGGCAGGCTTGTATTTCg                                                                                                                 >  8:142518/1‑104 (MQ=255)
                                                                                    aaCCAGAGCCTGCTCAAGAATGCCGGTGATATCGTCAAGGCCGGTGAAGCCATTTCCACCGTTGGCGACAGCGGTGGTCAGGACGCGGCAGGCTTGTATTTCGCTATCCGCCAGCAGGGCCGGCCTACCGACCCTTCGCAGTGGTGCCg                                                                   >  7:247635/1‑149 (MQ=255)
                                                                                       cAGAGCCTGCTCAAGAATGCCGGTGATATCGTCAAGGCCGGTGAAGCCATTTCCACCGTTGGCGACAGCGGTGGTCAGGACGCGGCAGGCTTGTATTTCGCTATCCGCCAGCAGGGCCGGCCTAc                                                                                        >  4:271849/1‑125 (MQ=255)
                                                                                       cAGAGCCTGCTCAAGAATGCCGGTGATATCGTCAAGGCCGGTGAAGCCATTTCCACCGTTGGCGACAGCGGTGGTCAGGACGCGGCAGGCTTGTATTTCGCTATCCGCCAGCAGGGCCGGCCTAc                                                                                        <  3:271849/125‑1 (MQ=255)
                                                                                          agCCTGCTCAAGAATGCCGGTGATATCGTCAAGGCCGGTGAAGCCATTTCCACCGTTGGCGACAGCGGTGGTCAGGACGCGGCAGGCTTGTATTTCGCTATCCGCCAGCAGGGCCGGCCTACCGACCCTTCGCAGTGGTGCCGGGGCTa                                                             <  4:48063/149‑1 (MQ=255)
                                                                                            ccTGCTCAAGAATGCCGGTGATATCGTCAAGGCCGGTGAAGCCATTTCCACCGTTGGCGACAGCGGTGGTCAGGACGCGGCAGGCTTGTATTTCGCTATCCGCCAGCAGGGCCGGCCTACCGACCCTTCGCAGTGGTGCCGGGGCTAGt                                                           <  8:63147/149‑1 (MQ=255)
                                                                                              tGCTCAAGAATGCCGGTGATATCGTCAAGGCCGGTGAACCCATTTCCACCGTTGGCGACAGCGGTGGTCAGGACGCGGCAGGCTTGTATTTCGCTATCCGCCAGCAGGGCCGGCCTACCGACCCTTCGCAGTGGTGCCGGGGCTAGTc                                                          >  4:98054/1‑148 (MQ=255)
                                                                                                 tCAAGAATGCCGGTGATATCGTCAAGGCCGGTGAAGCCATTTCCACCGTTGGCGACAGCGGTGGTCAGGACGCGGCAGGCTTGTATTTCGCTATCCGCCAGCAGGGCCGGCCTACCGCCCCTTCGAAGTGGTGCCGGGGCTAGTCAAAt                                                      >  4:370918/1‑149 (MQ=255)
                                                                                                 tCAAGAATGCCGGTGATATCGTCAAGGCCGGTGAAGCCATTTCCACCGTTGGCGACAGCGGTGGTCAGGACGCGGCAGGCTTGTATTTCGCTATCCGCCAGCAGGGCCGGCCTACCGACCCTTCGCAGTGGTGCCGGGGCTAGTCAAAt                                                      >  2:294310/1‑149 (MQ=255)
                                                                                                 tCAAGAATGCCGGTGATATCGTCAAGGCCGGTGAAGCCATTTCCACCGTTGGCGACAGCGGTGGTCAGGACGCGGCAGGCTTGTATTTCGCTATCCGCCAGCAGGGCCGGCCTACCGACCCTTCGCAGTGGTGCCGGGGCTAGTCAAAt                                                      <  6:63647/149‑1 (MQ=255)
                                                                                                      aaTGCCGGTGATATCGTCAAGGCCGGTGAAGCCATTTCCACCGTTGGCGACAGCGGTGGTCAGGACGCGGCAGGCTTGTATTTCGCTATCCGCCAGCAGGGCCGGCCTACCGACCCTTCGCAGTGGTGCCGGGGCTAGTCAAATTTCTa                                                 >  3:376200/1‑149 (MQ=255)
                                                                                                          ccGGTGATATCGTCAAGGCCGGTGAAGCCATTTCCACCGTTGGCGACAGCGGTGGTCAGGACGCGGCAGGCTTGTATTTCGCTATCCGCCAGCAGGGCCGGCCTACCGACCCTTCGCAGTGGTGCCGGGGCTAGTCAAATTTCTACGgc                                             >  4:98495/1‑149 (MQ=255)
                                                                                                          ccGGTGATATCGTCAAGGCCGGTGAAGCCATTTCCACCGTTGGCGACAGCGGTGGTCAGGACGCGGCAGGCTTGTATTTCGCTATCCGCCAGCAGGGCCGGCCTACCGACCCTTCGCAGTGGTGCCGGGGCTAGTCAAATTTCTACGgc                                             <  6:77723/149‑1 (MQ=255)
                                                                                                            ggTGATATCGTCAAGGCCGGTGAAGCCATTTCCACCGTTGGCGACAGCGGTGGTCAGGACGCGGCAGGCTTGTATTCCGCTATCCGCCAGCAGGGCCGGCCTACCGACCCTTCGCAGTGGTGCCGGGGCTAGTCAAATTTCTACGgcg                                            >  5:26590/1‑148 (MQ=255)
                                                                                                                 tatCGTCAAGGCCGGTGAAGCCATTTCCACCGTTGGCGACAGCGGTGGTCAGGACGCGGCAGGCTTGTATTTCGCTATCCGCCAGCAGGGCCGGCCTACCGACCCTTCGCAGTGGTGCCGGGGCTAGTCAAATTTCTACGGCGCAGC‑Ct                                     <  8:247635/149‑1 (MQ=255)
                                                                                                                 tatCGTCAAGGCCGGTGAAGCCATTTCCACCCTTGGCGACAGCGGAGGTCAGGACGCGGCAGGCTTCTATTTCGCTATCCGCCAGCAGGGCCGGCCTACCGACCCTTCGCAGTGGTGCCGGGGCTAGTCAAATTTCTACGGCGCAGC‑Ct                                     >  1:210810/1‑149 (MQ=255)
                                                                                                                     gTCAAGGCCGGTGAAGCCATTTCCACCGTTGGCGACAGCGGTGGTCAGGACGCGGCAGGCTTGTATTTCGCTATCCGCCAGCAGGGCCGGCCTACCGACCCTTCGCAGTGGTGCCGGGGCTAGTCAAATTTCTACGGCGCAGC‑CTGAGt                                 <  3:61097/149‑1 (MQ=255)
                                                                                                                     gTCAAGGCCGGTGAAGCCATTTCCACCGTTGGCGACAGCGGTGGTCAGGACGCGGCAGGCTTGTATTTCGCTATCCGCCAGCAGGGCCGGCCTACCGACCCTTCGCAGTGGTGCCGGGGCTAGTCAAATTTCTACGGCGCAGC‑CTGAGt                                 >  2:323646/1‑149 (MQ=255)
                                                                                                                         aGGCCGGTGAAGCCATTTCCACCGTTGGCGACAGCGGTGGTCAGGACGCGGCAGGCTTGTATTTCGCTATCCGCCAGCAGGGCCGGCCTACCGACCCTTCGCAGTGGTGCCGGGGCTAGTCAAATTTCTACGGCGCAGC‑CTGAGtgctg                             >  1:34958/1‑149 (MQ=255)
                                                                                                                          ggCCGGTGAAGCCATTTCCACCGTTGGCGACAGCGGTGGTCAGGACGCGGCAGGCTTGTATTTCGCTATCCGCCAGCAGGGCCGGCCTACCGACCCTTCGCAGTGGTGCCGGGGCTAGTCAAATTTCTACGGCGCAGC‑CTGAGtgctgc                            >  1:394393/1‑149 (MQ=255)
                                                                                                                           gCCGGTGAAGCCATTTCCACCGTTGGCGACAGCGGTGGTCAGGACGCGGCAGGCTTGTATTTCGCTATCCGCCAGCAGGGCCGGCCTACCGACCCTTCGCAGTGGTGCCGGGGCTAGTCAAATTTCTACGGCGCAGC‑CTGAGTGCTgcg                           >  7:383678/1‑149 (MQ=255)
                                                                                                                             cGGTGAAGCCATTTCCACCGTTGGCGACAGCGGTGGTCAGGACGCGGCAGGCTTGTATTTCGCTATCCGCCAGCAGGGCCGGCCTACCGACCCTTCGCAGTGGTGCCGGGGCTAGTCAAATTTCTACGGCGCAGC‑CTGAGTGCTGCGcc                         <  2:394393/149‑1 (MQ=255)
                                                                                                                              ggTGAAGCCATTTCCACCGTTGGCGACAGCGGTGGTCAGGACGCGGCAGGCTTGTATTTCGCTATCCGCCAGCAGGGCCGGCCTACCGACCCTTCGCAGTGGTGCCGGGGCTAGTCAAATTTCTACGGCGCAGCACTCAG‑GCTGCGCCg                        <  4:376200/149‑1 (MQ=255)
                                                                                                                                      cATTTCCACCGTTGGCGACAGCGGTGGTCAGGACGCGGCAGGCTTGTATTTCGCTATCCGCCAGCAGGGCCGGCCTACCGACCCTTCGCAGTGGTGCCGGGGCTAGTCAAATTTCTACGGCGCAGC‑CTGAGTGCTGCGCCGATGCTccc                >  2:223283/1‑149 (MQ=255)
                                                                                                                                       aTTTCCACCGTTGGCGACAGCGGTGGTCAGGACGCGGCAGGCTTGTATTTCGCTATCCGCCAGCAGGGCCGGCCTACCGACCCTTCGCAGTGGTGCCGGGGCTAGTCAAATTTCTACGGCGCAGC‑CTGAGTGCTGCGCCGATGCTccc                >  2:48437/1‑148 (MQ=255)
                                                                                                                                       aTTTCCACCGTTGGCGACAGCGGTGGTCAGGACGCGGCAGGCTTGTATTTCGCTATCCGCCAGCAGGGCCGGCCTACCGACCCTTCGCAGTGGTGCCGGGGCTAGTCAAATTTCTACGGCGCAGC‑CTGAGTGCTGCGCCGATGCTccc                >  8:166724/1‑148 (MQ=255)
                                                                                                                                       aTTTCCACCGTTGGCGACAGCGGTGGTCAGGACGCGGCAGGCTTGTATTTCGCTATCCGCCAGCAGGGCCGGCCTACCGACCCTTCGCAGTGGTGCCGGGGCTAGTCAAATTTCTACGGCGCAGC‑CTGAGTGCTGCGCCGATGCTccc                >  5:218934/1‑148 (MQ=255)
                                                                                                                                        tttCCACCGTTGGCGACAGCGGTGGTCAGGACGCGGCAGGCTTGTATTTCGCTATCCGCCAGCAGGGCCGGCCTACCGACCCTTCGCAGTGGTGCCGGGGCTAGTCAAATTTCTACGGCGCAGC‑CTGAGTGCTGCGCCGATGCTCCCCa              <  3:176633/149‑1 (MQ=255)
                                                                                                                                         gaccACCGCTGTCGCCAACGGTGGTCAGGACGCGGCAGGCTTGTATTTCGCTATCCGCCAGCAGGGCCGGCCTACCGACCCTTCGCAGTGGTGCCGGGGCTAGTCAAATTTCTACGGCGCAGC‑CTGAGTGCTGCGCCGATGCTCCCCAg             >  2:112584/3‑149 (MQ=255)
                                                                                                                                          accACCGCTGTCGCCAACGGTGGTCAGGACGCGGCAGGCTTGTATTTCGCTATCCGCCAGCAGGGCCGGCCTACCGACCCTTCGCAGTGGTGCCGGGGCTAGTCAAATTTCTACGGCGCAGc‑c                                      >  6:212549/2‑123 (MQ=255)
                                                                                                                                          accACCGCTGTCGCCAACGGTGGTCAGGACGCGGCAGGCTTGTATTTCGCTATCCGCCAGCAGGGCCGGCCTACCGACCCTTCGCAGTGGTGCCGGGGCTAGTCAAATTTCTACGGCGCAGc‑c                                      <  5:212549/122‑1 (MQ=255)
                                                                                                                                           ccACCGTTGGCGACAGCGGTGGTCAGGACGCGGCAGGCTTGTATTTCGCTATCCGCCAGCAGGGCCGGCCTACCGACCCTTCGCAGTGGTGCC‑GGGCTAGTCAAATTTCTACGGCGCAGC‑CTGAGTGCTGCGCCGATGCTCCCCAgg            >  5:361281/1‑147 (MQ=255)
                                                                                                                                                    gCGACAGCGGTGGTCAGGACGCGGCAGGCTTGTATTTCGCTATCCGCCAGCAGGGCCGGCCTACCGACCCTTCGCAGTGGTGCCGGGGCTAGTCAAATTTCTACGGCGCAGC‑CTGAGTGCTGCGCCGATGCTCCCCAGGGAGCGCTTAc  <  1:57169/149‑1 (MQ=255)
                                                                                                                                                    gCGACAGCGGTGGTCAGGACGCGGCAGGCTTGTATTTCGCTATCCGCCAGCAGGGCCGGCCTACCGACCCTTCGCAGTGGTGCCGGGGCTAGTCAAATTTCTACGGCGCAGC‑CTGAGTGCTGCGCCGATGCTCCCCAGGGAGCGCTTAc  <  3:259677/149‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                     
GTGGTATTCGCGGACTGGTTGCGCGGTGCCGGGCTTCTGGTCATTCTCGACCATGGCAACGGTTACCTGAGCCTGTACGGCCACAACCAGAGCCTGCTCAAGAATGCCGGTGATATCGTCAAGGCCGGTGAAGCCATTTCCACCGTTGGCGACAGCGGTGGTCAGGACGCGGCAGGCTTGTATTTCGCTATCCGCCAGCAGGGCCGGCCTACCGACCCTTCGCAGTGGTGCCGGGGCTAGTCAAATTTCTACGGCGCAGC‑CTGAGTGCTGCGCCGATGCTCCCCAGGGAGCGCTTAC  >  NC_002947/5763702‑5763998

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: ‑