Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,184,158:1 +C 6.6% coding (227/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,1581.C 6.6% 242.5 / 9.7 61coding (227/1458 nt)mltFmembrane‑bound lytic murein transglycosylase F
Reads supporting (aligned to +/- strand):  ref base . (30/27);  new base C (2/2);  total (32/29)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TAGGGATGGGAGTATCTCACCTGGGTTTTGCGCCGCTCACTGCTGACCAGGCCGGCGGCAGCCAGTACCGGGCCGGAGGGCTTGCCGAGGGCGTCGTACAGCTCGTCCAGGTTGTCGGCGGTCTCGATCTGCAGCTTTACGCCGAGATCGTCGGCGAAATGCTGGACCAGTTCGTACTCAAAGCCGGTTTCGCCGTTACGGTCCTGGAAGTAGGTTGCCGGGCTGTTGCGTGTGATCACGCGCAGCACGCCATCCTCCTTCACGCGCTCGAGGGTGCTGGGTTTTTCAACG  >  NC_002947/1184014‑1184304
                                                                                                                                                 |                                                                                                                                                   
tAGGTATGGGAGTATCTCACCTGGTTGTTGCGCCGCTCACTGCTGACCAGGCGGGCTGCAGCCAGTACCGGTCCGGAGGGCTTTCCGAGTGCGTCTTACAGCTCGTCCTGGTTGTCGGCGGTCTCTATCTGCAGCTTTACGCCGAGAtc                                                                                                                                                 <  6:169066/149‑1 (MQ=255)
         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                                                                                                                                  >  3:305543/1‑124 (MQ=255)
                                        tgctgACCAGGCCGGCGGCAGCCAGTACCGGGCCGGAGGGCTTGCCGAGGGCGTCGTACAGCTCGTCCAGGTTGTCGGCGGTCTCGATCTGCAGCTTTACGCCGAGATCGTCGGCGAAATGCTg                                                                                                                                  <  4:305543/124‑1 (MQ=255)
                                             aCCAGGCCGGCGGCAGCCAGTACCGGGCCGGAGGGCTTGCCGAGGGCGTCGTACAGCTCGTCCAGGTTGTCGGCGGTCTCGATCTGCAGCTTTACGCCGAGATCGTCGGCGAAATGCTGGACCAGTTCGTACTCAAAGCCGGTTTCGcc                                                                                                   >  3:252493/1‑149 (MQ=255)
                                             aCCAGGCCGGCGGCAGCCAGTACCGGGCCGGAGGGCTTGCCGAGGGCGTCGTACAGCTCGTCCAGGTTGTCGGCGGTCTCGATCTGCAGCTTTACGCCGAGATCGTCGGCGAAATGCTGGACCAGTTCGTACTCAAAGCCGGTTTCGc                                                                                                    <  5:286258/148‑1 (MQ=255)
                                                         gcAGCCAGTACCGGGCCGGAGGGCTTGCCGAGGGCGTCGTACAGCTCGTCCAGGTTGTCGGCGGTCTCGATCTGCAGCTTTACGCCGAGATCGTCGGCGAAATGCTGGACCAGTTCGTACTCAAAGCCGGTTTCGCCGTTACGGTCCTg                                                                                       <  8:161746/149‑1 (MQ=255)
                                                           aGCCAGTACCGGGCCGGAGGGCTTGCCGAGGGCGTCGTACAGCTCGTCCAGGTTGTCGGCAGTCTCGATCTGCAGCTTTACGCCGAGATCGTCGGCGAAATGCTGGACCAGTTCGTACTCAAAGCCGGTTTCGCCGTTACGGTCCTgg                                                                                      >  7:150811/1‑148 (MQ=255)
                                                           aGCCAGTACCGGGCCGGAGGGCTTGCCGAGGGCGTCGTACAGCTCGTCCAGGTTGTCGGCAGTCTCGATCTGCAGCTTTACGCCGAGATCGTCGGCGAAATGCTGGACCAGTTCGTACTCAAAGCCGGTTTCGCCGTTACGGTCCTgg                                                                                      <  8:150811/148‑1 (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                                                            <  4:252493/149‑1 (MQ=255)
                                                                                    ccGAGGGCGTCGTACAGCTCGTCCAGGTTGTCGGCGGTCTCGATCTGCAGCTTTACGCCGAGATCGTCGGCGAAATGCTGGACCAGTTCGTACTCAAAGCCGGTTTCGCCGTTACGGTCCTGGAAGTAGGTTGCCGGGCTGTTGCgtgt                                                            <  1:203900/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                   <  8:10271/143‑1 (MQ=255)
                                                                                                                                   cAGCTTTACGCCGAGATCGTCGGCGAAATGCTGGACCAGTTCGTACTCAAAGCCGGTTTCGCCGTTACGGTCCTGGAAGTAGGTTGCCGGGCTGTTGCGTGTGATCACGCGCAGCACGCCATCCTCCTTCACGCGCTCGAggg                   <  4:98096/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)
                                                                                                                                                 |                                                                                                                                                   
TAGGGATGGGAGTATCTCACCTGGGTTTTGCGCCGCTCACTGCTGACCAGGCCGGCGGCAGCCAGTACCGGGCCGGAGGGCTTGCCGAGGGCGTCGTACAGCTCGTCCAGGTTGTCGGCGGTCTCGATCTGCAGCTTTACGCCGAGATCGTCGGCGAAATGCTGGACCAGTTCGTACTCAAAGCCGGTTTCGCCGTTACGGTCCTGGAAGTAGGTTGCCGGGCTGTTGCGTGTGATCACGCGCAGCACGCCATCCTCCTTCACGCGCTCGAGGGTGCTGGGTTTTTCAACG  >  NC_002947/1184014‑1184304

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: