Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 2,238,940 C→A 8.4% intergenic (+70/+58) uvrB → / ← PP_1975 uvrABC excinuclease protein B/EmrB/QacA family drug resistance transporter

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029472,238,9400CA 8.4% 438.0 / 23.8 143intergenic (+70/+58)uvrB/PP_1975uvrABC excinuclease protein B/EmrB/QacA family drug resistance transporter
Reads supporting (aligned to +/- strand):  ref base C (68/63);  new base A (3/9);  total (71/72)
Fisher's exact test for biased strand distribution p-value = 1.29e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.21e-01

Only 100 of 150 total aligned reads displayed.

CCCGCGACCTGGAGTTCGAAGCTGCGGCGCAGTTGCGTGACGAGATTGCGCAGCTGCGGGAGCGGTTGATCAGCAGCTGAAGCCGGGCGCGGGCCCTGTGGGAGCGGGCGTGCCCGCGAAGAATTCAACGCGGTGGCTGGCACCGGCTGCGCCGGTGTTCGCGGGCATGCCCGCTCCCACAGGGACCGCATCTGATCTGTAATAGCCTCAATGGGCTTCCCCAGCCTTCAACCCTTCCGGCAACTTGCGGGTCAGCAGCACCGCCACCATGCTCACCGCCAACCCGATCCCGACAAAG  >  NC_002947/2238791‑2239088
                                                                                                                                                     |                                                                                                                                                    
cccGCGACCTGGAGTTCGAAGCTGCGGCGCAGTTGCGTGACGAGATTGCGCAGCTGCGGGAGCGGTTGATCAGCAGCTGAAGCCGGGCGCGGGCCCTGTGGGAGCGGGCGTGCCCGCGAAGAATTCAACGCGGTGGCTGGCACCGGCTgc                                                                                                                                                      <  1:5002027/150‑1 (MQ=255)
cccGCGACCTGGAGTTCGAAGCTGCGGCGCAGTTGCGTGACGAGATTGCGCAGCTGCGGGAGCGGTTGATCAGCAGCTGAAGCCGGGCGCGGGCCCTGTGGGAGCGGGCGTGCCCGCGAAGAATTCAACGCGGTGGCTGGCACCGGCTgc                                                                                                                                                      >  2:1019174/1‑150 (MQ=255)
  cgcgACCTGGAGTTCGAAGCTGCGGCGCAGTTGCGTGACGAGATTGCGCAGCTGCGGGAGCGGTTGATCAGCAGCTGAAGCCGGGCGCGGGCCCTGTGGGAGCGGGCGTGCCCGCGAAGAATTCAACGCGGTGGCTGGCACCGGCTgcgc                                                                                                                                                    <  2:8697087/150‑1 (MQ=255)
  cgcgACCTGGAGTTCGAAGCTGCGGCGCAGTTGCGTGACGAGATTGCGCAGCTGCGGGAGCGGTTGATCAGCAGCTGAAGCCGGGCGCGGGCCCTGTGGGAGCGGGCGTGCCCGCGAAGAATTCAACGCGGTGGCTGGCACCGGCTgcgc                                                                                                                                                    <  1:3699186/150‑1 (MQ=255)
  cgcgACCTGGAGTTCGAAGCTGCGGCGCAGTTGCGTGACGAGATTGCGCAGCTGCGGGAGCGGTTGATCAGCAGCTGAAGCCGGGCGCGGGCCCTGTGGGAGCGGGCGTGCCCGCGAAGAATTCAACGCGGTGGCTGGCACCGGCTgcgc                                                                                                                                                    <  2:3225353/150‑1 (MQ=255)
  cgcgACCTGGAGTTCGAAGCTGCGGCGCAGTTGCGTGACGAGATTGCGCAGCTGCGGGAGCGGTTGATCAGCAGCTGAAGCCGGGCGCGGGCCCTGTGGGAGCGGGCGTGCCCGCGAAGAATTCAACGCGGTGGCTGGCACCGGCTgcgc                                                                                                                                                    <  2:6380825/150‑1 (MQ=255)
      aCCTGGAGTTCGAAGCTGCGGCGCAGTTGCGTGACGAGATTGCGCAGCTGCGGGAGCGGTTGATCAGCAGCTGAAGCCGGGCGCGGGCCCTGTGGGAGCGGGCGTGCCCGCGAAGAATTCAACGCGGTGGCTGGCACCGGCTGCGCCGgt                                                                                                                                                <  1:2224505/150‑1 (MQ=255)
      aCCTGGAGTTCGAAGCTGCGGCGCAGTTGCGTGACGAGATTGCGCAGCTGCGGGAGCGGTTGATCAGCAGCTGAAGCCGGGCGCGGGCCCTGTGGGAGCGGGCGTGCCCGCGAAGAATTCAACGCGGTGGCTGGCACCGGCTGCGCCGgt                                                                                                                                                <  2:703249/150‑1 (MQ=255)
      aCCTGGAGTTCGAAGCTGCGGCGCAGTTGCGTGACGAGATTGCGCAGCTGCGGGAGCGGTTGATCAGCAGCTGAAGCCGGGCGCGGGCCCTGTGGGAGCGGGCGTGCCCGCGAAGAATTCAACGCGGTGGCGAACACCGGCGCAGCCGgt                                                                                                                                                <  2:3776501/150‑1 (MQ=39)
        cTGGAGTTCGAAGCTGCGGCGCAGTTGCGTGACGAGATTGCGCAGCTGCGGGAGCGGTTGATCAGCAGCTGAAGCCGGGCGCGGGCCCTGTGGGAGCGGGCGTGCCCGCGAAGAATTCAACGCGGTGGCTGGCACCGGCTGCGCCGgtgt                                                                                                                                              >  1:8883391/1‑150 (MQ=255)
        cTGGAGTTCGAAGCTGCGGCGCAGTTGCGTGACGAGATTGCGCAGCTGCGGGAGCGGTTGATCAGCAGCTGAAGCCGGGCGCGGGCCCTGTGGGAGCGGGCGTGCCCGCGAAGAATTCAACGCGGTGGCTGGCACCGGCGCAGCCGgtgc                                                                                                                                              <  2:2968787/150‑2 (MQ=255)
          ggAGTTCGACGCTGCGGCGCAGTTGCGTGACGAGATTGCGCAGCTGCGGGAGCGGTTGATCAGCAGCTGAAGCCGGGCGCGGGCCCTGTGGGAGCGGGCGTGCCCGCGAAGAATTCAACGCGGTGGCTGGCACCGGCGCAGCCGgtgcca                                                                                                                                            <  1:2220157/150‑4 (MQ=38)
          ggAGTTCGAAGCTGCGGCGCAGTTGCGTGACGAGATTGCGCAGCTGCGGGAGCGGTTGATCAGCAGCTGAAGCCGGGCGCGGGCCCTGTGGGAGCGGGCGTGCCCGCGAAGAATTCAACGCGGTGGCTGGCACCGGCGCAGCCGgtgcca                                                                                                                                            <  1:2220207/150‑4 (MQ=255)
              ttCGAAGCTGCGGCGCAGTTGCGTGACGAGATTGCGCAGCTGCGGGAGCGGTTGATCAGCAGCTGAAGCCGGGCGCGGGCCCTGTGGGAGCGGGCGTGCCCGCGAAGAATTCAACGCGGTGGCTGGCACCGGCTGCGCCGGTGTTcgcg                                                                                                                                         <  1:7129129/149‑1 (MQ=255)
              ttCGAAGCTGCGGCGCAGTTGCGTGACGAGATTGCGCAGCTGCGGGAGCGGTTGATCAGCAGCTGAAGCCGGGCGCGGGCCCTGTGGGAGCGGGCGTGCCCGCGAAGAATTCAACGCGGTGGCTGGCACCGGCTGCGCCGGTGTTcgcg                                                                                                                                         >  2:7129129/1‑149 (MQ=255)
              ttCGAAGCTGCGGCGCAGTTGCGTGACGAGATTGCGCAGCTGCGGGAGCGGTTGATCAGCAGCTGAAGCCGGGCGCGGGCCCTGTGGGAGCGGGCGTGCCCGCGAAGAATTCAACGCGGTGGCTGGCACCGGCTGCGCCGGTGTTCGCgg                                                                                                                                        <  2:452601/150‑1 (MQ=255)
               tCGAAGCTGCGGCGCAGTTGCGTGACGAGATTGCGCAGCTGCGGGAGCGGTTGATCGGCAGCTGAAGCCGGGCGCGGGCCCTGTGGGAGCGGGCGTGCCCGCGAAGAATTCAACGCGGTGGCTGGCACCGGCTGCGCCGGTGTTCGCggg                                                                                                                                       <  1:4673794/150‑1 (MQ=255)
                 gAAGCTGCGGCGCAGTTGCGTGACGAGATTGCGCAGCTGCGGGAGCGGTTGATCAGCAGCTGAAGCCGGGCGCGGGCCCTGTGGGAGCGGGCGTGCCCGCGAAGAATTCAACGCGGTGGCTGGCACCGGCTGCGCCGGTGTTCGCGGGCa                                                                                                                                     <  2:6873384/150‑1 (MQ=255)
                   aGCTGCGGCGCAGTTGCGTGACGAGATTGCGCAGCTGCGGGAGCGGTTGATCAGCAGCTGAAGCCGGGCGCGGGCCCTGTGGGAGCGGGCGTGCCCGCGAAGAATTCAACGCGGTGGCTGGCACCGGCTGCGCCGGTGTTCGCGGGCATg                                                                                                                                   <  2:2701308/150‑1 (MQ=255)
                    gCTGCGGCGCAGTTGCGTGACGAGATTGCGCAGCTGCGGGAGCGGTTGATCAGCAGCTGAAGCCGGGCGCGGGCCCTGTGGGAGCGGGCGTGCCCGCGAAGAATTCAACGCGGTGGCTGGCACCGGCTGCGCCGGTGTTCGCGGGCATGc                                                                                                                                  <  2:7350022/150‑1 (MQ=255)
                     cTGCGGCGCAGTTGCGTGACGAGATTGCGCAGCTGCGGGAGCGGTTGATCAGCAGCTGAAGCCGGGCGCGGGCCCTGTGGGAGCGGGCGTGCCCGCGAAGAATTCAACGCGGTGGCTGGCACCGGCTGCGCCGGTGTTCGCGGGCATGcc                                                                                                                                 <  2:7783556/150‑1 (MQ=255)
                      tGCGGCGCAGTTGCGTGACGAGATTGCGCAGCTGCGGGAGCGGTTGATCAGCAGCTGAAGCCGGGCGCGGGCCCTGTGGGAGCGGGCGTGCCCGCGAAGAATTCAACGCGGTGGCTGGCACCGGCTGCGCCGGTGTTCGCGGGCATGccc                                                                                                                                <  1:5761711/150‑1 (MQ=255)
                      tGCGGCGCAGTTGCGTGACGAGATTGCGCAGCTGCGGGAGCGGTTGATCAGCAGCTGAAGCCGGGCGCGGGCCCTGTGGGAGCGGGCGTGCCCGCGAAGAATTCAACGCGGTGGCTGGCACCGGCTGCGCCGGTGTTCGCGGGCATGccc                                                                                                                                <  1:5761760/150‑1 (MQ=255)
                       gcggcgCAGTTGCGTGACGAGATTGCGCAGCTGCGGGAGCGGTTGATCAGCAGCTGAAGCCGGGCGCGGGCCCTGTGGGAGCGGGCGTGCCCGCGAAGAATTCAACGCGGTGGCTGGCACCGGCTGCGCCGGTGTTCGCGGGCATGCCCg                                                                                                                               >  2:3950594/1‑150 (MQ=255)
                          gcgcAGTTGCGTGACGAGATTGCGCAGCTGCGGGAGCGGTTGATCAGCAGCTGAAGCCGGGCGCGGGCCCTGTGGGAGCGGGCGTGCCCGCGAAGAATTCAACGCGGTGGCTGGCACCGGCTGCGCCGGTGTTCGCGGGCATGCCCGCTc                                                                                                                            >  1:4562027/1‑150 (MQ=255)
                          gcgcAGTTGCGTGACGAGATTGCGCAGCTGCGGGAGCGGTTGATCAGCAGCTGAAGCCGGGCGCGGGCCCTGTGGGAGCGGGCGTGCCCGCGAAGAATTCAACGCGGTGGCTGGCACCGGCTGCGCCGGTGTTCGCGGGCATGCCCGCTc                                                                                                                            >  1:3636220/1‑150 (MQ=255)
                            gcAGTTGCGTGACGAGATTGCGCAGCTGCGGGAGCGGTTGATCAGCAGCTGAAGCCGGGCGCGGGCCCTGTGGGAGCGGGCGTGCCCGCGAAGAATTCAACGCGGTGGCTGGCACCGGCTGCGCCGGTGTTCGCGGGCATGCCCGCTccc                                                                                                                          <  1:4177574/150‑1 (MQ=255)
                              aGTTGCGTGACGAGATTGCGCAGCTGCGGGAGCGGTTGATCAGCAGCTGAAGCCGGGCGCGGGCCCTGTGGGGGCGGGCGTGCCCGCGAAGAAATAAAACCGGTGGCGGGCACCCGCTGCGCCCGTGTTAGCGGGCAAGCCCCCGCCcac                                                                                                                        >  1:4935419/1‑150 (MQ=255)
                                   cGTGACGAGATTGCGCAGCTGCGGGAGCGGTTGATCCGCAGCTGAAGCCGGGCGCGGGCCCTGTGGGAGCGGGCGTGCCCGCGAAGAATTCAACGCGGTGGCTGGCACCGGCTGCGCCGGTGTTCGCGGGCATGCCCGCTCCCACAGGGa                                                                                                                   <  1:485133/150‑1 (MQ=255)
                                    gTGACGCGTTTGCGCAGCTGCGGGCGCGGTTGATCAGCAGCTGAAGCCGGGCGCGGGCCCTGTGGGAGCGGGCGTGCCCGCGAAGAATTCAACGCGGTGGCTGGCACCGGCTGCGCCGGTGTTCGCGGGCATGCCCGCTCCCACAGGGAc                                                                                                                  <  1:5951409/150‑1 (MQ=255)
                                     tGACGAGATTGCGCAGCTGCGGGAGCGGTTGATCAGCAGCTGAAGCCGGGCGCGGGCCCTGTGGGAGCGGGCGTGCCCGCGAAGAATTCAACGCGGTGGCTGGCACCGGCTGCGCCGGTGTTCGCGGGCATGc                                                                                                                                  >  1:7133210/1‑133 (MQ=255)
                                      gACGAGATTGCGCAGCTGCGGGAGCGGTTGATCAGCAGCTGAAGCCGGGCGCGGGCCCTGTGGGAGCGGGCGTGCCCGCGAAGAATTCAAAGCCGTGGCTGGCACCGGCTGCGCCGGTGTTCGCGGGCATGCCCGCTCCCACAGGGACCg                                                                                                                >  1:8069883/1‑150 (MQ=255)
                                         gagaTTGCGCAGCTGCGGGAGCGGTTGATCAGCAGCTGAAGCCGGGCGCGGGCCCTGTGGGGGCGGGCGGGCCCCCGAAGAATTCAAACCCGGGGCGGGCACCCGCTGCGCCGGTGGTGGCGGGGCGGCCCCCGCCCCCCGGGAccccca                                                                                                             >  1:1755719/1‑146 (MQ=255)
                                         gagaTTGCGCAGCTGCGGGAGCGGTTGATCAGCAGCTGAAGCCGGGCGCGGGCCCTGTGGGAGCGGGCGTGCCCGCGAAGAATTCAACGCGGTGGCTGGCACCGGCTGCGCCGGTGTTCGCGGGCATGCCCGCTCCCACAgg                                                                                                                     <  2:1328377/142‑1 (MQ=255)
                                         gagaTTGCGCAGCTGCGGGAGCGGTTGATCAGCAGCTGAAGCCGGGCGCGGGCCCTGTGGGAGCGGGCGTGCCCGCGAAGAATTCAACGCGGTGGCTGGCACCGGCTGCGCCGGTGTTCGCGGGCATGCCCGCTCCCACAgg                                                                                                                     >  1:1328377/1‑142 (MQ=255)
                                            aTTGCGCAGCTGCGGGAGCGGTTGATCAGCAGCTGAAGCCGGGCGCGGGCCCTGTTGGAGCGGGCGTGCCCGCGAAGAATTAAAAGAGGTGGCTGGCACCGGCTGCGCCGGTGTTCGCGGGCATGACCGCTCCCACACGGCACGCatcat                                                                                                          >  1:8966259/1‑148 (MQ=25)
                                                  cAGCTGCGGGAGCGGTTGATCAGCAGCTGAAGCCGGGCGCGGGCCCTGTGGGAGCGGGCGTGCCCGCGAAGAATTCAACGCGGTGGCTGGCACCGGCTGCGCCGGTGTTCGCGGGCATGCCCGCTCCCACAGGGACCGCATCTGATCTGt                                                                                                    <  2:4774030/150‑1 (MQ=255)
                                                  cAGCTGCGGGAGCGGTTGATCAGCAGCTGAAGCCGGGCGCGGGCCCTGTGGGAGCGGGCGTGCCCGCGAAGAATTCAACGCGGTGGCGGGCACCGGCTGCGCCGGTGTTCGCGGGCATGCCCGCGCCCACAGGGACCGCATCTGATCTGt                                                                                                    <  2:4774092/150‑1 (MQ=255)
                                                   aGCTGCGGGAGCGGTTGATCAGCAGCTGAAGCCGGGCGCGGGCCCTGTGGGAGCGGGCGTGCCCGCGAAGAATTCAACGCGGTGGCTGGCACCGGCTGCGCCGGTGTTCGCGGGCATGCCCGCTCCCACAGGGACCGCATCTGATCTGTa                                                                                                   >  2:6387964/1‑150 (MQ=255)
                                                      ggggggCGCGGTTGACCCGCCGCTGACGCCGGGCGGGGGCCCTGGGGGCGGGGGCGTGCCCGCCAAGAATTCAACGCGGTGGCTGGCCCCGGCGGCGCGGGTGTGCGCGGGCCGGCCCGCTCCCACAGGGACCGCATCTGATCTGTAATa                                                                                                <  1:3370159/147‑1 (MQ=255)
                                                         gggAGCGGTTGATCAGCAGCTGAAGCCGGGCGCGGGCCCTGTGGGAGCGGGCGTGCCCGCGACGAATTCAACGCGGTGGCTGGCACCGGCTGCGCCGGTGTTCGCGGGCATGcc                                                                                                                                 >  2:4691485/1‑114 (MQ=255)
                                                          ggAGCGGTTGATCAGCAGCTGAAGCCGGGCGCGGGCCCTGTGGGAGCGGGCGTGCCCCCGAAGAATTCAACGCGGTGGCTGGCACCCGCGGCGCCCGTGTTAGCGGGCAATCCCCCTCCCACAGGGACCGCAACAGATCTGTAATAGcca                                                                                            >  1:8756207/1‑149 (MQ=25)
                                                                 ttGATCAGCAGCTGAAGCCGGGCGCGGGCCCTGTGGGAGCGGGCCTGCCCCCCAAGAATTCAACCCGGTGGCTGGCACCCGCTGCGCCCGTGTTCGCGGGCATGCCCCCACCCCCACGGACCCCAACTGAACAGTAAAAACCACAAAggg                                                                                     >  1:8985996/1‑150 (MQ=21)
                                                                   gATCAGCAGCTGAAGCCGGGCGCGGGCCCTGTGGGAGCGGGCGTGCCCGCGAAGAATTAAAAGCGGTGGCTGGCACCGGCTGCGCCCGTGTTAGCCGGCATGCCCCCTCCCACAGGGACCGCATCTGATCTGTAATAGCCTCAAAGGGCt                                                                                   >  1:6660363/1‑150 (MQ=37)
                                                                     accgccgcTGACGCGGGGCGGGGGCCGTGTGGGCGGGGGCGGGCCCGCGAATAATTCACCGCGGGGGCTGGCACCGGCGGCGCGGGTTTTCGCGGGCTGGCCCGCTCCCACAGGGACCGCATCTGATCTGTAATAGCCTCAATGGGCTTc                                                                                 <  1:3060655/144‑1 (MQ=255)
                                                                       agcagcTGAAGCCGGGCGCGGGCCCTGTGGGAGCGGGCGTGCCCGCGAAGAATTCAACGCGGTGGCTGGCACCGGCTGCGCCGGTGTTCGCGGGCAAGCCCGCTCCCACAGGGACCGCATCTGATCGGTAATAGCCTCAAGGGGCGTccc                                                                               >  2:3050710/1‑150 (MQ=255)
                                                                         cagcTGAAGCCGGGCGCGGGCCCTGTGGGTGCGGGCGGGCCCGCGAATAATTCAACGCGGTGGCTGGCACCGGCTGCGCCGGTGTTCGCGGGCATGCCCGCTCCCACAGGGACCGCATCTGATCTGTAATAGCCTCAATGGGCTTCCCCa                                                                             <  1:3418007/150‑1 (MQ=255)
                                                                          agcTGAAGCCGGGCGCGGGCCCTGTGGGAGCGGGCGTGCCCGCGAAGAATTCAACGCGGTGGCTGGCACCGGCTGCGCCGGTGTTCGCGGGCATGCCCGCTCCCACAGGGACCGCATCTGATCTGTAATAGCCTCAATGGGCTTCCCCAg                                                                            <  2:7248822/150‑1 (MQ=255)
                                                                            gcataGCCGGGCGCGGGCCCTGTGGGAGCGGGCGTGCCCGCGAAGAATTCAACGCGGTGGCTGGCACCGGCTGCGCCGGTGTTCGCGGGCATGCCCGCTCCCACAGGGACCGCATCTGATCTGTAATAGCCTCAATGGGCTTCCCCAGcc                                                                          <  2:2465478/146‑1 (MQ=255)
                                                                            cTGAAGCCGGGCGCGGGCCCTGTGGGAGCGGGCGTGCCCGCGAAGAATTCAACGCGGTGGCTGGCACCGGCTGCGCCGGTGTTCGCGGGCATGCCCGCTCCCACAGGGACCGCATCTGATCTGTAATAGCCTCAATGGGCTTCCCCAGcc                                                                          <  2:1758798/150‑1 (MQ=255)
                                                                               aaGCCGGGCGCGGGCCCTGTGGGAGCGGGCGTGCCCGCGAAGAATTCAACGCGGTGGCTGGCACCGGCTGCGCCGGTGTTCGCGGGCATGCCCGCTCCCACAGGGACCGCATCTGATCTGTAATAGCCTCAATGGGCTTCCCCAGCCTTc                                                                       <  2:2735215/150‑1 (MQ=255)
                                                                                aGCGGGGCGGGGGCCCTGTGGGAGCGGGCGTGCCCGCTAAGAATTCAACGCGGTGGCTGGCACGGGCTGCGCCGGTTTTCGCGGGATTGCCCGCTCCAACAGGGACCGCATCTGATCTGTAATAGCCTCAATGGGCTTCCCCATCCTTCa                                                                      <  1:2502508/150‑1 (MQ=25)
                                                                                 gCCGGGGGCGGGCCGTGGGGGAGCGGGCGTGCCCGCGACTATTTCAAGGGGGGGGCTGGCACCGGCTGCGCGGGTGTGCGCGGGCTTGCCCGCTCCCACAGGGACCGCATCTGATCTGTAATAGCCTCAATGGGCTTCCCCAGCCTTCaa                                                                     <  1:3856841/150‑1 (MQ=25)
                                                                                       cgcgGGCCCTGTGGGAGCGGGCGTGCCCGCGTAGAATTCAACGCGGTGGCTGGCTCCGGGTGCGCGGGTGTTCGCGGGGATGCCCGCTCCCACAGGGACCGCATCTGATCTGTAATAGCCTCAATGGGCTTCCCCAGCCTTCAACCCTTc                                                               <  2:4836035/150‑1 (MQ=37)
                                                                                       cgcgGGCCCTGTGGGAGCGGGCGTGCCCGCGAAGAATTCAACGCGGTGGCTGGCACCGGCTGCGCCGGTGTTCGCGGGCATGCCCGCTCCCACAGGGACCGCATCTGATCTGTAATAGCCTCAATGGGCTTCCCCAGCCTTCAACCCTTc                                                               <  2:8628621/150‑1 (MQ=255)
                                                                                          gggCCCTGTGGGAGCGGGCGTGCCCGCGAAGAATTCAACGCGGTGGCTGGCACCGGCTGCGCCGGTGTTCGCGGGCATGCCCGCTCCCACAGGGACCGCATCTGATCTGTAATAGCCTCAATGGGCTTCCCCAGCCTTCAACCCGTCCgg                                                            >  2:7221821/1‑150 (MQ=255)
                                                                                           ggCCCTGTGGGCGCGGGCGGGCCCGCGAAGAATTCACCGCGGTGGCTGGCACGGGGGGCGCCGGTGTTCGCGGGCTTGCCCGCTCCCACAGGGACCGCATCTGATCTGTAATAGCCTCAATGGGCTTCCCCAGCCTTCAACCCTTCCGGc                                                           <  1:2486580/150‑1 (MQ=255)
                                                                                                tgtgGGAGCGGGCGTGCCCGCGAAGACTTCAACGCGGTGGCTGGCACCGGCTGCGCCGGTGTTCGCGGGCATGCCCGCTCCCACAGGGACCGCATCTGAGCTGTAATAGCCTCAATGGGCTTCCCCAGCCTTCAACCCGTCCGGcccctt                                                      >  1:4080792/1‑145 (MQ=255)
                                                                                                 gtgGGAGCGGGCGTGCCCGCGAAGAATTCAACGCGGTGGCTGGCACCGGCTGCGCCGGTGTTCGCGGGCATGCCCGCTCCCACAGGGACCGCATCTGATCTGTAATAGCCTCAATGGGCTTCCCCAGCCTTCAACCCTTCCGGCAACTTg                                                     >  2:6902993/1‑150 (MQ=255)
                                                                                                   gggAGCGGGCGGGCCCGCGAAGAATTCAACGCGGTGGCTGGCACCGGCTGCGCCGGTGTTCGCGGGCATGCCCGCTCCCACAGGGACCGCATCTGATCTGTAATAGCCTCAATGGGCTTCCCCAGCCTTCAACCCTTCCGGCAACTTGCg                                                   <  2:7467247/150‑1 (MQ=255)
                                                                                                    ttcgCGGGCGTGCCCGCGAAGAATTCAACGCGGTGGCTGGCACCGGCTGCGCCGGTGTTCGCGGGCATGCCCGCTCCCACAGGGACCGCATCTGATCTGTAATAGCCTCAATGGGCTTCCCCAGCCTTCAACCCTTCCGGCAACTTGCgg                                                  >  1:5052987/4‑150 (MQ=255)
                                                                                                    ggAGGGGGCGTGCCCGCGAATAATTCAACGGGGTGGCGGGCACGGGCTGCGCCGGTGTTCGCGGGCATGCCCGCTCCCACAGGGACCGCATCTGATCTGTAATAGCCTCAATGGGCTTCCCCAGCCTTCAACCCTTCCGGCAACTTGCgg                                                  <  1:209174/150‑1 (MQ=38)
                                                                                                    ggAGCGGGCGTGCCCGCGAAGAATTCAACGCGGTGGCTGGCACCGGCTGCGCCGGTGTTCGCGGGCATGCCCGCTCCCACGGGGACCGCATCTGATCTGTAATAGCCTCAATGGGCTTCCCCAGCCTTCAACCCTTCCGGCAACTTGCgg                                                  <  1:4771600/150‑1 (MQ=255)
                                                                                                    ggAGCGGGCGTGCCCGCGAAGAATTCAACGCGGTGGCTGGCACCGGCTGCGCCGGTGTTCGCGGGCATGCCCGCTCCCACAGGGACCGCATCTGATCTGTAATAGCCTCAATGGGCTTCCCCAGCCTTCAACCCTTCCGGCAACTTGCgg                                                  >  1:7198686/1‑150 (MQ=255)
                                                                                                          ggCGTGCCCGCGAAGAATTCAACGCGGTGGCTGGCACCGGCTGCGCCGGTGTTCGCGGGCATGCCCGCTCCCACAGGGACCGCATCTGATCTGTAATAGCCTCAATGGGCTTCCCCAGCCTTCAACCCTTCCGGCAACTTGCGGGTcagc                                            >  1:3789916/1‑150 (MQ=255)
                                                                                                           gCGTGCCCGCGAAGAATTCAACGCGGTGGCTGGCACCGGCTGCGCCGGTGTTCGCGGGCATGCCCGCTCCCACAGGGACCGCATCTGATCTGTAATAGCCTCAATGGGCTTCCCCAGCCTTCAACCCTTCCGGCAACTTGCGGGTcagca                                           >  2:5001716/1‑150 (MQ=255)
                                                                                                           gCGTGCCCGCGAAGAATTCAACGCGGTGGCTGGCACCGGCTGCGCCGGTGTTCGCGGGCATGCCCGCTCCCACAGGGACCGCATCTGATCTGTAATAGCCTCAATGGGCTTCCCCAGCCTTCAACCCTTCCGGCAACTTGCGGGTcagca                                           >  1:7875248/1‑150 (MQ=255)
                                                                                                              tGCCCGCGAAGAATTCAACGCGGTGGCTGGCACCGGCTGCGCCGGTGTTCGCGGGCATGCCCGCTCCCACAGGGACCGCATCTGATCTGTAATAGCCTCAATGGGCTTCCCCAGCCTTCAACCCTTCCGGCAACTTGCGGGTCAgcagca                                        <  2:5367919/150‑1 (MQ=255)
                                                                                                               gCCCGCTAATAATTCAACGCGGTGGCTGGCACCGGCTGCGCCGGTTTTCGCGGGCATGCCAGCTCCCACAGGGACCGCATCTGATCTGTAATAGCCTCAATGGGCTTCCCCAGCATTCAACCCTTCCGGCAACTTGCGGGTCAGCAGCAc                                       <  2:6478386/150‑1 (MQ=255)
                                                                                                               gCCCGCGAAGAATTCAACGCGGTGGCTGGCACCGGCTGCGCCGGTGTTCGCGGGCATGCCCGCTCCCACAGGGACCGCATCTGATCTGTAATAGCCTCAATGGGCTTCCCCAGCCTTCAACCCTTCCGGCAACTTGCGGGTCAGCAGCAc                                       <  2:6479206/150‑1 (MQ=255)
                                                                                                                 ccGCGAAGAATTCAACGCGGTGGCTGGCACCGGCTGCGCCGGTGTTCGCGGGCATGCCCGCTCCCACAGGGACCGCATCTGATCTGTAATAGCCTCAATGGGCTTCCCCAGCCTTCAACCCTTCCGGCAACTTGCGGGTCAGCAGCACCg                                     <  1:764629/150‑1 (MQ=255)
                                                                                                                 ccGCGAAGAATTCAACGCGGTGGCTGGCACCGGCTGCGCCGGTGTTCGCGGGCATGCCCGCTCCCACAGGGACCGCATCTGATCTGTAATAGCCTCAATGGGCTTCCCCAGCCTTCAACCCTTCCGGCAACTTGCGGGTCAGCAGCACCg                                     >  2:1708124/1‑150 (MQ=255)
                                                                                                                 ccGCGAAGAATTCAACGCGGTGGCTGGCACCGGCTGCGCCGGTGTTCGCGGGCATGCCCGCTCCCACAGGGACCGCATCTGATCTGTAATAGCCTCAATGGGCTTCCCCAGCCTTCAACCCTTCCGGCAACTTGCGGGTCAGCAGCACCg                                     >  2:3352253/1‑150 (MQ=255)
                                                                                                                 ccGCGAAGAATTCAACGCGGTGGCTGGCACCGGCTGCGCCGGTGTTCGCGGGCATGCCCGCTCCCACAGGGACCGCATCTGATCTGTAATAGCCTCAATGGGCTTCCCCAGCCTTCAACCCTTCCGGCAACTTGCGGGTCAGCAGCACCg                                     >  2:582544/1‑150 (MQ=255)
                                                                                                                 ccGCGAAGAATTCAACGCGGTGGCTGGCACCGGCTGCGCCGGTGTTCGCGGGCATGCCCGCTCCCACAGGGACCGCATCTGATCTGTAATAGCCTCAATGGGCTTCCCCAGCCTTCAACCCTTCCGGCAACTTGCGGGTCAGCAGCACCg                                     >  2:582548/1‑150 (MQ=255)
                                                                                                                           ttCAACGCGGTGGCTGGCACCGGCTGCGCCGGTGTTCGCGGGCATGCCCGCTCCCACAGGGACCGCATCTGATCTGTAATAGCCTCAATGGGCTTCCCCAGCCTTCAACCCTTCCGGCAACTTGCGGGTCAGCAGCACCGCCACCATGCt                           >  1:2302475/1‑150 (MQ=255)
                                                                                                                            tCAACGCGGTGGCTGGCACCGGCTGCGCCGGTGTTCGCGGGCATGCCCGCTCCCACAGGGACCGCATCTGATCTGTAATAGCCTCAATGGGCTTCCCCAGCCTTCAACCCTTCCGGCAACTTGCGGGTCAGCAGCACCGCCACCATGCTc                          >  2:5363161/1‑150 (MQ=255)
                                                                                                                             cAACGGGGTGGCTGGCACCGGCTGCGCCGGTGTTCGCGGGCATGCCCGCTCCCACAGGGACCGCATCTGATCTGTAATAGCCTCAATGGGCTTCCCCAGCCTTCAACCCTTCCGGCAACTTGCGGGTCAGCAGCACCGCCACCATGCTCa                         <  1:1019174/150‑1 (MQ=255)
                                                                                                                              aaCGCGGTGGCTGGCACCGGCTGCGCCGGTGTTCGCGGGCATGCCCGCTCCCACAGGGACCGCATCTGATCTGTAATAGCCTCAATGGGCTTCCCCAGCCTTCAACCCTTCCGGCAACTTGCGGGTCAGCAGCACCGCCACCATGCTCAc                        >  2:1993921/1‑150 (MQ=255)
                                                                                                                              aaCGCGGTGGCTGGCACCGGCTGCGCCGGTGTTCGCGGGCATGCCCGCTCCCACAGGGACCGCATCTGATCTGTAATAGCCTCAATGGGCTTCCCCAGCCTTCAACCCTTCCGGCAACTTGCGGGTCAGCAGCACCGCCACCATGCTCAc                        >  1:1405323/1‑150 (MQ=255)
                                                                                                                               aCGCGGTGGCTGGCACCGGCTGCGCCGGTGTTCGCGGGCATGCCCGCTCCCACAGGGACCGCATCTGATCTGTAATAGCCTCAATGGGCTTCCCCAGCCTTCAACCCTTCCGGCAACTTGCGGGTCAGCAGCACCGCCACCATGCTCAcc                       >  2:779048/1‑150 (MQ=255)
                                                                                                                                 gcgGTGGCTGGCACCGGCTGCGCCGGTGTTCGCGGGCATGCCCGCTCCCACAGGGACCGCATCTGATGTGTAATAGACTCAATGGGCTTCCCCAGCCTTCAACCCTTCCGGCAACTTGCGGGTCAGCAGCACCGCCACCATGCTCACCGc                     >  1:7279036/1‑150 (MQ=255)
                                                                                                                                 gcgGTGGCTGGCACCGGCTGCGCCGGTGTTCGCGGGCATGCCCGCTCCCACAGGGACCGCATCTGATCTGTAATAGCCTCAATGGGCTTCCCCAGCCTTCAACCCTTCCGGCAACTTGCGGGTCAGCAGCACCGCCACCATGCTCACCGc                     >  2:6198597/1‑150 (MQ=255)
                                                                                                                                 gcgGTGGCTGGCACCGGCTGCGCCGGTATTCGCGGGCATGCCCGCTCCCACAGGGACCGCATCTGATCTGTAATAGCCTCAATGGGCTTCCCCAGCCTTCAACCCTTCCGGCAAc                                                        <  1:2915291/115‑1 (MQ=255)
                                                                                                                                     tggctggcCCCGGCTGCGCCGGTGTTCGCGGGCATGCCCGCTCCCACGGGGTCCGCCTCTGATCTGTAATAGCCTCAGTGGGCTTCCGCAGCCTTCAACCCTTCCGGCAACTTGCGGGTCAGCTGCACCGCCACCATGCTCTCCGCCAAc                 >  1:8352030/1‑150 (MQ=38)
                                                                                                                                      cgcgaacACCGGCGCAGCCGGTGCTCGCGGGCATGCCCGCTCCCACAGGGACCGCATCTGATCTGTAATAGCCTCAATGGGCTTCCCCAGCCTTCAACCCTTCCGGCAACTTGCGGGTCAGCAGCACCGCCACCATGCTCACCGCCAAcc                >  1:6536194/7‑150 (MQ=255)
                                                                                                                                       gctggcACCGGCTGCGCCGGTGTTCGCGGGCATGCCCGCTCCCACAGGGACCGCATCTGATCTGTAATAGCCTCAATGGGCTTCCCCAGCCTTCAACCCTTCCGGCAACTTGCGGGTCAGCAGCACCGCCACCATGCTCACCGCCAAccc               <  2:8985966/150‑1 (MQ=255)
                                                                                                                                       gctggcACCGGCTGCGCCGGTGTTCGCGGGCATGCCCGCTCCCACAGGGACCGCATCTGATCTGTAATAGCCTCAATGGGCTTCCCCAGCCTTCAACCCTTCCGGCAACTTGCGGGTCAGCAGCACCGCCACCATGCTCACCGCCAAccc               <  2:8985996/150‑1 (MQ=255)
                                                                                                                                         tggcACCGGCTGCGCCGGTGTTCGCGGGCATGCCCGCTCCCACAGGGACCGCATCTGATCTGTAATAGCCTCAATGGGCTTCCCCAGCCTTCAACCCTTCCGGCACCTTGCGGGTCAGCAGCACCGCCACCATGCTCACCGCCAACCCGa             >  2:5222290/1‑150 (MQ=255)
                                                                                                                                         tggcACCGGCTGCGCCGGTGTTCGCGGGCATGCCCGCTCCCACAGGGACCGCATCTGATCTGTAATAGCCTCAATGGGCTTCCCCAGCCTTCAACCCTTCCGGCAACTTGCGGGTCAGCAGCACCGCCACCATGCTCACCGCCAACCCGa             >  2:686075/1‑150 (MQ=255)
                                                                                                                                         tggcACCGGCTGCGCCGGTGTTCGCGGGCATGCCCGCTCCCACAGGGACCGCATCTGATCTGTAATAGCCTCAATGGGCTTCCCCAGCCTTCAACCCTTCCGGCAACTTGCGGGTCAGCAGCACCGCCACCATGCTCACCGCCAACCCGa             >  1:9114757/1‑150 (MQ=255)
                                                                                                                                         tggcACCGGCTGCGCCGGTGTTCGCGGGCATGCCCGCTCCCACAGGGACCGCATCTGATCTGTAATAGCCTCAATGGGCTTCCCCAGCCTTCAACCCTTCCGGCAACTTGCGGGTCAGCAGCACCGCCACCATGCTCACCGCCAACCCGa             >  1:8000451/1‑150 (MQ=255)
                                                                                                                                          aacACCGGCGCAGCCGGTGTTCGCGGGCATGCCCGCTCCCACAGGGACCGCATCTGATCTGTAATAGCCTCAATGGGCTTCCCCAGCCTTCAACCCTTCCGGCAACTTGCGGGTCAGCAGCACCGCCACCATGCTCACCGCCAACCCGAt            >  1:1536074/3‑150 (MQ=255)
                                                                                                                                          aacACCGGCGCAGCCGGTGTTCGCGGGCATGCCCGCTCCCACAGGGACCGCATCTGATCTGTAATAGCCTCAATGGGCTTCCCCAGCCTTCAACCCTTCCGGCAACTTGCGGGTCAGCAGCACCGCCACCATGCTCACCGCCAACCCGAt            >  2:3179564/3‑150 (MQ=255)
                                                                                                                                            cACCGGCTGCGCCGGTGTTCGCGGGCATGCCCGCTCCCACAGGGACCGCATCTGATCTGTAATAGCCTCAATGGGCTTCCCCAGCCTTCAACCCTTCCGGCAACTTGCGGGTCAGCAGCACCGCCACCATGCTCACCGCCAACCCGATcc          <  2:3367627/150‑1 (MQ=255)
                                                                                                                                            cACCGGCTGCGCCGGTGTTCGCGGGCATGCCCGCTCCCACAGGGACCGCATCTGATCTGTAATAGCCTCAATGGGCTTCCCCAGCCTTCAACCCTTCCGGCAACTTGCGGGTCAGCAGCACCGCCACCATGCTCACCGCCAACCCGATcc          <  2:8004797/150‑1 (MQ=255)
                                                                                                                                              ccGGCTGCGCCGGTGTTCGCGGGCATGCCCGCTCCCACAGGGACCGCATCTGATCTGTAATAGCCTCAATGGGCTTCCCCAGCCTTCAACCCTTCCGGCAACTTGCGGGTCAGCAGCACCGCCACCATGCTCACCGCCAACCCGATCCCg        >  1:4793555/1‑150 (MQ=255)
                                                                                                                                              ccGGCTGCGCCGGTGTTCGCGGGCATGCCCGCTCCCACAGGGACCGCATCTGATCTGTAATAGCCTCAATGGGCTTCCCCAGCCTTCAACCCTTCCGGCAACTTGCGGGTCAGCAGCACCGCCACCATGCTCACCGCCAACCCGATCCCg        <  2:8069883/150‑1 (MQ=255)
                                                                                                                                               cGGCTGCGCCGGTGTTCGCGGGCATGCCCGCTCCCACAGGGACCGCATCTGATCTGTAATAGCCTCAATGGGCTTCCCCAGCCAACAACCAATCCGGCAACTTGAGGGTCAGCAGCACCGCCACCATGCTCACCGCCAACCCGATCCCGa       >  1:6155260/1‑150 (MQ=39)
                                                                                                                                                    gcgcCGGTGTTCGCGGGCATGCCCGCTCCCACAGGGACCGCATCTGATCTGTAATAGCCTCAATGGGCTTCCCCAGCCTTCAACCCTTCCGGCAACTTGCGGGTCAGCAGCACCGCCACCATGCTCACCGCCAACCCGATCCCGACAAAg  <  2:4886739/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                    
CCCGCGACCTGGAGTTCGAAGCTGCGGCGCAGTTGCGTGACGAGATTGCGCAGCTGCGGGAGCGGTTGATCAGCAGCTGAAGCCGGGCGCGGGCCCTGTGGGAGCGGGCGTGCCCGCGAAGAATTCAACGCGGTGGCTGGCACCGGCTGCGCCGGTGTTCGCGGGCATGCCCGCTCCCACAGGGACCGCATCTGATCTGTAATAGCCTCAATGGGCTTCCCCAGCCTTCAACCCTTCCGGCAACTTGCGGGTCAGCAGCACCGCCACCATGCTCACCGCCAACCCGATCCCGACAAAG  >  NC_002947/2238791‑2239088

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 11 ≤ ATCG/ATCG < 25 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: