Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,551,996 T→C 24.3% S409P (TCA→CCA)  PP_4882 → iron ABC transporter permease

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,551,9960TC24.3% 117.7 / 37.6 70S409P (TCA→CCA) PP_4882iron ABC transporter permease
Reads supporting (aligned to +/- strand):  ref base T (29/24);  new base C (9/8);  total (38/32)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TCTGCCCTATCTGCTGCACGCCCTGCCCGGCCTGGTGATTGCCCTGACGTTGGTGTTCTTCGCCCTGCACTACGTGCCGGCGCTGTACCAGACCACGGCCTTGCTGCTGTTGGCCTATGCGCTGCTGTTCCTGCCGTTGGCGCAATCACCGGTGCGCACCGCGCTGAACAAGGCTTCGCCGACGCTGGAAGAAGCCGCGCGCACCCTTGGCGCGAGCAGCTTTGCCGCGTTCTGCCGGGTGACCTTGCCGATCATCTTCCCGGCCATGGCGGCAGCCTTTGCGCTGGTGTTC  >  NC_002947/5551851‑5552142
                                                                                                                                                 |                                                                                                                                                  
tCTGCCCTATCTGCTGCACGCCCTGCCCGGCCTGGTGATTGCCCTGACGTTGGTGTTCTTCGCACTGCACTACGTGCCGGCGCTGTACCAGACCACGGCCTTGCTGCTGTTGGCCTATGCGCTGCTGTTCCTGCCGTTGGCGCAATCAc                                                                                                                                                 >  5:330675/1‑149 (MQ=255)
    cccTATCTGCTGCACGCCCTGCCCGGCCTGGTGATTGCCCTGACGTTGGTGTTCTTCGCCCTGCACTACGTGCCGGCGCTGTACCAGACCACGGCCTTGCTGCTGTTGGCCTATGCGCTGCTGTTCCTGCCGTTGGCGCAATCACCg                                                                                                                                               <  2:117076/147‑1 (MQ=255)
    cccTATCTGCTGCACGCCCTGCCCGGCCTGGTGATTGCCCTGACGTTGGTGTTCTTCGCCCTGCACTACGTGCCGGCGCTGTACCAGACCACGGCCTTGCTGCTGTTGGCCTATGCGCTGCTGTTCCTGCCGTTGGCGCAATCACCGGt                                                                                                                                             <  4:206609/149‑1 (MQ=255)
    cccTATCTGCTGCACGCCCTGCCCGGCCTGGTGATTGCCCTGACGTTGGTGTTCTTCGCCCTGCACTACGTGCCGGCGCTGTACCAGACCACGGCCTTGCTGCTGTTGGCCTATGCGCTGCTGTTCCTGCCGTTGGCGCAACCACCg                                                                                                                                               <  2:103450/147‑1 (MQ=255)
    cccTATCTGCTGCACGCCCTGCCCGGCCTGGTGATTGCCCTGACGTTGGTGTTCTTCGCCCTGCACTACGTGCCGGCGCTGTACCAGACCACGGCCTTGCTGCTGTTGGCCTATGCGCTGCTGTTCCTGCCGTTGGCGCAACCACCGGt                                                                                                                                             >  7:146969/1‑149 (MQ=255)
       tATCTGCTGCACGCCCTGCCCGGCCTGGTGATTGCCCTGACGTTGGTGTTCTTCGCCCTGCACTACGTGCCGGCGCTGTACCAGACCACGGCCTTGCTGCTGTTGGCCTATGCGCTGCTGTTCCTGCCGTTGGCGCAATCACCGGTgcg                                                                                                                                          >  4:6630/1‑149 (MQ=255)
             ctgcACGCCCTGCCCGGCCTAGTGATTGCCCTGACGTTGGTGTTCTTCGCCCTGCACTACGTGCCGGCGCTGTACCAGACCACGGCCTTGCTGCTGTTGGCCTATGCGCTGCTGTTCCTGCCGTTGGCGCAATCACCGGTGCGCa                                                                                                                                        <  5:69393/145‑1 (MQ=255)
                cACGCCCTGCCCGGCCTGGTGATTGCCCTGACGTTGGTGTTCTTCGCCCTGCACTACGTGCCGGCGCTGTACCAGACCACGGCCTTGCTGCTGTTGGCCTATGCGCTGCTGTTCCTGCCGTTGGCGCAATCACCGGTGCGCACCACGCt                                                                                                                                 >  6:345011/1‑149 (MQ=255)
                cACGCCCTGCCCGGCCTGGTGATTGCCCTGACGTTGGTGTTCTTCGCCCTGCACTACGTGCCGGCGCTGTACCAGACCACGGCCTTGCTGCTGTTGGCCTATGCGCTGCTGTTCCTGCCGTTGGCGCAACCACCGGTGCGCACCGCGCt                                                                                                                                 <  8:289076/149‑1 (MQ=255)
                    cccTGCCCGGCCTGGTGATTGCCCTGACGTTGGTGTTCTTCGCCCTGCACTACGTGCCGGCGCTGTACCAGACCACGGCCTTGCTGCTGTTGGCCTATGCGCTGCTGTTCCTGCCGTTGGCGCAATCACCGGTGCGCACCGCGCTGAAc                                                                                                                             <  5:307512/149‑1 (MQ=255)
                     ccTGCCCGGCCTGGTGATTGCCCTGACGTTGGTGTTCTTCGCCCTGCACTACGTGCCGGCGCTGTACCAGACCACGGCCTTGCTGCTGTTGGCCTATGCGCTGCTGTTCCTGCCGTTGGCGCAATCACCGGTGCGCACCGCGCTGAAc                                                                                                                             >  4:14499/1‑148 (MQ=255)
                      cTGCCCGGCCTGGTGATTGCCCTGACGTTGGTGTTCTTCGCCCTGCACTACGTGCCGGCGCTGTACCAGACCACGGCCTTGCTGCTGTTGGCCTATGCGCTGCTGTTCCTGCCGTTGGCGCAATCACCGGTGCGCACCTCGCTGAACaa                                                                                                                           <  5:345011/149‑1 (MQ=255)
                        gCCCGGCCTGGTGATTGCCCTGACGTTGGTGTTCTTCGCCCTGCACTACGTGCCGGCGCTGTACCAGACCACGGCCTTGCTGCTGTTGGCCTATGCGCTGCTGTTCCTGCCGTTGGCGCAATCACCGGTGCGCACCGCGCTGAACAAgg                                                                                                                         >  8:265015/1‑149 (MQ=255)
                                tGGTGATTGCCCTGACGTTGGTGTTCTTCGCCCTGCACTACGTGCCGGCGCTGTACCAGACCACGGCCTTGCTGCTGTTGGCCTATGCGCTGCTGTTCCTGCCGTTGGCGCAATCACCGGTGCGCACCGCGCTGAACAAGGCTTCGcc                                                                                                                  >  7:29399/1‑148 (MQ=255)
                                    gATTGCCCTGACGTTGGTGTTCTTCGCCCTGCACTACGTGCCGGCGCTGTACCAGACCACGGCCTTGCTGCTGTTGGCCTATGCGCTGCTGTTCCTGCCGTTGGCGCAACCACCGGTGCGCACCGCGCTGAACAAGGCTTCGCCGACGc                                                                                                             >  6:146872/1‑149 (MQ=255)
                                     aTTGCCCTGACGTTGGTGTTCTTCGCCCTGCACTACGTGCCGGCGCTGTACCAGACCACGGCCTTGCTGCTGTTGGCCTATGCGCTGCTGTTCCTGCCGTTGGCGCAATCACCGGTGCGCACCGCGCTGAACAAGGCTTCGCCGACGCt                                                                                                            >  8:11832/1‑149 (MQ=255)
                                        gCCCTGACGTTGGTGTTCTTCGCCCTGCACTACGTGCCGGCGCTGTACCAGACCACGGCCTTGCTGCTGTTGGCCTATGCGCTGCTGTTCCTGCCGTTGGCGCAATCACCGGTGCGCACCGCGCTGAACAAGGCTTCGCCGACGCTgg                                                                                                          <  2:168187/148‑1 (MQ=255)
                                          ccTGACGTTGGTGTTCTTCGCCCTGCACTACGTGCCGGCGCTGTACCAGACCACGGCCTTGCTGCTGTTGGCCTATGCGCTGCTGTTACTGCCGTTGGCGCAATCACCGGTGCGCACCGCGCAGAACAAGGCTTCGCCGACGCTgg                                                                                                          >  5:27815/1‑146 (MQ=255)
                                           cTGACGTTGGTGTTCTTCGCCCTGCACTACGTGCCGGCGCTGTACCAGACCACGGCCTTGCTGCTGTTGGCCTATGCGCTGCTGTTCCTGCCGTTGGCGCAATCACCGGTGCGCACCGCGCTGAACAAGGCTTCGCCGACGCTGgaag                                                                                                       >  7:124092/1‑148 (MQ=255)
                                                    gtgtTCTTCGCCCTGCACTACGTGCCGGCGCTGTACCAGACCACGGCCTTGCTGCTGTTGGCCTATGCGCTGCTGTTCCTGCCGTTGGCGCAATCACCGGTGCGCACCGCGCTGAACAAGGCTTCGCCGACGCTGGAAGAAGCcgcgc                                                                                              <  3:256403/148‑1 (MQ=255)
                                                     tgtTCTTCGCCCTGCACTACGTGCCGGCGCTGTACCAGACCACGGCCTTGCTGCTGTTGGCCTATGCGCTGCTGTTCCTGCCGTTGGCGCAATCACCGGTGCGCACCGCGCTGAACAAGGCTTCGCCGACGCTGGAAGAAGCCgcgcgc                                                                                            >  5:35322/1‑149 (MQ=255)
                                                         cttcGCCCTGCACTACGTGCCGGCGCTGTACCAGACCACGGCCTTGCTGCTGTTGGCCTATGCGCTGCTGTTCCTGCCGTTGGCGCAATCACCGGTGCGCACCGCGCTGAACAAGGCTTCGCCGACGCTGgaagaa                                                                                                     <  7:187367/136‑1 (MQ=255)
                                                         cttcGCCCTGCACTACGTGCCGGCGCTGTACCAGACCACGGCCTTGCTGCTGTTGGCCTATGCGCTGCTGTTCCTGCCGTTGGCGCAATCACCGGTGCGCACCGCGCTGAACAAGGCTTCGCCGACGCTGgaagaa                                                                                                     >  8:187367/1‑136 (MQ=255)
                                                         cttcGCCCTGCACTACGTGCCGGCGCTGTACCAGACCACGGCCTTGCTGCTGTTGGCCTATGCGCTGCTGTTCCTGCCGTTGGCGCAATCACCGGTGCGCACCGCGCTGAACAAGGCTTCGCCGACGCTGGAAGAAGCCGCGCGCAccc                                                                                        >  2:274903/1‑149 (MQ=255)
                                                               ccTGCACTACGTGCCGGCGCTGTACCAGACCACGGCCTTGCTGCTGTTGGCCTATGCGCTGCTGTTCCTGCCGTTGGCGCAATCACCGGTGCGCACCGCGCTGAACAAGGCTTCGCCGACGCTGGAAGAAGCCGCGCGCACCCTTGgcg                                                                                  <  7:11832/149‑1 (MQ=255)
                                                                cTGCACTACGTGCCGGCGCTGTACCAGACCACGGCCTTGCTGCTGTTGGCCTATGCGCTGCTGTTCCTGCCGTTGGCGCAACCACCGGTGCGCACCGCGCTGAACAAGGCTTCGCCGACGCTGGAAGAAGCCGCGCGCACCCTTGgcgc                                                                                 <  5:146872/149‑1 (MQ=255)
                                                                      tACGTGCCGGCGCTGTACCAGACCACGGCCTTGCTGCTGTTGGCCTATGCGCTGCTGTTCCTGCCGTTGGCGCAATCACCGGTGCGCACCGCGCTGAACAAGGCTTCGCCGACGCTGGAAGAAGCCGCGCGCACCCTTGGCGCGagcag                                                                           >  3:116013/1‑149 (MQ=255)
                                                                          tGCCGGCGCTGTACCAGACCACGGCCTTGCTGCTGTTGGCCTATGCGCTGCTGTTCCTGCCGTTGGCGCAATCACCGGTGCGCACCGCGCTGAACAAGGCTTCGCCGACGCTGGAAGAAGCCGCGCGCACCCTTGGCGCGAGCAGCttt                                                                       <  3:6630/149‑1 (MQ=255)
                                                                          tGCCGGCGCTGTACCAGACCACGGCCTTGCTGCTGTTGGCCTATGCGCTGCTGTTCCTGCCGTTGGCGCAACCACCGGTGCGCACCGCGCTGAACAAGGCTTCGCCGACGCTGGAAGAAGCCGCGCGCACCCTTGGCGCGAGCAGCttt                                                                       <  6:269394/149‑1 (MQ=255)
                                                                             cGGCGCTGTACCAGACCACGGCCTTGCTGCTGTTGGCCTATGCGCTGCTGTTCCTGCCGTTGGCGCAACCACCGGTGCGCACCGCGCTGAACAAGGCTTCGCCGACGCTGGAAGAAGCCGCGCGCACCCTTGGCGCGAGCAGCTTTGcc                                                                    >  5:278055/1‑149 (MQ=255)
                                                                             cGGCGCTGTACCAGACCACGGCCTTGCTGCTGTTGGCCTATGCGCTGCTGTTCCTGCCGTTGGCGCAACCACCGGTGCGCACCGCGCTGAACAAGGCTTCGCCGACGCTGGAAGAAGCCGCGCGCACCCTTGGCGCGAGCAGCTTTGcc                                                                    >  4:106464/1‑149 (MQ=255)
                                                                             cGGCGCTGTACCAGACCACGGCCTTGCTGCTGTTGGCCTATGCGCTGCTGTGCCTGCCGTTGGCGCAATCACCGGTGCGCACCGCGCTGAACAAGGCTTCGCCGACGCTGGAAGAAGCCGCGCGCACCCTTGGCGCGAGCAGCTTTGcc                                                                    <  3:14499/149‑1 (MQ=255)
                                                                              ggCGCTGTACCAGACCACGGCCTTGCTGCTGTTGGCCTATGCGCTGCTGTTCCTGCCGTTGGCGCAACCACCGGTGCGCACCGCGCTGAACAAGGCTTCGCCGACGCTGGAAGAAGCCGCGCGCACCCTTGGCGCGagc                                                                             >  2:61754/1‑139 (MQ=255)
                                                                                   tGTACCAGACCACGGCCTTGCTGCTGTTGGCCTATGCGCTGCTGTTCCTGCCGTTGGCGCAACCACCGGTGCGCACCGCGCTGAACAAGGCTTCGCCGACGCTGGAAGAAGCCGCGCGCACCCTTGGCGCGAGCAGCTTTGCCGCGTTc                                                              >  2:202418/1‑149 (MQ=255)
                                                                                    gTACCAGACCACGGCCTTGCTGCTGTTGGCCTATGCGCTGCTGTTCCTGCCGTTGGCGCAATCACCGGTGCGCACCGCGCTGAACAAGGCTTCGCCGACGCTGGAAGAAGCCGCGCGCACCCTTGGCGCGAGCAGCTTTGCCGCGTTCt                                                             >  8:66838/1‑149 (MQ=255)
                                                                                       ccAGACCACGGCCTTGCTGCTGTTGGCCTATGCGCTGCTGTTCCTGCCGTTGGCGCAATCACCGGTGCGCACCGCGCTGAACAAGGCTTCGCCGACGCTGGAAGAAGCCGCGCGCACCCTTGGCGCGAGCAGCTTTGCCGCGTTCTGc                                                           <  6:46194/148‑1 (MQ=255)
                                                                                        cAGACCACGGCCTTGCTGCTGTTGGCCTATGCGCTGCTGTTCCTGCCGTTGGCGCAATCACCGGTGCGCACCGCGCTGAACAAGGCTTCGCCGACGCTGGAAGAAGCCGCGCGCACCCTTGGCGCGAGCAGCTTTGCCGCGTTCTGCCg                                                         >  2:119684/1‑149 (MQ=255)
                                                                                        cAGACCACGGCCTTGCTGCTGTTGGCCTATGCGCTGCTGTTCCTGCCGTTGGCGCAATCACCGGTGCGCACCGCGCTGAACAAGGCTTCGCCGACGCTGGAAGAAGCCGCGCGCACCCTTGGCGCGAGCAGCTTTGCCGCGTTCTGCCg                                                         >  8:37375/1‑149 (MQ=255)
                                                                                          gACCACGGCCTTGCTGCTGTTGGCCTATGCGCTGCTGTTCCTGCCGTTGGCGCAACCACCGGTGCGCACCGCGCTGAACAAGGCTTCGCCGACGCTGGAAGAAGCCGCGCGCACCCTTGGCGCGAGCAGCTTTGCCGCGTTCTGCCggg                                                       >  6:346801/1‑149 (MQ=255)
                                                                                           aCCACGGCCTTGCTGCTGTTGGCCTATGCGCTGCTGTTCCTGCCGTTGGCGCAACCACCGGTGCGCACCGCGCTGAACAAGGCTTCGCCGACGCTGGAAGAAGCCGCGCGCACCCTTGGCGCGAGCAGCTTTGCCGCGTTCTGCCGGGt                                                      <  8:146969/149‑1 (MQ=255)
                                                                                             cACGGCCTTGCTGCTGTTGGCCTATGCGCTGCTGTTCCTGCCGTTGGCGCAATCACCGGTGCGCACCGCGCTGAACAAGGCTTCGCCGACGCTGGAAGAAGCCGCGCGCACCCTTGGCGCGAGCAGCTTTGCCGCGTTCTGCCGGGTGa                                                    >  3:188010/1‑149 (MQ=255)
                                                                                               cGGCCTTGCTGCTGTTGGCCTATGCGCTGCTGTTCCTGCCGTTGGCGCAATCACCGGTGCGCACCGCGCTGAACAAGGCTTCGCCGACGCTGGAAGAAGCCGCGCGCACCCTTGGCGCGAGCAGCTTTGCCGCGTTCTGCCGGGTGAcc                                                  >  7:202925/1‑149 (MQ=255)
                                                                                                 gCCTTGCTGCTGTTGGCCTATGCGCTGCTGTTCCTGCCGTTGGCGCAACCACCGGTGCGCACCGCGCTGAACAAGGCTTCGCCGACGCTGGAAGAAGCCGCGCGCACCCTTGGCGCGAGCAGCTTTGCCGCGTTCTGCCGGGTGACCtt                                                <  5:346801/149‑1 (MQ=255)
                                                                                                    ttGCTGCTGTTGGCCTATGCGCTGCTGTTCCTGCCGTTGGCGCAATCACCGGTGCGCACCGCGCTGAACAAGGCTTCGCCGACGCTGGAAGAAGCCGCGCGCACCCTTGGCGCGAGCAGCTTTGCCGCGTTCTGCCGGGTGACCTTGcc                                             >  5:100678/1‑149 (MQ=255)
                                                                                                       ctgctgTTGGCCTATGCGCTGCTGTTCCTGCCGTTGGCGCAATCACCGGTGCGCACCGCGCTGAACAAGGCTTCGCCGACGCTGGAAGAAGCCGCGCGCACCCTTGGCGCGAGCAGCTTTGCCGCGTTCTGCCGGGTGACCTTGCCGat                                          >  7:42635/1‑149 (MQ=255)
                                                                                                       ctgctgTTGGCCTATGCGCTGCTGTTCCTGCCGTTGGCGCAATCACCGGTGCGCACCGCGCTGAACAAGGCTTCGCCGACGCTGGAAGAAGCCGCGCGCACCCTTGGCGCGAGCAGCTTTGCCGCGTTCTGCCGGGTGACCTTGCCGa                                           <  4:116013/148‑1 (MQ=255)
                                                                                                          ctgTTGGCCTATGCGCTGCTGTTCCTGCCGTTGGCGCAATCACCGGTGCGCACCGCGCTGAACAAGGCTTCGCCGACGCTGGAAGAAGCCGCGCGCACCCTTGGCGCGAGCAGCTTTGCCGCGTTCTGCCGGGTGACCTTGCCGatc                                         <  6:330675/147‑1 (MQ=255)
                                                                                                          ctgTTGGCCTATGCGCTGCTGTTCCTGCCGTTGGCGCAATCACCGGTGCGCACCGCGCTGAACAAGGCTTCGCCGACGCTGGAAGAAGCCGCGCGCACCCTTGGCGCGAGCAGCTTTGCCGCGTTCTGCCGGGTGACCTTGCCGatc                                         <  8:42635/147‑1 (MQ=255)
                                                                                                          ctgTTGGCCTATGCGCTGCTGTTCCTGCCGTTGGCGCAATCACCGGTGCGCACCGCGCTGAACAAGGCTTCGCCGACGCTGGAAGAAGCCGCGCGCACCCTTGGCGCGAGCAGCTTTGCCGCGTTCTGCCGGGTGACCTTGCCGatc                                         <  8:202925/147‑1 (MQ=255)
                                                                                                             ttGGCCTATGCGCTGCTGTTCCTGCCGTTGGCGCAACCACCGGTGCGCACCGCGCTGAACAAGGCTTCGCCGACGCTGGAAGAAGCCGCGCGCACCCTTGGCGCGAGCAGCTTTGCCGCGTTCTGCCGGGTGACCTTGCCGatcatc                                      <  1:49457/147‑1 (MQ=255)
                                                                                                             ttGGCCTATGCGCTGCTGTTCCTGCCGTTGGCGCAACCACCGGTGCGCACCGCGCTGAACAAGGCTTCGCCGACGCTGGAAGAAGCCGCGCGCACCCTTGGCGCGAGCAGCTTTGCCGCGTTCTGCCGGGTGACCTTGCCGATCATCtt                                    >  8:35947/1‑149 (MQ=255)
                                                                                                               ggCCTATGCGCTGCTGTTCCTGCCGTTGGCGCAATCACCGGTGCGCACCGCGCTGAACAAGGCTTCGCCGACGCTGGAAGAAGCCGCACGCACCCTTGGCGCGAGCAGCTTTGCCGCGTTCTGCCGGGTGACCTTGCCGATCATCTTc                                   <  2:176837/148‑1 (MQ=255)
                                                                                                                 ccTATGCGCTGCTGTTCCTGCCGTTGGCGCAATCACCGGTGCGCACCGCGCTGAACAAGGCTTCGCCGACGCTGGAAGAAGCCGCGCGCACCCTTGGCGCGAGCAGCTTTGCCGCGTTCTGCCGGGTGACCTTGCCGATCATCTTCCCg                                <  1:274903/149‑1 (MQ=255)
                                                                                                                  cTATGCGCTGCTGTTCCTGCCGTTGGCGCAACCACCGGTGCGCACCGCGCTGAACAAGGCTTCGCCGACGCTGGAAGAAGCCGCGCGCACCCTTGGCGCGAGCAGCTTTGCCGCGTTCTGCCGGGTGACCTTGCCGATCATCTTCCCgg                               >  6:104491/1‑149 (MQ=255)
                                                                                                                  cTATGCGCTGCTGTTCCTGCCGTTCGCGCAATCACCGGTGCGCACCGCGCTGAACAAGGCTTCGCCGACGCTGGAAGAAGCCGCGCGCACCCTTGGCGCGAGCAGCATTGGCGCGTTCTGCCGGGTGACCTTGCCGATCTTCTTCCCgg                               >  8:304683/1‑149 (MQ=255)
                                                                                                                      gcgcTGCTGTTCCTGCCGTTGGCGCAATCACCGGTGCGCACCGCGCTGAACAAGGCTTCGCCGACGCTGGAAGAAGCCGCGCGCACCCTTGGCGCGAGCAGCTTTGCCGCGTTCTGCCGGGTGACCTTGCCGATCATCTTCCCGGCCa                            <  7:265015/148‑1 (MQ=255)
                                                                                                                       cgcTGCTGTTCCTGCCTTTGGCGCAATCACCGGTGCGCACCGCGCTGAACAAGGCTTCGCCGACGCTGGAAGAAGCCGCGCGCACCCTTGGCGCGAGCAGCTTTGCCGCGTTCTGCCGGGTGACCTTGCCGATCATCTTCCCGGCCAt                           <  1:119684/148‑1 (MQ=255)
                                                                                                                        gctgctGTTCCTGCCGTTGGCGCAATCACCGGTGCGCACCGCGCTGAACAAGGCTTCGCCGACGCTGGAAGAAGCCGCGCGCACCCTTGGCGCGAGCAGCTTTGCCGCGTTCTGCCGGGTGACCTTGCCGATCATCTTCCCGGCCATgg                         >  4:165654/1‑149 (MQ=255)
                                                                                                                         ctgctgTTCCTGCCGTTGGCGCAATCACCGGTGCGCACCGCGCTGAACAAGGCTTCGCCGACGCTGGAAGAAGCCGCGCGCACCCTTGGCGCGAGCAGCTTTGCCGCGTTCTGCCGGGTGACCTTGCCGATCATCTTCCCGGCCATggc                        <  6:35322/149‑1 (MQ=255)
                                                                                                                                  cTGCCGTTGGCGCAATCACCGGTGCGCACCGCGCTGAACAAGGCTTCGCCGACGCTGGAAGAAGCCGCGCGCACCCTTGGCGCGAGCAGCTTTGCCGCGTTCTGCCGGGTGACCTTGCCGATCATCTTCCCGGCCATGGCGGCAGCCtt               >  7:270074/1‑149 (MQ=255)
                                                                                                                                  cTGCCGTTGGCGCAATCACCGGTGCGCACCGCGCTGAACAAGGCTTCGCCGACGCTGGAAGAAGCCGCGCGCACCCTTGGCGCGAGCAGCTTTGCCGCGTTCTGCCGGGTGACCTTGCCGATCATCTTCCCGGCCATGGCGGCAGCCtt               >  4:33758/1‑149 (MQ=255)
                                                                                                                                  cTGCCGTTGGCGCAATCACCGGTGCGCACCGCGCTGAACAAGGCTTCGCCGACGCTGGAAGAAGCCGCGCGCACCCTTGGCGCGAGCAGCTTTGCCGCGTTCTGCCGGGTGACCTTGCCGATCATCTTCCCGGCCATGGCGGCAGCCtt               >  4:210183/1‑149 (MQ=255)
                                                                                                                                   tGCCGTTGGCGCAATCACCGGTGCGCACCGCGCTGAACAAGGCTTCGCCGACGCTGGAAGAAGCCGCGCGCACCCTTGGCGCGAGCAGCTTTGCCGCGTTCTGCCGGGTGACCTTGCCGATCATCTTCCCGGCCATGGCGGCAGCCttt              >  4:79899/1‑149 (MQ=255)
                                                                                                                                   tGCCGTTGGCGCAATCACCGGTGCGCACCGCGCTGAACAAGGCTTCGCCGACGCTGGAAGAAGCCGCGCGCACCCTTGGCGCGAGCAGCTTTGCCGCGTTCTGCCGGGTGACCTTGCCGATCATCTTCCCGGCCATGGCGGCAGCCttt              >  5:259461/1‑149 (MQ=255)
                                                                                                                                     ccGTTGGCGCAACCACCGGTGCGCACCGCGCTGAACAAGGCTTCGCCGACGCTGGAAGAAGCCGCGCGCACCCTTGGCGCGAGCAGCTTTGCCGCGTTCTGCCGGGTGACCTTGCCGATCATCTTCCCGGCCATGGCGGCAGCCttt              <  5:104491/147‑1 (MQ=255)
                                                                                                                                       gTTGGCGCAATCACCGGTGCGCACCGCGCTGAACAAGGCTTCGCCGACGCTGGAAGAAGCCGCGCGCACCCTTGGCGCGAGCAGCTTTGCCGCGTTCTGCCGGGTGACCTTGCCGATCATCTTCCCGGCCATGGCGGCAGCCTTTgcgc          <  6:100678/149‑1 (MQ=255)
                                                                                                                                          ggCGCAATCACCGGTGCGCACCGCGCAGAACAAGGCTTCGCCGACGCTGGAAGAAGCCGCGCGCACCCTTGGCGCGAGCAGCTTTGCCGCGTTCTGCCGGGTGACCTTGCCGATCATCTTCCCGGCCATGGCGGCAGCCTTTGCGCTg        <  6:27815/148‑1 (MQ=255)
                                                                                                                                           gcgcAATCACCGGTGCGCACCGCGCTGAACAAGGCTTCGCCGACGCTGGAAGAAGCCGCGCGCACCCTTGGCGCGAGCAGCTTTGCCGCGTTCTGCCGGGTGACCTTGCCGATCATCTTCCCGGCCATGGCGGCAGCCTTTGCGCTgg       <  7:66838/148‑1 (MQ=255)
                                                                                                                                             gcAATCACCGGTGCGCACCGCGCTGAACAAGGCTTCGCCGACGCTGGAAGAAGCCGCGCGCACCCTTGGCGCGAGCAGCTTTGCCGCGTTCTGCCGGGTGACCTTGCCGATCATCTTCCCGGCCATGGCGGCAGCCTTTGCGCTGgtgt    >  2:74945/1‑149 (MQ=255)
                                                                                                                                               aaTCACCGGTGCGCACCGAGCTGAACAAGGCTTCGCCGACGCTGGAAGAAGCCGCGCGCACCCTTGGAGCGAGCAGCTTTGCCGCGTTCTGCCGGGTGACCTTGCCGATCATCTTCCCGGCCATGGCGGCAGCCTTTGCGCTGGTGTTc  >  7:56395/1‑149 (MQ=255)
                                                                                                                                                 |                                                                                                                                                  
TCTGCCCTATCTGCTGCACGCCCTGCCCGGCCTGGTGATTGCCCTGACGTTGGTGTTCTTCGCCCTGCACTACGTGCCGGCGCTGTACCAGACCACGGCCTTGCTGCTGTTGGCCTATGCGCTGCTGTTCCTGCCGTTGGCGCAATCACCGGTGCGCACCGCGCTGAACAAGGCTTCGCCGACGCTGGAAGAAGCCGCGCGCACCCTTGGCGCGAGCAGCTTTGCCGCGTTCTGCCGGGTGACCTTGCCGATCATCTTCCCGGCCATGGCGGCAGCCTTTGCGCTGGTGTTC  >  NC_002947/5551851‑5552142

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