Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,763,847 T→C 6.1% V408A (GTT→GCT)  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,8470TC 6.1% 243.8 / 4.4 82V408A (GTT→GCT) PP_5057M23/M37 family peptidase
Reads supporting (aligned to +/- strand):  ref base T (45/32);  new base C (3/2);  total (48/34)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.99e-01

GTGGTATTCGCGGACTGGTTGCGCGGTGCCGGGCTTCTGGTCATTCTCGACCATGGCAACGGTTACCTGAGCCTGTACGGCCACAACCAGAGCCTGCTCAAGAATGCCGGTGATATCGTCAAGGCCGGTGAAGCCATTTCCACCGTTGGCGACAGCGGTGGTCAGGACGCGGCAGGCTTGTATTTCGCTATCCGCCAGCAGGGCCGGCCTACCGACCCTTCGCAGTGGTGCCGGGGCTAGTCAAATTTCTACGGCGCAGC‑CTGAGTGCTGCGCCGATGCTCCCCAGG  >  NC_002947/5763702‑5763988
                                                                                                                                                 |                                                                                                                                              
gtgGTATTCGCGGACTGGTTGCGCGGTGCCGGGCTTCTGGTCATTCTCGACCATGGCAACGGTTACCTGAGCCTGTACGGCCACAACCAGAGCCTGCTCAAGAATGCCGGTGATATCGTCAAGGCCGGTGAAGCCATTTCCACCGTTgg                                                                                                                                             <  8:3993/149‑1 (MQ=255)
 tgGTATTCGCGGACTGGTTGCGCGGTGCCGGGCTTCTGGTCATTCTCGACCATGGCAACGGTTACCTGAGCCTGTACGGCCACAACCAGAGCCTGCTCAAGAATGCCGGTGATATCGTCAAGGCCGGTGAAGCCATTTCCACCGTTGGc                                                                                                                                            <  3:160062/149‑1 (MQ=255)
 tgGTATTCGCGGACTGGTTGCGCGGTGCCGGGCTTCTGGTCATTCTCGACCATGGCAACGGTTACCTGAGCCTGTACGGCCACAACCAGAGCCTGCTCAAGAATGCCGGTGATATCGTCAAGGCCGGTGAAGCCATTTCCACCGTTGGc                                                                                                                                            >  2:88357/1‑149 (MQ=255)
 tgGTATTCGCGGACTGGTTGCGCGGTGCCGGGCTTCTGGTCATTCTCGACCATGGCAACGGTTACCTGAGCCTGTACGGCCACAACCAGAGCCTGCTCAAGAATGCCGGTGATATCGTCAAGGCCGGTGAAGCCATTTCCACCGTTGGc                                                                                                                                            >  6:232075/1‑149 (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                                                                                                      >  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                                                                                                     >  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                                                                                                       >  8:142518/1‑104 (MQ=255)
                                                                                   caACCAGAGCCTGCTCAAGAATGCCGGTGATATCGTCAAGGCCGGTGAAGCCATTTCCACCGTTGGCGACAGCGGTGGTCAGGACGCGGCAGGCTTGTATTTCg                                                                                                       <  7:142518/104‑1 (MQ=255)
                                                                                    aaCCAGAGCCTGCTCAAGAATGCCGGTGATATCGTCAAGGCCGGTGAAGCCATTTCCACCGTTGGCGACAGCGGTGGTCAGGACGCGGCAGGCTTGTATTTCGCTATCCGCCAGCAGGGCCGGCCTACCGACCCTTCGCAGTGGTGCCg                                                         >  7:247635/1‑149 (MQ=255)
                                                                                       cAGAGCCTGCTCAAGAATGCCGGTGATATCGTCAAGGCCGGTGAAGCCATTTCCACCGTTGGCGACAGCGGTGGTCAGGACGCGGCAGGCTTGTATTTCGCTATCCGCCAGCAGGGCCGGCCTAc                                                                              <  3:271849/125‑1 (MQ=255)
                                                                                       cAGAGCCTGCTCAAGAATGCCGGTGATATCGTCAAGGCCGGTGAAGCCATTTCCACCGTTGGCGACAGCGGTGGTCAGGACGCGGCAGGCTTGTATTTCGCTATCCGCCAGCAGGGCCGGCCTAc                                                                              >  4:271849/1‑125 (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                       >  2:323646/1‑149 (MQ=255)
                                                                                                                     gTCAAGGCCGGTGAAGCCATTTCCACCGTTGGCGACAGCGGTGGTCAGGACGCGGCAGGCTTGTATTTCGCTATCCGCCAGCAGGGCCGGCCTACCGACCCTTCGCAGTGGTGCCGGGGCTAGTCAAATTTCTACGGCGCAGC‑CTGAGt                       <  3:61097/149‑1 (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      >  8:166724/1‑148 (MQ=255)
                                                                                                                                       aTTTCCACCGTTGGCGACAGCGGTGGTCAGGACGCGGCAGGCTTGTATTTCGCTATCCGCCAGCAGGGCCGGCCTACCGACCCTTCGCAGTGGTGCCGGGGCTAGTCAAATTTCTACGGCGCAGC‑CTGAGTGCTGCGCCGATGCTccc      >  2:48437/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                            <  5:212549/122‑1 (MQ=255)
                                                                                                                                          accACCGCTGTCGCCAACGGTGGTCAGGACGCGGCAGGCTTGTATTTCGCTATCCGCCAGCAGGGCCGGCCTACCGACCCTTCGCAGTGGTGCCGGGGCTAGTCAAATTTCTACGGCGCAGc‑c                            >  6:212549/2‑123 (MQ=255)
                                                                                                                                           ccACCGTTGGCGACAGCGGTGGTCAGGACGCGGCAGGCTTGTATTTCGCTATCCGCCAGCAGGGCCGGCCTACCGACCCTTCGCAGTGGTGCC‑GGGCTAGTCAAATTTCTACGGCGCAGC‑CTGAGTGCTGCGCCGATGCTCCCCAgg  >  5:361281/1‑147 (MQ=255)
                                                                                                                                                 |                                                                                                                                              
GTGGTATTCGCGGACTGGTTGCGCGGTGCCGGGCTTCTGGTCATTCTCGACCATGGCAACGGTTACCTGAGCCTGTACGGCCACAACCAGAGCCTGCTCAAGAATGCCGGTGATATCGTCAAGGCCGGTGAAGCCATTTCCACCGTTGGCGACAGCGGTGGTCAGGACGCGGCAGGCTTGTATTTCGCTATCCGCCAGCAGGGCCGGCCTACCGACCCTTCGCAGTGGTGCCGGGGCTAGTCAAATTTCTACGGCGCAGC‑CTGAGTGCTGCGCCGATGCTCCCCAGG  >  NC_002947/5763702‑5763988

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