Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,184,163 Δ1 bp 6.6% coding (222/1458 nt) 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,184,1630G. 6.6% 198.6 / 7.2 61coding (222/1458 nt)mltFmembrane‑bound lytic murein transglycosylase F
Reads supporting (aligned to +/- strand):  ref base G (31/26);  new base . (2/2);  total (33/28)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.63e-01

GAGTATCTCACCTGGGTTTTGCGCCGCTCACTGCTGACCAGGCCGGCGGCAGCCAGTACCGGGCCGGAGGGCTTGCCGAGGGCGTCGTACAGCTCGTCCAGGTTGTCGGCGGTCTCGATCTGCAGCTTTACGCCGAGATCGTCGGCGAAATGCTGGACCAGTTCGTACTCAAAGCCGGTTTCGCCGTTACGGTCCTGGAAGTAGGTTGCCGGGCTGTTGCGTGTGATCACGCGCAGCACGCCATCCTCCTTCACGCGCTCGAGGGTGCTGGGTTTTTCAACGCAG  >  NC_002947/1184023‑1184307
                                                                                                                                             |                                                                                                                                                 
gAGTATCTCACCTGGGTTTTGCGCCGCTCACTGCTGACCAGGCCGGCGGCAGCCAGTACCGGGCCGGAGGGCTTGCCGAGGGCGTCGTACAGCTCGTCCAGGTTGTCGGCGGTCTCGATCTGCAGCTTTACGCCGAGATCGTCGGCGa                                                                                                                                            >  8:84580/1‑148 (MQ=255)
 aGTATCTCACCTGGGTTTTGCGCCGCTCACTGCTGACCAGGCCGGCGGCAGCCAGTACCGGGCCGGAGGGCTTGCCGAGGGCGTCGTACAGCTCGTCCAGGTTGTCGGCGGTCTCGATCTGCAGCTTTACGCCGAGATCGTCGGCGaaa                                                                                                                                          >  8:31601/1‑149 (MQ=255)
       tcACCTGGGTTTTGCGCCGCTCACTGCTGACCAGGCCGGCGGCAGCCAGTACCGGGCCGGAGGGCTTGCCGAGGGCGTCGTACAGCTCGTCCAGGTTGTCGGCGGTCTCGATCTGCAGCTTTACGCCGAGATCGTCGGCGAAATGGTg                                                                                                                                     >  2:301737/1‑148 (MQ=255)
       tcACCTGGGTTTTGCGCCGCTCACTGCTGACCAGGCCGGCGGCAGCCAGTACCGGGCCGGAGGGCTTGCCGAGGGCGTCGTACAGCTCGTCCAGGTTGTCGGCGGTCTCGATCTGCAGCATTTCGCCGACGATCTCGGCGTAAAGCTgc                                                                                                                                    >  4:47320/1‑148 (MQ=255)
         aCCTGGGTTTTGCGCCGCTCACTGCTGACCAGGCCGGCGGCAGCCAGTACCGGGCCGGAGGGCTTGCCGAGGGCGTCGTACAGCTCGTCCAGGTTGTCGGCGGTCTCGATCTGCAGCTTTACGCCGAGATCGTCGGCGAAATGCTg                                                                                                                                     <  7:84580/146‑1 (MQ=255)
         aCCTGGGTTTTGCGCCGCTCACTGCTGACCAGGCCGGCGGCAGCCAGTACCGGGCCGGAGGGCTTGCCGAGGGCGTCGTACAGCTCGTCCAGGTTGTCGGCGGTCTCGATCTGCAGCATTTCGCCGACGATCTCGGCGTAAAGCTgcag                                                                                                                                  <  3:47320/149‑4 (MQ=38)
         aCCTGGGTTTTGCGCCGCTCACTGCTGACCAGGCCGGCGGCAGCCAGTACCGGGCCGGAGGGCTTGCCGAGGGCGTCGTACAGCTCGTCCAGGTTGTCGGCGGTCTCGATCTGCAGCATTTCGCCGACGATCTCGGCGTAAAGCTgc                                                                                                                                    <  2:55713/147‑2 (MQ=255)
           cTGGGTTTTGCGCCGCTCACTGCTGACCAGGCCGGCGTCAGCCAGTACCGGGACGGAGGGCTTGCCGAGGGCGTCGTACAGCTCGTCCAGGTTGTCGGCGGTCTCGATCTGCAGATTTACGCCGAGATCGTCGGCGAAATGCTGGAc                                                                                                                                  >  2:203900/1‑147 (MQ=255)
           cTGGGTTTTGCGCCGCTCACTGCTGACCAGGCCGGCGGCAGCCAGTACCGGGCCGGAGGGCTTGCCGAGGGCGTCGTACAGCTCGTCCAGGTTGTCGGCGGTCTCGATCTGCAGCTTTACGCCGAGATCGTCGGCGaaa                                                                                                                                          <  4:134000/139‑1 (MQ=255)
           cTGGGTTTTGCGCCGCTCACTGCTGACCAGGCCGGCGGCAGCCAGTACCGGGCCGGAGGGCTTGCCGAGGGCGTCGTACAGCTCGTCCAGGTTGTCGGCGGTCTCGATCTGCAGCTTTACGCCGAGATCGTCGGCGaaa                                                                                                                                          >  3:134000/1‑139 (MQ=255)
           cTGGGTTTTGCGCCGCTCACTGCTGACCAGGCCGGCGGCAGCCAGTACCGGGCCGGAGGGCTTGCCGAGGGCGTCGTACAGCTCGTCCAGGTTGTCGGCGGTCTCGATCTGCAGCTTTACGCCGAGATCGTCGGCGAAATGCTGGAc                                                                                                                                  >  1:110513/1‑147 (MQ=255)
              ggTTTTGCGCCGCTCACTGCTGACCAGGCCGGCGGCAGCCAGTACCGGGCCGGAGGGCTTGCCGAGGGCGTCGTACAGGTCGTCCAGGTTGTCGGCGGTCTCGATCTGCAGCTTTACGCCGAGATCGTCGGCGAAATGCTGGACCAGtt                                                                                                                             <  8:257024/149‑1 (MQ=255)
                             aCTGCTGACCAGGCCGGCGGCAGCCAGTACCGGGCCGGAGGGCTTGCCGAGGGCGTCGTACTGCTCGTCCAGGTTGTCGGCGGTCTCGATCTGCAGCTTTACGCCGAGATCGTCGGCGAAATGCTGGACCAGTTCGTACTCAAAGCCgg                                                                                                             >  1:89708/1‑149 (MQ=255)
                             aCTGCTGACCAGGCCGGCGGCAGCCAGTACCGGGCCGGAGGGCTTGCCGAGGGCGTCGTACAGCTCGTCCAGGTTGTCGGCGGTCTCGATCTGCAGCTTTACGCCGAGATCGTCGGCGAAATGCTg                                                                                                                                     <  6:90200/126‑1 (MQ=255)
                             aCTGCTGACCAGGCCGGCGGCAGCCAGTACCGGGCCGGAGGGCTTGCCGAGGGCGTCGTACAGCTCGTCCAGGTTGTCGGCGGTCTCGATCTGCAGCTTTACGCCGAGATCGTCGGCGAAATGCTg                                                                                                                                     >  5:90200/1‑126 (MQ=255)
                             aCTGCTGACCAGGCCGGCGGCAGCCAGTACCGGGCCGGAGGGCTTGCCGAGGGCGTCGTACAGCTCGTCCAGGTTGTCGGCGGTCTCGATCTGCAGCTTTACGCCGAGATCGTCGGCGAAATGCTGGACCAGTTCGTACTCAAAGCCg                                                                                                              <  7:31601/148‑1 (MQ=255)
                               tgctgACCAGGCCGGCGGCAGCCAGTACCGGGCCGGAGGGCTTGCCGAGGGCGTCGTACAGCTCGTCCAGGTTGTCGGCGGTCTCGATCTGCAGCTTTACGCCGAGATCGTCGGCGAAATGCTg                                                                                                                                     <  4:305543/124‑1 (MQ=255)
                               tgctgACCAGGCCGGCGGCAGCCAGTACCGGGCCGGAGGGCTTGCCGAGGGCGTCGTACAGCTCGTCCAGGTTGTCGGCGGTCTCGATCTGCAGCTTTACGCCGAGATCGTCGGCGAAATGCTg                                                                                                                                     >  3:305543/1‑124 (MQ=255)
                                    aCCAGGCCGGCGGCAGCCAGTACCGGGCCGGAGGGCTTGCCGAGGGCGTCGTACAGCTCGTCCAGGTTGTCGGCGGTCTCGATCTGCAGCTTTACGCCGAGATCGTCGGCGAAATGCTGGACCAGTTCGTACTCAAAGCCGGTTTCGcc                                                                                                      >  3:252493/1‑149 (MQ=255)
                                    aCCAGGCCGGCGGCAGCCAGTACCGGGCCGGAGGGCTTGCCGAGGGCGTCGTACAGCTCGTCCAGGTTGTCGGCGGTCTCGATCTGCAGCTTTACGCCGAGATCGTCGGCGAAATGCTGGACCAGTTCGTACTCAAAGCCGGTTTCGc                                                                                                       <  5:286258/148‑1 (MQ=255)
                                                gcAGCCAGTACCGGGCCGGAGGGCTTGCCGAGGGCGTCGTACAGCTCGTCCAGGTTGTCGGCGGTCTCGATCTGCAGCTTTACGCCGAGATCGTCGGCGAAATGCTGGACCAGTTCGTACTCAAAGCCGGTTTCGCCGTTACGGTCCTg                                                                                          <  8:161746/149‑1 (MQ=255)
                                                  aGCCAGTACCGGGCCGGAGGGCTTGCCGAGGGCGTCGTACAGCTCGTCCAGGTTGTCGGCAGTCTCGATCTGCAGCTTTACGCCGAGATCGTCGGCGAAATGCTGGACCAGTTCGTACTCAAAGCCGGTTTCGCCGTTACGGTCCTgg                                                                                         <  8:150811/148‑1 (MQ=255)
                                                  aGCCAGTACCGGGCCGGAGGGCTTGCCGAGGGCGTCGTACAGCTCGTCCAGGTTGTCGGCAGTCTCGATCTGCAGCTTTACGCCGAGATCGTCGGCGAAATGCTGGACCAGTTCGTACTCAAAGCCGGTTTCGCCGTTACGGTCCTgg                                                                                         >  7:150811/1‑148 (MQ=255)
                                                   gCCAGTACCGGGCCGGAGGGCTTGCCGAGGGCGTCGTACAGCTCGTCCAGGTTGTCGGCGGTCTCGATCTGCAGCTTTACGCCGAGATCGTCGGCGAAATGCTGGACCAGTTCGTACTCAAAGCCGGTTTCGCCGTTACGGTCCTGGaa                                                                                       <  7:28952/149‑1 (MQ=255)
                                                        tACCGGGCCGGAGGGCTTGCCGAGGGCGTCGTACAGCTCGTCCAGGTTGTCGGCGGTCTCGATCTGCAGCTTTACGCCGAGATCGTCGGCGAAATGCTGGACCAGTTCGTACTCAAAGCCGGTTTCg                                                                                                        <  3:188666/127‑1 (MQ=255)
                                                        tACCGGGCCGGAGGGCTTGCCGAGGGCGTCGTACAGCTCGTCCAGGTTGTCGGCGGTCTCGATCTGCAGCTTTACGCCGAGATCGTCGGCGAAATGCTGGACCAGTTCGTACTCAAAGCCGGTTTCg                                                                                                        >  4:188666/1‑127 (MQ=255)
                                                                 ggAGGGCTTGCCGAGGGCGTCGTACAGCTCGTCCAGGTTGTCGGCGGTCTCGATCTGCAGCTTTACGCCGAGATCGTCGGCGAAATGCTGGACCAGTTCGTACTCAAAGCCGGTTTCGCCGTTACGGTCCTGGAAGTAGGTTGCCGGGc                                                                         <  5:262401/149‑1 (MQ=255)
                                                                       cTTGCCGAGGGCGTCGTACAGCTCGTCCAGGTTGTCGGCGGTCTCGATCTGCAGCTTTACGCCGAGATCGTCGGCGAAATGCTGGACCAGTTCGTACTCAAAGCCGGTTTCGCCGTTACGGTCCTGGAAGTAGGTTGCCGGGCTGTTGc                                                                   <  2:110513/149‑1 (MQ=255)
                                                                        ttGCCGAGGGCGTCGTACAGCTCGTCCAGGTTGTCGGCGGTCTCGATCTGCAGCTTTACGCCGAGATCGTCGGCGAAATGCTGGACCAGTTCGTACTCAAAGCCGGTTTCGCCGTTACGGTCCTGGAAGTAGGTTGCCGGGCTGTTGCg                                                                  >  3:138760/1‑149 (MQ=255)
                                                                           ccGAGGGCGTCGTACAGCTCGTCCAGGTTGTCGGCGGTCTCGATCTGCAGCTTTACGCCGAGATCGTCGGCGAAATGCTGGACCAGTTCGTACTCAAAGCCGGTTTCGCCGTTACGGTCCTGGAAGTAGGTTGCCGGGCTGTTGCgtgt                                                               <  1:203900/149‑1 (MQ=255)
                                                                           ccGAGGGCGTCGTACAGCTCGTCCAGGTTGTCGGCGGTCTCGATCTGCAGCTTTACGCCGAGATCGTCGGCGAAATGCTGGACCAGTTCGTACTCAAAGCCGGTTTCGCCGTTACGGTCCTGGAAGTAGGTTGCCGGGCTGTTGCgtgt                                                               <  4:252493/149‑1 (MQ=255)
                                                                             gAGGGCGTCGTACTGCTCGTCCAGGTTGTCGGCGGTCTCGATCTGCAGCTTTACGCCGAGATCGTCGGCGAAATGCTGGACCAGTTCGTACTCAAAGCCGGTTTCGCCGTTACGGTCCTGGAAGTAGGTTGCCGGGCTGTTGCGTGTGa                                                             <  2:89708/149‑1 (MQ=255)
                                                                               gggCGTCGTACAGCTCGTCCAGGTTGTCGGCGGTCTCGATCTGCAGCTTTACGCCGAGATCGTCGGCGAAATGCTGGACCAGTTCGTACTCAAAGCCGGTTTCGCCGTTACGGTCCTGGAAGTAGGTTGCCGGGCTGTTGCGTGTGATc                                                           <  8:40531/149‑1 (MQ=255)
                                                                                     cgtACAGTTCGTCTAGGTTGTAGGCTGTCTCGATCTGCAGCTTTAGGCCTAGATCGTCGGCGAAATGCTTTTCCTGTTCGTACTCAAAGCCTTTTTCGCCGTTACGTTCCTGTAACTATGTTGCCGGGCTGTTGCGTGTGATCActcg                                                      <  4:174973/148‑1 (MQ=255)
                                                                                               gTCCAGGTTGTCGGCGGTCTCGATCTGCAGCTTTACGCCGAGATCGTCGGCGAAATGCTGGACCAGTTCGTACTCAAAGCCGGTTTCGCCGTTACGGTCCTGGAAGTAGGTTGCCGGGCTGTTGCGTGTGATCACGCGCAGCACGCCAt                                           >  3:237291/1‑149 (MQ=255)
                                                                                                    ggTTGTCGGCGGTCTCGATCTGCAGCTTTACGCCGAGATCGTCGGCGAAATGCTGGACCAGTTCGTACTCAAAGCCGGTTTCGCCGTTACGGTCCTGGAAGTAGGTTGCCGGGCTGTTGCGTGTGATCACGCGCAGCACGCCAtcctcc                                      >  8:206642/1‑149 (MQ=255)
                                                                                                     gTTGTCGGCGGTCTCGATCTGCAGCTTTACGCCGAGATCGTCGGCGAAATGCTGGACCAGTTCGTACTCAAAGCCGGTTTCGCCGTTACGGTCCTGGAAGTAGGTTGCCGGGCTGTTGCGTGTGATCACGCGCAGCACGCCATcctcct                                     >  3:98096/1‑149 (MQ=255)
                                                                                                        gTCGGCGGTCTCGATCTGCAGCTTTACGCCGAGATCGTCGGCGAAATGCTGGACCAGTTCGTACTCAAAGCCGGTTTCGCCGTTACGGTCCTGGAAGTAGGTTGCCGGGCTGTTGCGTGTGATCACGCGCAGCACGCCa                                            >  1:235437/1‑139 (MQ=255)
                                                                                                        gTCGGCGGTCTCGATCTGCAGCTTTACGCCGAGATCGTCGGCGAAATGCTGGACCAGTTCGTACTCAAAGCCGGTTTCGCCGTTACGGTCCTGGAAGTAGGTTGCCGGGCTGTTGCGTGTGATCACGCGCAGCACGCCa                                            <  2:235437/139‑1 (MQ=255)
                                                                                                            gcggTCTCGATCTGCAGCTTTACGCCGAGATCGTCGGCGAAATGCTGGACCAGTTCGTACTCAAAGCCGGTTTc                                                                                                         >  7:215009/1‑74 (MQ=255)
                                                                                                            gcggTCTCGATCTGCAGCTTTACGCCGAGATCGTCGGCGAAATGCTGGACCAGTTCGTACTCAAAGCCGGTTTc                                                                                                         <  8:215009/74‑1 (MQ=255)
                                                                                                               gTCTCGATCTTCAGCTTTACGCCGAGATCGTCGGCGAAATGCTGGATCAGTTCGTACTCAAAGCCGGTTTCGCCGTTACGGTCCTgg                                                                                         <  1:145745/87‑1 (MQ=255)
                                                                                                               gTCTCGATCTGCAGCTTTACGCCGAGATCGTCGGCGAAATGCTGGATCAGTTCGTACTCAAAGCCGGTTTCGCCGTTACGGTCCTgg                                                                                         >  2:145745/1‑87 (MQ=255)
                                                                                                                tctcGATCTGCAGCTTTACGCCGAGATCGTCGGCGAAATGCTGGACCAGTTCGTACTCAAAGCCGGTTTCGCCGTTACGGTCCTGGAAGTAGGTTGCCGGGCTGTTGCGTGTGATCACGCGCAGCACGCCATCCTCCTTCACGCGCTc                           >  2:79899/1‑148 (MQ=255)
                                                                                                                 ctcGATCTGCAGCTTTACGCCGAGATCGTCGGCGAAATGCTGGACCAGTTCGTACTCAAAGCCGGTTTCGCCGTTACGGTCCTGGAAGTAGGTTGCCGGGCTGTTGCGTGTGATCACGCGCAGCACGCCATCCTCCTTCACGCGCTc                           >  3:105770/1‑147 (MQ=255)
                                                                                                                  tcGATCTGCAGCTTTACGCCGAGATCGTCGGCGAAATGCTGGACCAGTTCGTACTCAAAGCCGGTTTCGCCGTTACGGTCCTGGAAGTAGGTTGCCGGGCTGTTGCGTGTGATCACGCGCAGCACGCCATCCTCCTTCACGCGCTCGAg                        >  7:10271/1‑149 (MQ=255)
                                                                                                                    gATCTGCAGCTTTACGCCGAGATCGTCGGCGAAATGCTGGACCAGTTCTTACTCAAAGCCGGTTTCGCCGTTACGGTCCTGGAAGTAGGTTGCCGGGCTGTTGCGTGTGATCACGCGCAGCACGCCATCCTCCTTCACGCGCTCGAgg                       <  7:267640/148‑1 (MQ=255)
                                                                                                                    gATCTGCAGCTTTACGCCGAGATCGTCGGCGAAATGCTGGACCAGTTCGTACTCAAAGCCGGTTTCGCCGTTACGGTCCTGGAAGTAGGTTGCCGGGCTGTTGCGTGTGATCACGCGCAGCACGCCATCCTCCTTCACGCGCTCGAggg                      >  4:263508/1‑149 (MQ=255)
                                                                                                                    gATCTGCAGCTTTACGCCGAGATCGTCGGCGAAATGCTGGACCAGTTCGTACTCAAAGCCGGTTTCGCCGTTACGGTCCTGGAAGTAGGTTGCCGGGCTGTTGCGTGTGATCACGCGCAGCACGCCATCCTCCTTCACGCGCTCGAggg                      >  4:259130/1‑149 (MQ=255)
                                                                                                                    gATCTGCAGCTTTACGCCGAGATCGTCGGCGAAATGCTGGACCAGTTCGTACTCAAAGCCGGTTTCGCCGTTACGGTCCTGGAAGTAGGTTGCCGGGCTGTTGCGTGTGATCACGCGCAGCACGCCATCCTCCTTCACGCGCTCGAgg                       <  7:206642/148‑1 (MQ=255)
                                                                                                                          cAGCTTTACGCCGAGATCGTCGGCGAAATGCTGGACCAGTTCGTACTCAAAGCCGGTTTCGCCGTTACGGTCCTGGAAGTAGGTTGCCGGGCTGTTGCGTGTGATCACGCGCAGCACGCCATCCTCCTTCACGCGCTCGAggg                      <  4:98096/143‑1 (MQ=255)
                                                                                                                          cAGCTTTACGCCGAGATCGTCGGCGAAATGCTGGACCAGTTCGTACTCAAAGCCGGTTTCGCCGTTACGGTCCTGGAAGTAGGTTGCCGGGCTGTTGCGTGTGATCACGCGCAGCACGCCATCCTCCTTCACGCGCTCGAggg                      <  8:10271/143‑1 (MQ=255)
                                                                                                                          cAGCTTTACGCCGAGATCGTCGGCGAAATGCTGGACCAGTTCGTACTCAAAGCCGGTTTCGCCGTTACGGTCCTGGAAGTAGGTTGCCGGGCTGTTGCGTGTGATCACGCGCAGCACGCCATCCTCCTTCACGCGCTCGAGGGTGCTgg                >  3:157399/1‑149 (MQ=255)
                                                                                                                           aGCATTTCGCCGACGATCTCGGCGTAAAGCTGGACCAGTTCGTACTCAAAGCCTAGTTTCGCCGTTACGGTCATGGAAGTAGGTTGCCGGCCTGTTGCGTGTGATCACGCGCAGCACGCCATCCTCCTTCACGCGCTCGAGGGTGCTgg                >  1:241875/1‑149 (MQ=255)
                                                                                                                                   gCCGAGATCGTCGGCGAAATGCTGGACCAGTTCGTACTCAAAGCCGGTTTCGCCGTTACGGTCCTGGAAGTAGGTTGCCGGGCTGTTGCGTGTGATCACGCGCAGCACGCCATCCTCCTTCACGCGCTCGAGGGTGCTGGGTTTTTCaa       >  6:288508/1‑149 (MQ=255)
                                                                                                                                   gCCGAGATCGTCGGCGAAATGCTGGACCAGTTCGTACTCAAAGCCGGTTTCGCCGTTACGGTCCTGGAAGTAGGTTGCCGGGCTGTTGCGTGTGATCACGCGCAGCACGCCATCCTCCTTCACGCGCTCGAGGGTGCTGGGTTTTTCaa       >  2:1254/1‑149 (MQ=255)
                                                                                                                                     cGAGATCGTCGGCGAAATGCTGGACCAGTTCGTACTCAAAGCCGGTTTCGCCGTTACGGTCCTGGAAGTAGGTTGCCGGGCTGTTGCGTGTGATCACGCGCAGCACGCGATCCTCCTTCACGCGCTCGAGGGTGCTGGGTTTTTc         <  1:212317/145‑1 (MQ=255)
                                                                                                                                     cGAGATCGTCGGCGAAATGCTGGACCAGTTCGTACTCAAAGCCGGTTTCGCCGTTACGGTCCTGGAAGTAGGTTGCCGGGCTGTTGCGTGTGATCACGCGCAGCACGCCATCCTCCTTCACGCGCTCGAGGGTGCTGGGTTTTTc         >  2:212317/1‑145 (MQ=255)
                                                                                                                                     cGAGATCGTCGGCGAAATGCTGGACCAGTTCGTACTCAAAGCCGGTTTCGCCGTTACGGTCCTGGAAGTAGGTTGCCGGGCTGTTGCGTGTGATCACGCGCAGCACGCCATCCTCCTTCACGCGCTCGAGGGTGCTGGGTTTTTCAACg     >  4:187899/1‑149 (MQ=255)
                                                                                                                                     cGAGATCGTCGGCGAAATGCTGGACCAGTTCGTACTCAAAGCCGGTTTCGCCGTTACGGTCCTGGAAGTAGGTTGCAGGGCTGTTGCGTGTGATCACGCGCAGCACGCCATCCTCCTTCACGCGCTCGAGGGTGCTGGGTTTTTCAACg     >  2:102045/1‑149 (MQ=255)
                                                                                                                                         gaTCGTCGGAGAAATGCTGGACCAGTTCGTACTCAAAGCCGGTTTCGCCGTTACGGTCCTGGAAGTAGGTTGACGGGCTGTTGCGTGTGATCACGCGCAGAACGCCATCCTCCTTCACGCGCTCGAGGGTGATGGGTTTTTCAACGCAg  >  6:258734/1‑149 (MQ=255)
                                                                                                                                             |                                                                                                                                                 
GAGTATCTCACCTGGGTTTTGCGCCGCTCACTGCTGACCAGGCCGGCGGCAGCCAGTACCGGGCCGGAGGGCTTGCCGAGGGCGTCGTACAGCTCGTCCAGGTTGTCGGCGGTCTCGATCTGCAGCTTTACGCCGAGATCGTCGGCGAAATGCTGGACCAGTTCGTACTCAAAGCCGGTTTCGCCGTTACGGTCCTGGAAGTAGGTTGCCGGGCTGTTGCGTGTGATCACGCGCAGCACGCCATCCTCCTTCACGCGCTCGAGGGTGCTGGGTTTTTCAACGCAG  >  NC_002947/1184023‑1184307

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