Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,183,456 A→T 7.4% 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,4560AT 7.4% 193.4 / 5.5 68I310N (ATT→AAT) mltFmembrane‑bound lytic murein transglycosylase F
Reads supporting (aligned to +/- strand):  ref base A (29/34);  new base T (3/2);  total (32/36)
Fisher's exact test for biased strand distribution p-value = 6.60e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.15e-01

GTACTTGGCGCCGCCCTGGATGCTCTGCCTCGGGTCCAGTCGGTTCGACACGCCCATGGCCTGGGCGGTGCGCTGGGTCAGCATCATCAGGCCGCGCACGCCGGTCTTGGAGGTGACTTCCGGCTGCCACATCGATTCCTGGTAGCCAATGGCTGCCAGCAGGCGCCAGTCGACCTGCTCGACCTTGGCGTAGCTCTTGAAGTGCTTCTCGTACTTGGGCAGGCGCTGTTGCAGGTGCTGGGCGAAGGTGTAGGCGCCGACATAGCCAAGTACGTCGACATGGCCGTAGTAGCGG  >  NC_002947/1183308‑1183602
                                                                                                                                                    |                                                                                                                                                  
gTACTTGGCGCCGCCCTGGATGCTCTGCCTCGGGTCCAGTCGGTTCGACACGCCCATGGCCTGGGCGGTGCGCTGGGTCAGCATCATCAGGCCGCGCACGCCGGTCTTGGAGGTGACTTCCGGCTGCCACATCGATTCCTGGTAGCCaa                                                                                                                                                    <  6:248156/149‑1 (MQ=255)
      ggCGCCGCCCTGGATGCTCTGCCTCGGGTCCAGTCGGTTCGACACGCCCATGGCCTGGGCGGTGCGCTGGGTCAGCATCATCAGGCCGCGCACGCCGGTCTTGGAGGTGACTTCCGGCTGCCACATCGATTCCTGGTAGCCAATGGc                                                                                                                                                <  6:36926/147‑1 (MQ=255)
               cTGGATGCTCTGCCTCGGGTCCAGTCGGTTCGACACGCCCATGGCCTGGGCGGTGCGCTGGGTCAGCATCATCAGGCCGCGCACGCCGGTCTTGGAGGTGACTTCCGGCTGCCACATCGATTCCTGGTAGCCAATGGCTGCCAGCAGgc                                                                                                                                     >  5:195435/1‑149 (MQ=255)
                 ggATGCTCTGCCTCGGGTCCAGTCGGTTCGACACGCCCATGGCCTGGGCGGTGCGCTGGGTCAGCATCATCAGGCCGCGCACGCCGGTCTTGGAGGTGACTTCCGGCTGCCACATCGATTCCTGGTAGCCAATGGCTGCCAGCAGgcgc                                                                                                                                   <  6:155724/149‑1 (MQ=255)
                 ggATGCTCTGCCTCGGGTCCAGTCGGTTCGACAAGCCCATGGCCTGGGAGGTGCGCTGGGTCAGCATCATCAGGCCGCGCACGCCGGTCTTGGAGGTGTCTTCCGGCTGCCACATCGATTCCTGGTAGCCAATGGCTGCCAGCAGgcgc                                                                                                                                   >  2:207362/1‑149 (MQ=255)
                  gATGCTGTGCCTCGGGTCCAGTCGGTTCGACACGCCCATGGCCTGGGCGGTGCGCTGGGTCAGCATCATCAGGCCGCGCACGCCGGTCTTGGAGGTGACTTCCGGCTGCCACATCGATTCCTGGTAGCCAATGGCTGCCAGCAGgcgc                                                                                                                                   <  1:147214/148‑1 (MQ=255)
                  gATGCTCTGCCTCGGGTCCAGTCGGTTCGACACGCCCATGGCCTGGGCGGTGCGCTGGGTCAGCATCATCAGGCCGCGCACGCCGGTCTTGGAGGTGACTTCCGGCTGCCACATCGATTCCTGGTAGCCAATGGCTGCCAGCAGgcgc                                                                                                                                   <  6:5713/148‑1 (MQ=255)
                  gATGCTCTGCCTCGGGTCCAGTCGGTTCGACACGCCCATGGCCTGGGCGGTGCGCTGGGTCAGCATCATCAGGCCGCGCACGCCGGTCTTGGAGGTGACTTCCGGCTGCCACATCGATTCCTGGTAGCCAATGGCTGCCAGCAGgcgc                                                                                                                                   <  1:28159/148‑1 (MQ=255)
                   aTGCTCTGCCTCGGGTCCAGTCGGTTCGACACGCCCATGGCCTGGGCGGTGCGCTGGGTCAGCATCATCAGGCCGCGCACGCCGGTCTTGGAGGTGACTTCCGGCTGCCACATCGATTCCTGGTAGCCAATGGCTGCCAGCAGGCGCCa                                                                                                                                 <  1:220990/149‑1 (MQ=255)
                      ctctGCCTCGGGTCCAGTCGGTTCGACACGCCCATGGCCTGGGCGGTGCGCTGGGTCAGCATCATCAGGCCGCGCACGCCGGTCTTGGAGGTGACTTCCGGCTGCCACATCGATTCCTGGTAGCCAATGGCTGCCAGCAGGCGCCAGt                                                                                                                               <  3:27383/148‑1 (MQ=255)
                           ccTCGGGTCCAGTCGGTTCGACACGCCCATGGCCTGGGCGGTGCGCTGGGTCAGCATCATCAGGCCGCGCACGCCGGTCTTGGAGGTGACTTCCGGCTGCCACATCGATTCCTGGTAGCCAATGGCTGCCAGCAGGCGCCAGTCGACCt                                                                                                                         >  8:371224/1‑149 (MQ=255)
                            gtCGGGGCCTGTCGGTTCGACACGCCCATGGCGTGGGCGGTGGGCTGGGTCAGCATCATCAGGCCGCGCACGCCGGTCTTGGAGGTGACTTCCGGCTGCCACATCGATTCCTGGTAGCCAATGGCTGCCAGCAGGCGCCAGTCGagcag                                                                                                                        <  6:263760/148‑5 (MQ=255)
                                    cAGTCGGTTCGACACGCCCATGGCCTGGGCGGTGCGCTGGGTCAGCATCATCAGGCCGCGCACGCCGGTCTTGGAGGTGACTTCCGGCTGCCACATCGATTCCTGGTAGCCAATGGCTGCCAGCAGGCGCCAGTCGACCTGCTCGACCt                                                                                                                >  8:202525/1‑149 (MQ=255)
                                    cAGTCGGTTCGACACGCCCATGGCCTGGGCGGTGCGCTGGGTCAGCATCATCAGGCCGCGCACGCCGGTCTTGGAGGTGACTTCCGGCTGCCACATCGATTCCTGGTAGCCAATGGCTGCCAGCAGGCGCCAGTCGACCTGCTCGACCt                                                                                                                >  8:208069/1‑149 (MQ=255)
                                    cAGTCGGTTCGACACGCCCATGGCCTGGGCGCTGCGCTGGGTCAGCATCATCAGGCCGCCCACGCCGGTCTTGCAGGCCACTTCCGGCTGCCACATCGATTCCTGGTAGCCAATGACTGCCAGCAGGCGCCATTCGACCTGCTCGACCt                                                                                                                >  8:254211/1‑149 (MQ=255)
                                     aGTCGGTTCGACACGCCCATGGCCTGGGCGGTGCGCTGGGTCAGCATCATCAGGCCGCGCACGCCGGTCTTGGAGGTGACTTCCGGCTGCCACATCGATTCCTGGTAGCCAATGGCTGCCAGCAGGCGCCAGTCGACCTGCTCGACCtt                                                                                                               >  8:192377/1‑149 (MQ=255)
                                      ggcGGTTCGACACGCCCATGGCCTGGGCGGTGCGCTGGGTCAGCATCATCAGGCCGCGCACGCCGGTCTTGGAGGTGACTTCCGGCTGCCACATCGATTCCTGGTAGCCAATGGCTGCCAGCAGGCGCCAGTCGACCTGCTCGACCTTg                                                                                                              <  8:126453/147‑1 (MQ=255)
                                      gTCGGTTCGACACGCCCATGGCCTGGGCGGTGCGCTGGGTCAGCATCATCAGGCCGCGCACGCCGGTCTTGGAGGTGACTTCCGGCTGCCACATCGATTCCTGGTAGCCAATGGCTGCCAGCAGGCGCCAGTCGACCTGCTCGAGCTTg                                                                                                              <  2:147196/149‑1 (MQ=255)
                                      gTCGGTTCGACACGCCCATGGCCTGGGCGGTGCGCTGGGTCAGCATCATCAGGCCGCGCACGCCGGTCTTGGAGGTGACTTCCGGCTGCCACATCGATTCCTGGTAGCCAATGGCTGCCAGCAGGCGCCAGTCGACCTGCTCGACCTTg                                                                                                              <  7:202525/149‑1 (MQ=255)
                                      gTCGGTTCGACACGCCCATGGCCTGGGCGGTGCGCTGGGTCAGCATCATCAGGCCGCGCACGCCGGTCTTGGAGGTGACTTCCGGCTGCCACATCGATTCCTGGTAGCCAATGGCTGCCAGCAGGCGCCAGTCGACCTGCTCGACCTTg                                                                                                              <  7:208069/149‑1 (MQ=255)
                                       tCGGTTCGACACGCCCATGGCCTGGGCGGTGCGCTGGGTCAGCATCATCAGGCCGCGCACGCCGGTCTTGGAGGTGACTTCCGGCTGCCACATCGATTCCTGGCAGCCATTGGCTACca                                                                                                                                           <  7:258792/119‑1 (MQ=255)
                                       tCGGTTCGACACGCCCATGGCCTGGGCGGTGCGCTGGGTCAGCATCATCAGGCCGCGCACGCCGGTCTTGGAGGTGACTTCCGGCTGCCACATCGATTCCTGGCAGCCATTGGCTACca                                                                                                                                           >  8:258792/1‑119 (MQ=255)
                                         ggTTCGACACGCCCATGGCCTGGGCGGTGCGCTGGGTCAGCATCATCAGGCCGCGCACGCCGGTCTTGGAGGTGACTTCCGGCTGCCACATCGATTCCTGGTAGCCAATGGCTGCCAGCAGGCGCCAGTCGACCTGCTCGACCTTGGc                                                                                                            <  7:192377/148‑1 (MQ=255)
                                           ttCGACACGCCCATGGCCTGGGCGGTGCGCTGGTTCAGCATCATCAGGCCGCGCACGCCGGTCTTGGAGGTGACTTCCGGCTGCCACATCGATTCCTGGTAGCCAATGGCTGCCAGCAGGCGCGAGTCGACCTGCTCGACCTTGGCGt                                                                                                          <  5:165285/148‑1 (MQ=255)
                                              gACACGCCCATGGCCTGGGCGGTGCGCTGGGTCAGCATCATCAGGCCGCGCACGCCGGTCTTGGAGGTGACTTCCGGCTGCCACATCGATTCCTGGTAGCCAATGGCTGCCAGCAGGCGCCAGTCGAc                                                                                                                           >  7:232555/1‑128 (MQ=255)
                                              gACACGCCCATGGCCTGGGCGGTGCGCTGGGTCAGCATCATCAGGCCGCGCACGCCGGTCTTGGAGGTGACTTCCGGCTGCCACATCGATTCCTGGTAGCCAATGGCTGCCAGCAGGCGCCAGTCGAc                                                                                                                           <  8:232555/128‑1 (MQ=255)
                                              gACACGCCCATGGCCTGGGCGGTGCGCTGGGTCAGCATCATCAGGCCGCGCACGCCGGTCTTGGAGGTGACTTCCGGCTGCCACATCGATTCCTGGTAGCCAATGGCTGCCAGCAGGCGCCAGTCGACCTGCTCGACCTTGGCGTAGct                                                                                                      >  5:313946/1‑149 (MQ=255)
                                                 acGCCCATGGCCTGGGCGGTGCGCTGGGTCAGCATCATCAGGCCGCGCACGCCGGTCTTGGAGGTGACTTCCGGCTGCCACATCGATTCCTGGTAGCCAATGGCTGCCAGCAGGCGCCAGTCGACCTGCTCGACCTTGGCGTAGCTCtt                                                                                                   >  4:150902/1‑149 (MQ=255)
                                                   gCCCATGGCCTGGGCGGTGCGCTGGGTCAGCATCATCAGGCCGCGCACGCCGGTCTTGGAGGTGACTTCCGGCTGCCACATCGATTCCTGGTAGCCAATGGCTGCCAGCAGGCGCCAGTCGACCTGCTCGACCTTGGCGTAGCTCTTGa                                                                                                 <  7:371224/149‑1 (MQ=255)
                                                    cccATGGACTGGGCGGTGCGCTGGGTCAGCATCATCAGGCCGCGCACGCCGGTCTTGGAGGTGACTTCCGGCTGCCACATCGATTCCTGGTAGCCAATGGCTGCCAGCAGGCGCCAGTCGACCTGCTCGACCTTGGCGTAGCTCTTGaa                                                                                                <  6:55229/149‑1 (MQ=255)
                                                      cATGGCCTGGGCGGTGCGCTGGGTCAGCATCATCAGGCCGCGCACGCCGGTCTTGGAGGTGACTTCCGGCTGCCACATCGATTCCTGGTAGCCAATGGCTGCCAGCAGGCGCCAGTCGACATGATCGACCATGGCGTAGATATTGaatt                                                                                              >  6:269423/1‑147 (MQ=255)
                                                        tGGCCTGGGCGGTGCGCTGGGTCAGCATCATCAGGCCGCGCACGCCGGTCTTGGAGGTGACTTCCGGCTGCCACATCGATTCCTGGTAGCCAATGGCTGCCAGCAGGCGCCAGTCGACCTGCTCGACCTTGGCGTAGCTCTTGAAGTGc                                                                                            <  6:195435/149‑1 (MQ=255)
                                                                   gTGCGCTGGGTCAGCATCATCAGGCCGCGCACGCCGGTCTTGGAGGTGACTTCCGGCTGCCACATCGATTCCTGGTAGCCAATGGCTGCCAGCAGGCGCCAGTCGACCTGCTCGACCTTGGCGTAGCTCTTGAAGTGCTTCTCGTACtt                                                                                 <  7:340617/149‑1 (MQ=255)
                                                                           ggTCAGCATCATCAGGCCGCGCACGCCGGTCTTGGAGGTGACTTCCGGCTGCCACATCGATTCCTGGTAGCCAATGGCTGCCAGCAGGCGCCAGTCGACCTGCTCGACCTTGGCGTAGCTCTTGAAGTGCTTCTCGTACTTGGGCAGgc                                                                         >  3:96402/1‑149 (MQ=255)
                                                                                 catcatCAGGCCGCGCACGCCGGTCTTGGAGGTGACTTCCGGCTGCCACATCGATTCCTGGTAGCCAATGGCTGCCAGCAGGCGCCAGTCGACCTGCACGACCTTGGCGTAGCTCTTGAAGTGCTTCTCGTACTTGGGCAGGCGCTGtt                                                                   >  1:128417/1‑149 (MQ=255)
                                                                                   tcatcaGGCCGCGCACGCCGGTCTTGGAGGTGACTTCCGGCTGCCACATCGATTCCTGGTAGCCAATGGCTGCCAGCAGGCGCCAGTCGACCTGCTCGACCTTGGCGTAGCTCTTGAAGTGCTTCTCGTACTTGGGCAGGCGCTGTTGc                                                                 <  6:313946/149‑1 (MQ=255)
                                                                                     atcaGGCCGCGCACGCCGGTCTTGGAGGTGACTTCCGGCTGCCACATGGATTCCTGGTAGCCAATGGCTGCCAGCAGGCGCCAGTCGACCTGCTCGACCTTGGCGTAGCTCTTGAAGTGCTTCTCGTACTTGGGCAGGCGCTGTTGCAg                                                               >  5:301518/1‑149 (MQ=255)
                                                                                     atcaGGCCGCGCACGCCGGTCTTGGAGGTGACTTCCGGCTGCCACATCGATTCCTGGTAGCCAATGGCTGCCAGCAGGCGCCAGTCGACCTGCTCGACCTTGGCGTAGCTCTTGAAGTGCTTCTCGTACTTGGGCAGGCGCTGTTGCAg                                                               >  2:242575/1‑149 (MQ=255)
                                                                                                  cGCCGGTCTTGGAGGTGACTTCCGGCTGCCACATCGATTCCTGGTAGCCAATGGCTGCCAGCAGGCGCCAGTCGACCTGCTCGACCTTGGCGTAGCTCTTGAAGTGCTTCTCGTACTTGGGCAGGCGCTGTTGCAGGTGCTGGGCGAAg                                                  <  8:82931/149‑1 (MQ=255)
                                                                                                        tCTTGGAGGTGACTTCCGGCTGCCACATCGATTCCTGGTAGCCAATGGCTGCCAGCAGGCGCCAGTCGACCTGCTCGACCTTGGCGTAGCTCTTGAAGTGCTTCTCGTACTTGGGCAGGCGCTGTTGCAGGTGCTGGGCGAAGGTGTAg                                            <  4:345546/149‑1 (MQ=255)
                                                                                                                   aCTTCCGGCTGCCACATCTATTCCTGGTAGCCAATGGCTGCCAGCAGGCGCCAGTCGACCTGCTCTACCTTGGCGTAGCTCTTGAAGTGCTTCTCGTACTTGGGCAGGCGCTGTTGCAGGTGCTGGGCGAAGGTGTAGGCGCCGACATa                                 <  7:254211/149‑1 (MQ=255)
                                                                                                                     ttCCGGCTGCCACATCGATTCCTGGTAGCCAATGGCTGCCAGCAGGCGCCAGTCGACCTGCTCGACCTTGGCGTAGCTCTTGAAGTGCTTCTCGTACTTGGGCAGGCGCTGTTGCAGGTGCTGGGCGAAGGTGTAGGCGCCGACATAGc                               >  1:128752/1‑149 (MQ=255)
                                                                                                                        cGGCTGCCACATCGATTCCTGGTAGCCAATGGCTGCCAGCAGGCGCCAGTCGACCTGCTCGACCTTGGCGTAGCTCTTGAAGTGCTTCTCGTACTTGGGCAGGCGCTGTTGCAGGTGCTGGGCGAAGGTGTAGGCGCCGACATAGCCaa                            <  6:301518/149‑1 (MQ=255)
                                                                                                                        cGGCTGCCACATCGATTCCTGGTAGACAATGGCTGCCAGCAGGCGCCAGTCGCCCTGCTCGACCTTGGCGTAGCTCTTGAAGTGCTTCTCGTACTTGGGCAGGCGCTGTTGCAGGTGCTGGGCGATGGTGTAGGCGGCGACATAGcc                              >  4:142748/1‑147 (MQ=255)
                                                                                                                              gcacaTCGATTACTGGTAGAAAATGGCTGCCAGCAGGCGCCAGTCGAACTGCTCGACCTTGGCGTAGCTCTTAAAGTGCTTCTCGTACTTTGGCAGGCGCTGTTGCTGGTGCTGGGCGAAGGTGTAGGCGCCGACATAGCCAAGTACga                      >  4:210/2‑148 (MQ=255)
                                                                                                                               cacaTCGATTCCTGGTAGCCAATGGCTGCCAGCAGGCGCCAGTCGACCTGCTCGACCTTGGCGTAGCTCTTGAAGTGCTTCTCGTACTTGGGCAGGCGCTGTTGCAGGTGCTGGGCGAAGGTGTAGGCGCCGACATAGCCAAGCACGTc                     >  6:255277/1‑149 (MQ=255)
                                                                                                                                 caTCGATTCCTGGTAGCCAATGGCTGCCAGCAGGCGCCAGTCGACCTGCTCGACCTTGGCGTAGCTCTTGAAGTGCTTCTCGTACTTGGGCAGGCGCTGTTGCAGGTGCTGGGCGAAGGTGTAGGCGCCGACATAGCCAAGTACGTCGa                   >  3:221763/1‑149 (MQ=255)
                                                                                                                                  aTCGATTCCTGGTAGCCAATGGCTGCCAGCAGGCGCCAGTCGACCTGCTCGACCTTGGCGTAGCTCTTGAAGTGCTTCTCGTACTTGGGCAGGCGCTGTTGCAGGTGCTGGGCGAAGGTGTAGGCGCCGACATAGCCAAGTACGTCGa                   >  4:214389/1‑148 (MQ=255)
                                                                                                                                    cGATTCCTGGTAGCCAATGGCTGCCAGCAGGCGCCAGTCGACCTGCTCGACCTTGGCGTAGCTCTTGAAGTGCTTCTCGTACTTGGGCAGGCGCTGTTGCAGGTGCTGGGCGAAGgtg                                               >  7:261249/1‑118 (MQ=255)
                                                                                                                                    cGATTCCTGGTAGCCAATGGCTGCCAGCAGGCGCCAGTCGACCTGCTCGACCTTGGCGTAGCTCTTGAAGTGCTTCTCGTACTTGGGCAGGCGCTGTTGCAGGTGCTGGGCGAAGgtg                                               <  8:261249/118‑1 (MQ=255)
                                                                                                                                    cGATTCCTGGTAGCCAATGGCTGCCAGCAGGCGCCAGTCGACCTGCTCGACCTTGGCGTAGCTCTTGAAGTGCTTCTCGTACTTGGGCAGGCGCTGTTGCAGGTGCTGGGCGAAGGTGTAGGCGCCGACATAGCCAAGTACGTCGACAt                >  2:196497/1‑149 (MQ=255)
                                                                                                                                    cGATTCCTGGTAGCCAATGGCTGCCAGCAGGCGCCAGTCGACCTGCTCGACCTTGGCGTAGCTCTTGAAGTGCTTCTCGTACTTGGGCAGGCGCTGTTGCAGGTGCTGGGCGAAGGTGTAGGCGCCGACATAGCCAAGTACGTCGACAt                <  1:253869/149‑1 (MQ=255)
                                                                                                                                     gcctgcTGGCAGCCATTGGCTACCAGCAGGCGCCAGTCGACCTGCTCGACCTTGGCGTAGCTCTTGAAGTGCTTCTCGTACTTGGGCAGGCGCTGTTGCAGGTGCTGGGCGAAGGTGTAGGCGCCGACATAGCCAAGTACGTCGACATg               >  5:42652/6‑149 (MQ=255)
                                                                                                                                      aTTCCTGGTAGCCAATGGCTGCCAGCAGGCGCCAGTCGACCTGCTCGACCTTGGCGTAGCTCTTGAAGTGCTTCTCGTACTTGGGCAGGCGCTGTTGCAGGTGCTGGGCGAAGGTGTAGGCGCCGACATAGCCAAGTACGTCGGCATgg              <  1:207362/149‑1 (MQ=255)
                                                                                                                                       ctgcTGGCAGCCATTGGCTACCAGCAGGCGCCAGTCGACCTGCTCGACCTTGGCGTAGCTCGTGAAGTGCTTCTCGTACTTGGGCAGGCGCTGTTGCAGGTGCTGGGCGAAGGTGTAGGCGCCGACATAGCCAAGTACGTCGACATGGc             >  4:11085/4‑149 (MQ=255)
                                                                                                                                        tCCTGGTAGCCAATGGCTGCCAGCAGGCGCCAGTCGACCTGCTCGACCTTGGCGTAGCTCTTGAAGTGCTTCTCGTACTTGGGCAGGCGCTGTTGCAGGTGCTGGGCg                                                     >  7:344219/1‑108 (MQ=255)
                                                                                                                                        tCCTGGTAGCCAATGGCTGCCAGCAGGCGCCAGTCGACCTGCTCGACCTTGGCGTAGCTCTTGAAGTGCTTCTCGTACTTGGGCAGGCGCTGTTGCAGGTGCTGGGCg                                                     <  8:344219/108‑1 (MQ=255)
                                                                                                                                          cTGGTAGCCAATGGCTTCCAGCAGGCGCCAGTCGACCTGCTCGACCTTGGCGTAGCTCTTGAAGTGCTTCTCGTACTTGGGCAGGCGCTGTTGCAGGTGCTGGGCGAAGGTGTAGGCGCCGACATAGCCAAGTACGTCGACATGGctg           <  3:150902/148‑3 (MQ=255)
                                                                                                                                          cTGGTAGCCAATGGCTGCCAGCAGGCGCCAGTCGACCTGCTCGACCTTGGCGTAGCTCTTGAAGTGCTTCTCGTACTTGGGCAGGCGCTGTTGCAGGTGCTGGGCGAAGGTGTAGGCGCCGACATAGCCAAGTACGTCGACATGGCCg           >  5:288863/1‑148 (MQ=255)
                                                                                                                                           tGGTAGCCAATGGCTGCCAGCAGGCGCCAGTCGACCTGCTCGACCTTGGCGTAGCTCTTGAAGTGCTTCTCGTACTTGGGCAGGCGCTGTTGCAGGTGCTGGGCGAAGGTGTAGGCGCCGACATAGCCAAGTACGTCGACATGGCCg           >  3:33040/1‑147 (MQ=255)
                                                                                                                                             gTAGCCAATGGCTGCCAGCAGGCGCCAGTCGACCTGCTCGACCTTGGCGTAGCTCTTGAAGTGCTTCTCGTACTTGGGCAGGCGCTGTTGCAGGTGCTGGGCGAAGGTGTAGGCGCCGACATAGCCAAGTACGTCGACATGGCCgtag        >  8:53626/1‑148 (MQ=255)
                                                                                                                                              tAGCCAATGGCTGCCAGCAGGCGCCAGTCGACCTGCTCGACCTTGGCGTAGCTCTTGAAGTGCTTCTCGTACTTGGGCAGGCGCTGTTGCAGGTGCTGGGCGAAGGTGTAGGCGCCGACATAGCCAAGTACGTCGACATGGCCgtagta      <  1:242575/149‑1 (MQ=255)
                                                                                                                                              tAGCCAATGGCTGACAGCAGGCGCCAGTCGACCTGCTCGACCTTGGCGTAGCTCTTGAAGTGCTTCTCGTACTTGGGCAGGCGCTGTTGCAGGTGCTTGGCGAAGGTGTAGGCGCCGACATAGCCAAGTACGTCGACATGGCCgtagta      <  5:269423/149‑1 (MQ=255)
                                                                                                                                              aaGCCAATGGCTTCCAGCAGGCGCCAGTCGACCTGCTCGACCTTGGCGTAGCTCTTGAATTTCTTCTCGTACTTGGGCAGGCTCTGTTGCAGGTGGTGGGCGAAGGTGTAGGCGCCGACATAGCCAAGTACGTCGACATGGCCgtagta      <  3:210/148‑1 (MQ=255)
                                                                                                                                               aGCCATTGGCTACCAGCAGGCGCCAGTCGACCTGCTCGACCTTGGCGTAGCTCTTGAAGTGCTTCTCGTACTTGGGCAGGCGCTGTTGCAGGTGCTGGGCGAAGGTGTAGGCGCCGACATAGCCAAGTACGTCGACATGGCCGtagtag     <  3:11085/149‑1 (MQ=255)
                                                                                                                                               aGCCAATGGCTGCCAGCAGGCGCCAGTCGACCTGCTCGACCTTGGCGTAGCTCTTGAAGTGCTTCTCGTACTTGGGCAGGCGCTGTTGCAGGTGCTGGGCGAAGGTGTAGgcgc                                        <  6:276483/114‑1 (MQ=255)
                                                                                                                                               aGCCAATGGCTGCCAGCAGGCGCCAGTCGACCTGCTCGACCTTGGCGTAGCTCTTGAAGTGCTTCTCGTACTTGGGCAGGCGCTGTTGCAGGTGCTGGGCGAAGGTGTAGgcgc                                        >  5:276483/1‑114 (MQ=255)
                                                                                                                                                 ccAATGGCTGCCAGCAGGCGCCAGTCGACCTGCTCGACCTTGGCGTAGCTCTTGAAGTGCTTCTCGTACTTGGGCAGGCGCTGTTGCAGGTGCTGGGCGAAGGTGTAGGCGCCGACATAGCCAAGTACGTCGACATGGCCGTAGTAGCg   <  2:128417/149‑1 (MQ=255)
                                                                                                                                                  cAATGGCTGCCAGCAGGCGCCAGTCGACCTGCTCGACCTTGGCGTAGCTCTTGAAGTGCTTCTCGTACTTGGGCAGGCGCTGTTGCAGGTGCTGGGCGAAGGTGTAGGCGCCGACATAGCCAAGTACGTCGACATGGCCGTAGTAGCgg  >  1:267366/1‑149 (MQ=255)
                                                                                                                                                    aTGGCTGCCAGCAGGCGCCAGTCGACCTGCTCGACCTTGGCGTAGCTCTTGAAGTGCTTCTCGTACTTGGGCAGGCGCTGTTGCa                                                                <  8:265473/85‑1 (MQ=255)
                                                                                                                                                    aTGGCTGCCAGCAGGCGCCAGTCGACCTGCTCGACCTTGGCGTAGCTCTTGAAGTGCTTCTCGTACTTGGGCAGGCGCTGTTGCa                                                                >  7:265473/1‑85 (MQ=255)
                                                                                                                                                    |                                                                                                                                                  
GTACTTGGCGCCGCCCTGGATGCTCTGCCTCGGGTCCAGTCGGTTCGACACGCCCATGGCCTGGGCGGTGCGCTGGGTCAGCATCATCAGGCCGCGCACGCCGGTCTTGGAGGTGACTTCCGGCTGCCACATCGATTCCTGGTAGCCAATGGCTGCCAGCAGGCGCCAGTCGACCTGCTCGACCTTGGCGTAGCTCTTGAAGTGCTTCTCGTACTTGGGCAGGCGCTGTTGCAGGTGCTGGGCGAAGGTGTAGGCGCCGACATAGCCAAGTACGTCGACATGGCCGTAGTAGCGG  >  NC_002947/1183308‑1183602

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