Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,183,456 T→C 9.4% Y312C (TAC→TGC)  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,4560TC 9.4% 209.8 / 8.4 74Y312C (TAC→TGC) mltFmembrane‑bound lytic murein transglycosylase F
Reads supporting (aligned to +/- strand):  ref base T (38/29);  new base C (3/4);  total (41/33)
Fisher's exact test for biased strand distribution p-value = 6.93e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.97e-01

ATGAAGTACTTGGCGCCGCCCTGGATGCTCTGCCTCGGGTCCAGTCGGTTCGACACGCCCATGGCCTGGGCGGTGCGCTGGGTCAGCATCATCAGGCCGCGCACGCCGGTCTTGGAGGTGACTTCCGGCTGCCACATCGATTCCTGGTAGCCAATGGCTGCCAGCAGGCGCCAGTCGACCTGCTCGACCTTGGCGTAGCTCTTGAAGTGCTTCTCGTACTTGGGCAGGCGCTGTTGCAGGTGCTGGGCGAAGGTGTAGGCGCCGACATAGCCAAGTACGTCGACATGGCCGTAGTA  >  NC_002947/1183309‑1183604
                                                                                                                                                   |                                                                                                                                                    
aTGAAGTACTTGGCGCCGCCCTGGATGCTCTGCCTCGGGTCCAGTCGGTTCGACACGCCCATGGCCTGGGCGGTGCGCTGGGTCAGCATCATCAGGCCGCGCACGCCGGTCTTGGAGGTGACTTCCGGCTGCCACATCGATTCCTGGTAg                                                                                                                                                    <  2:1128767/150‑1 (MQ=255)
 tGAAGTACTTGGCGCCGCCCTGGATGCTCTGCCTCGGGTCCAGTCGGTTCGACACGCCCATGGCCTGGGCGGTGCGCTGGGTCAGCATCATCAGGCCGCGCACGCCGGTCTTGGAGGTGACTTCCGGCTGCCACATCGATTCCTGGTAGc                                                                                                                                                   >  1:111140/1‑150 (MQ=255)
 tGAAGTACTTGGCGCCGCCCTGGATGCTCTGCCTCGGGTCCAGTCGGTTCGACACGCCCATGGCCTGGGCGGTGCGCTGGGTCAGCATCATCAGGCCGCGCACGCCGGTCTTGGAGGTGACTTCCGGCTGCCACATCGATTCCTGGCAGc                                                                                                                                                   >  2:1132774/1‑150 (MQ=255)
     gTACTTGGCGCCGCCCTGGATGCTCTGCCTCGGGTCCAGTCGGTTCGACACGCCCATGGCCTGGGCGGTGCGCTGGGTCAGCATCATCAGGCCGCGCACGCCGGTCTTGGAGGTGACTTCCGGCTGCCACATCGATTCCTGGTAGCCAAt                                                                                                                                               >  1:1395296/1‑150 (MQ=255)
      tACTTGGCGCCGCCCTGGATGCTCTGCCTCGGGTCCAGTCGGTTCGACACGCCCATGGCCTGGGCGGTGCGCTGGGTCAGCATCATCAGGCCGCGCACGCCGGTCTTGGAGGTGACTTCCGGCTGCCACATCGATTCCTGGTAGCCAATg                                                                                                                                              <  2:1240271/150‑1 (MQ=255)
          tGGCGCCGCCCTGGATGCTCTGCCTCGGGTCCAGTCGGTTCGACACGCCCATGGCCTGGGCGGTGCGCTGGGTCAGCATCATCAGGCCGCGCACGCCGGTCTTGGAGGTGACTTCCGGCTGCCACATCGATTCCTGGTCGCCAACGCCTg                                                                                                                                          <  2:19231/150‑1 (MQ=255)
             cgccgcCCTGGATGCTCTGCCTCGGGTCCAGTCGGTTCGACACGTCAATGGCATGGGGGCGGCGCTGACGCAGCATCATCAGGCCGCGCACGCCGGTCTTGGAGGTGACTTCCGGCTGCCACATCGATTCCTGGCAGCCATTGGCTACca                                                                                                                                       <  1:525972/150‑1 (MQ=255)
             cgccgcCCTGGATGCTCTGCCTCGGGTCCAGTCGGTTCGACACGCCCATGGCCTGGGCGGTGCGCTGGGTCAGCATCATCAGGCCGCGCACGCCGGTCTTGGAGGTGACTTCCGGCTGCCACATCGATTCCTGGCAGCCATTGGCTACca                                                                                                                                       <  1:1132774/150‑1 (MQ=255)
                            tctGCCTCGGGTCCAGTCGGTTCGACACGCCCATGGCCTGGGCGGTGCGCTGGGTCAGCATCATCAGGCCGCGCACGCCGGTCTTGGAGGTGACTTCCGGCTGCCACATCGATTCCTGGCAGCCATTGGCTACca                                                                                                                                       >  1:632493/1‑135 (MQ=255)
                            tctGCCTCGGGTCCAGTCGGTTCGACACGCCCATGGCCTGGGCGGTGCGCTGGGTCAGCATCATCAGGCCGCGCACGCCGGTCTTGGAGGTGACTTCCGGCTGCCACATCGATTCCTGGCAGCCATTGGCTACca                                                                                                                                       <  2:632493/135‑1 (MQ=255)
                            tctGCCTCGGGTCCAGTCGGTTCGACACGCCCATGGCCTGGGCGGTGCGCTGGGTCAGCATCATCAGGCCGCGCACGCCGGTCTTGGAGGTGACTTCCGGCTGCAACATCGATTCCTGGTAGCCAATGGCTGCCAGCAGGCGCCAGTCGa                                                                                                                        >  1:1373266/1‑150 (MQ=255)
                             ctGCCTCGGGTCCAGTCGGTTCGACACGCCCATGGCCTGGGCGGTGCGCTGGGTCAGCATCATCAGGCCGCGCACGCCGGTCTTGGAGGTGACTCCCGCATTTTACATCGATTCCTGGTAGCCAATGGCTGCCAGCAGGCGCCAGTCGAc                                                                                                                       <  1:1327532/150‑1 (MQ=255)
                                ccTCGGGTCCAGTCGGTTCGACACGCCCATGGCCTGGGCGGTGCGCTGGGTCAGCATCATCAGGCCGCGCACGCCGGTCTTGGAGGTGACTTCCGGCTGCCACATCGATTCCTGGTAg                                                                                                                                                    <  1:117089/118‑1 (MQ=255)
                                ccTCGGGTCCAGTCGGTTCGACACGCCCATGGCCTGGGCGGTGCGCTGGGTCAGCATCATCAGGCCGCGCACGCCGGTCTTGGAGGTGACTTCCGGCTGCCACATCGATTCCTGGTAg                                                                                                                                                    >  2:117089/1‑118 (MQ=255)
                                ccTCGGGTCCAGTCGGTTCGACACGCCCATGGCCTGGGCGGTGCGCTGGGTCAGCATCATCAGGCCGCGCACGCCGGTCTTGGAGGTGACTTCCGGCTGCCACATCGATTCCTGGTAGCCAATg                                                                                                                                              <  2:1186398/124‑1 (MQ=255)
                                ccTCGGGTCCAGTCGGTTCGACACGCCCATGGCCTGGGCGGTGCGCTGGGTCAGCATCATCAGGCCGCGCACGCCGGTCTTGGAGGTGACTTCCGGCTGCCACATCGATTCCTGGTAGCCAATg                                                                                                                                              >  1:1186398/1‑124 (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                                                                                                                                         <  1:696756/89‑1 (MQ=255)
                                                                        gTGCGCTGGGTCAGCATCATCAGGCCGCGCACGCCGGTCTTGGAGGTGACTTCCGGCTGCCACATCGATTCCTGGTAGCCAATGGCTGc                                                                                                                                         >  2:696756/1‑89 (MQ=255)
                                                                         tGCGCTGGGTCAGCATCATCAGGCCGCGCACGCCGGTCTTGGAGGTGACTTCCGGCTGCCACATCGATTCCTGGTAGCCAATGGCTGCCAGCAGGCGCCAGTCGACCTGCTCGACCTTg                                                                                                          >  1:586022/1‑119 (MQ=255)
                                                                         tGCGCTGGGTCAGCATCATCAGGCCGCGCACGCCGGTCTTGGAGGTGACTTCCGGCTGCCACATCGATTCCTGGTAGCCAATGGCTGCCAGCAGGCGCCAGTCGACCTGCTCGACCTTg                                                                                                          <  2:586022/119‑1 (MQ=255)
                                                                                 gTCAGCATCATAAGGCCGCGCACGCCGGTCTTGGAGGTGACTTCCGGCTGCCACATCGATTCCTGGTAGCCAATGGCTGCCAGCAGGCGCCAGTCGACCTGCTCGACCTTGGCGTAGCTCTTGAAGTGCTTCTCGTACTTGGGCAGgcgc                                                                   >  2:54428/1‑150 (MQ=255)
                                                                                  tCAGCATCATCAGGCCGCGCACGCCGGTCTTGGAGGTGACTTCCGGCTGCCACATCGATTCCTGGCAGCCATTGGCTACcaggaat                                                                                                                                  <  2:1206850/86‑5 (MQ=38)
                                                                                  tCAGCATCATCAGGCCGCGCACGCCGGTCTTGGAGGTGACTTCCGGCTGCCACATCGATTCCTGGCAGCCATTGGCTACcaggaat                                                                                                                                  >  1:1206850/1‑82 (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)
                                                                                               gCCGCGCACGCCGGTCTTGGAGGTGACTTCCGGCTGCCACATCGATTCCTGGTAGCCAATGGCt                                                                                                                                           >  1:265422/1‑64 (MQ=255)
                                                                                               gCCGCGCACGCCGGTCTTGGAGGTGACTTCCGGCTGCCACATCGATTCCTGGTAGCCAATGGCt                                                                                                                                           <  2:265422/64‑1 (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                                                                                                          >  2:1196737/1‑88 (MQ=255)
                                                                                                        gCCGGTCTTGGAGGTGACTTCCGGCTGCCACATCGATTCCTGGTAGCCAATGGCTGCCAGCAGGCGCCAGTCGACCTGCTCGACCTTg                                                                                                          <  1:1196737/88‑1 (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                              >  2:59676/1‑146 (MQ=255)
                                                                                                                          ttCCGGCTGCCACATCGATTCCTGGTAGCCAATGGCTGCCAGCAGGCGCCAGTCGACCTGCTCGACCTTGGCGTAGCTCTTGAAGTGCTTCTCGTACTTGGGCAGGCGCTGTTGCAGGTGCTGGGCGAAGGTGTAGGCGCCGACAt                              <  1:59676/146‑1 (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)
                                                                                                                                                   |                                                                                                                                                    
ATGAAGTACTTGGCGCCGCCCTGGATGCTCTGCCTCGGGTCCAGTCGGTTCGACACGCCCATGGCCTGGGCGGTGCGCTGGGTCAGCATCATCAGGCCGCGCACGCCGGTCTTGGAGGTGACTTCCGGCTGCCACATCGATTCCTGGTAGCCAATGGCTGCCAGCAGGCGCCAGTCGACCTGCTCGACCTTGGCGTAGCTCTTGAAGTGCTTCTCGTACTTGGGCAGGCGCTGTTGCAGGTGCTGGGCGAAGGTGTAGGCGCCGACATAGCCAAGTACGTCGACATGGCCGTAGTA  >  NC_002947/1183309‑1183604

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: ‑