Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 2,776,107 A→C 7.8% K152Q (AAG→CAG)  sotB → sugar efflux transporter

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029472,776,1070AC 7.8% 214.7 / 7.5 75K152Q (AAG→CAG) sotBsugar efflux transporter
Reads supporting (aligned to +/- strand):  ref base A (30/39);  new base C (2/4);  total (32/43)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.92e-01

TTCATTGTCAGCCACTTGCTGTCCTGGCTTTCGCAAAGCTTCGCCATGCTGCTGGTCAGCCGCATCGGCATTGCCCTGGCGCATGCGGTGTTCTGGTCCATCACCGCCTCGCTGGCGGTGCGTGTGGCGCCGCCGGGCCAGCAGGCCAAGGCGCTGGGCCTGCTGGCCACCGGTACTACCCTGGCCATGGTCCTGGGCATCCCGCTGGGCCGTGTGGTGGGCGAGGCACTGGGCTGGCGTATCACGTTCCTGAGCATTGCCGGGGTAGCGCTAGCCACCATGCTTTGCCTGAT  >  NC_002947/2775960‑2776252
                                                                                                                                                   |                                                                                                                                                 
ttCATTGTCAGCCACTTGCTGTCCTGGCTTTCGCAAAGCTTCGCCATGCTGCTGGTCAGCCGCATCGGCATTGCCCTGGCGCATGCGGTGTTCTGGTCCATCACCGCCTCGCTGGCGGTGCGTGTGGCGCCGCCGGGCCAGCAGGCCa                                                                                                                                                   <  8:26395/148‑1 (MQ=255)
   aTTGTCAGCCACTTGCTGTCCTGGCTTTCGCAAAGCTTCGCCATGCTGCTGGTCAGCCGCATCGGCATTGCCCTGGCGCATGCGGTGTTCTGGTCCATCACCGCCTCGCTGGCGGTGCGTGTGGCGCCGCCGGGCCAGCAGGCCAAggg                                                                                                                                               >  3:177174/1‑148 (MQ=255)
   aTTGTCAGCCACTTGCTGTCCTGGCTTTCGCAAAGCTTCGCCATGCTGCTGGTCAGCCGCATCGGCATTGCCCTGGCGCATGCGGTGTTCTGGTCCATCACCGCCTCGCTGGCGGTGCGTGTGGCGCCGCCGGGCCAGCAGGCCAAgg                                                                                                                                                <  7:97256/148‑1 (MQ=255)
   aTTGTCAGCCACTTGCTGTCCTGGCTTTCGCAAAGCTTCGCCATGCTGCTGGTCAGCCGCATCGGCATTGCCCTGGCGCATGCGGTGTTCTGGTCCATCACCGCCTCGCTGGCGGTGCGTGTGGCGCCGCCGGGCCAGCAGGCCAAgg                                                                                                                                                <  4:236935/148‑1 (MQ=255)
       tCAGCCACTTGCTGTCCTGGCTTTCGCAAAGCTTCGCCATGCTGCTGGTCAGCCGCATCGGCATTGCCCTGGCGCATGCGGTGTTCTGGTCCATCACCGCCTCGCTGGCGGTGCGTGTGGCGCCGCCGGGCCAGCAGGCCAAGGCGCt                                                                                                                                            <  6:278576/148‑1 (MQ=255)
         aGCCACTTGCTGTCCTGGCTTTCGCAAAGCTTCGCCATGCTGCTGGTCAGCCGCATCGGCATTGCCCTGGCGCATGCGGTGTTCTGGTCCATCACCGCCTCGCTGGCGGTGCGTGTGGCGCCGCCGGGCCAGCAGGCCAAGGCGCTgg                                                                                                                                          <  4:201010/148‑1 (MQ=255)
         aGCCACTTGCTGTCCTGGCTTTCGCAAAGCTTCGCCATGCTGCTGGTCAGCCGCATCGGCATTGCCCTGGCGCATGCGGTGTTCTGGTCCATCACCGCCTCGCTGGCGGTGCGTGTGGCGCCGCCGGGCCAGCAGGCCAAGGCGCTgg                                                                                                                                          <  4:31161/148‑1 (MQ=255)
          gCCACTTGCTGTCCTGGCTTTCGCAAAGCTTCGCCATGCTCCTGGTCAGCCGCATCGGCATTGCCCTGGCGCATGCGGTGTTCTGGTCCATCACCGCCTCGCTGGCGGTGCGTGTGGCGCCGCCGGGCCAGCAGGCCAAGGCGCTGGGc                                                                                                                                        >  3:42128/1‑149 (MQ=255)
            cACTTGCTGTCCTGGCTTTCGCAAAGCTTCGCCATGCTGCTGGTCAGCCGCATCGGCATTGCCCTGGCGCATGCGGTGTTCTGGTCCATCACCGCCTCGCTGGCGGTGCGTGTGGCGCCGCCGGGCCAGCAGGCCAAGGCGCTGGGCct                                                                                                                                      >  6:377024/1‑149 (MQ=255)
             aCTTGCTGTCCTGGCTTTCGCAAAGCTTCGCCATGCTGCTGGTCAGCCGCATCGGCATTGCCCTGGCGCATGCGGTGTTCTGGTCCATCACCGCCTCGCTGGCGGTGCGTGTGGCGCCGCCGGGCCAGCAGGCCAAGGCGCTGGGCctg                                                                                                                                     <  7:40819/149‑1 (MQ=255)
                 gCTGTCCTGGCTTTCGCAAAGCTTCGCCATGCTGCTGGTCAGCCGCATCGGCATTGCCCTGGCGCATGCGGTGTTCTGGTCCATCACCGCCTCGCTGGCGGTGCGTGTGGCGCCGCCGGGCCAGCAGGCCCAGCGCCTTGGCctgctg                                                                                                                                  <  1:419803/148‑1 (MQ=255)
                 gCTGTCCTGGCTTTCGCAAAGCTTCGCCATGCTGCTGGTCAGCCGCATCGGCATTGCCCTGGCGCATGCGGTGTTCTGGTCCATCACCGCCTCGCTGGCGGTGCGTGTGGCGCCGCCGGGCCAGCAGGCCCAGCGCCTTGGCctgct                                                                                                                                   <  8:201774/147‑1 (MQ=255)
                 gCTGTCCTGGCTTTCGCAAAGCTTCGCCATGCTGCTGGTCAGCCGCATCGGCATTGCCCTGGCGCATGCGGTGTTCTGGTCCATCACCGCCTCGCTGGCGGTGCGTGTGGCGCCGCCGGGCCAGCAGGCCAAGGCGCTGGGCctgctg                                                                                                                                  <  5:377024/148‑1 (MQ=255)
                 gCTGTCCTGGCTTTCGCAAAGCTTCGCCATGCTGCTGGTCAGCCGCATCGGCATTGCCCTGGCGCATGCGGTGTTCTGGTCCATCACCGCCTCGCTGGCGGTGCGTGTGGCGCCGCCGGGCCAGCAGGCCAAGGCGCTGGGCctgct                                                                                                                                   <  1:158951/147‑1 (MQ=255)
                 gCTGTCCTGGCTTTCGCAAAGCTTCGCCATGCTGCTGGTCAGCCGCATCGGCATTGCCCTGGCGCATGCGGTGTTCTGGTCCATCACCGCCTCGCTGGCGGTGCGTGTGGCGCCGCCGGGCCAGCAGGCCAAGGCGCTGGGCctgct                                                                                                                                   <  4:42128/147‑1 (MQ=255)
                 gCTGTCCTGGCTTTCGCAAAGCTTCGCCATGCTGCTGGTCAGCCGCATCGGCATTGCCCTGGCGCATGCGGTGTTCTGGTCCATCACCGCCTCGCTGGCGGTGCGTGTGGCGCCGCCGGGCCAGCAGGCCAAGGCGCTGAGCcttctg                                                                                                                                  <  6:252779/148‑1 (MQ=255)
                  cTGTCCTGGCTTTCGCAAAGCTTCGCCATGCTGCTGGTCAGCCGCATCGGCATTGCCCTGGCGCATGCGGTGTTCTGGTCCATCACCGCCTCGCTGGCGGTGCGTGTGGCGCCGCCGGGCCAGCAGGCCCAGCGCCTTGGCCTGCTgg                                                                                                                                 <  5:384305/148‑1 (MQ=255)
                  cTGTCCTGGCTTTCGCAAAGCTTCGCCATGCTGCTGGTCAGCCGCATCGGCATTGCCCTGGCGCATGCGGTGTTCTGGTCCATCACCGCCTCGCTGGCGGTGCGTGTGGCGCCGCCGGGCCAGCAGGCCCAGCGCCTTGGCCTGCTgg                                                                                                                                 <  8:222392/148‑1 (MQ=255)
                    gTCCTGGCTTTCGCAAAGCTTCGCCATGCTGCTGGTCAGCCGCATCGGCATTGCCCTGGCGCATGCGGTGTTCTGGTCCATCACCGCCTCGCTGGCGGTGCGTGTGGCGCCGCCGGGCCAGCAGGCCAAGGCGCTGGGCCTGCTGGCCa                                                                                                                              >  7:156471/1‑149 (MQ=255)
                     tCCTGGCTTTCGCAAAGCTTCGCCATGCTGCTGGTCAGCCGCATCGGCATTGCCCTGGCGCATGCGGTGTTCTGGTCCATCACCGCCTCGCTGGCGGTGCGTGTGGCGCCGCCGGGCCAGCAGGCCAAGGCGCTGGGCCTGCTGGCCAc                                                                                                                             >  8:56657/1‑149 (MQ=255)
                     tCCTGGCTTTCGCAAAGCTTCGCCATGCTGCTGGTCAGCCGCATCGGCATTGCCCTGGCGCATGCGGTGTTCTGGTCCATCACCGCCTCGCTGGCGGTGCGTGTGGCGCCGCCGGGCCAGCAGGCCAAGGCGCTGGGCCTGCTGGCCAc                                                                                                                             <  1:150636/149‑1 (MQ=255)
                       cTGGCTTTCGCAAAGCTTCGCCATGCTGCTGGTCAGCCGCATCGGCATTGCCCTGGCGCATGCGGTGTTCTGGTCCATCACCGCCTCGCTGGCGGTGCGTGTGGCGCCGCCGGGCCAGCAGGCCAAGGCGCTGGGCCTGCTGGCCACCg                                                                                                                           >  5:82768/1‑149 (MQ=255)
                       cTGGCTTTCGCAAAGCTTCGCCATGCTGCTGGTCAGCCGCATCGGCATTGCCCTGGCGCATGCGGTGTTCTGGTCCATCACCGCCTCGCTGGCGGTGCGTGTGGCGCCGCCGGGCCAGCAGGCCAAGGCGCTGGGCCTGCTGGCCACCg                                                                                                                           >  6:407069/1‑149 (MQ=255)
                        tGGCTTTCGCAAAGCTTCGCCATGCTGCTGGTCAGCCGCATCGGCATTGCCCTGGCGCATGCGGTGTTCTGGTCCATCACCGCCTCGCTGGCGGTGCGTGTGGCGCCGCCGGGCCAGCAGGCCAAGGCGCTGGGCCTGCTGGCCACCgg                                                                                                                          >  7:197664/1‑149 (MQ=255)
                           cTTTCGCAAAGCTTCGCCATGCTGCTGGTCAGCCGCATCGGCATTGCCCTGGCGCATGCGGTGTTCTGGTCCATCACCGCCTCGCTGGCGGTGCGTGTGGCGCCGCCGGGCCAGCAGGCCAAGGCGCTGGGCCTGCTGGCCACCGGtac                                                                                                                       >  2:37249/1‑149 (MQ=255)
                            tttCGCAAAGCTTCGCCATGCTGCTGGTCAGCCGCATCGGCATTGCCCTGGCGCATGCGGTGTTCTGGTCCATCACCGCCTCGCTGGCGGTGCGTGTGGCGCCGCCGGGCCAGCAGGCCAAGGCGCTGGGCCTGCTGGCCACCgg                                                                                                                          >  5:159580/1‑145 (MQ=255)
                            tttCGCAAAGCTTCGCCATGCTGCTGGTCAGCCGCATCGGCATTGCCCTGGCGCATGCGGTGTTCTGGTCCATCACCGCCTCGCTGGCGGTGCGTGTGGCGCCGCCGGGCCAGCAGGCCAAGGCGCTGGGCCTGCTGGCCACCgg                                                                                                                          <  6:159580/145‑1 (MQ=255)
                              tCGCAAAGCTTCGCCATGCTGCTGGTCAGCCGCATCGGCATTGCCCTGGCGCATGCGGTGTTCTGGTCCATCACCGCCTCGCTGGCGGTGCGTGTGGCGCCGCCGGGCCAGCAGGCCAAGGCGCTGGGCCTGCTGGCCACCGGtacta                                                                                                                     <  3:337029/148‑1 (MQ=255)
                                 cAAAGCTTCGCCATGCTGCTGGTCAGCCGCATCGGCATTGCCCTGGCGCATGCGGTGTTCTGGTCCATCACCGCCTCGCTGGCGGTGCGTGTGGCGCCGCCGGGCCAGCAGGCCAAGGCGCTGGGCCTGCTGGCCACCGGTACTACCCt                                                                                                                 >  5:361258/1‑149 (MQ=255)
                                     gCTTCGCCATGCTGCTGGTCAGCCGCATCGGCATTGCCCTGGCGCATGCGGTGTTCTGGTCCATCACCGCCTCGCTGGCGGTGCGTGTGGCGCCGCCGGGCCAGCAGGCCAAGGCGCTGGGCCTGCTGGCCACCGGTACTACCCTGGcc                                                                                                             <  8:156471/149‑1 (MQ=255)
                                      cTTCGCCATGCTGCTGGTCAGCCGCATCGGCATTGCCCTGGCGCATGCGGTGTTCTGGTCCATCACCGCCTCGCTGGCGGTGCGTGTGGCGCCGCCGGGCCAGCAGGCCAAGGCGCTGGGCCTGCTGGCCACCGGTACTACCCTGGCCa                                                                                                            <  3:255969/149‑1 (MQ=255)
                                         cGCCATGCTGCTGGTCAGCCGCATCGGCATTGCCCTGGCGCATGCGGTGTTCTGGTCCATCACCGCCTCGCTGGCGGTGCGTGTGGCGCCGCCGGGCCAGCAGGCCAAGGCGCTGGGCCTGCTGGCCACCGGTACTACCCTGGc                                                                                                              <  4:47439/144‑1 (MQ=255)
                                         cGCCATGCTGCTGGTCAGCCGCATCGGCATTGCCCTGGCGCATGCGGTGTTCTGGTCCATCACCGCCTCGCTGGCGGTGCGTGTGGCGCCGCCGGGCCAGCAGGCCAAGGCGCTGGGCCTGCTGGCCACCGGTACTACCCTGGCCagat                                                                                                         >  3:47439/1‑146 (MQ=255)
                                         cGCCATGCTGCTGGTCAGCCGCATCGGCATTGCCCTGGCGCATGCGGTGTTCTAGTCCATCACCGCCTCGCTGGCGGTGCGTGTGGCGCCGCCGGGCCAGCAGGCCAAGGCGCTGGGCCTGCTGGCCACCg                                                                                                                           <  3:21682/131‑1 (MQ=255)
                                            cATGCTGCTGGTCAGCCGCATCGGCATTGCCCTGGCGCATGCGGTGTTCTGGTCCATCACCGCCTCGCTGGCGGTGCGTGTGGCGCCGCCGGGCCAGCAGGCCAAGGCGCTGGGCCTGCTGGCCACCGGTACTACCCTGGCCAGGGTcc                                                                                                      >  6:148564/1‑149 (MQ=255)
                                             aTGCTGCTGGTCAGCCGCATCGGCATTGCCCTGGCGCATGCGGTGTTCTGGTCCATCACCGCCTCGCTGGCGGTGCGTGTGGCGCCGCCGGGCCAGCAGGCCAAGGCGCTGGGCCTGCTGGCCACCGGTACTACCCTGGCCATGGTCCt                                                                                                     <  8:112003/149‑1 (MQ=255)
                                             aTGCTGCTGGTCAGCCGCATCGGCATTGCCCTGGCGCATGCGGTGTTCTGGTCCATCACCGCCTCGCTGGCGGTGCGTGTGGCGCCGCCGGGCCACCAGGCCAACGCGCTGGGGCAGCTGGGCCCCGGTACTACCCTGGCCATGGtgat                                                                                                     >  8:147367/1‑146 (MQ=255)
                                                ctgctgGTCAGCCGCATCGGCATTGCCCTGGCGCATGCGGTGTTCTGGTCCATCACCGCCTCGCTGGCGGTGCGTGTGGCGCCGCCGGGCCAGCAGGCCAAGGCGCTGGGCCTGCTGGCCACCGGTACTACCCTGGCCATGGTCCTgg                                                                                                   <  6:361258/148‑1 (MQ=255)
                                                ctgctgGTCAGCCGCATCGGCATTGCCCTGGCGCATGCGGTGTTCTGGTCCATCACCGCCTCGCTGGCGGTGCGTGTGGCGCCGCCGGGCCAGCAGGCCAAGGCGCTGGGCCTGCTGGCCACCGGTACTACCCTGGCCATGGTCCTgg                                                                                                   <  3:372983/148‑1 (MQ=255)
                                                 tgctgGTCAGCCGCATCGGCATTGCCCTGGCGCATGCGGTGTTCTGGTCAATCACCGCCTCGCTGGAGGTGCGTGTGGCGCCGCCGGGCCAGCAGGCACAGGCGCTGGGCCTGCTGGCCACCGGTAATACCCTGGCCATGGTCCTggg                                                                                                  >  3:85682/1‑148 (MQ=255)
                                                  gctgGTCAGCCGCATCGGCATTGCCCTGGCGCATGCGGTGTTCTGGTCCATCACCGCCTCGCTGGCGGTGCGTGTGGCGCCGCCGGGCCAGCAGGCCAAGGCGCTGGGCCTGCTGGCCa                                                                                                                              <  1:12996/119‑1 (MQ=255)
                                                  gctgGTCAGCCGCATCGGCATTGCCCTGGCGCATGCGGTGTTCTGGTCCATCACCGCCTCGCTGGCGGTGCGTGTGGCGCCGCCGGGCCAGCAGGCCAAGGCGCTGGGCCTGCTGGCCa                                                                                                                              >  2:12996/1‑119 (MQ=255)
                                                      gTCAGCCGCATCGGCATTGCCCTGGCGCATGCGGTGTTCTGGTCCATCACCGCCTCGCTGGCGGTGCGTGTGGCGCCGCCGGGCCAGCAGGCCAAGGCGCTGGGCCTGCTGGCCACCGGTACTACCCTGGCCATGGTCCTGGGCATcc                                                                                             <  7:56657/148‑1 (MQ=255)
                                                           ccGCATCGGCATTGCCCTGGCGCATGCGGTGTTCTGGTCCATCACCGCCTCGCTGGCGGTGCGTGTGGCGCCGCCGGGCCAGCAGGCCAAGGCGCTGGGCCTGCTGGCCACCGGTACTACCCTGGCCATGGTCCTGGGCATCCCGCTgg                                                                                       <  6:371565/149‑1 (MQ=255)
                                                            cGCATCGGCATTGCCCTGGCGCCTGCGGTGTTCTGGTCCATCACCGCCTCGCTGGCGGTGCGTGTGGCGCCGCCGGGCCAGCAGGCCAAGGCGCTGGGCCTGCTGGCCACCGGTACTACCCTGGCCATGGGCCTGGGCATCCCGCTggg                                                                                      <  3:242002/149‑1 (MQ=255)
                                                                    cATTGCCCTGGCGCATGCGGTGTTCTGGTCCATCACCGCCTCGCTGGCGGTGCGTGTGGCGCCGCCGGGCCAGCAGGCCAAGGCGCTGGGCCTGCTGGCCACCGGTACTACCCTGGCCATGGTCCTGGGCATCCCGCTGGGCCGTgtgg                                                                              <  5:3159/149‑1 (MQ=255)
                                                                        gCCCTGGGGGATGCGGTGTTCTGGTCCATCCCCGCCTCGCTGGCGGTGCGTGTGGCGCCGCCGGGCCAGCAGGCCAAGGCGCTGGGCCTGCTGGCCACCGGTACTACCCTGGCCATGGTCCTGGGCATCCCGCTGGGCCGTGtggtgg                                                                           <  2:393620/148‑1 (MQ=255)
                                                                         cccTGGCGCATGCGGTGTTCTGGTCCATCACCGCCTCGCTGGCGGTGCGTGTGGCGCCGCCGGGCCAGCAGGCCAAGGCGCTGGGCctgct                                                                                                                                   >  8:335710/1‑91 (MQ=255)
                                                                         cccTGGCGCATGCGGTGTTCTGGTCCATCACCGCCTCGCTGGCGGTGCGTGTGGCGCCGCCGGGCCAGCAGGCCAAGGCGCTGGGCctgct                                                                                                                                   <  7:335710/91‑1 (MQ=255)
                                                                           cTGGCGCATGCGGTGTTCTGGTCCATCACCGCCTCGCTGGCGGTGCGTGTGGCGCCGCCGGGCCAGCAGGCCAAGGCGCTGGGCCTGCTGGCCACCGGTACTACCCTGGCCATGGTCCTGGGCATCCCGCTGGGCCGTGTGGTGGGCGa                                                                       <  1:279369/149‑1 (MQ=255)
                                                                            tggggtaTGGGGTGTTTTGGTCTTTCACGGGGTCGGTGGCGGTGGGTGTGGGGCGGCCGGGCCAGCAGGCCAAGGCGCTGGGCCTGCTGGCCACCGGTACTACCCTGGCCATGGTCCTGGGCATCCCGCTGGGCCGTGTGGTGGGCGAg                                                                      <  2:433777/143‑1 (MQ=255)
                                                                                gcATGCGGTGTTCTGGTCCATCACCGCCTCGCTGGCGGTGCGTGTGGCGCCGCCGGGCCAGCAGGCCAAGGCGCTGGGCCTGCTGGCCACCGGTACTACCCTGGCCATGGTCCTGGGCATCCCGCTGGGCCGTGTGGTGGGCGAGGCAc                                                                  >  3:52461/1‑149 (MQ=255)
                                                                                   tGCGGTGTTCTGGTCCATCACCGCCTCGCTGGCGGTGCGTGTGGCGCCGCCGGGCCAGCAGGCCAAGGCGCTGGGCCTGCTGGCCACCGGTACTACCCTGGCCATGGTCCTGGGCATCCCGCTGGGCCGTGTGGTGGGCGAGGCACTgg                                                               <  5:407069/149‑1 (MQ=255)
                                                                                     cGGTGTTCTGGTCCATCACCGCCTCGCTGGCGGTGCGTGTGGCGCCGCCGGGCCAGCAGGCCAAGGCGCTGGGCCTGCTGGCCACCGGTACTACCCTGGCCATGGTCCTGGGCATCCCGCGGGGCCGTGTGGTGGGCGAGGCACTGGGc                                                             >  5:333834/1‑149 (MQ=255)
                                                                                                       ccGCCTCGCTGGCGGTGCGTGTGGCGCCGCCGGGCCAGCAGGCCAAGGCGCTGGGCCTGCTGGCCACCGGTACTACCCTGGCCATGGTCCTGGGCATCCCGCTGGGCCGTGTGGTGGGCGAGGCACTGGGCTGGCGTATCACGTTCCTg                                           >  7:122275/1‑149 (MQ=255)
                                                                                                        cGCCTCGCTGGCGGTGCGTGTGGCGCCGCCGGGCCAGCAGGCCAAGGCGCTGGGCCTGCTGGCCACCGGTACTACCCTGGCCATGGTCCTGGGCATCCCGCTGGGCCGTGTGGTGGGCGAGGCACTGGGCTGGCGTATCACGTTCCt                                            >  1:62640/1‑147 (MQ=255)
                                                                                                        cGCCTCGCTGGCGGTGCGTGTGGCGCCGCCGGGCCAGCAGGCCAAGGCGCTGGGCCTGCTGGCCACCGGTACTACCCTGGCCATGGTCCTGGGCATCCCGCTGGGCCGTGTGGTGGGCCAGGCACTGGGCTGGCGTATCACGTTCCTGa                                          <  1:37249/149‑1 (MQ=255)
                                                                                                         gCCTCGCTGGCGGTGCGTGTGGCGCCGCCGGGCCAGCAGGCCAAGGCGCTGGGCCTGCTGGCCACCGGTACTACCCTGGCCATGGTCCTGGGCATCCCGCTGGGCCGTGTGGTGGGCGAGGCACTGGGCTGGCGTATCACGTTCCTGAg                                         >  3:130067/1‑149 (MQ=255)
                                                                                                              gCTGGCGGTGCGTGTGGCGCCGCCGGGCCAGCAGGCCAAGGCGCTGGGCCTGCTGGCCACCGGTACTACCCAGGCCATGGTCCTGGGCATCCCGCTGGGCCGTGTGGTGGGCGCGGCACTGGGGTGGCGGATAACGGTGCTGAGCAttt                                    >  8:779/1‑148 (MQ=255)
                                                                                                                  tcGGTGCGTGTGGCGCCGCCGGGCCAGCAGGCCAAGGCGCGGGGCCTGCTGGCCACCGGTACTACCCTGGCCATGGTCCTGGGCATCCCGCTGGGCCGTGTGGTGGGCGAGGCACTGGGCTGGCGTATCACGTTCCTGAGCATTGCCg                                 <  1:290309/147‑1 (MQ=255)
                                                                                                                    ggTGCGTGTGGCGCCGCCGGGCCAGCAGGCCAAGGCGCTGGGCCTGCTGGCCCCCCGTACTACCCTGGCCATGGTACTGGGCATCCCGca                                                                                         >  8:34228/1‑89 (MQ=255)
                                                                                                                    ggTGCGTGTGGCGCCGCCGGGCCAGCAGGCCAAGGCGCTGGGCCTGCTGGCCACCGGTACTACCCTGGCCATGGTCCTGGGCATCCCGCt                                                                                         <  7:34228/90‑1 (MQ=255)
                                                                                                                        cGTGTGGCGCCGCCGGGCCAGCAGGCCAAGGCGCTGGGCCTGCTGGCCACCGGTACTACCCTGGCCATGGTCCTGGGCATCCCGCTGGGCCGTGTGGTGGGCGAGGCACTGGGCTGGCGTATCACGTTCCTGAGCATTGCCGGGGTAgc                          <  4:143769/149‑1 (MQ=255)
                                                                                                                          tgtgGCGCCGCCGGGCCAGCAGGCCAGGGCGCTGGGCCTGCTGGCCCCCGGTACTAACCTGGCCATGGTCCTGGGCATCCCGCTGGGCCGTGTGGTGGGCGAGGCACTGGGCTGGCGTATCACGTTCCTGAGCATTGCCGGGGTAgcgc                        >  1:12932/1‑149 (MQ=255)
                                                                                                                          tgtgGCGCCGCCGGGCCAGCAGGCCAAGGCGCTGGGCCTGCTGGCCACCGGTACTACCCTGGGCATGGTCCTGGGCATCCCGCTGGGCCATGTGGTGGGCGAGGCACTGGGCTGGCGTATCACGTTCCTGAGCATTGCCGGGGTAgcgc                        >  5:184477/1‑149 (MQ=255)
                                                                                                                          tgtgGCGCCGCCGGGCCAGCAGGCCAAGGCGCTGGGCCTGCTGGCCACCGGTACTACCCTGGCCATGGTCCTGGGCATCCCGCTGGGCCGTGTGGTGGGCGAGGCACTGGGCTGGCGTATCACGTTCCTGAGCATTGCCGGGGTAgcgc                        <  8:197664/149‑1 (MQ=255)
                                                                                                                              gcgcCGCCGGGCCAGCAGGCCAAGGCGCTGGGCCTGCTGGCCACCGGTACTACCCTGGCCATGGTCCTGGGCATCCCGCTGGGCCGTGTGGTGGGCGAGGCACTGGGCTGGCGTATCACGTTCCTGAGCATTGCCGGGGTAGCGCTAGc                    >  1:257108/1‑149 (MQ=255)
                                                                                                                                taccggtggCCAGCAGGCCCAGCGCCTTGGCCTGCTGGCCCCCGGTACTACCCTGGCCATGGTCCTGGGCATCCCGCTGGGCCGTGTGGTGGGCGAGGCACTGGGCTGGCGTATCACGTTCCTGAGCATTGCCGGGGTAGCGCTAGaca                  >  8:150773/8‑149 (MQ=38)
                                                                                                                                gccgccGGGCCAGCAGGCCAAGGCGCTGGGCCTGCTGGCCACCGGTACTACCCTGGCCATGGTCCTGGGCATCCCGCTGGGCCGTGTGGTGGGCGAGGCACTGGGCTGGCGTATCACCTTCCTGAGCATTGCCGGGGTAGCGCTAGcca                  >  6:350580/1‑149 (MQ=255)
                                                                                                                                   gccgGGCCAGCAGGCCAAGGCGCTGGGCCTGCTGGCCACCGGTACTACCCTGGCCATGGTCCTGGGCATCCCGCTGGGCCGTGTGGTGGGCGAGGCACTGGGCTGGCGTATCACGTTCCTGAGCATTGCCGGGGTAGCGCTAGccacca               <  4:177174/149‑1 (MQ=255)
                                                                                                                                   gccgGGCCAGCAGGCCAAGGCGCTGGGCCTGCTGGCCACCGGTACTACCCTGGCCATGGTCCTGGGCATCCCGCTGGGCCGTGTGGTGGGCGAGGCACTGGGCTGGCGTATCACGTTCCTGAGCATTGCCGGGGTAGCGCTAGccacca               >  7:303176/1‑149 (MQ=255)
                                                                                                                                           agcatGCCAAGGCGCTGGGCCTGCTGGCCACCGGTACTACCCTGGCCCTGGTCCTGGGCATCCCGCTGGGCCTGGTGGTGGGCGAGGCACTGGGCTGGCGTATCACGTTCCTGAGCATTGCCGGGGTAGCGCTAGCCACCATGCTTTGc       <  4:85682/149‑1 (MQ=255)
                                                                                                                                           agcagGCCAAGGCGCTGGGCCTGCTGGCCACCGGTACTACCCTGGCCATGGTCCTGGGCATCCCGCTGGGCCGTGTGGTGGGCGAGGCACTGGGCTGGCGTATCACGTTCCTGAGCATTGCCGGGGTAGCGCTAGCCACCATGCTTTGc       >  2:291322/1‑149 (MQ=255)
                                                                                                                                            gcagGCCAAGGCGCTGGGCCTGCTGGCCACCGGTACTACCCTGGCCATGGTCCTGGGCATCCCGCTGGGCCGTGTGGTGGGCGAGGCACTGGGCTGGCGTATCACGTTCCTGAGCATTGCCGGGGTAGCGCTAGCCACCATGCTTTGcc      <  7:779/149‑1 (MQ=255)
                                                                                                                                             cagGCCAAGGCGCTGGGCCTGCTGGCCACCGGTACTACCCTGGCCATGGTCCTGGGCATCCCGCTGGGCCGTGTGGTGGGCGAGGCACTGGGCTGGCGTATCACGTTCCTGAGCATTGCCGGGGTAGCGCTAGCCACCATGCTTTGCCt     <  4:130067/149‑1 (MQ=255)
                                                                                                                                                gCCAAGGCGCTGGGCCTGCTGGCCACCGGTACTACCCTGGCCATGGTCCTGGGCATCCCGCTGGGCCGTGTGGTGGGCGAGGCACTGGGCTGGCGTATCACGTTCCTGAGCATTGCCGGGGTAGCGCTAGCCACCATGCTTTGCCTGAt  >  5:74724/1‑149 (MQ=255)
                                                                                                                                                   |                                                                                                                                                 
TTCATTGTCAGCCACTTGCTGTCCTGGCTTTCGCAAAGCTTCGCCATGCTGCTGGTCAGCCGCATCGGCATTGCCCTGGCGCATGCGGTGTTCTGGTCCATCACCGCCTCGCTGGCGGTGCGTGTGGCGCCGCCGGGCCAGCAGGCCAAGGCGCTGGGCCTGCTGGCCACCGGTACTACCCTGGCCATGGTCCTGGGCATCCCGCTGGGCCGTGTGGTGGGCGAGGCACTGGGCTGGCGTATCACGTTCCTGAGCATTGCCGGGGTAGCGCTAGCCACCATGCTTTGCCTGAT  >  NC_002947/2775960‑2776252

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 21 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑