Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,183,468 G→A 8.4% A308V (GCA→GTA)  mltF ← membrane‑bound lytic murein transglycosylase F

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029471,183,4680GA 8.4% 200.4 / 6.5 71A308V (GCA→GTA) mltFmembrane‑bound lytic murein transglycosylase F
Reads supporting (aligned to +/- strand):  ref base G (38/27);  new base A (2/4);  total (40/31)
Fisher's exact test for biased strand distribution p-value = 3.93e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.77e-01

TGGCGCCGCCCTGGATGCTCTGCCTCGGGTCCAGTCGGTTCGACACGCCCATGGCCTGGGCGGTGCGCTGGGTCAGCATCATCAGGCCGCGCACGCCGGTCTTGGAGGTGACTTCCGGCTGCCACATCGATTCCTGGTAGCCAATGGCTGCCAGCAGGCGCCAGTCGACCTGCTCGACCTTGGCGTAGCTCTTGAAGTGCTTCTCGTACTTGGGCAGGCGCTGTTGCAGGTGCTGGGCGAAGGTGTAGGCGCCGACATAGCCAAGTACGTCGACATGGCCGTAGTAGCGGTCCTTCAAT  >  NC_002947/1183319‑1183617
                                                                                                                                                     |                                                                                                                                                     
tGGCGCCGCCCTGGATGCTCTGCCTCGGGTCCAGTCGGTTCGACACGCCCATGGCCTGGGCGGTGCGCTGGGTCAGCATCATCAGGCCGCGCACGCCGGTCTTGGAGGTGACTTCCGGCTGCCACATCGATTCCTGGTCGCCAACGCCTg                                                                                                                                                       <  2:19231/150‑1 (MQ=255)
   cgccgcCCTGGATGCTCTGCCTCGGGTCCAGTCGGTTCGACACGTCAATGGCATGGGGGCGGCGCTGACGCAGCATCATCAGGCCGCGCACGCCGGTCTTGGAGGTGACTTCCGGCTGCCACATCGATTCCTGGCAGCCATTGGCTACca                                                                                                                                                    <  1:525972/150‑1 (MQ=255)
   cgccgcCCTGGATGCTCTGCCTCGGGTCCAGTCGGTTCGACACGCCCATGGCCTGGGCGGTGCGCTGGGTCAGCATCATCAGGCCGCGCACGCCGGTCTTGGAGGTGACTTCCGGCTGCCACATCGATTCCTGGCAGCCATTGGCTACca                                                                                                                                                    <  1:1132774/150‑1 (MQ=255)
                  tctGCCTCGGGTCCAGTCGGTTCGACACGCCCATGGCCTGGGCGGTGCGCTGGGTCAGCATCATCAGGCCGCGCACGCCGGTCTTGGAGGTGACTTCCGGCTGCCACATCGATTCCTGGCAGCCATTGGCTACca                                                                                                                                                    <  2:632493/135‑1 (MQ=255)
                  tctGCCTCGGGTCCAGTCGGTTCGACACGCCCATGGCCTGGGCGGTGCGCTGGGTCAGCATCATCAGGCCGCGCACGCCGGTCTTGGAGGTGACTTCCGGCTGCCACATCGATTCCTGGCAGCCATTGGCTACca                                                                                                                                                    >  1:632493/1‑135 (MQ=255)
                  tctGCCTCGGGTCCAGTCGGTTCGACACGCCCATGGCCTGGGCGGTGCGCTGGGTCAGCATCATCAGGCCGCGCACGCCGGTCTTGGAGGTGACTTCCGGCTGCAACATCGATTCCTGGTAGCCAATGGCTGCCAGCAGGCGCCAGTCGa                                                                                                                                     >  1:1373266/1‑150 (MQ=255)
                   ctGCCTCGGGTCCAGTCGGTTCGACACGCCCATGGCCTGGGCGGTGCGCTGGGTCAGCATCATCAGGCCGCGCACGCCGGTCTTGGAGGTGACTCCCGCATTTTACATCGATTCCTGGTAGCCAATGGCTGCCAGCAGGCGCCAGTCGAc                                                                                                                                    <  1:1327532/150‑1 (MQ=255)
                      ccTCGGGTCCAGTCGGTTCGACACGCCCATGGCCTGGGCGGTGCGCTGGGTCAGCATCATCAGGCCGCGCACGCCGGTCTTGGAGGTGACTTCCGGCTGCCACATCGATTCCTGGTAGCCAATGGCTGCCAGCAGGCGCCAGTCGACCTg                                                                                                                                 >  2:1157098/1‑150 (MQ=255)
                        tCGGGTCCAGTCGGTTCGACACGCCCATGGCCTGGGCGGTGCGCTGGGTCAGCATCATCAGGCCGCGCACGCCGGTCTTGGAGGTGACTTCCGGCGGCCACATCGATTCCTGGTAGCCAATGGCTGCCAGCAGGCGCCAGTCGACCTGCt                                                                                                                               >  1:96816/1‑150 (MQ=255)
                        tCGGGTCCAGTCGGTTCGACACGCCCATGGCCTGGGCGGGGCGCTGGGTCAGCATCATCAGGCCGCGCACGCCGGTCTTGGAGGTGACTTCCGGCTGCCACATCGATTCCTGGTAGCCAATGGCTGCCAGCAGGCGCCAGTCGACCTGCt                                                                                                                               >  1:114688/1‑150 (MQ=255)
                                aGTCGGTTCGACACGCCCATGGCCTGGGCGGTGCGCTGGGTCAGCATCATCAGGCCGCGCACGCCGGTCTTGGAGGTGACTTCCGGCTGCCACATCGATTCCTGGTAGCCAATGGCTGCCAGCAGGCGCCAGTCGACCTGCTCGACCTTg                                                                                                                       <  2:199213/150‑1 (MQ=255)
                                aGTCGGTTCGACACGCCCATGGCCTGGGCGGGGCGCTGGGTCAGCATTATCAGGCGGCGGACGCCGGTCTTGGAGGTGACTTCCGGCTGCCACATCGATTCCTGGTAGCCAATGGCTGCCAGCAGGCGCCAGTCGACCTGCTCGACCTTg                                                                                                                       <  2:685820/150‑1 (MQ=255)
                                 gTCGGTTCGACACGCCCATGGCCTGGGCGGTGCGCTGGGTCAGCATCATCAGGCCGCGCACGCCGGTCTTGGAGGTGACTTCCGGCTGCCACATCGATTCCTGGTAGCCAATGGCTGCCAGCAGGCGCCAGTCGACCTGCTCGACCTTgg                                                                                                                      >  1:1322081/1‑150 (MQ=255)
                                 gTCGGTTCGACACGCCCATGGCCTGGGCGGTGCGCTGGGTCAGCATCATCAGGCCGCGCACGCCGGTCTTGGAGGTGACTTCCGGCTGCCACATCGATTCCTGGTAGCCAATGGCTGCCAGCAGGCGCCAGTCGACCTGCTCGACCTTgg                                                                                                                      <  1:919626/150‑1 (MQ=255)
                                  tCGGTTCGACACGCCCATGGCCTGGGCGGTGCGCTGGGTCAGCATCATCAGGCCGCGCACGCCGGTCTTGGAGGTGACTTCCGGCTGCCACATCGATTCCTGGTAGCCAATGGCTGCCAGCAGGCGCCAGTCGACCTGCTCGACCTTGGc                                                                                                                     <  2:1404823/150‑1 (MQ=255)
                                        cGACACGCCCATGGCCTGGGCGGTGCGCTGGGTCAGCATCATCAGGCCGCGCACGCCGGTCTTGGAGGTGACTTCCGGCTGCCACATCGATTCCTGGTAGCCAATGGCTGCCAGCAGGCGCCAGTCGACCTGCTCGACCTTGGCGTAGct                                                                                                               >  2:1340206/1‑150 (MQ=255)
                                                      ccTGGGCGGTGCGCTGGGTCAGCATCATCAGGCCGCGCACGCCGGTCTTGGAGGTGACTTCCGGCTGCCACATCGATTCCTGGTAGCCAATGGCTGCCAGCAGGCGCCAGTCGACCTGCTCGACCTTGGCGTAGCTCTTGAAGTGCTtct                                                                                                 <  1:1018316/150‑1 (MQ=255)
                                                       cTGGGCGGTGCGCTGGGTCAGCATCATCAGGCCGCGCACGCCGGTCTTGGAGGTGACTTCCGGCTGCCACATCGATTCCTGGTAGCCAATGGCTGCCAGCAGGCGCCAGTCGACCTGCTCGACCTTGGCGTAGCTCTTGAAGTGCTtctc                                                                                                >  2:1054724/1‑150 (MQ=255)
                                                          ggCGGTGCGCTGGGTCAGCATCATCAGGCCGCGCACGCCGGTCTTGGAGGTGACTTCCGGCTGCCACATCGATTCCTGGTAGCCAATGGCTGCCAGCAGGCGCCAGTCGACCTGCTCGACCTTGGCGTAGCTCTTGAAGTGCTTCTCGTa                                                                                             >  1:1017761/1‑150 (MQ=255)
                                                              gTGCGCTGGGTCAGCATCATCAGGCCGCGCACGCCGGTCTTGGAGGTGACTTCCGGCTGCCACATCGATTCCTGGTAGCCAATGGCTGc                                                                                                                                                      >  2:696756/1‑89 (MQ=255)
                                                              gTGCGCTGGGTCAGCATCATCAGGCCGCGCACGCCGGTCTTGGAGGTGACTTCCGGCTGCCACATCGATTCCTGGTAGCCAATGGCTGc                                                                                                                                                      <  1:696756/89‑1 (MQ=255)
                                                               tGCGCTGGGTCAGCATCATCAGGCCGCGCACGCCGGTCTTGGAGGTGACTTCCGGCTGCCACATCGATTCCTGGTAGCCAATGGCTGCCAGCAGGCGCCAGTCGACCTGCTCGACCTTg                                                                                                                       >  1:586022/1‑119 (MQ=255)
                                                               tGCGCTGGGTCAGCATCATCAGGCCGCGCACGCCGGTCTTGGAGGTGACTTCCGGCTGCCACATCGATTCCTGGTAGCCAATGGCTGCCAGCAGGCGCCAGTCGACCTGCTCGACCTTg                                                                                                                       <  2:586022/119‑1 (MQ=255)
                                                                       gTCAGCATCATAAGGCCGCGCACGCCGGTCTTGGAGGTGACTTCCGGCTGCCACATCGATTCCTGGTAGCCAATGGCTGCCAGCAGGCGCCAGTCGACCTGCTCGACCTTGGCGTAGCTCTTGAAGTGCTTCTCGTACTTGGGCAGgcgc                                                                                >  2:54428/1‑150 (MQ=255)
                                                                        tCAGCATCATCAGGCCGCGCACGCCGGTCTTGGAGGTGACTTCCGGCTGCCACATCGATTCCTGGCAGCCATTGGCTACcaggaat                                                                                                                                               >  1:1206850/1‑82 (MQ=38)
                                                                        tCAGCATCATCAGGCCGCGCACGCCGGTCTTGGAGGTGACTTCCGGCTGCCACATCGATTCCTGGCAGCCATTGGCTACcaggaat                                                                                                                                               <  2:1206850/86‑5 (MQ=38)
                                                                              tcatcaGGCCGCGCACGCCGGTCTTGGAGGTGACTTCCGGCTGCCACATCGATTCCTGGTAGCCCATGTGTGCCAGCAGGCGCCAGTCGACCTGCTCGACCTTGGCGTAGCTCTTGAAGTGCTTCTCGTACTTGGGCAGGCGCTGTTGCa                                                                         <  1:1157098/150‑1 (MQ=255)
                                                                                 tcaGGCCGCGCACGCCGGTCTTGGAGGTGACTTTCTGCGTCCACATCGATTCCTGGTAGCCAATGGCTGCCAGCAGGCGCCAGTCGACCTGCTCGACCTTGGCGTAGCTCTTGAAGTGCTTCTCGTACTTGGGCAGGCGCTGTTGCAGGt                                                                      <  2:114688/150‑1 (MQ=255)
                                                                                   aGGCCGCGCACGCCGGTCTTGGAGGTGACTTCCGGCTTCCACATCGATTCCTGGTAGCCAATGGCTGCCAGCAGGCGCCAGTCGACCTGCTCGACCTTGGCGTAGCTCTTGAAGTGCTTCTCGTACTTGGGCAGGCGCTGTTGCAGGTGc                                                                    >  1:1154145/1‑150 (MQ=255)
                                                                                   aGGCCGCGCACGCCGGTCTTGGAGGTGACTTCCGGCTGCCACATCGATTCCTGGTAGCCAATGGCTGCCAGCAGGCGCCAGTCGACCTGCTCGACCTTGGCGTAGCTCTTGAAGTGCTTCTCGTACTTGGGCAGGCGCTGTTGCAGGTGc                                                                    <  1:54428/150‑1 (MQ=255)
                                                                                    ggCCGCGCACGCCGGTCTTGGAGGTGACTTCCGGCTGCCACATCGATTCCTGGTAGCCAATGGCTGCCAGCAGGCGCCAGTCGACCTGCTCGACCTTGGCGTAGCTCTTGAAGTGCTTCTCGTACTTGGGCAGGCGCTGTTGCAGGTGCt                                                                   >  2:1125107/1‑150 (MQ=255)
                                                                                    ggCCGCGCACGCCGGTCTTGGAGGTGACTTCCGGCTGCCACATCGATTCCTGGTAGCCAATGGCTGCCAGCAGGCGCCAGTCGACCTGCTCGACCTTGGCGTAGCTCTTGAAGTGCTTCTCGTACTTGGGCAGGCGCTGTTGCAGGTGCt                                                                   >  2:845030/1‑150 (MQ=255)
                                                                                     gCCGCGCACGCCGGTCTTGGAGGTGACTTCCGGCTGCCACATCGATTCCTGGTAGCCAATGGCTGCCAGCAGGCGCCAGTCGACCTGCTCGACCTTGGCGTAGCTCTTGAAGTGCTTCTCGTACTTGGGCAGGCGCTGTTGCAGGTGCTg                                                                  >  1:659825/1‑150 (MQ=255)
                                                                                       cgcgCCCGCCGGTCTTGGAGGTGACTTCCGGCTGCCACATCGATTCCTGGTAGCCAATGGCTGCCAGCAGGCGCCAGTCGCCCTGCTCGACCTTGGCGTAGCTCTTGAAGTGCTTCTCGTACTTGGGCAGGCGCTGTTGCAGGTGCTggg                                                                <  2:6890/150‑1 (MQ=255)
                                                                                           cACGCCGGTCTTGGAGGTGACTTCCGGCTGCCACATCGATTCCTGGTAGCCAATGGCTGCCAGCAGGCGCCAGTCGACCTGCTCGACCTTGGCGTAGCTCTTGAAGTGCTTCTCGTACTTGGGCAGGCGCTGTTGCAGGTGCTGGGCGaa                                                            <  2:1017761/150‑1 (MQ=255)
                                                                                            aCGCCGGTCTTGGAGGTGACTTCCGGCTGCCACATCGATTCCTGGTAGCCAATGGCTGCCAGCAGGCGCCAGTCGACCTGCTCGACCTTGGCGTAGCTCTTGAAGTGCTTCTCGTACTTGGGCAGGCGCTGTTGCAGGTGCTGGGCGAAg                                                           <  1:845030/150‑1 (MQ=255)
                                                                                              gCCGGTCTTGGAGGTGACTTCCGGCTGCCACATCGATTCCTGGTAGCCAATGGCTGCCAGCAGGCGCCAGTCGACCTGCTCGACCTTg                                                                                                                       <  1:1196737/88‑1 (MQ=255)
                                                                                              gCCGGTCTTGGAGGTGACTTCCGGCTGCCACATCGATTCCTGGTAGCCAATGGCTGCCAGCAGGCGCCAGTCGACCTGCTCGACCTTg                                                                                                                       >  2:1196737/1‑88 (MQ=255)
                                                                                                cGGTCTTGGAGGTGACTTCCGGCTGCCACATCGATTCCTGGTAGCCAATGGCTGCCAGCAGGCGCCAGTCGACCTGCTCGACCTTGGCGTAGCTCTTGAAGTGCTTCTCGTACTTGGGCAGGCGCTGTTGCAGGTGCTGGGCGAAGgtgt                                                       >  2:1073411/1‑150 (MQ=255)
                                                                                                 ggTCTTGGAGGTGACTTCCGGCTGCCACATCGATTCCTGGTAGCCAATGGCTGCCAGCAGgcgc                                                                                                                                            <  1:684215/64‑1 (MQ=255)
                                                                                                 ggTCTTGGAGGTGACTTCCGGCTGCCACATCGATTCCTGGTAGCCAATGGCTGCCAGCAGgcgc                                                                                                                                            >  2:684215/1‑64 (MQ=255)
                                                                                                 ggTCTTGGAGGTGACTTCCGGCTGCCACATCGATTCCTGGTAGCCAATGGCTGCCAGCAGGCGCCAGTCGACCTGCTCGACCTTGGCGTAGCTCTTGAAGTGCTTCTCGTACTTGGGCAGGCGCTGTTGCAGGTGCTGGGCGAAGGTGTa                                                      <  2:408133/150‑1 (MQ=255)
                                                                                                  gTCTTGGAGGTGACTTCCGGCTGCCACATCGATTCCTGGTAGCCAATGGCTGCCAGCAGGCGCCAGTCGACCTGCTCGACCTTGGCGTAGCTCTTGAAGTGCTTCTCGTACTTGGGCAGGCGCTGTTGCAGGTGCTGGGCGAAGGTGTAg                                                     <  1:1340206/150‑1 (MQ=255)
                                                                                                     ttGGAGGTGACTTCCGGCTGCCACATCGATTCCTGGTAGCCAATGGCTGCCAGCAGGCGCCAGTCGACCTGCTCGACCTTGGCGTAGCTCTTGAAGTGCTTCTCGTACTTGGGCAGGCGCTGTTGCAGGTGCTGGGCGAAGGTGTAGgcg                                                  >  2:845282/1‑150 (MQ=255)
                                                                                                     ttGGAGGTGACTTCCGGCTGCCACATCGATTCCTGGTAGCCAATGGCTGCCAGCAGGCGCCAGTCGACCTGCTCGACCTTGGCGTAGCTCTTGAAGTGCTTCTCGTACTTGGGCAGGCGCTGTTGCAGGTGCTGGGCGAAGGTGTAGgcg                                                  >  2:1340416/1‑150 (MQ=255)
                                                                                                         aGGTGACTTCCGGCTGCCACATCGATTCCTGGTAGCCAATGGCTGCCAGCAGGCGCCAGTCGACCTGCTCGACCTTGGCGTAGCTCTTGAAGTGCTTCTCGTACTTGGGCAGGCGCTGTTGCAGGTGCTGGGCGAAGGTGTAGGCGCCGa                                              >  1:1346693/1‑150 (MQ=255)
                                                                                                           gTGACTTCCGGCTGCCACATCGATTCCTGGTAGCCAATGGCTGCCAGCAGGCGCCAGTCGACCTGCTCGCCCTTGGCGTAGCTCTTGAAGTGCTTCTCGTACTTGGGCAGGCGCTGTTGCAGGTGCTGGGCGAAGGTGTAGACGCCGACa                                            >  2:1002493/1‑150 (MQ=255)
                                                                                                              aCTTCCGGCTGCCACATCGATTCCTGGTAGCCAATGGCTGCCAGCAGGCGCCAGTCGACCTGCTCGACCTTGGCGTAGCTCTTGAAGTGCTTCTCGTACTTGGGCAGGCGCTGTTGCAGGTGCTGGGCGAAGGTGTAGGCGCCGACATAg                                         <  2:1322081/150‑1 (MQ=255)
                                                                                                                ttCCGGCTGCCACATCGATTCCTGGTAGCCAATGGCTGCCAGCAGGCGCCAGTCGACCTGCTCGACCTTGGCGTAGCTCTTGAAGTGCTTCTCGTACTTGGGCAGGCGCTGTTGCAGGTGCTGGGCGAAGGTGTAGGCGCCGACAt                                           <  1:59676/146‑1 (MQ=255)
                                                                                                                ttCCGGCTGCCACATCGATTCCTGGTAGCCAATGGCTGCCAGCAGGCGCCAGTCGACCTGCTCGACCTTGGCGTAGCTCTTGAAGTGCTTCTCGTACTTGGGCAGGCGCTGTTGCAGGTGCTGGGCGAAGGTGTAGGCGCCGACAt                                           >  2:59676/1‑146 (MQ=255)
                                                                                                                ttCCGGCTGCCACATCGATTCCTGGTAGCCAATGGCTGCCAGCAGGCGCCAGTCGACCTGCTCGACCTTGGCGTAGCTCTTGAAGTGCTTCTCGTACTTGGGCAGGCGCTGTTGCAGGTGCTGGGCGAAGGTGTAGGCGCCGACATAGcc                                       >  2:505076/1‑150 (MQ=255)
                                                                                                                   cGGCTGCCACATCGATTCCTGGTAGCCAATGGCTGCCAGCAGGCGCCAGTCGACCTGCTCGACCTTGGCGTAGCTCTTGAAGTGCTTCTCGTACTTGGGCAGGCGCTGTTGCAGGTGCTGGGCGAAGGTGTAGGCGCCGACATAGCCAAg                                    >  1:145419/1‑150 (MQ=255)
                                                                                                                         ccACATCGATTCCTGGTAGCCAATGGCTGCCAGCAGGCGCCAGTCGACCTGCTCGCCCTTGGCGTAGCTCTTGAAGTGCTTCTCGTACTTGGGCAGGCGCTGTTGCAGGTGCTGGGCTTCTGTGTTGTCGCCGACATAGCCTAGTACGtt                              >  1:653711/1‑149 (MQ=255)
                                                                                                                         ccACATCGATTCCTGGTAGCCAATGGCTGCCAGCAGGCGCCAGTCGACCTGCTCGACCTTGGCGTAGCTCTTGAAGTGCTTCTCGTACTTGGGCAGGCGCTGTTGCAGGTGCTGGGCGAAGGTGTAGGCGCCGACATAGCCAAGTACGTc                              >  2:654582/1‑150 (MQ=255)
                                                                                                                           acaTCGATTCCTGGTAGCCAATGGCTGCCAGCAGGCGCCAGTCGACCTGCTCGACCTTGGCGTAGCTCTTGAAGTGCTTCTCGTACTTGGGCAGGCGCTGTTGCAGGTGCTGGGCGAAGGTGTAGGCGCCGACATAGCCAAGTACGTCGa                            >  1:226384/1‑150 (MQ=255)
                                                                                                                             aTCGATTCCTGGTAGCCAATGGCTGCCAGCAGGCGCCAGTCGACCTGCTCGACCTTGGCGTAGCTCTTGAAGTGCTTCTCGTACTTGGGCAGGCGCTGTTGCAGGTGCTGGGCGAAGGTGTa                                                      >  2:1088197/1‑122 (MQ=255)
                                                                                                                             aTCGATTCCTGGTAGCCAATGGCTGCCAGCAGGCGCCAGTCGACCTGCTCGACCTTGGCGTAGCTCTTGAAGTGCTTCTCGTACTTGGGCAGGCGCTGTTGCAGGTGCTGGGCGAAGGTGTa                                                      <  1:1088197/122‑1 (MQ=255)
                                                                                                                              tCGATTCCTGGTAGCCAATGGCTGCCAGCAGGCGCCAGTCGACCTGCTCGACCTTGGCGTAGCTCTTGAAGTGCTTCTCGTACTTGGGCAGGCGCTGTTGCAGGTGCTGGGCGAAGGTGTAGGCGCCGACATAGCCAAGTACGTCGACAt                         >  1:763861/1‑150 (MQ=255)
                                                                                                                                  ttCCTGGTAGCCAATGGCTGCCAGCAGGCGCCAGTCGACCTGCTCGACCTTg                                                                                                                       <  1:1159617/52‑1 (MQ=255)
                                                                                                                                  ttCCTGGTAGCCAATGGCTGCCAGCAGGCGCCAGTCGACCTGCTCGACCTTg                                                                                                                       >  2:1159617/1‑52 (MQ=255)
                                                                                                                                     cTGGTAGCCAATGGCTGCCAGCAGGCGCCAGTCGACCTGCTCGACCTTGGCGTAGCTCTTGAAGTGCTTCTCGTACTTGGGCAGGCGCTGTTGCAGGTGCTGGGCGAAGGTGTAGGCGCCGACATAGCCAAGTACGTCGACATGGCCgta                  <  1:1340416/150‑1 (MQ=255)
                                                                                                                                      tGGTAGCCAATGGCTGCCAGCAGGCGCCAGTCGACCTGCTCGACCTTGGCGTAGCTCTTGAAGTGCTTCTCGTACTTGGGCAGGCGCTGTTGCAGGTGCTGGGCGAAGGTGTAGGCGCCGACATAGCCAAGTACGTCGACATGGCCgtag                 >  1:1221754/1‑150 (MQ=255)
                                                                                                                                        gTAGCCAATGGCTGCCAGCAGGCGCCAGTCGACCTGCTCGACCTTGGCGTAGCTCTTGAAGTGCTTCTCGTACTTGGGCAGGCGCTGTTGCAGGTGCTGGGCGAAGGTGTAGGCGCCGACATAGCCAAGTACGTCGACATGGCCgtagta               >  1:1001925/1‑150 (MQ=255)
                                                                                                                                            ccAATGGCTGCCAGCAGGCGCCAGTCGACCTGCTCGACCTTGGCGTAGCTCTTGAAGTGCTTCTCGTACTTGGGCAGGCGCTGTTGCAGGTGCTGGGCgca                                                            >  2:54129/1‑99 (MQ=255)
                                                                                                                                            ccAATGGCTGCCAGCAGGCGCCAGTCGACCTGCTCGACCTTGGCGTAGCTCTTGAAGTGCTTCTCGTACTTGGGCAGGCGCTGTTGCAGGTGCTGGGCGaa                                                            <  1:54129/101‑1 (MQ=255)
                                                                                                                                            ccAATGGCTGCCAGCAGGCGCCAGTCGACCTGCTCGACCTTGGCGTAGCTCTTGAAGTGCTTCTCGTACTTGGGCAGGCGCTGTTGCAGGTGCTGGGCGAAGGTGTAGGCGCCGACATAGCCAAGTACGTCGACATCGCAGTAGTAGAgg           >  2:80292/1‑150 (MQ=255)
                                                                                                                                            ccAATGGCTGCCAGCAGGCGCCAGTCGACCTGCTCGACATTGGCGTAGCTCTTGAAGTGCTTCTCGTACTTGGGCAgg                                                                                   >  1:661898/1‑78 (MQ=255)
                                                                                                                                            ccAATGGCTGCCAGCAGGCGCCAGTCGACCTGCTCGACATTGGCGTAGCTCTTGAAGTGCTTCTCGTACTTGGGCAgg                                                                                   <  2:661898/78‑1 (MQ=255)
                                                                                                                                             cAATGGCTGCCAGCAGGCGCCAGTCGACCTGCTCGACCTTGGCGTAGCTCTTGAAGTGCTTCTCGTACTTGGGCAGGCGCTGTTGCAGGTGCTGGGCGAAGGTGTAGGCGCCGACATAGCCAAGTACGTCGACATGGCCGTAGTAGCGGt          <  1:845282/150‑1 (MQ=255)
                                                                                                                                              aaTGGCTGCCAGCAGGCGCCAGTCGACCTGCTCGACCTTGGCGTAGCTCTTGAAGTGCTTCTCGTACTTGGGCAGGCGCTGTTGCAGGTGCTGGGCGa                                                             <  1:703306/98‑1 (MQ=255)
                                                                                                                                              aaTGGCTGCCAGCAGGCGCCAGTCGACCTGCTCGACCTTGGCGTAGCTCTTGAAGTGCTTCTCGTACTTGGGCAGGCGCTGTTGCAGGTGCTGGGCGAAGGTGTAGGCGCCGACATAGCCAAGTACGTCGACATGGCCGTCGTAGCGGTc         >  1:18592/1‑150 (MQ=255)
                                                                                                                                              aaTGGCTGCCAGAAGGAGCCAGTCGACCTGCTCGACCTTGGCGTAGCTCTTGAAGTGCTTCTCGTACTTGGGCAGGCGCTGTTGCAGGTGCTGGGCGa                                                             >  2:703306/1‑98 (MQ=255)
                                                                                                                                                     gCCAGCAGGCGCCAGTCGACCTGCTCGACCTTGGCGTAGCTCTTGAAGTGCTTCTCGTACTTGGGCAGGCGCTGTTGCAGGTGCTGGGCGAAGGTGTAGGCGCCGACATAGCCAAGTACGTCGACATGGCCGTAGTAGCGGTCCTTCAAt  <  2:1154145/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
TGGCGCCGCCCTGGATGCTCTGCCTCGGGTCCAGTCGGTTCGACACGCCCATGGCCTGGGCGGTGCGCTGGGTCAGCATCATCAGGCCGCGCACGCCGGTCTTGGAGGTGACTTCCGGCTGCCACATCGATTCCTGGTAGCCAATGGCTGCCAGCAGGCGCCAGTCGACCTGCTCGACCTTGGCGTAGCTCTTGAAGTGCTTCTCGTACTTGGGCAGGCGCTGTTGCAGGTGCTGGGCGAAGGTGTAGGCGCCGACATAGCCAAGTACGTCGACATGGCCGTAGTAGCGGTCCTTCAAT  >  NC_002947/1183319‑1183617

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