Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 2,351,983 T→A 5.7% S428C (AGT→TGT)  PP_2067 ← EmrB/QacA family drug resistance transporter

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029472,351,9830TA 5.7% 446.9 / 11.2 139S428C (AGT→TGT) PP_2067EmrB/QacA family drug resistance transporter
Reads supporting (aligned to +/- strand):  ref base T (67/64);  new base A (3/5);  total (70/69)
Fisher's exact test for biased strand distribution p-value = 4.93e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.99e-01

Only 100 of 142 total aligned reads displayed.

GGCCATCAGGTAAGCCTCTCGTCGAATGGTCGCCTCCAACGCTTCCAGCATTTGCTGCTGAGGCGGCAGCCAGCTTTGCACCGTCTGGTTGATACGCGCCTGCACGCTTGGCTCGAACAGCGTCACCTGCTCGTGAATGCGCTCACTGTGAAAGCGTTCACGGGTCGCCACCAGTTGCGTCAACCCGGCCGTACCGATGGCGCCACCCAGGTTGCGCAACATCGAGAACACCGCCGAGGCTGAGCCGGCTTCACGCTTGTCCAGCCCGGCCACCGCCAGCACCGATAACGCAAC  >  NC_002947/2351837‑2352130
                                                                                                                                                   |                                                                                                                                                   
ggCCATCAGGTAAGCCTCTCGTCGAATGGTCGCCTCCAACGCTTCCAGCATTTGCTGCTGAGGCGGCAGCCAGCTTTGCACCGTCTGGTTGATACGCGCCTGCACGCTTGGCTCGAACAGCGTCACCTGCTCGTGAATGCGCTCACtgt                                                                                                                                                   >  8:317013/1‑149 (MQ=255)
    aTCAGGTAAGCCTCTCGTCGAATGGTCGCCTCCAACGCTTCCAGCATTTGCTGCTGAGGCGGCAGCCAGCTTTGCACCGTCTGGTTGATACGCGCCTGCACGCTTGGCTCGAACAGCGTCACCTGCTCGTGAATGCGCTCACTGTGaaa                                                                                                                                               <  4:105250/149‑1 (MQ=255)
       aGGTAAGCCTCTCGTCGAATGGTCGCCTCCAACGCTTCCAGCATTTGCTGCTGAGGCGGCAGCCAGCTTTGCACCGTCTGGTTGATACGCGCCTGCACGCTTGGCTCGAACAGCGTCACCTGCTCGTGAACGCTTTCACAGTGagcgca                                                                                                                                            <  2:12272/149‑6 (MQ=255)
         gTAAGCCTCTCGTCGAATGGTCGCCTCCAACGCTTCCAGCATTTGCTGCTGAGGCGGCAGCCAGCTTTGCACCGTCTGGTTGATACGCGCCTGCACGCTTGGCTCGAACAGCGTCACCTGCTCGTGAATGCGCTCACTGTGaaa                                                                                                                                               <  8:19394/144‑1 (MQ=255)
         gTAAGCCTCTCGTCGAATGGTCGCCTCCAACGCTTCCAGCATTTGCTGCTGAGGCGGCAGCCAGCTTTGCACCGTCTGGTTGATACGCGCCTGCACGCTTGGCTCGAACAGCGTCACCTGCTCGTGAATGCGCTCACTGTGaaa                                                                                                                                               >  7:19394/1‑144 (MQ=255)
          tAAGCCTCTCGTCGAATGGTCGCCTCCAACGCTTCCAGCATTTGCTGCTGAGGCGGCAGCCAGCTTTGCACCGTCTGGTTGATACGCGCCTGCACGCTTGGCTCGAACAGCGTCACCTGCTCGTGAATGCGCTCACTGTGAAAGCGtt                                                                                                                                          >  1:76500/1‑148 (MQ=255)
            aGCCTCTCGTCGAATGGTCGCCTCCAACGCTTCCAGCATTTGCTGCTGAGGCGGCAGCCAGCTTTGCACCGTCTGGTTGATACGCGCCTGCACGCTTGGCTCGAACAGCGTCACCTGCTCGTGAATGCGCTCACTGTGAAAGCGTTCAc                                                                                                                                       <  8:149933/149‑1 (MQ=255)
            aGCCTCTCGTCGAATGGTCGCCTCCAACGCTTCCAGCATTTGCTGCTGAGGCGGCAGCCAGCTTTGCACCGTCTGGTTGATACGCGCCTGCACGCTTGGCTCGAACAGCGTCACCTGCTCGTGAACGCTTTCACAGTGAGCGCATTCAc                                                                                                                                       <  2:234699/149‑1 (MQ=255)
                tctcGTCGAATGGTCGCCTCCAACGCTTCCAGCATTTGCTGCTGAGGCGGCAGCCAGCTTTGCACCGTCTGGTTGATACGCGCCTGCACGCTTGGCTCGAACAGCGTCACCTGCTCGTGAATGCGCTCACTGTGAAAGCGTTCACGGGt                                                                                                                                   >  2:26015/1‑149 (MQ=255)
                 ctcGTCGAATGGTCGCCTCCAACGCTTCCAGCATTTGCTGCTGAGGCGGCAGCCAGCTTTGCACCGTCTGGTTGATACGCGCCTGCACGCTTGGCTCGAACAGCGTCACCTGCTCGTGAATGCGCTCACTGTGAAAGCGTTCACGGGTc                                                                                                                                  >  3:352218/1‑149 (MQ=255)
                    gtcgAATGGTCGCCTCCAACGCTTCCAGCATTTGCTGCTGAGGCGGCAGCCAGCTTTGCACCGTCTGGTTGATACGCGCCTGCACGCTTGGCTCGAACAGCGTCACCTGCTCGTGAATGCGCTCACTGTGAAAGCGTTCACGGGTc                                                                                                                                  <  3:258680/146‑1 (MQ=255)
                     tcgAATGGTCGCCTCCAACGCTTCCAGCATTTGCTGCTGAGGCGGCAGCCAGCTTTGCACCGTCTGGTTGATACGCGCCTGCACGCTTGGCTCGAACAGCGTCACCTGCTCGTGAATGCGCTCACTGTGAAAGCGTTCACGGGTCGcca                                                                                                                              >  3:176931/1‑149 (MQ=255)
                     tcgAATGGTCGCCTCCAACGCTTCCAGCATTTGCTGCTGAGGCGGCAGCCAGCTTTGCACCGTCTGGTTGATACGCGCCTGCACGCTTGGCTCGAACAGCGTCACCTGCTCGTGAACGCTTTCACAGTGAGCGCATTCACgag                                                                                                                                    >  4:105438/1‑141 (MQ=255)
                     tcgAATGGTCGCCTCCAACGCTTCCAGCATTTGCTGCTGAGGCGGCAGCCAGCTTTGCACCGTCTGGTTGATACGCGCCTGCACGCTTGGCTCGAACAGCGTCACCTGCTCGTGAACGCTTTCACAGTGAGCGCATTCACgag                                                                                                                                    <  3:105438/143‑3 (MQ=255)
                        aaTGGTCGCCTCCAACGCTTCCAGCATTTGCTGCTGAGGCGGCAGCCAGCTTTGCACCGTCTGGTTGATACGCGCCTGCACGCTTGGCTCGAACAGCGTCACCTGCTCGTGAATGCGCTCACTGTGAAAGCGTTCACGGGTCGccacca                                                                                                                           <  4:352218/149‑1 (MQ=255)
                             tCGCCTCCAACGCTTCCAGCATTTGCTGCTGAGGCGGCAGCCAGCTTTGCACCGTCTGGTTGATACGCGCCTGCACGCTTGGCTCGAACAGCGTCACCTGCTCGTGAATGCGCTCACTGTGAAAGCGt                                                                                                                                           <  5:79543/128‑1 (MQ=255)
                             tCGCCTCCAACGCTTCCAGCATTTGCTGCTGAGGCGGCAGCCAGCTTTGCACCGTCTGGTTGATACGCGCCTGCACGCTTGGCTCGAACAGCGTCACCTGCTCGTGAATGCGCTCACTGTGAAAGCGTTCACGGGTCGCCACCAGTTGc                                                                                                                      >  2:48183/1‑149 (MQ=255)
                               gCCTCCAACGCTTCCAGCATTTGCTGCTGAGGCGGCAGCCAGCTTTGCACCGTCTGGTTGATACGCGCCTGCACGCTTGGCTCGAACAGCGTCACCTGCTCGTGAATGCGCTCACTGTGAAAGCGTTCACGGGTCGCCACCAGTTGCGt                                                                                                                    >  7:333261/1‑149 (MQ=255)
                                   gcAACGCTTCCTGCATTTGCTGCTGAGGCGGCAGCCAGCTTTGCACCGTCTGGTTGATACGCGCCTGCACGCTTGGCTCGAAGAGCGTCACCTGCTCGTGAATGCGCTCACTGTGa                                                                                                                                                 <  6:102470/115‑1 (MQ=255)
                                   ccAACGCTTCCAGCATTTGCTGCTGAGGCGGCAGCCAGCTTTGCACCGTCTGGTTGATACGCGCCTGCACGCTTGGCTCGAACAGCGTCACCTGCTCGTGAATGCGCTCACTGTGAAAGCGTTCACGGGTCGCCACCAGTTGCGt                                                                                                                    >  5:287984/1‑145 (MQ=255)
                                   ccAACGCTTCCAGCATTTGCTGCTGAGGCGGCAGCCAGCTTTGCACCGTCTGGTTGATACGCGCCTGCACGCTTGGCTCGAACAGCGTCACCTGCTCGTGAATGCGCTCACTGTGAAAGCGTTCACGGGTCGCCACCAGTTGCGt                                                                                                                    <  6:287984/145‑1 (MQ=255)
                                   ccAACGCTTCCAGCATTTGCAGCTGAGGCGCCAGCCAGCTTTGCACCGTCTGGTTGATACGCGCCTGCACGCTTGGCTCGAACACCGTCAACTGCTCGTGAATACCCTCACTGTGa                                                                                                                                                 >  5:102470/1‑116 (MQ=255)
                                    cAACGCTTCCAGCATTTGCTGCTGAGGCGGCAGCCAGCTTTGCACCGTCTGGTTGATACGCGCCTGCACGCTTGGCTCGAACAGCGTCACCTGCTCGTGAATGCGCTCACTGTGAAAGCGTTCACGGGTCGCCACCAGTTGCGTCAAcc                                                                                                               >  5:306643/1‑149 (MQ=255)
                                    cAACGCTTCCAGCATTTGCTGCTGAGGCGGCAGCCAGCTTTGCACCGTCTGGTTGATACGCGCCTGCACGCTTGGCTCGAACAGCGTCACCTGCTCGTGAATGCGCTCACTGTGAAAGCGTTCACGGGTCGCCACCAGTTGCGTCAAcc                                                                                                               >  5:107108/1‑149 (MQ=255)
                                    cAACGCTTCCAGCATTTGCTGCTGAGGCGGCAGCCAGCTTTGCACCGTCTGGTTGATACGCGCCTGCACGCTTGGCTCGAACAGCGTCACCTGCTCGTGAATGCGCTCACTGTGAAAGCGTTCACGGGTCGCCACCAGTTGCGTCAAcc                                                                                                               >  8:77565/1‑149 (MQ=255)
                                       cGCTTCCAGCATTTGCTGCTGAGGCGGCAGCCAGCTTTGCACCGTCTGGTTGATACGCGCCTGCACGCTTGGCTCGAACAGCGTCACCTGCTCGTGAATGCGCTCACTGTGAAAGCGTTCACGGGTCGCCACCAGTTGCGTCAACCCgg                                                                                                            <  3:128647/149‑1 (MQ=255)
                                         cTTCCAGCATTTGCTGCTGAGGCGGCAGCCAGCTTTGCACCGTCTGGTTGATACGCGCCTGCACGCTTGGCTCGAACAGCGTCACCTGCTCGTGAATGCGCTCACTGTGAAAGCGTTCACGGGTCGCCACCAGTTGCGTCAACCCGGcc                                                                                                          <  8:239335/149‑1 (MQ=255)
                                          ttCCAGCATTTGCTGCTGAGGCGGCAGCCAGCTTTGCACCGTCTGGTTGATACGCGCCTGCACGCTTGGCTCGAACAGCGTCACCTGCTCGTGAATGCGCTCACTGTGAAAGCGTTCACGGGTCGCCACCAGTTGCGTCAACCCGGCCg                                                                                                         >  1:272486/1‑149 (MQ=255)
                                            ccAGCATTTGCTGCTGAGGCGGCAGCCAGCTTTGCACCGTCTGGTTGATACGCGCCTGCACGCTTGGCTCGAACAGCGTCACCTGCTCGTGAATGCGCTCACTGTGAAAGCGtt                                                                                                                                          >  1:152865/1‑114 (MQ=255)
                                            ccAGCATTTGCTGCTGAGGCGGCAGCCAGCTTTGCACCGTCTGGTTGATACGCGCCTGCACGCTTGGCTCGAACAGCGTCACCTGCTCGTGAATGCGCTCACTGTGAAAGCGtt                                                                                                                                          <  2:152865/114‑1 (MQ=255)
                                             cAGCATTTGCTGCTGAGGCGGCAGCCAGCTTTGCACCGTCTGGTTGATACGCGCCTGCACGCTTGGCTCGAACAGCGTCACCTGCTCGTGAATGCGCTCACTGTGAAAGCGTTCACGGGTCGCCACCAGTTGCGTCAACCCg                                                                                                             <  1:214443/142‑1 (MQ=255)
                                             cAGCATTTGCTGCTGAGGCGGCAGCCAGCTTTGCACCGTCTGGTTGATACGCGCCTGCACGCTTGGCTCGAACAGCGTCACCTGCTCGTGAATGCGCTCACTGTGAAAGCGTTCACGGGTCGCCACCAGTTGCGTCAACCCGGCCGTa                                                                                                       <  5:53076/148‑1 (MQ=255)
                                             cAGCATTTGCTGCTGAGGCGGCAGCCAGCTTTGCACCGTCTGGTTGATACGCGCCTGCACGCTTGGCTCGAACAGCGTCACCTGCTCGTGAATGCGCTCACTGTGAAAGCGTTCACGGGTCGCCACCAGTTGCGTCAACCCGGCCGTAc                                                                                                      <  5:124667/149‑1 (MQ=255)
                                             cAGCATTTGCTGCTGAGGCGGCAGCCAGCTCTGCACCGTCTGGTTGATACGCGCCTGCACGCTTGGCTCGAACAGCGTCACCTGCTCGTGAATGCGCTCACTGTGAAAGCGTTCACGGGTCGCCACCAGTTGCGTCAACCCGGCCGTa                                                                                                       <  8:188356/148‑1 (MQ=255)
                                              aGCATTTGCTGCTGAGGCGGCAGCCAGCTTTGCACCGTCTGGTTGATACGCGCCTGCACGCTTGGCTCGAACAGCGTCACCTGCTCGTGAATGCGCTCACTGTGAAAGCGTTCACGGGTCGCCACCAGTTGCGTCAACCCGGCCGTAcc                                                                                                     >  2:260387/1‑149 (MQ=255)
                                               gCATTTGCTGCTGAGGCGGCAGCCAGCTTTGCACCGTCTGGTTGATACGCGCCTGCACGCTTGGCTCGAACAGCGTCACCTGCTCGTGAATGCGCTCACTGTGAAAGCGTTCACGGGTCGCCACCAGTTGCGTCAACCCGGCCGTAcc                                                                                                     >  4:14254/1‑148 (MQ=255)
                                                  tttGCTGCTGAGGCGGCAGCCAGCTTTGCACCGTCTGGTTGATACGCGCCTGCACGCTTGGCTCGAACAGCGTCACCTGCTCGTGAATGCGCTCACTGTGa                                                                                                                                                 >  6:67503/1‑101 (MQ=255)
                                                  tttGCTGCTGAGGCGGCAGCCAGCTTTGCACCGTCTGGTTGATACGCGCCTGCACGCTTGGCTCGAACAGCGTCACCTGCTCGTGAATGCGCTCACTGTGAAAGCGTTCACGGGTCGCCACCAGTTGCGTCAACCCGGCCGTACCGATg                                                                                                 <  7:222969/149‑1 (MQ=255)
                                                    tgctgcTGAGGCGGCAGCCAGCTTTGCACCGTCTGGTTGATACGCGCCTGCACGCTTGGCTCGAACAGCGTCACCTGCTCGTGAATGCGCTCACTGTGAAAGCGTTCACGGGTCGCCACCAGTTGCGTCAACCCGGCCGTACCGTTgg                                                                                                <  6:33421/148‑1 (MQ=255)
                                                     gctgctGAGGCGGCAGCCAGCTTTGCACCGTCTGGTTGATACGCGCCTGCACGCTTGGCTCGAACAGCGTCACCTGCTCGTGAATGCGCTCACTGTGAAAGCGTTCACGGGTCGCCACCAGTTGCGTCAACCCGGCCGTACCGATGgc                                                                                               <  1:48183/148‑1 (MQ=255)
                                                            aGGCGGCAGCCAGCTTTGCACCGTCTGGTTGATACGCGCCTGCACGCTTGGCTCGAGACAGCGTCACCTGCTCGTGAATGCGCTCACTGTGAAAGCGTTCACGGGTCGCCACCAGTTGCGTCAACCCGGCCGTACCGATGGCGCCAcc                                                                                         >  1:280973/1‑148 (MQ=255)
                                                            aGGCGGCAGCCAGCTTTGCACCGTCTGGTTGATACGCGCCTGCACGCTTGGCTCGAACAGCGTCACCTGCTCGTGAATGCGCTCACTGTGAAAGCGTTCACGGGTCGCCACCAGTTGCGTCAACCCGGCCGTACCGATGGCGCCAcc                                                                                         >  3:85938/1‑147 (MQ=255)
                                                              gcggcAGCCAGCTTTGCACCGTCTGGTTGATACGCGCCTGCACGCTTGGCTCGAACAGCGTCACCTGCTCGTGAATGCGCTCACTGTGAAAGCGTTCACGGGTCGCCACCAGTTGCGTCAACCCGGCCGTACCGATGGCGCCACCCAgg                                                                                     >  5:213666/1‑149 (MQ=255)
                                                              gcggcAGCCAGCTTTGCACCGTCTGGTTGATACGCGCCTGCACGCTTGGCTCGAACAGCGTCACCTGCTCGTGAATGCGCTCACTGTGAAAGCGTTCACGGGTCGCCACCAGTTGCGTCAACCCGGCCGTACCGATGGCGCCACCCAgg                                                                                     <  7:317013/149‑1 (MQ=255)
                                                                ggcAGCCAGCTTTGCACCGTCTGGTTGATACGCGCCTGCACGCTTGGCTCGAACAGCGTCACCTGCTCGTGAATGCGCTCACTGTGAAAGCGTTCACGGGTCGCCACCAGTTGCGTCAACCCGGCCGTACCGATGGCGCCACCCAGGtt                                                                                   <  4:156515/149‑1 (MQ=255)
                                                                  cagccagcTTTGCACCGTCTGGTTGATACGCGCCTGCACGCTTGGCTCGAACAGCGTCACCTGCTCGTGAATGCGCTCACTGTGAAAGCGTTCACGGGTCGCCACCAGTTGc                                                                                                                      <  3:295841/112‑1 (MQ=255)
                                                                      cagcTTTGCACCGTCTGGTTGATACGCGCCTGCACGCTTGGCTCGAACAGCGTCACCTGCTCGTGAATGCGCTCACTGTGAAAGCGTTCACGGGTCGCCACCAGTTGCGTCACCCCGGCAGTACCGATGGCGCCACCCAGGTTGCGTaa                                                                             >  8:115357/1‑149 (MQ=255)
                                                                       agcTTTGCACCGTCTGGTTGATACGCGCCTGCACGCTTGGCTCGAACAGCGTCACCTGCTCGTGAATGCGCTCACTGTGAAAGCGTTCACGGGTCGCCACCAGTTGCGTCAACCCGGCCGTACCGATGGCGCCACCCAGGTTGCGCaa                                                                             >  3:160652/1‑148 (MQ=255)
                                                                          tttGCACCGTCTGGTTGATACGCGCCTGCACGCTTGGCTCGAACAGCGTCACCTGCTCGTGAATGCGCTCACTGTGAAAGCGTTCACGGGTCGCCACCAGTTGCGTCAACCCg                                                                                                             >  5:219381/1‑113 (MQ=255)
                                                                          tttGCACCGTCTGGTTGATACGCGCCTGCACGCTTGGCTCGAACAGCGTCACCTGCTCGTGAATGCGCTCACTGTGAAAGCGTTCACGGGTCGCCACCAGTTGCGTCAACCCg                                                                                                             <  6:219381/113‑1 (MQ=255)
                                                                             gCACCGTCTGGTTGATACGCGCCTGCACGCTTGGCTCGAACAGCGTCACCTGCTCGTGAATGCGCTCACTGTGAAAGCGTTCACGGGTCGCCACCAGTTGCGTCAACCCGGCCGTACCGATGGCGCCACCCAGGTTGCGCAACATCga                                                                       <  2:280973/148‑1 (MQ=255)
                                                                             gCACCGTCTGGTTGATACGCGCCTGCACGCTTGGCTCGAACAGCGTCACCTGCTCGTGAATGCGCTCACTGTGAAAGCGTTCACGGGTCGCCACCAGTTGCGTCAACCCGGCCGTACCGATGGCGCCACCCAGGTTGCGCAACATCga                                                                       <  1:260387/148‑1 (MQ=255)
                                                                              cACCGTCTGGTTGATACGCGCCTGCACGCTTGGCTCGAACAGCGTCACCTGCTCGTGAATGCGCTCACTGTGAAAGCGTTCACGGGTCGCCACCAGTTGCGTCAACCCGGCCGTACCGATGGCGCCACCCAGGTTGCGCAACATCgaga                                                                     <  4:160652/149‑1 (MQ=255)
                                                                              cACCGTCTGGTTGATACGCGCCTGCACGCTTGGCTCGAACAGCGTCACCTGCTCGTGAATGCGCTCACTGTGAAAGCGTTCACGGGTCGCCACCAGTTGCGTCAACCCGGCCGTACCGATGGCGCCACCCAGGTTGCGCAACATCgag                                                                      <  3:359968/148‑1 (MQ=255)
                                                                              cACCGTCTGGTTGATACGCGCCTGCACGCTTGGCTCGAACAGCGTCACCTGCTCGTGAATGCGCTCACTGTGAAAGCGTTCACGGGTCGCCACCAGTTGCGTCAACCCGGCCGTACCGATGGCGCCACCCAGGTTGCGCAACATCgag                                                                      >  4:13647/1‑148 (MQ=255)
                                                                               aCCGTCTGGTTGATACGCGCCTGCACGCTTGGCTCGAACAGCGTCACCTGCTCGTGAATGCGCTCACTGTGAAAGCGTTCACGGGTCGCCACCAGTTGCGTCAACCCGGCCGTACCGATGGCGCCACCCAGGTTGCGCAACATCGAGaa                                                                    <  6:107108/149‑1 (MQ=255)
                                                                                     tGGTTGATACGCGCCTGCACGCTTGGCTCGAACAGCGTCACCTGCTCGTGAATGCGCTCACTGTGAAAGCGTTCACGGGTCGCCACCAGTTGCGTCaa                                                                                                                 >  2:208861/1‑98 (MQ=255)
                                                                                     tGGTTGATACGCGCCTGCACGCTTGGCTCGAACAGCGTCACCTGCTCGTGAATGCGCTCACTGTGAAAGCGTTCACGGGTCGCCACCAGTTGCGTCAAACAGGCAGTACCGATGGCGCCACCCAGGTTGCGCAACATCGAGAacaactc                                                              >  1:95222/1‑145 (MQ=255)
                                                                                       gTTGATACGCGCCTGCAGGCTTGGCTCGAACAGCGTCACCTGCTCGTGAATGCGCTCACTGTGAAAGCGTTCACGGGTCGCCACCAGTTGCGTCAACCCGGCCGTACCGTTGGCGCCACCCAGGTTGCGCAACATCGAGAACAccgccg                                                            <  5:8492/149‑1 (MQ=255)
                                                                                       gTTGATACGCGCCTGCACGCTTGGCTCGAACAGCGTCACCTGCTCGTGAATGCGCTCACTGTGAAAGCGTTCACGGGTCGCCACCAGTTGCGTCAACCCGGCCGTACCGATGGCGCCACCCAGGTTgc                                                                                 <  6:95938/128‑1 (MQ=255)
                                                                                        ttGATACGCGCCTGCACGCTTGGCTCGAACAGCGTCACCTGCTCGTGAATGCGCTCACTGTGAAAGCGTTCACGGGTCGCCACCAGTTGCGTCAACCCGGCCGTACCGATGGCGCCACCCAGGTTGCGCAACATCGAGAACACCGCCGa                                                           <  2:76500/149‑1 (MQ=255)
                                                                                             aCGCGCCTGCACGCTTGGCTCGAACAGCGTCACCTGCTCGTGAATGCGCTCACTGTGAAAGCGTTCACGGGTCGCCACCAGTTGCGTCAACCCGGCCGTACCGATGGCGCCACCCAGGTTGCGCACCATCGAGCACACCGCCGCGGCTg                                                      >  8:246890/1‑149 (MQ=255)
                                                                                             aCGCGCCTGCACGCTTGGCTCGAACAGCGTCACCTGCTCGTGAATGCGCTCACTGTGAAAGCGTTCACGGGTCGCCACCAGTTGCGTCAACCCGGCCGTACCGATGGCGCCACCCAGGTTGCGCAACATCGAGAACACCGCCGAGGCt                                                       >  7:250894/1‑148 (MQ=255)
                                                                                               gcgcCTGCACGCTTGGCTCGAACAGCGTCACCTGCTCGTGAATGCGCTCACTGTGAAAGCGTTCACGGGTCGCCACCAGTTGCGTCAACCCGGCCGTACCGATGGCGCCACCCAGGTTGCGCa                                                                              >  2:71244/1‑123 (MQ=255)
                                                                                               gcgcCTGCACGCTTGGCTCGAACAGCGTCACCTGCTCGTGAATGCGCTCACTGTGAAAGCGTTCACGGGTCGCCACCAGTTGCGTCAACCCGGCCGTACCGATGGCGCCACCCAGGTTGCGCa                                                                              <  1:71244/123‑1 (MQ=255)
                                                                                                   cTGCACGCTTGGCTCGAACAGCGTCACCTGCTCGTGAATGCGCTCACTGTGAAAGCGTTCACGGGTCGCCACCAGTTGCGTCAACCCGGCCGTACCGATGGCGCCACCCAGGTTGCGCAACATCGAGAACACCGCCGAGGCTGAGCCgg                                                >  1:167292/1‑149 (MQ=255)
                                                                                                     gCACGCTTGGCTCGAACAGCGTCACCTGCTCGTGAATGCGCTCACTGTGAAAGCGTTCACGGGTCGCCACCAGtt                                                                                                                        <  4:69083/75‑1 (MQ=255)
                                                                                                     gCACGCTTGGCTCGAACAGCGTCACCTGCTCGTGAATGCGCTCACTGTGAAAGCGTTCACGGGTCGCCACCAGTTGCGTCAACCCGGCCGTACCGATGGCGCCACCCAGGtt                                                                                   <  1:119909/112‑1 (MQ=255)
                                                                                                      cACGCTTGGCTCGAACAGCGTCACCTGCTCGTGAATGCGCTCACTGTGAAAGCGTTCACGGGTCGCCACCAGTTGCGTCAACCCGGCCGTACCGATGGCGCCACCCAGGTTGCGCAACATCGAGAACACCGCCGAGGCTGAGCCGGCtt                                             >  1:219324/1‑149 (MQ=255)
                                                                                                         gCTTGGCTCGAACAGCGTCACCTGCTCGTGAATGCGCTCACTGTGAAAGCGTTCACGGGTCGCCAACAGTTGCGTCAACCAGGCCGTACCGATGGCGCCCCCCAGGTTGCGCAACATCGAGAACACCGCCGATGCTGAGCCGGCTTCAc                                          >  5:85017/1‑149 (MQ=255)
                                                                                                           ttGGCTCGAACAGCGTCACCTGCTCGTGAATGCGCTCACTGTGAAAGCGTTCACGGGTCGCCACCAGTTGCGTCAACCCGGCCGTACCGATGGCGCCACCCAGGTTGCGCAACATCGAGAACACCGCCGAGGCTGAGCCGGCTTCACGc                                        >  7:45930/1‑149 (MQ=255)
                                                                                                           ctGGCTCGAACAGCGTCACCTGCTCGTGAATGCGCTCACTGTGAAAGCGTTCACGGGTCGCCACCAGTTGCGTCAACCCGGCCGTACCGATGGCGCCACCCAGGTTGCGCAACATCGAGAACACCGCCGAGGCTGAGCCGGCTTCACGc                                        <  6:114041/148‑1 (MQ=255)
                                                                                                             ggCTCGAACAGCGTCACCTGCTCGTGAATGCGCTCACTGTGAAAGCGTTCACGGGTCGCCACCAGTTGCGTCAACCCGGCCGTACCGATGGCGCCACCCAGGTTGCGCAACATCGAGAACACCGCCGAGGCTGAGCCGGCTTCACGCtt                                      >  6:150268/1‑149 (MQ=255)
                                                                                                              gCTCGAACAGCGTCACCTGCTCGTGAATGCGCTCACTGTGAAAGCGTTCACGGGTCGCCACCAGTTGCGTCAACCCGGCCGTACCGATGGCGCCACCCAGGTTGCGCAACATCGAGAACACCGCCGAGGCTGAGCCGGCTTCACGCtt                                      <  4:176931/148‑1 (MQ=255)
                                                                                                              gCTCGAACAGCGTCACCTGCTCGTGAATGCGCTCACTGTGAAAGCGTTCACGGGTCGCCACCAGTTGCGTCAACCCGGCCGTACCGATGGCGCCACCCAGGTTGCGCAACATCGAGAACACCGCCGAGGCTGAGCCGGCTTCACGCtt                                      <  8:333261/148‑1 (MQ=255)
                                                                                                              gCTCGAACAGCGTCACCTGCTCGTGAATGCGCTCACTGTGAAAGCGTTCACGGGTCGCCACCAGTTGCGTCAACCCGGCCGTACCGATGGCGCCACCCAGGTTGCGCAACATCGAGAACACCGCCGAGGCTGAGCCGGCTTCACGCTTg                                     <  2:219324/149‑1 (MQ=255)
                                                                                                              gCTCGAACAGCGTCACCTGCTCGTGAATGCGCTCACTGTGAAAGCGTTCACGGGTCGCCACCAGTTGCGTCAACCCGGCCGTACCGATGGCGCCACCCAGGTTGCGCAACATCGAGAACACCGCCGAGGCTGAGCCGGCTTCACGCTTg                                     <  7:434/149‑1 (MQ=255)
                                                                                                               cTCGAACAGCGTCACCTGCTCGTGAATGCGCTCACTGTGAAAGCGTTCACGGGTCGCCACCAGTTGCGTCAACCCGGCCGTACCGATGGCGCCACCCAGGTTGCGCAACATCGAGAACACCGCCGAGGCTGAGCCGGCTTCACGCTTGt                                    >  2:309448/1‑149 (MQ=255)
                                                                                                                  gAACAGCGTCACCTGCTCGTGAATGCGCTCACTGTGAAAGCGTTCACGGGTCGCCACCAGTTGCGTCAACCCGGCCGTACCGATGGCGCCACCCAGGTTGCGCAACATCGAGAACACCGCCGAGGCTGAGCCGGCTTCACGCTTGTcc                                  <  7:77565/148‑1 (MQ=255)
                                                                                                                   aaCAGCGTCACCTGCTCGTGAATGCGCTCACTGTGAAAGCGTTCACGGGTCGCCACCAGTTGCGTCAACCCGGCCGTACCGATGGCGCCACCCAGGTTGCGCAACATCGAGAACACCGCCGAGGCTGAGCCGGCTTCACGCTTGTc                                   >  6:48759/1‑146 (MQ=255)
                                                                                                                   aaCAGCGTCACCTGCTCGTGAATGCGCTCACTGTGAAAGCGTTCACGGGTCGCCACCAGTTGCGTCAACCCGGCCGTACCGATGGCGCCACCCAGGTTGCGCAACATCGAGAACACCGCCGAGGCTGAGCCGGCTTCACGCTTGTCCa                                 <  5:48759/148‑1 (MQ=255)
                                                                                                                             aCCTGCTCGTGAATGCGCTCACTGTGAAAGCGTTCACGGGTCGCCACCAGTTGCGTCAACCCGGCCGTACCGATGGCGCCACCCAGGTTGCGCAACATCGAGAACACCGCCGAGGCTGAGCCGGCTTCACGCTTGTCCAGCCCGGCCAc                       >  3:104567/1‑149 (MQ=255)
                                                                                                                               cTGCTCGTGAATGCGCTCACTGTGAAAGCGTTCACGGGTCGCCACCAGTTGCGTCAACCCGGCCGTACCGATGGCGCCACCCAGGTTGCGCAACATCGAGAACACCGCCGAGGCTGAGCCGGCTTCACGc                                        <  5:146509/130‑1 (MQ=255)
                                                                                                                               cTGCTCGTGAATGCGCTCACTGTGAAAGCGTTCACGGGTCGCCACCAGTTGCGTCAACCCGGCCGTACCGATGGCGCCACCCAGGTTGCGCAACATCGAGAACACCGCCGAGGCTGAGCCGGCTTCACGc                                        >  6:146509/1‑130 (MQ=255)
                                                                                                                                   tCGTGAATGCGCTCACTGTGAAAGCGTTCACGGGTCGCCACCAGTTGCGt                                                                                                                    >  5:197990/1‑50 (MQ=255)
                                                                                                                                   tCGTGAATGCGCTCACTGTGAAAGCGTTCACGGGTCGCCACCAGTTGCGt                                                                                                                    <  6:197990/50‑1 (MQ=255)
                                                                                                                                   tCGTGAATGCGCTCACTGTGAAAGCGTTCACGGGTCGCCACCAGTTGCGTCAACCCGGCCGTACCGATGGCGCCACCAAGGTTGCGCAACATCGAGAACCCCGCCGAGGATGACCCGGCTTCACGCTTGTCCAGCCCGGCCACCGCCAg                 >  8:261453/1‑149 (MQ=255)
                                                                                                                                   ccGTGAACGCTTTCACAGTGAGCGCATTCACGGGTCGCCACCAGTTGCGTCAACCCGGCCGTACCGATGGCGCCACCCAGGTTGCGCAACATCGAGAACACCGCCGAGGCTGAGCCGGCTTCACGCTTGTCCAGCCCGGCCACCGCCAg                 >  6:360927/2‑149 (MQ=255)
                                                                                                                                     gTGAATGCGCTCACTGTGAAAGCGTTCACGGGTCGCCACCAGTTGCGTCAACCCGGCCGTACCGATGGCGCCACCCAGGTTGCGCAACATCGAGAACACCGCCGAGGCTGAGCCGGCTTCACGCTTGTCCAGCCCGGCCACCGCCAGCa               <  4:85938/149‑1 (MQ=255)
                                                                                                                                     gTGAACGCTTTCACAGTGAGCGCATTCACGGGTCGCCACCAGTTGCGTCAACCCGGCCGTACCGATGGCGCCACCCAGGTTGCGCAACATCGAGAACACCGCCGAGGCTGAGCCGGCTTCACGCTTGTCCAGCCCGGCCACCGCCAGCa               >  4:72527/1‑149 (MQ=255)
                                                                                                                                      tGAATGCGCTCACTGTGAAAGCGTTCACGGGTCGCCACCAGTTGCGTCAACCCGGCCGTACCGATGGCGCCACCCAGGTTGCGCAACATCGAGAACACCGCCGAGGCTGAGCCGGCTTCACGCTTGTCCAGCCCGGCCACCGCCAGCAc              <  3:13647/149‑1 (MQ=255)
                                                                                                                                          tGCGCTCACTGTGAAAGCGTTCACGGGTCGCCACCAGTTGCGTCAACCCGGCCGTACCGATGGCGCCACCCAGGTTGCGCAACATCGAGAACACCGCCGAGGCTGAGCCGGCTTCACGc                                        >  4:227019/1‑119 (MQ=255)
                                                                                                                                          tGCGCTCACTGTGAAAGCGTTCACGGGTCGCCACCAGTTGCGTCAACCCGGCCGTACCGATGGCGCCACCCAGGTTGCGCAACATCGAGAACACCGCCGAGGCTGAGCCGGCTTCACGc                                        <  3:227019/119‑1 (MQ=255)
                                                                                                                                           gcgcTCACTGTGAAAGCGTTCACGGGTCGCCACCAGTTGCGTCAACCCGGCCGTACCGATGGCTCCACCCAGGTTGCGCAACATCGAGAACACCCCCGAGGCTGAGCCGGCTTCACGCTTGTCCAGCCCGGCCACCGCCAGCCCCGAt          >  7:235414/1‑148 (MQ=255)
                                                                                                                                               tCACTGTGAAAGCGTTCACGGGTCGCCACCAGTTGCGTCAACCCGGCCGTACCGATGGCGCCACCCAGGTTGCGCAACATCGAGAACACCGCCGAGGCTGAGCCGGCTTc                                            >  3:81630/1‑110 (MQ=255)
                                                                                                                                               tCACTGTGAAAGCGTTCACGGGTCGCCACCAGTTGCGTCAACCCGGCCGTACCGATGGCGCCACCCAGGTTGCGCAACATCGAGAACACCGCCGAGGCTGAGCCGGCTTCACGCTTGTCCAGCCCGGCCACCGCCAGCACCGATAACg      >  2:251024/1‑148 (MQ=255)
                                                                                                                                               tCACTGTGAAAGCGTTCACGGGTCGCCACCAGTTGCGTCAACCCGGCCGTACCGATGGCGCCACCCAGGTTGCGCAACATCGAGAACACCGCCGAGGCTGAGCCGGCTTCACGCTTGTCCAGCCCGGCCACCGCCAGCACCGATAACg      >  3:9736/1‑148 (MQ=255)
                                                                                                                                               tCACTGTGAAAGCGTTCACGGGTCGCCACCAGTTGCGTCAACCCGGCCGTACCGATGGCGCCACCCAGGTTGCGCAACATCGAGAACACCGCCGAGGCTGAGCCGGCTTCACGCTTGTCCAGCCCGGCCACCGCCAGCACCGATAACg      >  6:366337/1‑148 (MQ=255)
                                                                                                                                                 aCTGTGAAAGCGTTCACGGGTCGCCACCAGTTGCGTCAACCCGGCCGTACCGATGGCGCCACCCAGGTTGCGCAACATCGAGAACACCGCCGAGGCTGAGCCGGCTTCACGCTTGTCCAGCCCGGCCACCGCCAGCACCGATAACGCaa   <  8:45930/149‑1 (MQ=255)
                                                                                                                                                   tgtgAAAGCGTTCACGGGTCGCCACCAGTTGCGTCAACCCGGCCGTACCGATGGCGCCACCCAGGTTGCGCAACATCGAGAACACCGCCGAGGCTGAGCCGGCTTCACGCTTGTCCAGCCCGGCCACCGCCAGCACCGATAACGCAAc  >  6:301138/1‑148 (MQ=255)
                                                                                                                                                   |                                                                                                                                                   
GGCCATCAGGTAAGCCTCTCGTCGAATGGTCGCCTCCAACGCTTCCAGCATTTGCTGCTGAGGCGGCAGCCAGCTTTGCACCGTCTGGTTGATACGCGCCTGCACGCTTGGCTCGAACAGCGTCACCTGCTCGTGAATGCGCTCACTGTGAAAGCGTTCACGGGTCGCCACCAGTTGCGTCAACCCGGCCGTACCGATGGCGCCACCCAGGTTGCGCAACATCGAGAACACCGCCGAGGCTGAGCCGGCTTCACGCTTGTCCAGCCCGGCCACCGCCAGCACCGATAACGCAAC  >  NC_002947/2351837‑2352130

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: