Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_002947 1,647,969 G→A Y265Y (TAC→TAT oprB‑II ← carbohydrate‑selective porin

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029471,647,9690GA95.6% 355.7 / 1.6 113Y265Y (TAC→TAToprB‑IIcarbohydrate‑selective porin
Reads supporting (aligned to +/- strand):  ref base G (4/1);  new base A (58/50);  total (62/51)
Fisher's exact test for biased strand distribution p-value = 3.76e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.84e-01

Only 100 of 116 total aligned reads displayed.

GACTGGTCGGACGCCAGCGAGGTGACCTGCTGCTGGGCGCCAATCCACAAGCCGTGCTTGCTCGAACTGCTGCGGTAGGCGGCGCCGCTGAGGGCGGCCGGCTGACCGTTGCTGTCCTTGAGAACATCTTGTGCCTTGGCATTACTGTAGTAGTAGCCGGCACGATATTCCCCTTTCAAGCCTTGGATACGTGGGGTCCATACCAGTTCGAACGGCATTACCGCACCCTGGGTGCCGCTGCCGCTGAGCTTGAAACCATTGCCGGATTCGAGGTTGGAAGGGTTCTGCTCGAAGA  >  NC_002947/1647820‑1648114
                                                                                                                                                     |                                                                                                                                                 
gACTGGTCGGACGCCAGCGAGGTGACCTGCTGCTGGGCGCCAATCCACAAGCCGTGCTTGCTCGAACTGCTGCGGTAGGCGGCGCCGCTGAGGGCGGCCGGCTGACCGTTGCTGTCCTTGAGAACATCTTGTGCCTTGGCATTACTgtag                                                                                                                                                   <  1:3326734/150‑1 (MQ=255)
     gTCGGACGCCAGCGAGGTGACCTGCTGCTGGGCGCCAATCCACAAGCCGTGCTTGCTCGAACTGCTGCGGTAGGCGGCGCCGCTGAGGGCGGCCGGCTGACCGTTGCTGTCCTTGAGAACATCTTGTGCCTTGGCATTACTgtaatagta                                                                                                                                              <  2:1155643/150‑1 (MQ=255)
          aCGCCAGCGAGGTGACCTGCTGCTGGGCGCCAATCCACAAGCCGTGCTTGCTCGAACTGCTGCGGTAGGCGGCGCCGCTGAGGGCGGCCGGCTGACCGTTGCTGTCCTTGAGAACATCTTGTGCCTTGGCATTACTGTAATAGTAGCCgg                                                                                                                                         <  2:469717/150‑1 (MQ=255)
            gCCAGCGAGGTGACCTGCTGCTGGGCGCCAATCCACAAGCCGTGCTTGCTCGAACTGCTGCGGTAGGCGGCGCCGCTGAGGGCGGCCGGCTGACCGTTGCTGTCCTTGAGAACATCTTGTGCCTTGGCATTACTGTAATAGTAGCCGGCa                                                                                                                                       <  1:2270977/150‑1 (MQ=255)
             ccAGCGAGGTGACCTGCTGCTGGGCGCCAATCCACAAGCCGTGCTTGCTCGAACCGCTGCGGTAGGCGGCGCCGCTGAGGGCGGCCGGCTGACCGTTGCTGTCCTTGAGAACATCTTGTGCCTTGGCATTACTGTAATAGTAGCCGGCAc                                                                                                                                      <  1:1239045/150‑1 (MQ=255)
              cAGCGAGGTGACCTGCTGCTGGGCGCCAATCCACAAGCCGTGCTTGCTCGAACTGCTGCGGTAGGCGGCGCCGCTGAGGGCGGCCGGCTGACCGTTGCTGTCCTTGAGAACATCTTGTGCCTTGGCATTACTGTAATAGTAGCCGGCACg                                                                                                                                     >  2:2347465/1‑150 (MQ=255)
                  gAGGTGACCTGCTGCTGGGCGCCAATCCACAAGCCGTGCTTGCTCGAACTGCTGCGGTAGGCGGCGCCGCTGAGGGCGGCCGGCTGACCGTTGCTGTCCTTGAGAACATCTTGTGCCTTGGCATTACTGTAATAGTAGCCGGCACGatat                                                                                                                                 >  1:56704/1‑150 (MQ=255)
                     gTGACCTGCTGCTGGGCGCCAATCCACAAGCCGTGCTTGCTCGAACTGCTGCGGTAGGCGGCGCTGCTGAGGGCGGCCGGCTGACCGTTGCTGTCCTTGAGAACATCTTGTGCCTTGGCATTACTGTAATAGTAGCCGGCACGATATTcc                                                                                                                              >  2:3600146/1‑150 (MQ=255)
                      tGACCTGCTGCTGGGCGCCAATCCACAAGCCGTGCTTGCTCGAACTGCTGCGGTAGGCGGCGCTGCTGAGGGCGGCCGGCTGACCGTTGCTGTCCTTGAGAACATCTTGTGCCTTGGCATTACTGTAATAGTAGCCGGCACGATATTccc                                                                                                                             <  1:3600146/150‑1 (MQ=255)
                          ctgctgctgGGCGCCAATCCACAAGCCGTGCTTGCTCGAACTGCTGCGGTAGGCGGCGCCGCTGAGGGCGGCCGGCTGACCGTTGCTGTCCTTGAGAACATCTTGTGCCTTGGACTTTCTGTAATGGTAGCCGGCACGATATGCCCCttt                                                                                                                         >  1:906179/1‑150 (MQ=255)
                           tgctgctgGGCGCCAATCCACAAGCCGTGCTTGCTCGAACTGCTGCGGTAGGCGGCGCCGCTGAGGGCGGCCGGCTGACCGTTGCTGTCCTTGAGAACATCTTGTGCCTTGGCATTACTGTAATAGTAGCCGGCACGATATTCCCCTTTc                                                                                                                        <  1:3048926/150‑1 (MQ=255)
                           tgctgctgGGCGCCAATCCACAAGCCGTGCTTGCTCGAACTGCTGCGGTAGGCGGCGCCGCTGAGGGCGGCCGGCTGACCGTTGCTGTCCTTGAGAACATCTTGTGCCTTGGCATTACTGTAATAGTAGCCGGCACGATATTCCCCTTTc                                                                                                                        >  2:1611499/1‑150 (MQ=255)
                               gctgGGCGCCAATCCACAAGCCGTGCTTGCTCGAACTGCTGCGGTAGGCGGCGCCGCTGAGGGCGGCCGGCTGCCCGTTGCTGTCCTTCAGATCATCTTGTGCCTTGGCATTACTGTAATAGTAGCCGGCACGATATTCCCCTTTCAAGc                                                                                                                    <  1:1314364/150‑1 (MQ=255)
                               gctgGGCGCCAATCCACAAGCCGTGCTTGCTCGAACTGCTGCGGTAGGCGGCGCCGCTGAGGGCGGCCGGCTGACCGTTGCTGTCCTTGAGAACATCTTGTGCCTTGGCATTACTGTAATAGTAGCCGGCACGATATTCCCCTTTCAAGc                                                                                                                    <  2:319358/150‑1 (MQ=255)
                                  gggNGCCAATCCACAAGCCGTGCTTGCTCGAACTGCTGCGGTAGGCGGCGCCGCTGAGGGCGGCCGGCTGACCGTTGCTGTCCTTGAGAACATCTTGTGCCTTNGNATTACTGTAATAGTANCCGGCACGATATTCCCCTTTCAAGCCtt                                                                                                                 <  1:2139598/150‑1 (MQ=255)
                                  gggCGCCAATCCACAAGCCGTGCTTGCTCGAACTGCTGCGGTAGGCGGCGCCGCTGAGGGCGGCCGGCTGACCGTTGCTGTCCTTGAGAACATCTTGTGCCTTGGCATTACTGTAATAGTAGCCGGCACGATATTCCCCTTTCAAGCCtt                                                                                                                 <  2:2177915/150‑1 (MQ=255)
                                      gcCAATCCACAAGCCGTGCTTGCTCGAACTGCTGCGGTAGGCGGCGCCGCTGAGGGCGGCCGGCTGACCGTTGCTGTCCTTGAGAACATCTTGTGCCTTGGCATTACTGTAATAGTAGCCGGCACGATATTCCCCTTTCAAGCCTTGGAt                                                                                                             >  1:417549/1‑150 (MQ=255)
                                        cAATCCACAAGCCGTGCTTGCTCGAACTGCTGCGGTAGGCGGCGCCGCTGAGGGCGGCCGGCTGACCGTTGCTGTCCTTGAGAACATCTTGTGCCTTGGCATTACTGTAATAGTAGCCGGCACGATATTCCCCTTTCAATCCTTGGATAc                                                                                                           >  2:79190/1‑150 (MQ=255)
                                        cAATCCACAAGCCGTGCTTGCTCGAACTGCTGCGGTAGGCGGCGCCGCTGAGGGCGGCCGGCTGACCGTTGCTGTCCTTGAGAACATCTTGTGCCTTGGCATTACTGTAATAGTAGCCGGCACGATATTCCCCTTTCAAGCCtt                                                                                                                 >  1:78538/1‑144 (MQ=255)
                                        cAATCCACAAGCCGTGCTTGCTCGAACTGCTGCGGTAGGCGGCGCCGCTGAGGGCGGCCGGCTGACCGTTGCTGTCCTTGAGAACATCTTGTGCCTTGGCATTACTGTAATAGTAGCCGGCACGATATTCCCCTTTCAAGCCtt                                                                                                                 <  2:78538/144‑1 (MQ=255)
                                         aaTCCACAAGCCGTGCTTGCTCGAACTGCTGCGGTAGGCGGCGCCGCTGAGGGCGGCCGGCTGACCGTTGCTGTCCTTGAGAACATCTTGTGCCTTGGCATTACTGTAATAGTAGCCGGCACGATATTCCCCTTTCAAGCCTTGGATACg                                                                                                          >  2:1323519/1‑150 (MQ=255)
                                          aTCCACAAGCCGTGCTTGCTCGAACTGCTGCGGTAGGCGGCGCCGCTGAGGGCGGCCGGCTGACCGTTGCTGTCCTTGAGAACATCTTGTGCCTTGGCATTACTGTAATAGTAGCCGGCACGATATTCCCCTTTCAAGCCTTGGATACGt                                                                                                         <  1:195580/150‑1 (MQ=255)
                                          aTCCACAAGCCGTGCTTGCTCGAACTGCTGCGGTAGGCGGCGCCGCTGAGGGCGGCCGGCTGACCGTTGCTGTCCTTGAGAACATCTTGTGCCTTGGCATTACTGTAATAGTAGCCGGCACGATATTCCCCTTTCAAGCCTTGGATACGt                                                                                                         <  2:3511955/150‑1 (MQ=255)
                                          aTCCACAAGCCGTGCTTGCTCGAACTGCTGCGGTAGGCGGCCCCGCTGAGGGCGGCCGGCTGAGCGTTGCTGTCCTTGAGAACATCTTGTGCCTTGGCATTACTGTAGTAGTAGCCGGCACGATATTCCGCTTTCAGGCCTTGGGTACGt                                                                                                         >  2:162390/1‑150 (MQ=255)
                                             cacaAGCCGTGCTTGCTCGAACTGCTGCGGTAGGCGGCGCCGCTGAGGGCGGCCGGCTGACCGTTGCTGTCCTTGAGAACATCTTGTGCCTTGGCATTACTGTAATAGTAGCCGGCACGATATTCCCCTTTCAAGCCTTGGATACGTggg                                                                                                      >  2:3027400/1‑150 (MQ=255)
                                             cacaAGCCGTGCTTGCTCGAACTGCTGCGGTAGGCGGCGCCGCTGAGGGCGGCCGGCTGACCGTTGCTGTCCTTGAGAACATCTTGTGCCTTGGCATTACTGTAATAGTAGCCGGCACGATATTCCCCTTTCAAGCCTTGGATACGTggg                                                                                                      <  2:2761065/150‑1 (MQ=255)
                                                  gCCGTGCTTGCTCGAACTGCTGCGGTAGGCGGCGCCGCGGAGGGCGGCCGGCTGACCGTTGCTGTCCTTGGGAACATCTTGTGCCTTGGCATTACTGTAATAGTAGCCGGCACGATATTCCCCTTTCAAGCCTTGGATACGTGGGGTCCa                                                                                                 >  2:723294/1‑150 (MQ=255)
                                                      tgcttgctCGAACTGCTGCGGTAGGCGGCGCCGCTGAGGGCGGCCGGCTGACCGTTGCTGTCCTTGAGAACATCTTGTGCCTTGGCATTACTGTAATAGTAGCCGGCACGATATTCCCCTTTCAAGCCTTGGATACGTggg                                                                                                      >  1:1796534/1‑141 (MQ=255)
                                                       gcttgctCGAACTGCTGCGGTAGGCGGCGCCGCTGAGGGCGGCCGGCTGACCGTTGCTGTCCTTGAGAACATCTTGTGCCTTGGCATTACTGTAATAGTAGCCGGCACGATATTCCCCTTTCAAGCCTTGGATACGTGGGGTCCATACCa                                                                                            <  2:417549/150‑1 (MQ=255)
                                                          tgctCGAACTGCTGCGGTAGGCGGCGCCGCTGAGGGCGGCCGGCTGACCGTTGCTGTCCTTGAGAACATCTTGTGCCTTGGCATTACTGTAATAGTAGCCGGCACGATATTCCCCTTTCAAGCCTTGGATACGTGGGGTCCATACCAGtt                                                                                         <  1:1611499/150‑1 (MQ=255)
                                                           gctCGAACTGCTGCGGTAGGCGGCGCCGCTGAGGGCGGCCGGCTGACCGTTGCTGTCCTTGAGAACATCTTGTGCCTTGGCATTACTGTAATAGTAGCCGGCACGATATTCCCCTTTCAAGCCTTGGATACGTggg                                                                                                      <  2:321176/136‑1 (MQ=255)
                                                           gctCGAACTGCTGCGGTAGGCGGCGCCGCTGAGGGCGGCCGGCTGACCGTTGCTGTCCTTGAGAACATCTTGTGCCTTGGCATTACTGTAATAGTAGCCGGCACGATATTCCCCTTTCAAGCCTTGGATACGTggg                                                                                                      <  2:2087279/136‑1 (MQ=255)
                                                           gctCGAACTGCTGCGGTAGGCGGCGCCGCTGAGGGCGGCCGGCTGACCGTTGCTGTCCTTGAGAACATCTTGTGCCTTGGCATTACTGTAATAGTAGCCGGCACGATATTCCCCTTTCAAGCCTTGGATACGTggg                                                                                                      >  1:321176/1‑136 (MQ=255)
                                                           gctCGAACTGCTGCGGTAGGCGGCGCCGCTGAGGGCGGCCGGCTGACCGTTGCTGTCCTTGAGAACATCTTGTGCCTTGGCATTACTGTAATAGTAGCCGGCACGATATTCCCCTTTCAAGCCTTGGATACGTggg                                                                                                      >  1:2087279/1‑136 (MQ=255)
                                                             tCGAACTGCTGCGGTAGGCGGCGCCGCTGAGGGCGGCCGGCTGACCGTTGCTGTCCTTGAGAACATCTTGTGCCTTGGCATTACTGTAATAGTAGCCGGCACGATATTCCCCTTTCAAGCCTTGGATACGTGGGGTCCATAc                                                                                              >  1:2595719/1‑142 (MQ=255)
                                                             tCGAACTGCTGCGGTAGGCGGCGCCGCTGAGGGCGGCCGGCTGACCGTTGCTGTCCTTGAGAACATCTTGTGCCTTGGCATTACTGTAATAGTAGCCGGCACGATATTCCCCTTTCAAGCCTTGGATACGTGGGGTCCATAc                                                                                              <  2:2595719/142‑1 (MQ=255)
                                                             tCGAACTGCTGCGGTAGGCGGCGCCGCTGAGGGCGGCCGGCTGACCGTTGCTGTCCTTGAGAACATCTTGTGCCTTGGCATTACTGTAATAGTAGCCGGCACGATATTCCCCTTTCAAGCCTTGGATACGTGGGGTCCATACCAGTTCGa                                                                                      >  1:3538260/1‑150 (MQ=255)
                                                               gAACTGCTGCGGTAGGCGGCGCCGCTGAGGGCGGCCGGCTGACCGTTGCTGTCCTTGAGAACATCTTGTGCCTTGGCATTACTGTAATAGTAGCCGGCACGATATTCCCCTTTCAAGCCTTGGATACGTGGGGTCCATACCAGTTCGAAc                                                                                    <  2:2972338/150‑1 (MQ=255)
                                                               gAACTGCTGCGGTAGGCGGCGCCGCTGAGGGCGGCCGGCTGACCGTTGCTGTCCTTGAGAACATCTTGTGCCTTGGCATTACTGTAATAGTAGCCGGCACGATATTCCCCTTTCAAGCCTTGGATACGTGGGGTCCATACCAGTTCGAAc                                                                                    >  1:3293988/1‑150 (MQ=255)
                                                                   tgctgcGGTAGGCGGCGCCGCTGAGGGCGGCCGGCTGACCGTTGCTGTCCTTGAGAACATCTTGTGCCTTGGCATTACTGTAGTAGTAGc                                                                                                                                            >  2:567763/1‑90 (MQ=255)
                                                                      tgcGGTAGGCGGCGCCGCTGAGGGCGGCCGGCTGACCGTTGCTGTCCTTGAGAACATCTTGTGCCTTGGCATTACTGTAATAGTAGCCGGCACGATATTCCCCTTTCAAGCCTTGGATACGTGGGGTCCATACCCGGTGGAACGGGATTa                                                                             >  1:1154307/1‑150 (MQ=255)
                                                                      tgcGGTAGGCGGCGCCGCTGAGGGCGGCCGGCTGACCGTTGCTGTCCTTGAGAACATCTTGTGCCTTGGCATTACTGTAATAGTAGCCGGCACGATATTCCCCTTTCAAGCCTTGGATACGTGGGGTCCATACCAGTTCGAACGGCATTa                                                                             >  1:1059702/1‑150 (MQ=255)
                                                                       gcGGTAGGCGGCGCCGCTGAGGGCGGCCGGCTGACCGTTGCTGTCCTTGAGAACATCTTGTGCCTTGGCATTACTGTAATAGTAGCCGGCACGATATTCCCCTTTCAAGCCTTGGATACGTGGGGTCCATACCAGTTCGAACGGCATTAc                                                                            >  2:3470091/1‑150 (MQ=255)
                                                                           tAGGCGGCGCCGCTGAGGGCGGCCGGCTGACCGTTGCTGTCCTTGAGAACATCTTGTGCCTTGGCATTACTGTAATAGTAGCCTGCACGATATTCCCCTTTCAAGCCTTCGATACGTGGGGCCCACCCCAGTTCGAACGGCATTACCgag                                                                        >  2:215640/1‑148 (MQ=255)
                                                                            aGGCGGCGCCGCTGAGGGCGGCCGGCTGACCGTTGCTGTCCTTGAGAACATCTTGTGCCTTGGCATTACTGTAATAGTAGCCGGCACGATATTCCCCTTTCAAGCCTTGGATACGTGGGGTCCATACCAGTTCGAACGGCATTACCGCAc                                                                       >  2:408233/1‑150 (MQ=255)
                                                                              gcggcgCCGCTGAGGGCGGCCGGCTGACCGTTGCTGTCCTTGAGAACATCTTGTGCCTTGGCATTGCTGTAATAGTAGCCGGCACGATATTCCCCTTTCAAGCCTTGGATACGTGGGGTCCATACCAGGTCGAACCGCATGACCGCAccc                                                                     >  2:1331403/1‑150 (MQ=255)
                                                                                ggcgCCGCTGAGGGCGGCCGGCTGACCGTTGCTGTCCTTGAGAACATCTTGTGCCTTGGCATTACTGTAATAGTAGCCGGCACGATATTCCCCTTTCAAGCCTTGGATACGTGGGGTCCATACCAGTTCGAACGGCATTACCGCAcc                                                                      >  1:161231/1‑147 (MQ=255)
                                                                                ggcgCCGCTGAGGGCGGCCGGCTGACCGTTGCTGTCCTTGAGAACATCTTGTGCCTTGGCATTACTGTAATAGTAGCCGGCACGATATTCCCCTTTCAAGCCTTGGATACGTGGGGTCCATACCAGTTCGAACGGCATTACCGCAcc                                                                      <  2:161231/147‑1 (MQ=255)
                                                                                 gcgcCGCTGAGGGCGGCCGGCTGACCGTTGCTGTCCTTGAGAACATCTTGTGCCTTGGCATTACTGTAATAGTAGCCGGCACGATATTCCCCTTTCAAGCCTTGGATACGTGGGGTCCATACCAGTTCGAACGGCATTACCGCACCCTgg                                                                  <  2:3324509/150‑1 (MQ=255)
                                                                                 gcgcCGCTGAGGGCGGCCGGCTGACCGTTGCTGTCCTTGAGAACATCTTGTGCCTTGGCATTACTGTAATAGTAGCCGGCACGATATTCCCCTTTCAAGCCTTGGATACGTGGGGTCCATACCAGTTCGAACGGCATTACCGCACCCTgg                                                                  >  1:1301736/1‑150 (MQ=255)
                                                                                                gccggcTGACCGTTGCTGTCCTTGAGAACATCTTGTGCCTTGGCATTACTGTAATAGTAGCCGGCACGATATTCCCCTTTCAAGCCTTGGATACGTGGGGTCCATAc                                                                                              >  1:3364789/1‑107 (MQ=255)
                                                                                                gccggcTGACCGTTGCTGTCCTTGAGAACATCTTGTGCCTTGGCATTACTGTAATAGTAGCCGGCACGATATTCCCCTTTCAAGCCTTGGATACGTGGGGTCCATAc                                                                                              <  2:3364789/107‑1 (MQ=255)
                                                                                                  cggcTGACCGTTGCTGTCCTTGAGAACATCTTGTGCCTTGGCATTACTGTAATAGTAGCCGGCACGATATTCTCCTTTCAAGCCTTGGATACATGGGGTCCCTACCAGTTCGCACGGCATTCCCGCACCCTGCATGCCGCCGCCGCTGAg                                                 >  2:424959/1‑150 (MQ=255)
                                                                                                    gcTGACCGTTGCTGTCCTTGAGAACATCTTGTGCCTTGGCATTACTGTAATAGTAGCCGGCACGATATTCCCCTTTCAAGCCTTGGATACGTGGGGTCCATACCAGTTCGAACGGc                                                                                 <  1:895572/116‑1 (MQ=255)
                                                                                                    gcTGACCGTTGCTGTCCTTGAGAACATCTTGTGCCTTGGCATTACTGTAATAGTAGCCGGCACGATATTCCCCTTTCAAGCCTTGGATACGTGGGGTCCATACCAGTTCGAACGGCATTACCGCACCCTGGGTGCCGCTGCCGCTGAGCt                                               >  2:1631216/1‑150 (MQ=255)
                                                                                                        aCCNTTGCTGTCCTTGAGAACATCTTGTGCCTTGGCATTACTGTAATAGTAGCCGGCACGATATTCCCCTTTCAAGCCTTGGATACGTGGGGTCCATACCAGTTCGAACGGCATTACCGCACCCTGGGCGCCGCTGCCGCTGACCTTgca                                           >  2:2243763/1‑148 (MQ=255)
                                                                                                        aCCGTTGCTGTCCTTGAGAACATCTTGTGCCTTGGCATTACTGTAATAGTAGCCGGCACGATATTCCCCTTTCCAGCCTTGGGTACGTGGGGGCCGTATCGGTTCGTAGGGTGTTTTTGCCGCCTGGGGGCCGCTGCCGCTGCGCTTggg                                           >  1:2370361/1‑148 (MQ=255)
                                                                                                         ccGTTGCTGTCCTTGAGAACATCTTGTGCCTTGGCATTACTGTAATAGTAGCCGGCACGATATTCCCCTTTCAAGcc                                                                                                                   <  1:1768958/77‑1 (MQ=255)
                                                                                                         ccGTTGCTGTCCTTGAGAACATCTTGTGCCTTGGCATTACTGTAATAGTAGCCGGCACGATATTCCCCTTTCAAGCCTTGGATACGTGGGGTCCATACCAGTTCGAACGGCATTACCGCACCCTgg                                                                  <  1:1579052/126‑1 (MQ=255)
                                                                                                         ccGTTGCTGTCCTTGAGAACATCGTGTGCCTTGGCATTACTGTAATAGTAGCCGGCACGATATTCCCCTTTCAAGcc                                                                                                                   >  2:1768958/1‑77 (MQ=255)
                                                                                                          cGTTGCTGTCCTTGAGAACATCTTGTGCCTTGGCATTACTGTAATAGTAGCCGGCACGATATTCCCCTTTCAAGCCTTGGATACGTGGGGTCCATACCAGTTCGAACGGCATTACCGCACCCTGGGTGCCGCTGCCGCTGAGCTTGAAAc                                         >  1:2523327/1‑150 (MQ=255)
                                                                                                           gTTGCTGTCCTTGAGAACATCTTGTGCCTTGGCATTACTGTAATAGTAGCCGGCACGATATTCCCCTTTCAAGCCTTGGATACGTGGGGTCCATACCAGTTCGAACGGCATTACCGCACCCTGGGTGCCGCTGCCGCTGAGc                                                >  2:143115/1‑142 (MQ=255)
                                                                                                           gTTGCTGTCCTTGAGAACATCTTGTGCCTTGGCATTACTGTAATAGTAGCCGGCACGATATTCCCCTTTCAAGCCTTGGATACGTGGGGTCCATACCAGTTCGAACGGCATTACCGCACCCTGGGTGCCGCTGCCGCTGAGc                                                <  1:143115/142‑1 (MQ=255)
                                                                                                            ttGCTGTCCTTGAGAACATCTTGTGCCTTGGCATTACTGTAATAGTAGCCGGCACGATATTCCCCTTTCAAGCCTTGGGTTCGTGGGGTCCCTACCAGTTCGAGCGGCATTGCCGCACCGCGGGGGAAGCGGCCGCCGAGATTGAAACCa                                       >  1:726801/1‑150 (MQ=255)
                                                                                                               cTGTCCTTGAGAACATCTTGTGCCTTGGCATTACTGTAATAGTAGCCGGCACGATATTCCCCTTTCAAGCCTTGGATACGTGGGGTCCATACCAGTTCGAACGGCATTACCGCACCCTGGGTGCCGCTGCCGCTGAGCTTGAAACCATTg                                    <  1:3027400/150‑1 (MQ=255)
                                                                                                                tGTCCTTGAGAACATCTTGTGCCTTGGCATTACTGTAATAGTAGCCGGCACGATATTCCCCTTTCAAGCCTTGGATACGTGGGGTCCATACCAGTTCGAACGGCATTACCGCACCCTGGGTGCCGCTGCCGCTGAGCTTGAAACCATTGc                                   <  2:3286977/150‑1 (MQ=255)
                                                                                                                  tCCTTGAGAACATCTTGTGCCTTGGCATTACTGTAATAGTAGCCGGCACGATATTCCCCTTTCAAGCCtt                                                                                                                 <  2:1283879/70‑1 (MQ=255)
                                                                                                                  tCCTTGAGAACATCTTGTGCCTTGGCATTACTGTAATAGTAGCCGGCACGATATTCCCCTTTCAAGCCtt                                                                                                                 >  1:1283879/1‑70 (MQ=255)
                                                                                                                    cTTGAGAACATCTTGTGCCTTGGCATTACTGTAGTAGTAGCCGGCACGATATTCCCCTTTCAAGCCTTGGATACGTGGGGTCCATACCAGTTCGAACGGCATTACCGCACGATGTGGGCCGCTGCCGCTGAGCTTGAAACCATCGCCGGa                               >  1:2882354/1‑150 (MQ=255)
                                                                                                                    cTTGAGAACATCTTGTGCCTTGGCATTACTGTAATAGTAGCCGGCACGATATTCCCCTTTCAAGCCTTGGATACGTGGGGTCCATACCAGTTCGAACGGCATTACCGCACCCTGGGTGCCGCTGCCGCTGAGCTTGAAACCATTGCCGGa                               >  1:3235537/1‑150 (MQ=255)
                                                                                                                     ttGAGAACATCTTGTGCCTTGGCATTACTGTAGTAGTAGCCGGCACGATATTCCCCTTTCAAGCCTTGGATACGTGGGGTCCATACCAGTTCGAACGGCATTACCGCACCCTGGGTGCCGCTGCCGCTGAGCTTGAAACCATTGCCGGAt                              >  1:2475246/1‑150 (MQ=255)
                                                                                                                         gaACATCTTGTGCCTTGGCATTACTGTAATAGTAGCCGGCACGATATTCCCCTTTCAAGCCTTGGATCCGTGGGAATCGTACCATGTCGAGCGGCATTACCGCACCCTGGGTGCCGCTGCCGCTGAGCTTGAAACCATTGCCGGATTCGa                          >  1:1308977/1‑150 (MQ=255)
                                                                                                                         gaACATCTTGTGCCTTGGCATTACTGTAATAGTAGCCGGCACGATATTCCCCTTTCAAGCCTTGGATACGTGGGGTCCATACCAGTTCGAACGGCATTACCGCACCCTGGGTGCCGCTGCCGCTGAGCTTGAAACCATTGCCGGATTCGa                          <  1:2347465/150‑1 (MQ=255)
                                                                                                                            cATCTTGTGCCTTGGCATTACTGTAATAGTAGCCGGCACGATATTCCCCTTTCAAGCCTTGGATACGTGGGGTCCATACCAGTTCGAACGGCATTACCGCACCCTGGGTGCCGCTGCCGCTGAGCTTGAAACCATTGCCGGATTCGAGGt                       >  1:2181265/1‑150 (MQ=255)
                                                                                                                             aTCTTGTGCCTTGGCATTACTGTAATAGTAGCCGGCACGATATTCCCCTTTCAAGCCTTGGATACGTGGGGTCCATACCAGTTCGAACGGCATTACCGCACCCTGGGTGCCGCTGCCGCTGAGCTTGAAACCATTGCCGGATTCGAGGtt                      >  2:3544575/1‑150 (MQ=255)
                                                                                                                                 tgtgCCTTGGCATTACTGTAATAGTAGCCGGCACGATATTCCCCTTTCAAGCCTTGGATACGTGGGGTCCATACCAGTTCGAACGGCATTACCGCACCCTGGGTGCCGCTGCCGCTGAGCTTGAAACCATTGCCGGATTCGAGGTTGGaa                  <  2:3538260/150‑1 (MQ=255)
                                                                                                                                   tgCCTTGGCATTACTGTAATAGTAGCCGGCACGATATTCCCCTTTCAAGCCTTGGATACGTGGGGTCCATACCAGTTCGAACGGCATTACCGCACCCTGGGTGCCGCTGCCGCTGAGCTTGAAACCATTGCCGGATTCGAGGTTGGAAgg                <  1:3071204/150‑1 (MQ=255)
                                                                                                                                     ccTTGGCATTACTGTAATAGTAGCCGGCACGATATTCCCCTTTCAAGCCTTGGATACGTGGGGTCCATACCAGTTCGAACGGCATTACCGCACCCTGGGTGCCGCTGCCGCTGAGCTTGAAACCATTGCCGGCTTCGAGGTGGGCAgggg              >  2:2321875/1‑149 (MQ=255)
                                                                                                                                      cTTGGCATTACTGTAATAGTAGCCGGCACGATATTCCCCTTTCAAGCCTTGGATACGTGGGGTCCATACCAGTTCGAAc                                                                                    <  2:860087/79‑1 (MQ=255)
                                                                                                                                      cTTGGCATTACTGTAATAGTAGCCGGCACGATATTCCCCTTTCAAGCCTTGGATACGTGGGGTCCATACCAGTTCGAAc                                                                                    >  1:860087/1‑79 (MQ=255)
                                                                                                                                      cTTGGCATTACTGTAATAGTAGCCGGCACGATATTCCCCTTTCAAGCCTTGGATACGTGGGGTCCATACCAGTTCGAACGGCATTACCGCACCCTGGGTGCCGCTGCCGCTGAGCTTGAAACCATTGCCGGATTCGAGGTTGGAAg                 <  1:1328800/146‑1 (MQ=255)
                                                                                                                                      cTTGGCATTACTGTAATAGTAGCCGGCACGATATTCCCCTTTCAAGCCTTGGATACGTGGGGTCCATACCAGTTCGAACGGCATTACCGCACCCTGGGTGCCGCTGCCGCTGAGCTTGAAACCATTGCCGGATTCGAGGTTGGAAg                 >  2:1328800/1‑146 (MQ=255)
                                                                                                                                      cTTGGCATTACTGTAATAGTAGCCGGCACGATATTCCCCTTTCAAGCCTTGGATACGTGGGGTCCATACCAGTTCGAACGGCATTACCGCACCCTGGGTGCCGCTGCCGCTGAGCTTGAAACCATTGCCGGATTCGAGGTTGGAAGGGtt             <  2:2506625/150‑1 (MQ=255)
                                                                                                                                      cTTGGCATTACTGTAATAGTAGCCGGCACGATATTCCCCTTTCAAGCCTTGGATACGTGGGGTCCATACCAGTTCGAACGGCATTACCGCACCCTGGGTGCCGCTGCCGCTGAGCTTGAAACCATTGCCGGATTCGAGGTTGGAAGGGtt             <  1:408233/150‑1 (MQ=255)
                                                                                                                                       ttGGCATTACTGTAATAGTAGCCGGCACGATATTCCCCTTTCAAGCCTTGGATACGTGGGGTCCATACCAGTTCGAACGGCATTACCGCACCCTGGGTGCCGCTGCCGCTGAGCTTGAAACCATTGCCGGATTCGAGGTTGGAAGGGTTc            >  2:1911215/1‑150 (MQ=255)
                                                                                                                                       ttGGCATTACTGTAATAGTAGCCGGCACGATATTCCCCTTTCAAGCCTTGGATACGTGGGGGCCATACCAGTTCGAACGGCATTACCGCACCCCGGGTGCCGCTGCCGATGGGCTTGAAACCATTGCCGCATCCGAGGTTGGAAGGGTTc            >  1:920201/1‑150 (MQ=255)
                                                                                                                                        tGGCATTACTGTAATAGTAGCCGGCACGATATTCCCCTTTCAAGCCTTGGATACGTGGGGTCCATACCAGTTCGAACGGCATTACCTCACCCTGGGTGCCGCTGCCGCTGAGCTTGAAACCATTGCCGGATTAGAGGTTGGAAGGGTTCt           <  2:1646184/150‑1 (MQ=255)
                                                                                                                                         ggCATTACTGTAATAGTAGCCGGCACGATATTCCCCTTTCAAGCCTTGGATACGTGGGGTCCATACCAGTGCGAACGGCATTACCGCACCCTGGGTGCCGCTGCCGCTGAGCTTGAAACCATTGCCGGATTCGCGGGTGGAAGGGTTCTg          >  1:2797834/1‑150 (MQ=255)
                                                                                                                                          gCATTACTGTAATAGTAGCCGGCACGATATTCCCCTTTCaa                                                                                                                      <  1:447085/41‑1 (MQ=255)
                                                                                                                                          gCATTACTGTAATAGTAGCCGGCACGATATTCCCCTTTCaa                                                                                                                      >  2:447085/1‑41 (MQ=255)
                                                                                                                                          gCATTACTGTAATAGTAGCCGGCACGATATTCCCCTTTCAAGCCTTGGATACGTGGGGTCCATACCAGTTCGAACGGCATTACCGCACCCTGGGTGCCGCTGCCGCTGAGCTTGAAACCATTGCCGGATTCGAGGTTGGAAGGGTTCTGc         >  2:2661742/1‑150 (MQ=255)
                                                                                                                                            aTTACTGTAATAGTAGCCGGCACGATATTCCCCTTTCAAGCCTTGGATACGTGGGGTCCATACCAGTTCGAACGGCATTACCGCACCCTGGGTGCCGCTGCCGCTGAGCTTGAAACCATTGCCGGATTCGAGGTTGGAAGGGTTCTGCTc       >  1:3312983/1‑150 (MQ=255)
                                                                                                                                             ttACTGTAATAGTAGCCGGCACGATATTCCCCTTTCAAGCCTTGGATACGTGGGGTCCATACCAGTTCGAACGGCATTACCGCACCCTGGGTGCCGCTGCCGCTGAGCTTGAAACCATTGCCGGATTCGAGGTTGGAAGGGTTCTGCTCg      >  2:3069879/1‑150 (MQ=255)
                                                                                                                                               aCTGTAATAGTAGCCGGCACGATATTCCCCTTTCAAGCCTTGGATACGTGGGGTCCATACCAGTTCGAACGGCATTACCGCACCCTGGGTGCCGCTGCCGCTGAGCTTGAAACCATTGCCGGATTCGAGGTTGGAAGGGTTCTGCTCGaa    >  1:2227662/1‑150 (MQ=255)
                                                                                                                                               aCTGTAATAGTAGCCGGCACGATATTCCCCTTTCAAGCCTTGGATACGTGGGGTCCATACCAGTTCGAACGGCATTACCGCACCCTGGGTGCCGCTGCCGCTGAGCTTGAAACCATTGCCGGATTCGAGGTTGGAAGGGTTCTGCTCGaa    <  2:56704/150‑1 (MQ=255)
                                                                                                                                                 tGTAATAGTAGCCGGCACGATATTCCCCTTTCAAGCCTTGGATACGTGGGGTCCATAc                                                                                              <  1:2139576/58‑1 (MQ=255)
                                                                                                                                                 tGTAATAGTAGCCGGCACGATATTCCCCTTTCAAGCCTTGGATACGTGGGGTCCATAc                                                                                              >  2:2139576/1‑58 (MQ=255)
                                                                                                                                                 tGTAATAGTAGCCGGCACGATATTCCCCTTTCAAGCCTTGGATACGTGGGGTCCATACCAGTTCGAACGGCATTACCGCACCCTGGGTGCCGCTGCCGCTGAGCTTGAAACCATTGCCGGATTCGAGGTTGGaa                  <  1:1347449/134‑1 (MQ=255)
                                                                                                                                                 tGTAATAGTAGCCGGCACGATATTCCCCTTTCAAGCCTTGGATACGTGGGGTCCATACCAGTTCGAACGGCATTACCGCACCCTGGGTGCCGCTGCCGCTGAGCTTGAAACCATTGCCGGATTCGAGGTTGGaa                  >  2:1347449/1‑134 (MQ=255)
                                                                                                                                                 tGTAATAGTAGCCGGCACGATATTCCCCTTTCAAGCCTTGGATACGTGGGGTCCATACCAGTTCGAACGGCATTACCGCACCCTGGGTGCCGCTGCCGCTGAGCTTGAAACCATTGCCGGATTCGAGGTTGGAAGGGTTCTGCTCGAAGa  >  2:20102/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                 
GACTGGTCGGACGCCAGCGAGGTGACCTGCTGCTGGGCGCCAATCCACAAGCCGTGCTTGCTCGAACTGCTGCGGTAGGCGGCGCCGCTGAGGGCGGCCGGCTGACCGTTGCTGTCCTTGAGAACATCTTGTGCCTTGGCATTACTGTAGTAGTAGCCGGCACGATATTCCCCTTTCAAGCCTTGGATACGTGGGGTCCATACCAGTTCGAACGGCATTACCGCACCCTGGGTGCCGCTGCCGCTGAGCTTGAAACCATTGCCGGATTCGAGGTTGGAAGGGTTCTGCTCGAAGA  >  NC_002947/1647820‑1648114

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