Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,517,283 T→G 8.0% V158G (GTG→GGG)  murE → UDP‑N‑acetylmuramoyl‑L‑alanyl‑D‑glutamate‑‑2, 6‑diaminopimelate ligase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029471,517,2830TG 8.0% 168.2 / 4.5 60V158G (GTG→GGG) murEUDP‑N‑acetylmuramoyl‑L‑alanyl‑D‑glutamate‑‑2, 6‑diaminopimelate ligase
Reads supporting (aligned to +/- strand):  ref base T (23/32);  new base G (3/2);  total (26/34)
Fisher's exact test for biased strand distribution p-value = 6.44e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

ACCAACGGCAAGACAAGTGTTACCCAACTGGTCGCCCAGGCGCTCGACCTGCTTGGCCAGCGCTGCGGCCTGATTGGCACCCTGGGTACCGGTTTCTACGGCGAGTTGCAGAGCGGCCGCCTGACCACGCCCGACCCGATTGCCGTGCAATCGACCCTCAACGACCTCAAGAAAGGCGGCGCCAGGGCCGTGGCCATGGAAGTGTCTTCCCATGC‑‑CCTGGAGCAAGGCCGGGTTGCCGCACTGGAATTCGACATCGCGGTGATGACCAACCTGTCCCGCGACC  >  NC_002947/1517138‑1517420
                                                                                                                                                 |                                                                                                                                           
aCCAACGGCAAGACAAGTGTTACCCAACTGGTCGCCCAGGCGCTCGACCTGCTTGGCCAGCGCTGCGGCCTGATTGGCACCCTGGGTACCGGTTTCTACGGCGAGTTGCAGAGCGGCCGCCTGACCACGCCCGACCCGATTGCCGTGc                                                                                                                                           >  2:109217/1‑148 (MQ=255)
 ccAACGGCAAGACAAGTGTTACCCAACTGGTCGCCCAGGCGCTCGACCTGCTTGGCCAGCGCTGCGGCCTGATTGGCACCCTGGGTACCGGTTTCTACGGCGAGTTGCAGAGCGGCCGCCTGACCACGCCCGACCCGATTGCCGTGCaa                                                                                                                                         <  3:319016/149‑1 (MQ=255)
         aaGACAAGTGTTACCCAACTGGTCGCCCAGGCGCTCGACCTGCTTGGCCAGCGCTGCGGCCTGATTGGCACCCTGGGTACCGGTTTCTACGGCGAGTTGCAGAGCGGCCGCCTGACCACGCCCGACCCGATTGCCGTGCAATCGACCCt                                                                                                                                 >  6:138453/1‑149 (MQ=255)
                gtgtTACCCAACTGGTCGCCCAGGCGCTCGACCTGCTTGGCCAGCGCTGCGGCCTGATTGGCACCCTGGGTACCGGTTTCTACGGCGAGTTGCAGAGCGGCCGCCTGACCACGCCCGACCCGATTGCCGGGCAATCGACCCTCAACGAc                                                                                                                          <  6:103722/149‑1 (MQ=255)
                gtgtTACCCAACTGGTCGCCCAGGCGCTCGACCTGCTTGGCCAGCGCTGCGGCCTGATTGGCACCCTGGGTACCGGTTTCTACGGCGAGTTGCAGAGCGGCCGCCTGACCACGCCCGACCCGATTGCACGGCAATCgggtc                                                                                                                                  >  8:258799/1‑137 (MQ=255)
                gtgtTACCCAACTGGTCGCCCAGGCGCTCGACCTGCTTGGCCAGCGCTGCGGCCTGATTGGCACCCTGGGTACCGGTTTCTACGGCGAGTTGCAGAGCGGCCGCCTGACCACGCCCGACCCGATTGCACGGCAATCgggtc                                                                                                                                  <  7:258799/141‑5 (MQ=255)
                 tgtTACCCAACTGGTCGCCCTGGCGCTCGACCTGCTTGGCCAGCGCTGCGGCCTGATTGGCACCCTGGGTACCGGTTTCTACGGCGAGTTGCAGAGCGGCCGCCTGACCACGCCCGACCCGATTGCCGTGCAATCGACCCTCAACGAcc                                                                                                                         <  5:45447/149‑1 (MQ=255)
                   ttACCCAACTGGTCGCCCAGGCGCTCGACCTGCTTGGCCAGCGCTGCGGCCTGATTGGCACCCTGGGTACCGGTTTCTACGGCGAGTTGCAGAGCGGCCGCCTGACCACGCCCGACCCGATTGCCGTGCAATCGACCCTCAACGACCTc                                                                                                                       <  5:155150/149‑1 (MQ=255)
                     aCCCAACTGGTCGCCCAGGCGCTCGACCTGCTTGGCCAGCGCTGCGGCCTGATTGGCACCCTGGGTACCGGTTTCTACGGCGAGTTGCAGAGCGGCCGCCTGACCACGCCCGACCCGATTGCCGTGCAATCGACCCTCAACGACCTc                                                                                                                       >  6:197415/1‑147 (MQ=255)
                      cccAACTGGTCGCCCAGGCGCTCGACCTGCTTGGCCAGCGCTGCGGCCTGATTGGCACCCTGGGTACCGGTTTCTACGGCGAGTTGCAGAGCGGCCGCCTGACCACGCCCGACCCGATTGCCGTGCAATCGACCCTCAACGACCTCaa                                                                                                                     <  1:331480/148‑1 (MQ=255)
                       ccAACTGGTCGCCCAGGCGCTCGACCTGCTTGGCCAGCGCTGCGGCCTGATTGGCACCCTGGGTACCGGTTTCTACGGCGAGTTGCAGAGCGGCCGCCTGACCACGCCCGACCCGATTGCCGTGCAATCGACCCTCAACGACCTCAAg                                                                                                                    >  1:77412/1‑148 (MQ=255)
                       ccAACTGGTCGCCCAGGCGCTCGACCTGCTTGGCCAGCGCTGCGGCCTGATTGGCACCCTGGGTACCGGTTTCTACGGCGAGTTGCAGAGCGGCCGCCTGACCACGCCCGACCCGATTGCCGTGCAATCGACCCTCAACGACCTCAAg                                                                                                                    >  8:183104/1‑148 (MQ=255)
                         aaCTGGTCGCCCAGGCGCTCGACCTGCTTGGCCAGCGCTGCGGCCTGATTGGCACCCTGGGTACCGGTTTCTACGGCGAGTTGCAGAGCGGCCGCCTGACCACGCCCGACCCGATTGCCGTGCAATCGACCCTCAACGACCTCAAGa                                                                                                                   >  7:161737/1‑147 (MQ=255)
                                   ccAGGCGCTCGACCTGCTTGGCCAGCGCTGCGGCCTGATTGGCACCCTGGGTACCGGTTTCTACGGCGAGTTGCAGAGCGGCCGCCTGACCACGCCCGACCCGATTGCCGTGCAATCGACCCTCaa                                                                                                                              <  2:155046/126‑1 (MQ=255)
                                   ccAGGCGCTCGACCTGCTTGGCCAGCGCTGCGGCCTGATTGGCACCCTGGGTACCGGTTTCTACGGCGAGTTGCAGAGCGGCCGCCTGACCACGCCCGACCCGATTGCCGTGCAATCGACCCTCaa                                                                                                                              >  1:155046/1‑126 (MQ=255)
                                   ccAGGCGCTCGACCTGCTTGGCCAGCGCTGCGGCCTGATTGGCACCCTGGGTACCGGTTTCTACGGCGAGTTGCAGAGCGGCCGCCTGACCACGCCCGACCCGATTGCCGTGCAATCGACCCTCAACGACCTCAAGAAAGGCGGCGCCa                                                                                                       <  3:345368/149‑1 (MQ=255)
                                     aGGCGCTCGACCTGCTTGGCCAGCGCTGCGGCCTGATTGGCACCCTGGGTACCGGTTTCTACGGCGAGTTGCAGAGCGGCCGCCTGACCACGCCCGACCCGATTGCCGTGCAATCGACCCTCAACGACCTCAAGAAAGGCGGCGCCAgg                                                                                                     >  2:111759/1‑149 (MQ=255)
                                        cgcTCGACCTGCTTGGCCAGCGCTGCGGCCTGATTGGCACCCTGGGTACCGGTTTCTACGGCGAGTTGCAGAGCGGCCGCCTGACCACGCCCGACCCGATTGCCGTGCAATCGAcc                                                                                                                                   >  8:214223/1‑116 (MQ=255)
                                        cgcTCGACCTGCTTGGCCAGCGCTGCGGCCTGATTGGCACCCTGGGTACCGGTTTCTACGGCGAGTTGCAGAGCGGCCGCCTGACCACGCCCGACCCGATTGCCGTGCAATCGAcc                                                                                                                                   <  7:214223/116‑1 (MQ=255)
                                         gcTCGACCTGCTTGGCCAGCGCTGCGGCCTGATTGGCACCCTGGGTACCGGTTTCTACGGCGAGTTGCAGAGCGGCCGCCTGACCACGCCCGACCCGATTGCCGTGCAATCGACCCTCAACGACCTCAAGAAAGGCGGCGCCAGGGCCg                                                                                                 >  5:25239/1‑149 (MQ=255)
                                             gACCTGCTTGGCCAGCGCTGCGGCCTGATTGGCACCCTGGGTACCGGTTTCTACGGCGAGTTGCAGAGCGGCCGCCTGACCACGCCCGACCCGATTGCCGTGCAATCGACCCTCAACGACCTCAAGAAAGGCGGCGCCAGGGCCGTGGc                                                                                             <  5:197415/149‑1 (MQ=255)
                                                      ggCCAGCGCTGCGGCCTGATTGGCACCCTGGGTACCGGTTTCTACGGCGAGTTGCAGAGCGGCCGCCTGACCACGCCCGACCCGATTGCCGTGCAATCGACCCTCAACGACCTCAAGAAAGGCGGCGCCAGGGCCGTGGCCATGGaa                                                                                      <  6:5700/147‑1 (MQ=255)
                                                      ggCCAGCGCTGCGGCCTGATTGGCACCCTGGGTACCGGTTTCTACGGCGAGTTGCAGAGCGGCCGCCTGACCACGCCCGACCCGATTGCCGTGCAATCGACCCTCAACGACCTCAAGAAAGGCGGCGCCAGGGCCGTGGCCATGGaa                                                                                      <  1:109217/147‑1 (MQ=255)
                                                       gCCAGCGCTGCGGCCTGATTGGCACCCTGGGTACCGGTTTCTACGGCGAGTTGCAGAGCGGCCGCCTGACCACGCCCGACCCGATTGCCGTGCAATCGACCCTCAACGACCTCAAGAAAGGCGGCGCCAGGGCCGTGGCCATGGAAgt                                                                                    <  2:77412/148‑1 (MQ=255)
                                                       gCCAGCGCTGCGGCCTGATTGGCACCCTGGGTACCGGTTTCTACGGCGAGTTGCAGAGCGGCCGCCTGACCACGCCCGACCCGATTGCCGTGCAATCGACCCTCAACGACCTCAAGAAAGGCGGCGCCAGGGCCGTGGCCATGGAAgt                                                                                    <  5:138453/148‑1 (MQ=255)
                                                           gcgcTGCGGCCTGATTGGCACCCTGGGTACCGGTTTCTACGGCGAGTTGCAGAGCGGCCGCCTGACCACGCCCGACCCGATTGCCGTGCAATCGACCCTCAACGACCTCAAGAAAGGCGGCGCCAGGGCCGTGGCCATGGAAGTGTCtt                                                                               <  4:343564/149‑1 (MQ=255)
                                                              cTGCGGCCTGATTGGCACCCTGGGTACCGGTTTCTACGGCGAGTTGCAGAGCGGCCGCCTGACCACGCCCGACCCGATTGCCGTGCAATCGACCCTCAACGACCTCAAGAAAGGCGGCGCCAGGGCCGTGGCCATGGAAGTGTCTTccc                                                                            <  5:250395/149‑1 (MQ=255)
                                                                    ccTGATTGGCACCCTGGGTACCGGTTTCTACGGCGAGTTGCAGAGCGGCCGCCTGACCACGCCCGACCCGATTGCCGTGCAATCGACCCTCAACGACCTCaa                                                                                                                     <  2:216048/102‑1 (MQ=255)
                                                                    ccTGATTGGCACCCTGGGTACCGGTTTCTACGGCGAGTTGCAGAGCGGCCGCCTGACCACGCCCGACCCGATTGCCGTGCAATCGACCCTCAACGACCTCaa                                                                                                                     >  1:216048/1‑102 (MQ=255)
                                                                       gATTGGCACCCTGGGTACCGGTTTCTACGGCGAGTTGCAGAGAGGCCGCCTGACCACGCCCGACCCGATTGCCGTGCAATCGACCCTCAACGACCTCAAGAAAGGCGCCGCCAGGGCCGTGGCCATGGAAGTGTCTTCCCATGC‑‑CCTg                                                                  >  4:159824/1‑148 (MQ=255)
                                                                             cACCCTGGGTACCGGTTTCTACGGCGAGTTGCAGAGCGGCCGCCTGACCACGCCCGACCCGATTGCCGTGCAATCGACCCTCAACGACCTCAAGAAAGGCGGCGCCAGGGCCGTGGCCATGGAAGTGTCTTCCCATGC‑‑CCTGGAGCAAg                                                           <  8:161737/149‑1 (MQ=255)
                                                                               cccTGGGTACCGGTTTCTACGGCGAGTTGCAGAGCGGCCGCCTGACCACGCCCGACCCGATTGCCGTGCAATCGACCCTCAACGACCTCAAGAAAGGCGGCGCCAGGGCCGTGGCCATGGAAGTGTCTTCCCATGC‑‑CCTGGAGCAAGGc                                                         >  8:160188/1‑149 (MQ=255)
                                                                                 cTGGGTACCGGTTTCTACGGCGAGTTGCAGAGCGGCCGCCTGACCACGCCCGACCCGATTGCCGTGCAATCGACCCTCAACGACCTCAAGAAAGGCGGCGCCAGGGCCGTGGCCATGGAAGTGTCTTCCCATGC‑‑CCTGGAGCAAGGCCg                                                       >  7:130402/1‑149 (MQ=255)
                                                                                    ggTACCGGTTTCTACGGCGAGTTGCAGAGCGGCCGCCTGACCACGCCCGACCCGATTGCCGTGCAATCGACCCTCAACGACCTCAAGAAAGGCGGCGCCAGGGCCGTGGCCATGGAAGTGTCTTCCCATGC‑‑CCTGGAGCAAGGCCgg                                                      <  4:158832/147‑1 (MQ=255)
                                                                                    ggTACCGGTTTCTACGGCGAGTTGCAGAGCGGCCGCCTGACCACGCCCGACCCGATTGCCGTGCAATCGACCCTCAACGACCTCAAGAAAGGCGGCGCCAGGGCCGTGGCCATGGAAGTGTCTTCCCATGC‑‑CCTGGAGCAAGGCCgg                                                      <  6:25239/147‑1 (MQ=255)
                                                                                        ccGGTTTCTACGGCGAGTTGTACAGCGGCCGCCTGACCACGCCCGACCCGATTCCCGTGCAATCGACCCTCAACGACCTCAAGAAAGGCGGCGCCAGGGCCGTGGCCATGGAAGTGTCTTCCCATGC‑‑CCTGGAGCAAGGCCGGGTTGcc                                                <  4:221745/149‑1 (MQ=255)
                                                                                         cGGTTTCTACGGCGAGTTGCAGAGCGGCCGCCTGCCCACGCCCGACCCGATTGCCGTGCAATCGACCCTCAACGCCCTCAAGAAAGGCGGCGCCAGGGCCGTGGCCATGGAAGTGTCGTCCCCTGC‑‑CCTGGAGCAAGGCCGCGTGGCCg                                               >  7:126635/1‑149 (MQ=255)
                                                                                         cGGTTTCTACGGCGAGTTGCAGAGCGGCCGCCTGACCACGCCCGACCCGATTGCCGTGCAATCGACCCTCAACGACCTCAAGAAAGGCGGCGCCAGGGCCGTGGCCATGGAAGTGTCTTCCCATGC‑‑CCTGGAGCAAg                                                           >  4:312888/1‑137 (MQ=255)
                                                                                         cGGTTTCTACGGCGAGTTGCAGAGCGGCCGCCTGACCACGCCCGACCCGATTGCCGTGCAATCGACCCTCAACGACCTCAAGAAAGGCGGCGCCAGGGCCGTGGCCATGGAAGTGTCTTCCCATGC‑‑CCTGGAGCAAg                                                           <  3:312888/137‑1 (MQ=255)
                                                                                         cGGTTTCTACGGCGAGTTGCAGAGCGGCCGCCTGACCACGCCCGACCCGATTGCCGTGCAATCGACCCTCAACGACCTCAAGAAAGGCGGCGCCAGGGCCGTGGCCATGGAAGTGTCTTCCCATGC‑‑CCTGGAGCAAGGCCGGGTTGCCg                                               <  3:312949/149‑1 (MQ=255)
                                                                                          ggTTTCTACGGCGAGTTGCAGAGCGGCCGCCTGACCACGCCCGACCCGATTGCCGTGCAATCGACCCTCAACGACCTCAAGAAAGGCGGCGCCAGGGCCGTGGCCATGGAAGTGTCTTCCCATGC‑‑CCTGGAGCAAGGCCGGGTTGCCGc                                              <  4:25452/149‑1 (MQ=255)
                                                                                           gTTTCTACGGCGAGTTGCAGAGCGGCCGCCTGACCACGCCCGACCCGATTGCCGTGCAATCGACCCTCAACGACCTCAAGAAAGGCGGCGCCAGGGCCGTGGCCATGGAAGTGTCt                                                                                >  6:214658/1‑116 (MQ=255)
                                                                                           gTTTCTACGGCGAGTTGCAGAGCGGCCGCCTGACCACGCCCGACCCGATTGCCGTGCAATCGACCCTCAACGACCTCAAGAAAGGCGGCGCCAGGGCCGTGGCCATGGAAGTGTCt                                                                                <  5:214658/116‑1 (MQ=255)
                                                                                                tACGGCGAGTTGCAGAGCGGCCGCCTGACCACGCCCGACCCGATTGCCGTGCAATCGACCCTCAACGACCTCAAGAAAGGCGGCGCCAGGGCCGTGGCCATGGAAGTGTCTTCCCATGC‑‑CCTGGAGCAAg                                                           >  1:8429/1‑130 (MQ=255)
                                                                                                tACGGCGAGTTGCAGAGCGGCCGCCTGACCACGCCCGACCCGATTGCCGTGCAATCGACCCTCAACGACCTCAAGAAAGGCGGCGCCAGGGCCGTGGCCATGGAAGTGTCTTCCCATGC‑‑CCTGGAGCAAg                                                           <  2:8429/130‑1 (MQ=255)
                                                                                                  cGGCGAGTTGCAGAGCGGCCGCCTGACCACGCCCGACCCGATTGCCGTGCAATCGACCCTCAACGACCTCAAGAAAGGCGGCGCCAGGGCCGTGGCCATGGAAGTGTCTTCCCATGC‑‑CCTGGAGCAAGGCCGGGTTGCCGCACTGGAAt                                      <  8:130402/149‑1 (MQ=255)
                                                                                                      gAGTTGCAGAGCGGCCGCCTGACCACGCCCGACCCGATTGCCGTGCAATCGACCCTCAACGACCTCAAGAAAGGCGGCGCCAGGGCCGTGGCCATGGAAGTGTCTTCCCATGC‑‑CCTGGAGCAAGGCCGGGTTGCCGCACTGGAATTCGa                                  <  7:160188/149‑1 (MQ=255)
                                                                                                         ttGCAGAGCGGCCGCCTGACCACGCCCGACCCGATTGCCGTGCAATCGACCCTCAACGACCTCAAGAAAGGCGGCGCCAGGGCCGTCGCCATGGAAGAGTCTTCCAATGC‑‑CCTGGAGAAAGGCCGGGTTGACGCACTGTAATTCGACAt                               >  5:237599/1‑149 (MQ=255)
                                                                                                             cgagCGGCCGCCTGACCACGCCCGCCCCGATTGCCGTGCAATCGACCCTCAACGACCTCCAGAAAGGCGGCGCCAGGGCCGTGGCCATGGAAGTGGCTTCCCATGG‑‑CCTGGAGCAAGGCCGGGTTGCCGCACTGGAATTCGACATcgcg                           <  8:126635/148‑1 (MQ=255)
                                                                                                             agagCGGCCGCCTGACCACGCCCGACCCGATTGCCGTGCAATCGACCCTCAACGACCTCAAGAAAGGCGGCGCCAGGGCCGTGGCCATGGAAGTGTCTTCCCATTC‑‑CCTGGAGCAACGCCGGGTTGCCGCAATGGAATTCGACATcgcg                           >  4:19934/1‑149 (MQ=255)
                                                                                                                  ggCCGCCTGACCACGCCCGACCCGATTGCCGTGCAATCGACCCTCAACGACCTCAAGAAAGGCGGCGCCAGGGCCGTGGCCATGGAAGTGTCTTCCCATGC‑‑CCTGGAGCAAGGCCGGGTTGCCGCACTGGAATTCGACATCGCGGtgat                      <  7:183104/149‑1 (MQ=255)
                                                                                                                  ggCCGCCTGACCACGCCCGACCCGATTGCCGTGCAATCGACCCTCAACGACCTCAAGAAAGGCGGCGCCAGGGCCGTGGCCATGGAAGTGTCTTCCCATGC‑‑CCTGGAGCAAGGCCGGGTTGCCGCACTGGAATTCGACATCGCGGtgat                      <  1:13136/149‑1 (MQ=255)
                                                                                                                      gccTGACCACGCCCGACCCGATTGCCGTGCAATCGACCCTCAACGACCTCAAGAAAGGCGGCGCCAGGGCCGTGGCCATGGAAGTGTCTTCCCATGC‑‑CCTGGAGCAAGGCCGGGTTGCCGCACTGGAATTCGACATCGCGGTGATGAc                   >  2:82043/1‑148 (MQ=255)
                                                                                                                      gccTGACCACGCCCGACCCAATTGCCGTGCAATCGACCCTCAACGACCTCAAGAAAGGCGGCGCCAGGGCCGTGGCCATGGAAGTGTCTTCCCATGC‑‑CCTGGAGCAAGGCCGGGTTGc                                                 >  8:221528/1‑118 (MQ=255)
                                                                                                                      gccTGACCACGCCCGACCCAATTGCCGTGCAATCGACCCTCAACGACCTCAAGAAAGGCGGCGCCAGGGCCGTGGCCATGGAAGTGTCTTCCCATGC‑‑CCTGGAGCAAGGCCGGGTTGc                                                 <  7:221528/118‑1 (MQ=255)
                                                                                                                                     ggtcgaTTGCACGGCAATCGACCCTCAACGACCTCAAGAAAGGCGGCGCCAGGGCCGTGGCCATGGAAGTGTCTTCCCATGC‑‑CCTGGAGCAAGGCCGGGTTGCCGCACTGGAATTCGACATCGCGGTGATGACCAACCTGTCCCGCGac   >  6:278087/4‑149 (MQ=255)
                                                                                                                                      gtcgaTTGCACGGCAATCGACCCTCAACGACCTCAAGAAAGGCGGCGCCAGGGCCGTGTCCATCGAAGTGTCTTACCATGC‑‑CCAGGAGCAAGGCCCGGTTCACGCACTTGAAATCGACATCGAGTAGATGACCAACCTGTCCCGCGacc  >  3:281060/3‑149 (MQ=255)
                                                                                                                                       ccgaTTGCCGTGCAAACGACCCACAAACACCTCAATAAATTGGGCGCCAGGGGCGTGGGCATTGTAGTGTTTTCCCATCCAACCTGGAGCAAGGCCGGGTTGCCGCACTGGAATTCGACATCGCGGTGATGACCAACCTGTCCCGCGac   <  1:111759/149‑1 (MQ=255)
                                                                                                                                              ccGTGCAATCGACCCTCAACGACCTCAAGAAAGGCGGCGCCAGGGCCGTGGCCATGGAAGTGTCTTCCCATGC‑‑CCTGGAGCAAGGCCGGGTTGCCGCACTGGAATTCGACATCGCGGTGATGAcc                  <  2:242398/125‑1 (MQ=255)
                                                                                                                                              ccGTGCAATCGACCCTCAACGACCTCAAGAAAGGCGGCGCCAGGGCCGTGGCCATGGAAGTGTCTTCCCATGC‑‑CCTGGAGCAAGGCCGGGTTGAAGCACTGGAATTCGACATCGCGGTGATGAcc                  >  1:242398/1‑125 (MQ=255)
                                                                                                                                                 |                                                                                                                                           
ACCAACGGCAAGACAAGTGTTACCCAACTGGTCGCCCAGGCGCTCGACCTGCTTGGCCAGCGCTGCGGCCTGATTGGCACCCTGGGTACCGGTTTCTACGGCGAGTTGCAGAGCGGCCGCCTGACCACGCCCGACCCGATTGCCGTGCAATCGACCCTCAACGACCTCAAGAAAGGCGGCGCCAGGGCCGTGGCCATGGAAGTGTCTTCCCATGC‑‑CCTGGAGCAAGGCCGGGTTGCCGCACTGGAATTCGACATCGCGGTGATGACCAACCTGTCCCGCGACC  >  NC_002947/1517138‑1517420

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