Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,329,091 T→C 6.1% A43A (GCA→GCG)  PP_4689 ← heme/hemin ABC transporter substrate‑binding protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,329,0910TC 6.1% 192.1 / 2.4 65A43A (GCA→GCG) PP_4689heme/hemin ABC transporter substrate‑binding protein
Reads supporting (aligned to +/- strand):  ref base T (36/25);  new base C (2/2);  total (38/27)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.35e-01

GCCCACCAGCACATCCGGGCGCAGGCTGAGGATGCCTTCGGCCGACAGTTGCCGCTGATAACCCACACTGGGCAAGGCTTTCAACGAGGCCGGGTGCTGACTGGTAGTGTCCACACCGACCAGGTGCGCTTCACCTCCCAATGCGGTTATCC‑ACTCGCTCAAGGCGCCACCCGCGCTGACCCAGCGTTGCGGCAGCTCGGCTGCCAAGGCTTGGGTAGAGAGCACAAGGCTTGCACACAGAGCGAGAAAAGCAGCTGGACGGCGCATCATCGGGTTCCTTCTGAAGGCG  >  NC_002947/5328950‑5329238
                                                                                                                                             |                                                                                                                                                    
gCCCACCAGCACATCCGGGCGCAGGCTGAGGATGCCTTCGGCCGACAGTTGCCGCTGATAACCCACACTGGGCAAGGCTTTCAACGAGGCCGGGTGCTGACTGGTAGTGTCCACACCGACCAGGTGCGCTTCACCTCCCAATGCGGTTa                                                                                                                                               >  8:198792/1‑149 (MQ=255)
  ccaccaGCACATCCGGGCGCAGGCTGAGGATGCCTTCGGCCGACAGTTGCCGCTGATAACCCACACTGGGCAAGGCTTTCAACGAGGCCGGGTGCTGACTGGTAGTGTCCACACCGACCAGGTGCGCTTCACCTCCCAATGCGGTTATc                                                                                                                                             <  7:200917/149‑1 (MQ=255)
    accaGCACATCCGGGCGCAGGCTGAGGATGCCTTCGGCCGACAGTTGCCGCTGATAACCCACACTGGGCAAGGCTTTCAACGAGGCCGGGTGCTGACTGGTAGTGTCCACACCGACCAGGTGCGCTTCACCTCCCAATGCg                                                                                                                                                   <  1:116333/141‑1 (MQ=255)
    accaGCACATCCGGGCGCAGGCTGAGGATGCCTTCGGCCGACAGTTGCCGCTGATAACCCACACTGGGCAAGGCTTTCAACGAGGCCGGGTGCTGACTGGTAGTGTCCACACCGACCAGGTGCGCTTCACCTCCCAATGCg                                                                                                                                                   >  2:116333/1‑141 (MQ=255)
    accaGCACATCCGGGCGCAGGCTGAGGATGCCTTCGGCCGACAGTTGCCGCTGATAACCCACACTGGGCAAGGCTTTCAACGAGGCCGGGTGCTGACTGGTAGTGTCCACACCGACCAGGTGCGCTTCACCTCCCAATGCGGTTATcc                                                                                                                                            >  3:130127/1‑148 (MQ=255)
        gCACATCCGGGCGCAGGCTGAGGATGCCTTCGGGCGACAGTTGCCGCTGATAACCCACACTGGGCAAGGCTTTCAACGAGGCCGGGTGCTGACTGGTAGTGTCCACACCGACCAGGTGCGCTTCACCTCCCAATGCGGTTATCC‑ACTCg                                                                                                                                      >  7:205524/1‑149 (MQ=255)
         cacaTCCGGGCGCAGGCTGAGGATGCCTTCGGCCGACAGTTGCCGCTGATAACCCACACTGGGCAAGGCTTTCAACGAGGCCGGGTGCTGACTGGTAGTGTCCACACCGACCAGGTGCGCTTCACCTCCCAATGCGGTTATCC‑ACTCGc                                                                                                                                     <  3:225787/149‑1 (MQ=255)
               cGGGCGCAGGCTTAGGATGCCTTCGGCCGACAGTTGCCGCTGATAACCCACACTGGGCAAGGCTTTCAACTAGGCCGGGTGCTGACTGGTAGTGTCAACACCGACCAGGTGCGCTTCACCTCCCAATGCGGTTATCC‑ACTCGCTCAAgg                                                                                                                               <  5:175124/149‑1 (MQ=255)
                      aGGCTGAGGATGCCTTCGGCCGACAGTTGCCGCTGTTAACCCACACTGGGCAAGGCTTTCAACGAGGCCGGGTGCTGACTGGTAGTGTCCACACCGACCAGGTGCGCTTCACCTCCCAATGCGGTTATCC‑ACTCGCTCAAGGCGCCAcc                                                                                                                        >  5:309198/1‑149 (MQ=255)
                            aGGATGCCTTCGGCCGACAGTTGCCGCTGATAACCCACACTGGGCAAGGCTTTCAAGGAGGCCGGGTGCTGACTGGTAGTGTCCACACCGACCAGGTGCGCTTCACCTCCCAATGCGGTTATCC‑ACTCGCTCAAGGCGCCACCcgcg                                                                                                                    <  3:332683/147‑1 (MQ=255)
                            aGGATGCCTTCGGCCGACAGTTGCCGCTGATAACCCACACTGGGCAAGGCTTTCAACGAGGCCGGGTGCTGACTGGTAGTGTCCACACCGACCAGGTGCGCTTCACCTCCCAATGCGGTTATCC‑ACTCGCTCAAGGCGCCACCCGCGCt                                                                                                                  <  8:250056/149‑1 (MQ=255)
                              gATGCCTTCGGCCGACAGTTGCCGCTGATAACCCACACTGGGCAAGGCTTTCAACGAGGCCGGGTGCTGACTGGTAGTGTCCACACCGACCAGGTGCGCTTCACCTCCCAATGCGGTTATCC‑a                                                                                                                                          >  2:156433/1‑123 (MQ=255)
                              gATGCCTTCGGCCGACAGTTGCCGCTGATAACCCACACTGGGCAAGGCTTTCAACGAGGCCGGGTGCTGACTGGTAGTGTCCACACCGACCAGGTGCGCTTCACCTCCCAATGCGGTTATCC‑a                                                                                                                                          <  1:156433/123‑1 (MQ=255)
                               aTGCCTTCGGCCGACAGTTGCCGCTGATAACCCACACTGGGCAAGGCTTTCAACGAGGCCGGGTGCTGACTGGTAGTGTCCACACAGACCAGGTGCGCTTCACCTCCCAATGCGGTTATCC‑ACTCGCTCAAGGCGCACCCCGCGCTGAc                                                                                                               >  4:252492/1‑149 (MQ=255)
                                   cTTCGGCCGACAGTTGCCGCTGATAACCCACACTGGGCAAGGCTTTCAACGAGGCCGGGTGCTGACTGGTAGTGTCCACACCGACCAGGTGCGCTTCACCTCCCAATGCGGTTATCC‑ACTCGCTCAAGGCGCCACCCGCGCTGACCCAg                                                                                                           >  8:306287/1‑149 (MQ=255)
                                   cTTCGGCCGACAGTTGCCGCTGATAACCCACACTGGGCAAGGCTTTCAACGAGGCCGGGTGCTGACTGGTAGTGTCCACACCGACCAGGTGCGCTTCACCTCCCAATGCGGTTATCC‑ACTCGCTCAAGGCGCCACCCGCGCTGACCCAg                                                                                                           >  5:173970/1‑149 (MQ=255)
                                    ttCGGCCGACAGTTGCCGCTGATAACCCACACTGGGCAAGGCTTTCAACGAGGCCGGGTGCTGACTGGTAGTGTCCACACCGACCAGGTGCGCTTCACCTCCCAATGCGGTTATCC‑ACTCGCTCAAGGCGGCCCCCGCGCTGCCCCAGc                                                                                                          >  8:258096/1‑149 (MQ=255)
                                       ggCCGACAGTTGCCGCTGATAACCCACACTGGGCAAGGCTTTCAACGAGGCCGGGTGCTGACTGGTAGTGTCCACACCGACCAGGTGCGCTTCACCTCCCAATGCGGTTATCC‑ACTCGCTCAAGGCGCCACCCGCGCTGACCCAGCGtt                                                                                                       <  4:130127/149‑1 (MQ=255)
                                        gCCGACAGTTGCCGCTGATAACCCACACTGGGCAAGGCTTTCAACGAGGCCGGGTGCTGACTGGTAGTGTCCACACCGACCAGGTGCGCTTCACCTCCCAATGCGGTTATCC‑ACTCGCTCAAGGCGCCACCCGCGCTGACCCAGCGTTg                                                                                                      >  4:229894/1‑149 (MQ=255)
                                         ccGACAGTTGCCGCTGATAACCCACACTGGGCAAGGCTTTCAACGAGGCCGGGTGCTGACTGGTAGTGTCCACACCGACCAGGTGCGCTTCACCTCCCAATGCGGTTATCC‑ACTCGCTCAAGGCGCCACCCGCGCTGACCCAGCGTTGc                                                                                                     <  3:196582/149‑1 (MQ=255)
                                           gACAGTTGCCGCTGTTAACCCACACTGGGCAAGGCTTTCAACGAGGCCGGGTGCTGACTGGTAGTGTCCACACCGACCAGGTGCGCTTCACCTCCCAATGCGGTTATCC‑ACTCGCTCAAGGCGCCACCCGCGCTGACCCAGCGTTGCgg                                                                                                   <  6:309198/149‑1 (MQ=255)
                                                          tAACCCACACTGGGCAAGGCTTTCAACGAGGCCGGGTGCTGACTGGTAGTGTCCACACCGACCAGGTGCGCTTCACCTCCCAATGCGGTTATCC‑ACTCGCTCAAGGCGCCACCCGCGCTGACCCAGCGTTGCGGCAGCTCGGCTGCCaa                                                                                    <  3:229894/149‑1 (MQ=255)
                                                              ccACACTGGGCAAGGCTTTCAACGAGGCCGGGTGCTGACTGGTAGTGTCCACACCGCCCAGGTGCGCTTCACCTCCCAATGCGGTTATCC‑ACTCGCTCAAGTCGCAACCCGAGCTGACCCAGCGTTGCGGCAGCTCGGCTCCCAAGGCt                                                                                >  4:67557/1‑149 (MQ=255)
                                                                 cacTGGGCAAGGCTTTCAACGAGGCCGGGTGCTGACTGGTAGTGTCCACACCGACCAGGTGCGCTTCACCTCCCAATGCGGTTATCC‑ACTCGCTCAAGGCGCCACCCGCGCTGACCCAGCGTTGCGGCAGCTCGGCTGCCAAGGCTTgg                                                                             <  5:212810/149‑1 (MQ=255)
                                                                     gggCAAGGCTTTCAACGAGGCCGGGTGCTGACTGGTAGTGTCCACACCGACCAGGTGCGCTTCACCTCCCAACGCGGTTATCC‑ACTCGCTCAAGGCGCCACCCGCGCTGACCCAGCGTTGCGGCAGCTCGGCTGCCAAGGCTTGGGt                                                                           <  2:96295/147‑1 (MQ=255)
                                                                     gggCAAGGCTTTCAACGAGGCCGGGTGCTGACTGGTAGTGTCCACACCGACCAGGTGCGCTTCACCTCCCAACGCGGTTATCC‑ACTCGCTCAAGGCGCCACCCGCGCTGACCCAGCGTTGCGGCAGCTCGGCTGCCAAGGCTTGGGt                                                                           >  1:96295/1‑147 (MQ=255)
                                                                      ggCAAGGCTTTCAACGAGGCCGGGTGCTGACTGGTAGTGTCCACACCGACCAGGTGCGCTTCACCTCCCAATGCGGTTATCC‑ACTCGCTCAAGGCGCCACCCGCGCTGACCCAGCGTTGCGGCAGCTCGGCTGCCAAGGCTTGGGTa                                                                          <  6:173970/147‑1 (MQ=255)
                                                                              tttCAACGAGGCCGGGTGCTGACTGGTAGTGTCCACACCGACCAGGTGCGCTTCACCTCCCAACGCGGTTATCC‑ACTCGCTCAAGGCGCCACCCGCGCTGACCCAGCGTTGCGGCAGCTCGGCt                                                                                         >  2:286028/1‑124 (MQ=255)
                                                                              tttCAACGAGGCCGGGTGCTGACTGGTAGTGTCCACACCGACCAGGTGCGCTTCACCTCCCAACGCGGTTATCC‑ACTCGCTCAAGGCGCCACCCGCGCTGACCCAGCGTTGCGGCAGCTCGGCt                                                                                         <  1:286028/124‑1 (MQ=255)
                                                                                tCAACGAGGCCGGGTGCTGACTGGTAGTGTCCACACCGACCAGGTGCGCTTCACCTCCCAATGCGGTTATCC‑ACTCGCTCAAGGCGCCACCCGCGCTGACCCAGCGTTGCGGCAGCTCGGCTGCCAAGGCTTGGGTAGAGAGCACAAgg                                                              >  3:326507/1‑149 (MQ=255)
                                                                                 cAACGAGGCCGGGTGCTGACTGGTAGTGTCCACACCGACCAGGTGCGCTTCACCTCCCAATGCGGTTATCC‑ACTCGCTCAAGGCGCCACCCGCGCTGACCCAGCGTTGCGGCAGCTCGGCTGCCAAGGCTTGGGTAGAGAGCACAAgg                                                              <  7:198792/148‑1 (MQ=255)
                                                                                  aaCGAGGCCGGGTGCTGACTGGTAGTGTCCACACCGACCAGGTGCGCTTCACCTCCCAATGCGGTTATCC‑ACTCGCTCAAGGCGCCACCCGCGCTGACCCAGCGTTGCGGCAGCTCGGCTGCCAAGGCTTGGGTAGAGAGCACAAGGCt                                                            >  2:261663/1‑149 (MQ=255)
                                                                                  aaCGAGGCCGGGTGCTGACTGGTAGTGTCCACACCGACCAGGTGCGCTTCACCTCCCAATGCGGTTATCC‑ACTCGCTCAAGGCGCCACCCGCGCTGACCCAGCGTTGCGGCAGCTCGGCTGCCAAGGCTTGGGTAGAGAGCACAAGGCt                                                            >  6:130056/1‑149 (MQ=255)
                                                                                        gCCGGGTGCTGACTGGTAGTGTCCACACCGACCAGGTGCGCTTCACCTCCCAATGCGGTTATCC‑ACTCGCTCAAGGCGCCACCCGCGCTGACCCAGCGTTGCGGCAGCTCGGCTGCCAAGGCTTGGGTAGAGAGCACAAGGCTTGcaca                                                      >  7:154310/1‑149 (MQ=255)
                                                                                        gCCGGGTGCTGACTGGTAGTGTCCACACCGACCAGGTGCGCTTCACCTCCCAATGCGGTTATCC‑ACTCGCTCAAGGCGCCACCCGCGCTGACCCAGCGTTGCGGCAGCTCGGCTGCCAAGGCTTGGGTAGAGAGCACAAGGCTTGcaca                                                      >  3:135115/1‑149 (MQ=255)
                                                                                           gggTGCTGACTGGTAGTGTCCACACCGACCAGGTGCGCTTCACCTCCCAATGCGGTTATCC‑ACTCGCTCAAGGCGCCACCCGCGCTGACCCAGCGTTGCGGCAGCTCGGCTGCCAAGGCTTGGGTAGAGAGCACAAGGCTTGcaca                                                      >  6:207451/1‑146 (MQ=255)
                                                                                             gTGCTGACTGGTAGTGTCCACACCGACCAGGTGCGCTTCACCTCCCAATGCGGTTATCC‑ACTCGCTCAAGGCGCCACCCGCGCTGACCCAGCGTTGCGGCAGCTCGGCTGCCAAGGCTTGGGTAGAGAGCACAAGGCTTGCACACagag                                                 >  6:17953/1‑149 (MQ=255)
                                                                                              tGCTGACTGGTAGTGTCCACACCGACCAGGTGCGCTTCACCTCCCAATGCGGTTATCC‑ACTCGCTCAAGGCGCCACCCGCGCTGACCCAGCGTTGCGGCAGCTCGGCTGCCAAGGCTTGGGTAGAGAGCACAAGGCTTGCACACAGAGc                                                <  3:89755/149‑1 (MQ=255)
                                                                                                cTGACTGGTAGTGTCCACACCGACCAGGTGCGCTTCACCTCCCAATGCGGTTATCC‑ACTCGCTCAAGGCGCCACCCGCGCTGACCCAGCGTTGCGGCAGCTCGGCTGCCAAGGCTTGGGTAGAGAGCACAAGGCTTGCACACAGAGCga                                              <  8:154310/149‑1 (MQ=255)
                                                                                                 tGACTTGTAGTGTCCACACAGACCAGGTGCGCTTCAACTCCCAATGCGGTTATCC‑ACTCGCTCAAGGCGCCACCCGCGCTGACCCAGCGTTGCGGCAGCTCGGCTGCCAAGGCTTGGGTAGAGAGCACAAGGCTTGCACACAGAGCgag                                             <  3:252492/149‑1 (MQ=255)
                                                                                                 tGACTGGTAGTGTCCACACCGACCAGGTGCGCTTCACCTCCCAATGCGGTTATCC‑ACTCGCTCAAGGCGCCACCCGCGCTGACCCAGCGTTGCGGCAGCTCGGCTGCCAAGGCTTGGGTAGAGAGCACAAGGCTTGCACACAGAGc                                                >  8:35592/1‑146 (MQ=255)
                                                                                                   aCTGGTAGTGTCCACACCGACCAGGTGCGCTTCACCTCCCAATGCGGTTATCC‑ACTCGCTCAAGGCGCCACCCGCGCTGACCCAGCGTTGCGGCAGCTCGGCTGCCCAGGCTTGGGTAGAGAGCACAAGGCTTGCACACAGAGCgaga                                            >  8:251210/1‑148 (MQ=255)
                                                                                                     tGGTAGTGTCCACACCGACCAGGTGCGCTTCACCTCCCAATGCGGTTATCCGACACGCTCAAGGCGCCACCCGCGATGACCCAGCGTTGCGGCAGCTCGGATGCCAAGGCTTGGGTAGAGAGCACAAGCCTTCCACACAGAGCGAGaa                                           >  1:85360/1‑148 (MQ=255)
                                                                                                       gTAGTGTCCACACCGACCAGGTGCGCTTCACCTCCCAATGCGGTTATCC‑ACTCGCTCAAGGCGCCACCCGCGCTGACCCAGCGTTGCGGCAGCTCGGCTGCCAAGGCTTGGGTAGAGAGCACAAGGCTTGCACACAGAGCGAGaaaa                                         <  5:17953/147‑1 (MQ=255)
                                                                                                       gTAGTGTCCACACCGACCAGGTGCGCTTCACCTCCCAATGCGGTTATCC‑ACTCGCTCAAGGCGCCACCCGCGCTGACCCAGCGTTGCGGCAGCTCGGCTGCCAAGGCTTGGGTAGAGAGCACAAGGCTTGCACACAGAGCGAGAAAag                                        >  4:328861/1‑148 (MQ=255)
                                                                                                        tAGTGTCCACACCGACCAGGTGCGCTTCACCTCCCAATGCGGTTATCC‑ACTCGCTCAAGGCGCCACCCGCGCTGACCCAGCGTTGCGGCAGCTCGGCTGCCAAGGCTTGGGTAGAGAGCACAAGGCTTGCACACAGAGCGAGAAAag                                        >  7:84323/1‑147 (MQ=255)
                                                                                                             tCCACACCGACCAGGTGCGCTTCACCTCCCAATGCGGTTATCC‑ACTCGCTCAAGGCGCCACCCGCGCTGACCCAGCGTTGCGGCAGCTCGGCTGCCAAGGCTTGGGTAGAGAGCACAAGGCTTGCACACAGAGCGAGAAAAGCAGCTgg                                 >  4:96444/1‑149 (MQ=255)
                                                                                                                acacCGACCAGGTGCGCTTCACCTCCCAATGCGGTTATCC‑ACTTGCTCAAGGCGCCACCCGCGCTGACCCAGCGTTGCGGCAGCTCGGCTGCCAAGGCTTGGGTAGAGAGCACAAGGCTTGCACACAGAGCGAGAAAAGCAGCTGGACg                              >  6:258722/1‑149 (MQ=255)
                                                                                                                acacCGACCAGGTGCGCTTCACCTCCCAATGCGGTTATCC‑ACTCGCTCAAGGCGCCACCCGCGCTGACCCAGCGTTGCGGCAGCTCGGCTGCCAAGGCTTGGGTAGAGAGCACAAGGCTTGCACACAGAGCGAGAAAAGCAGCTGGACg                              >  1:116504/1‑149 (MQ=255)
                                                                                                                acacCGACCAGGTGCGCTTCACCTCCCAATGCGGTTATCC‑ACTCGCTCAAGGCGCCACCCGCGCTGACCCAGCGTTGCGGCAGCTCGGCTGCCAAGGCTTGGGTAGAGAGCACAAGGCTTGCACACAGAGCGAGAAAAGCAGCTGGACg                              >  6:142129/1‑149 (MQ=255)
                                                                                                                 cacCGACCAGGTGCGCTTCACCTCCCAATGCGGTTATCC‑ACTCGCTCAAGGCGCCACCCGCGCTGACCCAGCGTTGCGGCAGCTCGGCTGCCAAGGCTTGGGTAGAGAGCACAAGGCTTGCACACagag                                                 <  7:172263/129‑1 (MQ=255)
                                                                                                                 cacCGACCAGGTGCGCTTCACCTCCCAATGCGGTTATCC‑ACTCGCTCAAGGCGCCACCCGCGCTGACCCAGCGTTGCGGCAGCTCGGCTGCCAAGGCTTGGGTAGAGAGCACAAGGCTTGCACACagag                                                 >  8:172263/1‑129 (MQ=255)
                                                                                                                    cGACCAGGTGCGCTTCACCTCCCAATGCGGTTATCC‑ACTCGCTCAAGGCGCCACCCGCGCTGACCCAGCGTTGCGGCAGCTCGGCTGCCAAGGCTTGGGTAGAGAGCACAAGGCTTGCACACAGAGCGAGAAAAGCAGCTGGACGgcgc                          <  2:85360/149‑1 (MQ=255)
                                                                                                                    cGACCAGGTGCGCTTCACCTCCCAATGCGGTAATCC‑ACTCGCTCAAGGCGCCACCCGCGCTGACCCAGCGTTGCGGCAGCTCGGCTGCCAAGGCTTGGGTAGAGAGCACAAGGCTTGCACACAGAGCgag                                             <  2:86631/130‑1 (MQ=255)
                                                                                                                    cGACCAGGTGCGCTTCACCTCCCAATGCGGTAATCC‑ACTCGCTCAAGGCGCCACCCGCGCTGACCCAGCGTTGCGGCAGCTCGGCTGCCAAGGCTTGGGTAGAGAGCACAAGGCTTGCACACAGAGCgag                                             >  1:86631/1‑130 (MQ=255)
                                                                                                                          ggTGCGCTTCACCTCCCAATGCGGTTATCC‑ACTCGCTCAAGGCGCCACCCGCGCTGACCCAGCGTTGCGGCAGCTCGGCTGCCAAGGCTTGGGTAGAGAGCACAAGGCTTGCACACAGAGCGAGAAAAGCAGCTGGACGGCGCatcatc                    >  1:8836/1‑149 (MQ=255)
                                                                                                                            tGCGCTTCACCTCCCAATGCGGTTATCC‑ACTCGCTCAAGGCGCCACCCGCGCTGACCCAGCGTTGCGGCAGCTCGGCTGCCAAGGCTTGGGTAGAGAGCACAAGGCTTGCACACAGAGCGAGAAAAGCAGCTGGACGGCGCATCATCgg                  <  8:205524/149‑1 (MQ=255)
                                                                                                                              cgcTTCACCTCCCAATGCGGTTATCC‑ACTCGCTCAAGGCGCCACCCGCGCTGACCCAGCGTTGCGGCAGCTCGGCTGCCAAGGCTTGGGTAGAGAGCACAAGGCTTGCACACAGAGCGAGAAAAGCAGCTGGACGGCGCATCATCggg                 <  4:326507/148‑1 (MQ=255)
                                                                                                                               gcTTCACCTCCCAATGCGGTTATCC‑ACTCGCTCAAGGCGCCACCCGCGCTGACCCAGCGTTGCGGCAGCTCGGCTGCCAAGGCTTGGGTAGAGAGCACAAGGCTTGCACACAGAGCGAGAAAAGCAGCTGGACGGCGCATCATCGGGtt               <  4:135115/149‑1 (MQ=255)
                                                                                                                                    aCCTCCCAATGCGGTTATCC‑ACTCGCTCAAGGCGCCACCCGCGCTGACCCAGCGTTGCGGCAGCTCGGCTGCCAAGGCTTGGGTAGAGAGCACAAGGCTTGCACACAGAGCGAGAAAAGCAGCTGGACGGCGCATCATCGGGTTCCTTc          <  3:67557/149‑1 (MQ=255)
                                                                                                                                      cTCCCAATGCGGTTATCC‑ACTCGCTCAAGGCGCCACCCGCGCTGACCCAGCGTTGCGGCAGCTCGGCTGCCAAGGCTTGGGTAGAGAGCACAAGGCTTGCACACAGAGCGAGAAAAGCAGCTGGACGGCGCATCATCGGGTTCCTTCt         >  2:198409/1‑148 (MQ=255)
                                                                                                                                       tCCCAATGCGGTTATCC‑ACTCGCTCAAGGCGCCACCCGCGCTGACCCAGCGTTGCGGCAGCTCGGCTGCCAAGGCTTGGGTAGAGAGCACAAGGCTTGCACACAGAGCGAGAAAAGCAGCTGGACGGCGCATCATCGGGTTCCTTCTg        >  7:170331/1‑148 (MQ=255)
                                                                                                                                         ccAATGCGGTTATCC‑ACTCGCTCAAGGCGCCACCCGCGCTGACCCAGCGTTGCGGCAGCTCGGCTGCCAAGGCTTGGGTAGAGAGCACAAGGCTTGCACACAGAGCGAGAAAAGCAGCTGGACGGCGCATCATCGGGTTCCTTCTGAAg     >  8:145510/1‑149 (MQ=255)
                                                                                                                                           aaTGCGGTTATCC‑ACTCGCTCAAGGCGCCACCCGCGCTGACCCAGCGTTGCGGCAGCTCGGCTGCCAAGGCTTGGGTAGAGAGCACAAGGCTTGCACACAGAGCGAGAAAAGCAGCTGGACGGCGCATCATCGGGTTCCTTCTGAAGGc   >  2:276307/1‑149 (MQ=255)
                                                                                                                                            aTGCGGTTATCC‑ACTCGCTCAAGGCGCCACCCGCGCTGACCCAGCGTTGCGGCAGCTCGGCTGCCAAGGCTTGGGTAGAGAGCACAAGGCTTGCACACAGAGCGAGAAAAGCAGCTGGACGGCGCATCATCGGGTTCCTTCTGAAGGCg  >  3:263010/1‑149 (MQ=255)
                                                                                                                                             |                                                                                                                                                    
GCCCACCAGCACATCCGGGCGCAGGCTGAGGATGCCTTCGGCCGACAGTTGCCGCTGATAACCCACACTGGGCAAGGCTTTCAACGAGGCCGGGTGCTGACTGGTAGTGTCCACACCGACCAGGTGCGCTTCACCTCCCAATGCGGTTATCC‑ACTCGCTCAAGGCGCCACCCGCGCTGACCCAGCGTTGCGGCAGCTCGGCTGCCAAGGCTTGGGTAGAGAGCACAAGGCTTGCACACAGAGCGAGAAAAGCAGCTGGACGGCGCATCATCGGGTTCCTTCTGAAGGCG  >  NC_002947/5328950‑5329238

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