Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,763,856 G→A 6.0% S411N (AGC→AAC)  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,8560GA 6.0% 253.4 / 3.7 83S411N (AGC→AAC) PP_5057M23/M37 family peptidase
Reads supporting (aligned to +/- strand):  ref base G (45/33);  new base A (3/2);  total (48/35)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.24e-01

TTCGCGGACTGGTTGCGCGGTGCCGGGCTTCTGGTCATTCTCGACCATGGCAACGGTTACCTGAGCCTGTACGGCCACAACCAGAGCCTGCTCAAGAATGCCGGTGATATCGTCAAGGCCGGTGAAGCCATTTCCACCGTTGGCGACAGCGGTGGTCAGGACGCGGCAGGCTTGTATTTCGCTATCCGCCAGCAGGGCCGGCCTACCGACCCTTCGCAGTGGTGCCGGGGCTAGTCAAATTTCTACGGCGCAGC‑CTGAGTGCTGCGCCGATGCTCCCCAGGGAGCGCTTACAGGA  >  NC_002947/5763708‑5764002
                                                                                                                                                    |                                                                                                                                                   
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                                                                                                                          >  3:332881/1‑143 (MQ=255)
                                 gTCATTCTCGACCATGGCAACGGTTACCTGAGCCTGTACGGCCACAACCAGAGCCTGCTCAAGAATGCCGGTGATATCGTCAAGGCCGGTGAAGCCATTTCCACCGTTGGCGACAGCGGTGGTCAGGACGCGGCAGGCTTGTa                                                                                                                          <  4:332881/143‑1 (MQ=255)
                                 gTCATTCTCGACCATGGCAACGGTTACCTGAGCCTGTACGGCCACAACCAGAGCCTGCTCAAGAATGCCGGTGATATCGTCAAGGCCGGTGAAGCCATTTCCACCGTTGGCGACAGCGGTGGTCAGGACGCGGCAGGCTTGTATTTCg                                                                                                                     <  8:347949/148‑1 (MQ=255)
                                 gTCATTCTCGACCATGGCAACGGTTACCTGAGCCTGTACGGCCACAACCAGAGCCTGCTCAAGAATGCCGGTGATATCGTCAAGGCCGGTGAAGCCATTTCCACCGTTGGCGACAGCGGTGGTCAGGACGCGGCAGGCTTGTATTTCGc                                                                                                                    >  7:168021/1‑149 (MQ=255)
                                 gTCATTCTCGACCATGGCAACGGTTACCTGAGCCTGTACGGCCACAACCAGAGCCTGCTCAAGAATGCCGGTGATATCGTCAAGGCCGGTGAAGCCATTTCCACCGTTGGCGACAGCGGTGGTCAGGACGCGGCAGGCTTGTATTTCGc                                                                                                                    >  5:180279/1‑149 (MQ=255)
                                 gTCATTCTCGACCATGGCAACGGTTACCTGAGCCTGTACGGCCACAACCAGAGCCTGCTCAAGAATGCCGGTGATATCGTCAAGGCCGGTGAAGCCATTTCCACCGTTGGCGACAGCGGTGGTCAGGACGCGGCAGGCTTGTATTTCGc                                                                                                                    >  8:73634/1‑149 (MQ=255)
                                   cATTCTCGACCATGGCAACGGTTACCTGAGCCTGTACGGCCACAACCAGAGCCTGCTCAAGAATGCCGGTGATATCGTCAAGGCCGGTGAAGCCATTTCCACCGCTGTCGCCAACggtg                                                                                                                                                >  4:261626/1‑119 (MQ=255)
                                   cATTCTCGACCATGGCAACGGTTACCTGAGCCTGTACGGCCACAACCAGAGCCTGCTCAAGAATGCCGGTGATATCGTCAAGGCCGGTGAAGCCATTTCCACCGCTGTCGCCAACggtg                                                                                                                                                <  3:261626/119‑1 (MQ=255)
                                    aTTCTCGACCATGGCAACGGTTACCTGAGCCTGTACGGCCACAACCAGAGCCTGCTCAAGAATGCCGGTGATATCGTCAAGGCCGGTGAAGCCATTTCCACCGTTGGCGACAGCGGTGGTCAGGACGCGGCAGGCTTGTATTTCGCt                                                                                                                   <  7:346996/147‑1 (MQ=255)
                                    aTTCTCGACCATGGCAACGGTTACCTGAGCCTGTACGGCCACAACCAGAGCCTGCTCAAGAATGCCGGTGATATCGTCAAGGCCGGTGAAGCCATTTCCACCGTTGGCGACAGCGGTGGTCAGGACGCGGCAGGCTTGTATTTCGCt                                                                                                                   >  8:346996/1‑147 (MQ=255)
                                      tctcGACCATGGCAACGGTTACCTGAGCCTGTACGGCCACAACCAGAGCCTGCTCAAGAATGCCGGTGATATCGTCAAGGCCGGTGAAGCCATTTCCACCGTTGGCGACAGCGGTGGTCAGGACGCGGCAGGCTTGTATTTCGCTATcc                                                                                                               >  2:192574/1‑149 (MQ=255)
                                          gACCATGGCAACGGTTACCTGAGCCTGTACGGCCACAACCAGAGCCTGCTCAAGAATGCCGGTGATATCGTCAAGGCCGGTGAAGCCATTTCCACCGTTGGCGACAGCGGTGGTCAGGACGCGGCAGGCt                                                                                                                              <  2:289541/130‑1 (MQ=255)
                                          gACCATGGCAACGGTTACCTGAGCCTGTACGGCCACAACCAGAGCCTGCTCAAGAATGCCGGTGATATCGTCAAGGCCGGTGAAGCCATTTCCACCGTTGGCGACAGCGGTGGTCAGGACGCGGCAGGCt                                                                                                                              >  1:289541/1‑130 (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                                                                                       >  3:48063/1‑149 (MQ=255)
                                                              gAGCCTGTACGGCCACAACCAGAGCCTGCTCAAGAATGCCGGTGATATCGTCAAGGCCGGTGAAGCCATTTCCACCGTTGGCGACAGCGGTGGTCAGGACGCGGCAGGCTTGTATTTCGCTATCCGCCAGCAGGGCCGGCCTACCGAcc                                                                                       >  4:259677/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                                                                             <  8:168021/149‑1 (MQ=255)
                                                                        ggCCACAACCAGAGCCTGCTCAAGAATGCCGGTGATATCGTCAAGGCCGGTGAAGCCATTTCCACCGTTGGCGACAGCGGTGGTCAGGACGCGGCAGGCTTGTATTTCGCTATCCGCCAGCAGGGCCGGCCTACCGACCCTTCGCAgtg                                                                             >  4:176633/1‑149 (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                                                                                                                     >  8:142518/1‑104 (MQ=255)
                                                                             caACCAGAGCCTGCTCAAGAATGCCGGTGATATCGTCAAGGCCGGTGAAGCCATTTCCACCGTTGGCGACAGCGGTGGTCAGGACGCGGCAGGCTTGTATTTCg                                                                                                                     <  7:142518/104‑1 (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                                                          <  6:63647/149‑1 (MQ=255)
                                                                                           tCAAGAATGCCGGTGATATCGTCAAGGCCGGTGAAGCCATTTCCACCGTTGGCGACAGCGGTGGTCAGGACGCGGCAGGCTTGTATTTCGCTATCCGCCAGCAGGGCCGGCCTACCGACCCTTCGCAGTGGTGCCGGGGCTAGTCAAAt                                                          >  2:294310/1‑149 (MQ=255)
                                                                                                aaTGCCGGTGATATCGTCAAGGCCGGTGAAGCCATTTCCACCGTTGGCGACAGCGGTGGTCAGGACGCGGCAGGCTTGTATTTCGCTATCCGCCAGCAGGGCCGGCCTACCGACCCTTCGCAGTGGTGCCGGGGCTAGTCAAATTTCTa                                                     >  3:376200/1‑149 (MQ=255)
                                                                                                    ccGGTGATATCGTCAAGGCCGGTGAAGCCATTTCCACCGTTGGCGACAGCGGTGGTCAGGACGCGGCAGGCTTGTATTTCGCTATCCGCCAGCAGGGCCGGCCTACCGACCCTTCGCAGTGGTGCCGGGGCTAGTCAAATTTCTACGgc                                                 <  6:77723/149‑1 (MQ=255)
                                                                                                    ccGGTGATATCGTCAAGGCCGGTGAAGCCATTTCCACCGTTGGCGACAGCGGTGGTCAGGACGCGGCAGGCTTGTATTTCGCTATCCGCCAGCAGGGCCGGCCTACCGACCCTTCGCAGTGGTGCCGGGGCTAGTCAAATTTCTACGgc                                                 >  4:98495/1‑149 (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                    >  5:218934/1‑148 (MQ=255)
                                                                                                                                 aTTTCCACCGTTGGCGACAGCGGTGGTCAGGACGCGGCAGGCTTGTATTTCGCTATCCGCCAGCAGGGCCGGCCTACCGACCCTTCGCAGTGGTGCCGGGGCTAGTCAAATTTCTACGGCGCAGC‑CTGAGTGCTGCGCCGATGCTccc                    >  8:166724/1‑148 (MQ=255)
                                                                                                                                 aTTTCCACCGTTGGCGACAGCGGTGGTCAGGACGCGGCAGGCTTGTATTTCGCTATCCGCCAGCAGGGCCGGCCTACCGACCCTTCGCAGTGGTGCCGGGGCTAGTCAAATTTCTACGGCGCAGC‑CTGAGTGCTGCGCCGATGCTccc                    >  2:48437/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)
                                                                                                                                               cGACAGCGGTGGTCAGGACGCGGCAGGCTTGTATTTCGCTATCCGCCAGCAGGGCCGGCCTACCGACCCTTCGCAGTGGTGCCGGGGCTAGTCAAATTTCTACGGCGCAGC‑CTGAGTGCTGCGCCGATGCTCCCCAg                 <  5:19546/137‑1 (MQ=255)
                                                                                                                                               cGACAGCGGTGGTCAGGACGCGGCAGGCTTGTATTTCGCTATCCGCCAGCAGGGCCGGCCTACCGACCCTTCGCAGTGGTGCCGGGGCTAGTCAAATTTCTACGGCGCAGC‑CTGAGTGCTGCGCCGATGCTCCCCAg                 >  6:19546/1‑137 (MQ=255)
                                                                                                                                               cGACAGCGGTGGTCAGGACGCGGCAGGCTTGTACTTCGCTATCCGCCAGCAGGGCCGGCCTACCGACCCTTCGCAGTGGTGCCGGGGCTAGTCAAATTTCTACGGCGCAGC‑CTGAGTGCTGCGCCGATGCTCCCCAGGGAGCGCTTACa     <  3:247517/149‑1 (MQ=255)
                                                                                                                                                 aCAGCTGTGGTCTGGACGCGGCATGGTTGTATTTGGCTATCCGCCAGGAGGGCCGGCCTACCGACCCTTCGCAGTTGTTCCGGGGCTAGTCAAATTTCTACGGCGCAGC‑CTGAGTGGTGGGCCTATGGTCCCCAGGGAGCGCTTACAgg   <  7:294120/149‑1 (MQ=255)
                                                                                                                                                  cAGCGGTGGTCAGGACGCGGCAGGCTTGTATTTCGCTATCCGCCAGCAGGGCCGGCCTACCGCCCCTTCGCAGTGGTGCCGGGGCTAGTCAAATTTCTACGGCGCAGC‑CTGAGTGCTGCGCCGATGCTCCCCAGGGAGCGCTTACAGGa  >  2:278955/1‑149 (MQ=255)
                                                                                                                                                    |                                                                                                                                                   
TTCGCGGACTGGTTGCGCGGTGCCGGGCTTCTGGTCATTCTCGACCATGGCAACGGTTACCTGAGCCTGTACGGCCACAACCAGAGCCTGCTCAAGAATGCCGGTGATATCGTCAAGGCCGGTGAAGCCATTTCCACCGTTGGCGACAGCGGTGGTCAGGACGCGGCAGGCTTGTATTTCGCTATCCGCCAGCAGGGCCGGCCTACCGACCCTTCGCAGTGGTGCCGGGGCTAGTCAAATTTCTACGGCGCAGC‑CTGAGTGCTGCGCCGATGCTCCCCAGGGAGCGCTTACAGGA  >  NC_002947/5763708‑5764002

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