Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,183,462 A→T 8.0% I310N (ATT→AAT)  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,4620AT 8.0% 225.8 / 9.5 75I310N (ATT→AAT) mltFmembrane‑bound lytic murein transglycosylase F
Reads supporting (aligned to +/- strand):  ref base A (39/30);  new base T (2/4);  total (41/34)
Fisher's exact test for biased strand distribution p-value = 4.01e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.74e-01

GTACTTGGCGCCGCCCTGGATGCTCTGCCTCGGGTCCAGTCGGTTCGACACGCCCATGGCCTGGGCGGTGCGCTGGGTCAGCATCATCAGGCCGCGCACGCCGGTCTTGGAGGTGACTTCCGGCTGCCACATCGATTCCTGGTAGCCAATGGCTGCCAGCAGGCGCCAGTCGACCTGCTCGACCTTGGCGTAGCTCTTGAAGTGCTTCTCGTACTTGGGCAGGCGCTGTTGCAGGTGCTGGGCGAAGGTGTAGGCGCCGACATAGCCAAGTACGTCGACATGGCCGTAGTAGCGGTC  >  NC_002947/1183314‑1183610
                                                                                                                                                    |                                                                                                                                                    
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                                                                                                                                             <  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)
                           ccTCGGGTCCAGTCGGTTCGACACGCCCATGGCCTGGGCGGTGCGCTGGGTCAGCATCATCAGGCCGCGCACGCCGGTCTTGGAGGTGACTTCCGGCTGCCACATCGATTCCTGGTAGCCAATg                                                                                                                                                    >  1:1186398/1‑124 (MQ=255)
                           ccTCGGGTCCAGTCGGTTCGACACGCCCATGGCCTGGGCGGTGCGCTGGGTCAGCATCATCAGGCCGCGCACGCCGGTCTTGGAGGTGACTTCCGGCTGCCACATCGATTCCTGGTAGCCAATg                                                                                                                                                    <  2:1186398/124‑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:845030/1‑150 (MQ=255)
                                                                                         ggCCGCGCACGCCGGTCTTGGAGGTGACTTCCGGCTGCCACATCGATTCCTGGTAGCCAATGGCTGCCAGCAGGCGCCAGTCGACCTGCTCGACCTTGGCGTAGCTCTTGAAGTGCTTCTCGTACTTGGGCAGGCGCTGTTGCAGGTGCt                                                            >  2:1125107/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                                                                                                                <  1:1196737/88‑1 (MQ=255)
                                                                                                   gCCGGTCTTGGAGGTGACTTCCGGCTGCCACATCGATTCCTGGTAGCCAATGGCTGCCAGCAGGCGCCAGTCGACCTGCTCGACCTTg                                                                                                                >  2:1196737/1‑88 (MQ=255)
                                                                                                     cGGTCTTGGAGGTGACTTCCGGCTGCCACATCGATTCCTGGTAGCCAATGGCTGCCAGCAGGCGCCAGTCGACCTGCTCGACCTTGGCGTAGCTCTTGAAGTGCTTCTCGTACTTGGGCAGGCGCTGTTGCAGGTGCTGGGCGAAGgtgt                                                >  2:1073411/1‑150 (MQ=255)
                                                                                                      ggTCTTGGAGGTGACTTCCGGCTGCCACATCGATTCCTGGTAGCCAATGGCTGCCAGCAGgcgc                                                                                                                                     >  2:684215/1‑64 (MQ=255)
                                                                                                      ggTCTTGGAGGTGACTTCCGGCTGCCACATCGATTCCTGGTAGCCAATGGCTGCCAGCAGgcgc                                                                                                                                     <  1:684215/64‑1 (MQ=255)
                                                                                                      ggTCTTGGAGGTGACTTCCGGCTGCCACATCGATTCCTGGTAGCCAATGGCTGCCAGCAGGCGCCAGTCGACCTGCTCGACCTTGGCGTAGCTCTTGAAGTGCTTCTCGTACTTGGGCAGGCGCTGTTGCAGGTGCTGGGCGAAGGTGTa                                               <  2:408133/150‑1 (MQ=255)
                                                                                                       gTCTTGGAGGTGACTTCCGGCTGCCACATCGATTCCTGGTAGCCAATGGCTGCCAGCAGGCGCCAGTCGACCTGCTCGACCTTGGCGTAGCTCTTGAAGTGCTTCTCGTACTTGGGCAGGCGCTGTTGCAGGTGCTGGGCGAAGGTGTAg                                              <  1:1340206/150‑1 (MQ=255)
                                                                                                          ttGGAGGTGACTTCCGGCTGCCACATCGATTCCTGGTAGCCAATGGCTGCCAGCAGGCGCCAGTCGACCTGCTCGACCTTGGCGTAGCTCTTGAAGTGCTTCTCGTACTTGGGCAGGCGCTGTTGCAGGTGCTGGGCGAAGGTGTAGgcg                                           >  2:1340416/1‑150 (MQ=255)
                                                                                                          ttGGAGGTGACTTCCGGCTGCCACATCGATTCCTGGTAGCCAATGGCTGCCAGCAGGCGCCAGTCGACCTGCTCGACCTTGGCGTAGCTCTTGAAGTGCTTCTCGTACTTGGGCAGGCGCTGTTGCAGGTGCTGGGCGAAGGTGTAGgcg                                           >  2:845282/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                                                                                                                >  2:1159617/1‑52 (MQ=255)
                                                                                                                                       ttCCTGGTAGCCAATGGCTGCCAGCAGGCGCCAGTCGACCTGCTCGACCTTg                                                                                                                <  1:1159617/52‑1 (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                                                                            <  2:661898/78‑1 (MQ=255)
                                                                                                                                                 ccAATGGCTGCCAGCAGGCGCCAGTCGACCTGCTCGACATTGGCGTAGCTCTTGAAGTGCTTCTCGTACTTGGGCAgg                                                                            >  1:661898/1‑78 (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)
                                                                                                                                                    |                                                                                                                                                    
GTACTTGGCGCCGCCCTGGATGCTCTGCCTCGGGTCCAGTCGGTTCGACACGCCCATGGCCTGGGCGGTGCGCTGGGTCAGCATCATCAGGCCGCGCACGCCGGTCTTGGAGGTGACTTCCGGCTGCCACATCGATTCCTGGTAGCCAATGGCTGCCAGCAGGCGCCAGTCGACCTGCTCGACCTTGGCGTAGCTCTTGAAGTGCTTCTCGTACTTGGGCAGGCGCTGTTGCAGGTGCTGGGCGAAGGTGTAGGCGCCGACATAGCCAAGTACGTCGACATGGCCGTAGTAGCGGTC  >  NC_002947/1183314‑1183610

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: