Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,521,290 C→A 11.7% Q169K (CAG→AAG)  murD → UDP‑N‑acetylmuramoylalanine‑‑D‑glutamate ligase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029471,521,2900CA11.7% 170.7 / 7.3 59Q169K (CAG→AAG) murDUDP‑N‑acetylmuramoylalanine‑‑D‑glutamate ligase
Reads supporting (aligned to +/- strand):  ref base C (27/25);  new base A (4/3);  total (31/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 = 9.95e-01

GCCAAGAGCACCGTCACGACGCTGGTCGGTGAAATGGCCGCCAAGGCCGGCAAACGCGTGGCCGTGGGCGGTAACCTCGGTACCCCGGCACTCGACCTGCTGGCCGATGACATCGAGCTGTACGTGCTGGAGCTTTCCAGCTTCCAGCTGGAAACCACTGACCAGCTCAACGCCGAAGTGGCGACCGTTCTGAACATCAGCGAAGACCACATGGACCGCTATAGCGGCCTTCCGGCCTACCACCTGGCCAAACACCGGATTTTCCGCGGTGCCCGCCAGGTGGTGGTCAACC  >  NC_002947/1521146‑1521437
                                                                                                                                                |                                                                                                                                                   
gCCAAGAGCACCGTCACGACGCTGGTCGGTGAAATGGCCGCCAAGGCCGGCAAACGCGTGGCCGTGGGCGGTAACCTCGGTACCCCGGCACTCGACCTGCTGGCCGATGACATCGAGCTGTACGTGCTGGAGCTTTCCCGCTTCCAGCTg                                                                                                                                                >  1:488044/1‑150 (MQ=255)
 ccAAGAGCACCGTCACGACGCTGGTCGGTGAAATGGCCGCCAAGGCCGGCAAACGCGTGGCCGTGGGCGGTAACCTCGGTACCCCGGCACTCGACCTGCTGGCCGATGACATCGAGCTGTACGTGCTGGAGCTTTCCAGGTTCCAGCTgg                                                                                                                                               <  2:452350/150‑1 (MQ=255)
  cAAGAGCACCGTCACGACGCTGGTCGGTGAAATGGCCGCCAAGGCCGGCAAACGCGTGGCCGTGGGCGGTAACCTCGGTACCCCGGCACTCGACCTGCTGGCCGATGACATCGAGCTGTACGTGCTGGAGCTTTCCAGCTTCCAGCTGGa                                                                                                                                              >  2:163658/1‑150 (MQ=255)
       gCACCGTCACGACGCTGGTCGGTGAAATGGCCGCCAAGGCCGGCAAACGCGTGGCCGTGGGCGGTAACCTCGGTACCCCGGCACTCGACCTGCTGGCCGATGACATCGAGCTGTACGTGCTGGAGCTTTCCAGCTTCCAGCTGGAAACCa                                                                                                                                         <  2:452188/150‑1 (MQ=255)
         aCCGTCACGACGCTGGTCGGTGAAATGGCCGCCAAGGCCGGCAAACGCGTGGCCGTGGGCGGTAACCTCGGTACCCCGGCACTCGACCTGCTGGCCGATGACATCGAGCTGTACGTGCTGGAGCTTTCCAGCTTCCAGCTGGAAACCACt                                                                                                                                       <  1:109984/150‑1 (MQ=255)
         aCCGTCACGACGCTGGTCGGTGAAATGGCCGCCAAGGCCGGCAAACGCGTGGCCGTGGGCGGTAACCTCGGTACCCCGGCACTCGACCTGCTGGCCGATGACATCGAGCTGTACGTGCTGGAGCTTTCCAGCTTCCAGCTGGAAACCACt                                                                                                                                       <  2:509467/150‑1 (MQ=255)
           cGTCACGACGCTGGTCGGTGAAATGGCCGCCAAGGCCGGCAAACGCGTGGCCGTGGGCGGTAACCTCGGTACCCCGGCACTCGACCTGCTGGCCGATGACGTCGAGCTGTACGTGCTGGAGCTTTCCAGCTTCCAGCTGGAAACCACTGa                                                                                                                                     >  1:1245297/1‑150 (MQ=255)
           cGTCACGACGCTGGTCGGTGAAATGGCCGCCAAGGCCGGCAAACGCGTGGCCGTGGGCGGTAACCTCGGTACCCCGGCACTCGACCTGCTGGCCGATGACATCGAGCTGTACGTGCTGGAGCTTTCCAGCTTCCAGCTGGAAACCACTGa                                                                                                                                     >  2:1363956/1‑150 (MQ=255)
                cgacgCTGGTCGGTGAAATGGCCGCCAAGGCCGGCAAACGCGTGGCCGTGGGCGGTAACCTCGGTACCCCGGCACTCGACCTGCTGGCCGATGACATCGAGCTGTACGTGCTGGAGCTTTCCAGCTTCCAGCTGGAAACCACTGACCAGc                                                                                                                                <  1:250161/150‑1 (MQ=255)
                 gacgCTGGTCGGTGAAATGGCCGCCAAGGCCGGCAAACGCGTGGCCGTGGGCGGTAACCTCGGTACCCCGGCACTCGACCTGCTGGCCGATGACATCGAGCTGTACGTGCTGGAGCTTTCCAGCTGGAAGCTGGaaa                                                                                                                                            <  1:754295/137‑1 (MQ=255)
                 gacgCTGGTCGGTGAAATGGCCGCCAAGGCCGGCAAACGCGTGGCCGTGGGCGGTAACCTCGGTACCCCGGCACTCGACCTGCTGGCCGATGACATCGAGCTGTACGTGCTGGAGCTTTCCAGCTGGAAGCTGGaaa                                                                                                                                            >  2:754295/1‑137 (MQ=255)
                    gCTGGTCGGTGAAATGGCCGCCAAGGCCGGCAAACGCGTGGCCGTGGGCGGTAACCTCGGTACCCCGGCACTCGACCTGCTGGCCGATGACATCGAGCTGTACGTGCTGGAGCTTTCCAGCTTCCAGCTGGAAACCACTGACCAGCTCaa                                                                                                                            >  1:1128707/1‑150 (MQ=255)
                    gCTGGTCGGTGAAATGGCCGCCAAGGCCGGCAAACGCGTGGCCGTGGGCGGTAACCTCGGTACCCCGGCACTCGACCTGCTGGCCGATGACATCGAGCTGTACGTGCTGGAGCTTTCCAGCTTCCAGCTGGAAACCACTGACCAGCTCaa                                                                                                                            <  1:632812/150‑1 (MQ=255)
                      tGGTCGGTGAAATGGCCGCCAAGGCCGGCAAACGCGTGGCCGTGGGCGGTAACCTCGGTACCCCGGCACTCGACCTGCTGGCCGATGACATCGAGCTGTACGTGCTGGAGCTTTCCAGCTTCCAGCTGGAAACCACTGACCAGCTc                                                                                                                              <  1:802574/146‑1 (MQ=255)
                      tGGTCGGTGAAATGGCCGCCAAGGCCGGCAAACGCGTGGCCGTGGGCGGTAACCTCGGTACCCCGGCACTCGACCTGCTGGCCGATGACATCGAGCTGTACGTGCTGGAGCTTTCCAGCTTCCAGCTGGAAACCACTGACCAGCTc                                                                                                                              >  2:802574/1‑146 (MQ=255)
                                aaTGGCCGCCAAGGCCGGCAAACGCGTGGCCGTGGGCGGTAACCTCGGTACCCCGGCACTCGACCTGCTGGCCGATGACATCGAGCTGTACGTGCTGGAGCTTTCCAGCTTCCAGCTGGAAACCACTGACCAGCTCAACGCCGAAGTgg                                                                                                                 >  2:863097/1‑149 (MQ=255)
                                aaTGGCCGCCAAGGCCGGCAAACGCGTGGCCGTGGGCGGTAACCTCGGTACCCCGGCACTCGACCTGCTGGCCGATGACATCGAGCTGTACGTGCTGGAGCTTTCCAGCTTCCAGCTGGAAACCACTGACCAGCTCAACGCCGAAGTgg                                                                                                                 <  1:863097/149‑1 (MQ=255)
                                aaTGGCCGCCAAGGCCGGCAAACGCGTGGCCGTGGGCGGTAACCTCGGTACCCCGGCACTCGACCTGCTGGCCGATGACATCGAGCTGTACGTGCTGGAGCTTTCCAGCTTCCAGCTGGAAACCACTGACCAGCTCAACGCCGAAGTGGc                                                                                                                <  1:106758/150‑1 (MQ=255)
                                       gccAAGGCCGGCAAACGCGTGGCCGTGGGCGGTAACCTCGGTACCCCGGCACTCGACCTGCTGGCCGATGACATCGAGCTGTACGTGCTGGAGCTTTCCAGCTGGAAGCTGGaaa                                                                                                                                            <  2:905648/115‑1 (MQ=255)
                                       gccAAGGCCGGCAAACGCGTGGCCGTGGGCGGTAACCTCGGTACCCCGGCACTCGACCTGCTGGCCGATGACATCGAGCTGTACGTGCTGGAGCTTTCCAGCTGGAAGCTGGaaa                                                                                                                                            >  1:905648/1‑115 (MQ=255)
                                         cAAGGCCGGCAAACGCGTGGCCGTGGGCGGTAACCTCGGTACCCCGGCACTCGACCTGCTGGCCGATGACATCGCGATGTACGTGCTGGAGCCTACCACGCTCACCCAATGTTCCTGGGATCAGCTCAACGCCGAAGTGGCGACCGTTCt                                                                                                       <  1:1315976/150‑1 (MQ=255)
                                         cAAGGCCGGCAAACGCGTGGCCGTGGGCGGTAACCTCGGTACCCCGGCACTCGACCTGCTGGCCGATGACATCGAGCTGTACGTGCTGGAGCTTTCCAGCTTCCAGCTGGAAACCACTGACCAGCTCAACGCCGAAGTGGCGACCGTTCt                                                                                                       <  1:687842/150‑1 (MQ=255)
                                                 gCAAACGCGTGGCCGTGGGCGGTAACCTCGGTACCCCGGCACTCGACCTGCTGGCCGATGACATCGAGCTGTACGTGCTGGAGCTTTCCAGCTTCCAGCTGGAAACCACTGACCAGCTCAACGCCGAAGTGGCGACCGTTCTGAACATCa                                                                                               >  2:721083/1‑150 (MQ=255)
                                                   aaaCGCGTGGCCGTGGGCGGTAACCTCGGTACCCCGGCACTCGACCTGCTGGCCGATGACATCGAGCTGTACGTGCTGGAGCTTTCCAGCTTCCAGCTGGAAACCACTGACCAGCTCAACGCCGAAGTGGCGACCGTTCTGAACATCAGc                                                                                             <  2:1245297/150‑1 (MQ=255)
                                                      cgcgTGGCCGTGGGCGGTAACCTCGGTACCCCGGCACTCGACCTGCTGGCCGATGACATCGAGCTGTACGTGCTGGAGCTTTCCAGCTTCCAGCTGGAAACCACTGACCAGCTCAACGCCGAAGTGGCGACCGTTCTGAACATCAGCGaa                                                                                          <  1:815868/150‑1 (MQ=255)
                                                           ggCCGTGGGCGGTAACCTCGGTACCCCGGCACTCGACCTGCTGGCCGATGACATCGAGCTGTACGTGCTGGAGCTTTCCAGCTTCCAGATGGAAACCACTGCCCCGCTCAACGCCGAAGTGGCGACCGTTCTGAACATCAGCGAAGACca                                                                                     <  2:488044/150‑1 (MQ=255)
                                                                           cTCGGTACCCCGGCACTCGACCTGCTGGCCGATGACATCGAGCTGTACGTGCTGGAGCTTTCCAGCTTCCAGCTGGAAACCACTGACCAGCTCAACGCCGAAGTGGCGACCGTTCTGAACATCAGCGaa                                                                                          <  2:817410/129‑1 (MQ=255)
                                                                           cTCGGTACCCCGGCACTCGACCTGCTGGCCGATGACATCGAGCTGTACGTGCTGGAGCTTTCCAGCTTCCAGCTGGAAACCACTGACCAGCTCAACGCCGAAGTGGCGACCGTTCTGAACATCAGCGaa                                                                                          >  1:817410/1‑129 (MQ=255)
                                                                           cTCGGTACCCCGGCACTCGACCTGCTGGCCGATGACATCGAGCTGTACGTGCTGGAGCTTTCCAGCTTCCAGCTGGAAACCACTGACCAGCTCAACGCCGAAGTGGCGACCGTTCTGAACATCAGCGAAGACCACATGGACCGCTATAGc                                                                     >  1:1290980/1‑150 (MQ=255)
                                                                            tCGGTACCCCGGCACTCGACCTGCTGGCCGATGACATCGAGCTGTACGTGCTGGAGCTTTCCAGCTTCCAGCTGGAAACCACTGACCAGCTCAACGCCGaa                                                                                                                     <  1:1176304/101‑1 (MQ=255)
                                                                            tCGGTACCCCGGCACTCGACCTGCTGGCCGATGACATCGAGCTGTACGTGCTGGAGCTTTCCAGCTTCCAGCTGGAAACCACTGACCAGCTCAACGCCGaa                                                                                                                     >  2:1176304/1‑101 (MQ=255)
                                                                            tCGGTACCCCGGCACTCGACCTGCTGGCCGATGACATCGAGCTGTACGTGCTGGAGCTTTCCAGCTTCCAGCTGGAAACCACTGACCAGCTCAACGCCGAAGTGGCGACCGTTCTGAACATCAGCGAAGACCACATGGACCGCTATAGCg                                                                    <  2:1277714/150‑1 (MQ=255)
                                                                              ggTACCCCGGCACTCGACCTGCTGGCCGATGACATCGAGCTGTACGTGCTGGAGCTTTCCAGCTTCCAGCTGGAAACCACTGACCAGCTCAACGCCGAAGTGGCGACCGTTCTGAACATCAGCGAAGACCACATGGACCGCTATAGCGGc                                                                  <  1:1363956/150‑1 (MQ=255)
                                                                                    ccGGCACTCGACCTGCTGGCCGATGACATCGAGCTGTACGTGCTGGAGCTTTCCAGCTTCCAGCTGGAAACCACTGACCAGCTCAACGCCGAAGTGGCGACCGTTCTGAACATCAGCGAAGACCACATGGACCGCTATAGCGGCCTTCCg                                                            >  2:319504/1‑150 (MQ=255)
                                                                                             gACCTACTGGCCGATGACATCGAGCTGTACGTGCTGGAGCTTTCCAGCTTCCAGCTGGAAACCACTGACCAGCTCAACGCCGAAGTGGCGACCGTTCTGAACATCAGCGAAGACCACATGGACCGCTATAGCGGCCTTCCGGCCTaccac                                                   <  2:1186115/150‑1 (MQ=255)
                                                                                                ctgctgGCCGATGACATCGAGCTGTACGTGCTGGAGCTTTCCAGCTTCCAGCTGGAAACCa                                                                                                                                         <  2:448739/61‑1 (MQ=255)
                                                                                                ctgctgGCCGATGACATCGAGCTGTACGTGCTGGAGCTTTCCAGCTTCCAGCTGGAAACCa                                                                                                                                         >  1:448739/1‑61 (MQ=255)
                                                                                                         gATGACATCGAGCTGTACGTGCTGGAGCTTTCCAGCTTCCAGCTGGAAACCACTGACCAGCTCAACGCCGAAGTGGCGACCGTTCTGAACATCAGCGAAGACCACATGGACCGCTATAGCGGCCTTCCGGCCTACCACCTGGCCAAacac                                       >  2:1113577/1‑150 (MQ=255)
                                                                                                          aTGACATCGAGCTGTACGTGCTGGAGCTTTCCAGCTTCCAGCTGGAAACCACTGACCAGCTCAACGCCGAAGTGGCGACCGTTCTGAACATCAGCGAAGACCACATGGACCGCTATAGCGGCCTTCCGGCCTACCACCTGGCCAAACAcc                                      >  1:159699/1‑150 (MQ=255)
                                                                                                             aCATCGAGCTGTACGTGCTGGAGCTTTCCAGCTTCCAGCTGGAAACCACTGACCAGCTCAACGCCGAAGTGGCGACCGTTCTGAACATCAGCGAAGACCACATGGACCGCTATAGCGGCCTTCCGGCCTACCACCTGGCCAAACACCGGa                                   <  1:721083/150‑1 (MQ=255)
                                                                                                             aCATCGAGCTGTACGTGCTGGAGCTTTCCAGCTTCCAGCTGGAAACCACTGACCAGCTCAACGCCGAAGTGGCGACCGTTCTGAACATCAGCGAAGACCACATGGACCGCTATAGCGGCCTTCCGGCCTACCACCTGGCCAAACACCGGa                                   <  2:1128707/150‑1 (MQ=255)
                                                                                                               aTCGAGCTGTACGTGCTGGAGCTTTCCAGCTTCCAGCTGGAAACCACTGACCAGCTCAACGCCGAAGTGGCGACCGTTCTGAACATCAGCGaa                                                                                          <  1:1345606/93‑1 (MQ=255)
                                                                                                               aTCGAGCTGTACGTGCTGGAGCTTTCCAGCTTCCAGCTGGAAACCACTGACCAGCTCAACGCCGAAGTGGCGACCGTTCTGAACATCAGCGaa                                                                                          >  2:1345606/1‑93 (MQ=255)
                                                                                                               aTCGAGCTGTACGTGCTGGAGCTTTCCAGCTTCCAGCTGGAAACCACTGACCAGCTCAACGCCGAAGTGGCGACCGTTCTGAACATCAGCGAAGACCACATGGACCGCTATAGCGGCCTTCCGGCCTACCACCTGGCCAAACACCGGAtt                                 >  1:1317590/1‑150 (MQ=255)
                                                                                                               aTCGAGCTGTACGTGCTGGAGCTTTCCAGCTTCCAGCTGGAAACCACTGACCAGCTCAACGCCGAAGTGGCGACCGTTCTGAACATCAGCGAAGACCACATGGACCGCTATAGCGGCCTTCCGGCCTACCACCTGGCCAAACACCGGAtt                                 >  2:948258/1‑150 (MQ=255)
                                                                                                               aTCGAGCTGTACGTGCTGGAGCTTTCCAGCTTCCAGCTGGAAACCACTGACCAGCTCAACGCCGAAGTGGCGACCGTTCTGAACATCAGCGAAGACCACATGGACCGCTATAGCGGCCTTCCGGCCTACCACCTGGCCAAACACCGGAtt                                 >  2:1207672/1‑150 (MQ=255)
                                                                                                                tCGAGCTGTACGTGCTGGAGCTTTCCAGCTTCCAGCTGGAAACCACTGACCAGCTCAACGCCGAAGTGGCGACCGTTCTGAACATCAGCGAAGACCACATGGACCGCTATAGCGGCCTTCCGGCCTACCACCTGGCCAAACACCGGAttt                                >  1:614787/1‑150 (MQ=255)
                                                                                                                  gAGCTGTACGTGCTGGAGCTTTCCAGCTTCCAGCTGGAAACCACTGACCAGCTCAACGCCGAAGTGGCGACCGTTTTGAACATCAGCGAAGACCACATGGACCGCTATAGCGGCCTTCCGGCCTACCACCTGGCCAAACACCGGATTTTc                              >  1:804970/1‑150 (MQ=255)
                                                                                                                        taggtggtgGCGCCTTGCAGCTTCCAGCGGTAGACCACTTACCAGCTCAACGACGAAGTGGCGACCGTCATGACCATCAGCGAAGACCACATGGACCGCTATAGCGGCCTTCCGGCCTACCACCTGGCCAAACACCGGATTTTCCGCGGt                        <  2:275092/147‑1 (MQ=255)
                                                                                                                        tACGTGCTGGAGCTTTCCAGCTTCCAGCTGGAAACCACTGACCAGCTCAACGCCGAAGTGGCGACCGTTCTGAACATCAGCGAAGACCACATGGACCGCTATAGCGGCCTTCCGGCCTACCACATGACCAAACACCGAATTTTCTGCGGt                        >  2:1150300/1‑150 (MQ=255)
                                                                                                                           gTGCTGGAGCTTTCCAGCTTCCAGCTGGAAACCACTGACCAGCTCAACGCCGAAGTGGCGACCGTTCTGAACATCAGCGAAGACCACATGGACCGCTATAGCGGCCTTCCGGCCTACCACCTGGCCAAACACCGGATTTTCCGCGGTGcc                     >  1:1333516/1‑150 (MQ=255)
                                                                                                                            tGCTGGAGCTTTCCAGCTTCCAGCTGGAAACCACTGACCAGCTCAACGCCGAAGTGGCGACCGTTCTGAACATCAGCGAAGACCACATGGACCGCTATAGCGGCCTTCCGGCCTACCACCTGGCCAAACACCGGATTTTCCGCGGTGccc                    <  2:804970/150‑1 (MQ=255)
                                                                                                                               tGGAGCTTTCCAGCTTCCAGCTGGAAACCACTGACCAGCGCAACGCCGAAGTGCCGACCGTTCTGCACATCAGCGAAGACCACATGGACCGCTATACCAGCCACCCGGCCTACCACCTGGCCAAACACCGGATTTTCCGCGGTGCCCGcc                 >  1:447312/1‑150 (MQ=255)
                                                                                                                                   gCTTTCCAGCTTCCAGCTGGAAACCACTGACCAGCTCAACGCCGAAGTGGCGACCGTTCTGAACATCAGCGAAGACCACATGGACCGCTATAGCGGCCTTCCGGCCTACCACCTGGCCAAACACCGGATTTTCCGCGGTGCCCGCCAggt             >  2:585629/1‑150 (MQ=255)
                                                                                                                                    gtttCCAGCTGGAAGCTGGAAACCACTGACCAGCTCAACGCCGAAGTGGCGACCGTTCTGAACATCAGCGAAGACCACATGGACCGCTATAGCGGCCTTCCGGCCTACCACCTGGCCAAACACCGGATTTTCCGCGGTGCCCGCCAggtg            >  2:1197114/2‑150 (MQ=255)
                                                                                                                                    cTTTCCAGCTTCCAGCTGGAAACCACTGACCAGCTCAACGCCGAAGTGGCGACCGTTCTGAACATCAGCGAAGACCACATGGACCGCTATAGCGGCCTTCCGGCCTACCCCCTGGCCAAACACCGGATTTTCCGCGGTGCCCGCCAggtg            >  1:1009146/1‑150 (MQ=255)
                                                                                                                                    cTTTCCAGCTTCCAGCTGGAAACCACTGACCAGCTCAACGCCGAAGTGGCGACCGTTCTGAACATCAGCGAAGACCACATGGACCGCTATAGCGGCCTTCCGGCCTACCACCTGGCCAAACACCGGATTTTCCGCGGTGCCCGCCAggtg            <  2:1067913/150‑1 (MQ=255)
                                                                                                                                    cTTTCCAGCTTCCAGCTGGAAACCACTGACCAGCTCAACGCCGAAGTGGCGACCGTTCTGAACATCAGCGAAGACCACATGGACCGCTATAGCGGCCTTCCGGCCTACCACCTGGCCAAACACCGGATTTTCCGCGGTGCCCGCCAggtg            >  1:932261/1‑150 (MQ=255)
                                                                                                                                     tttCCAGCTGGAAGCTGGAAACCACTGACCAGCTCAACGCCGAAGTGGCGACCGTTCTGAACATCAGCGAAGACCACATGGACCGCTATAGAGGCCTTCCGCCCTACCACCTGGCCAAACACCGGATTTTCCGCGGTGCCCGCCAggtgg           >  2:403214/1‑150 (MQ=255)
                                                                                                                                      ttccagctggaagcTGGAAACCACTGACCAGCTCAACGCCGAAGTGGCGACCGTTCTGAACACCAGCGAAGACCACATGGACCGCTATAGCGGCCTTCCGGCCTACCACCTGGCCAAACACCGGCTTTTCCGCGGTGCCCGCCAggtggt          >  1:528602/1‑150 (MQ=255)
                                                                                                                                       tccagcttccagctGGAAACCACTGACCAGCTCAACGCCGAAGTGGCGACCGTTCTGAACATCAGCGAAGACCACATGGACCGCTATAGCGGCCTTCCGGCCTACCACCTGGCCAAACACCGGATTTTCCGCGGTGCCCGCCAggtggtg         >  2:629389/1‑150 (MQ=255)
                                                                                                                                       tccagcttccagctGGAAACCACTGACCAGCTCAACGCCGAAGTGGCGACCGTTCTGAACATCAGCGAAGACCACATGGACCGCTATAGCGGCCTTCCGGCCTACCACCTGGCCAAACACCGGATTTTCCGCGGTGCCCGCCAggtggtg         <  2:614787/150‑1 (MQ=255)
                                                                                                                                          agcttccagctGGAAACCACTGACCAGCTCAACGCCGAAGTGGCGACCGTTCTGAACATCAGCGAAGACCACATGGACCGCTATAGCGGCCTTCCGGCCTACCACCTGGCCAAACACCGGATTTTCCGCGGTGCCCGCCAGGTGGTGGTc      <  2:159699/150‑1 (MQ=255)
                                                                                                                                           gcttccagctGGAAACCACTGACCATCTCAACGCCGAAGTGGCGACCGTTCTGAACATCAGCGAAGACCACATGGACCGCTATAGCGGCCTTCCGGCCTACCACCTGGCCAAACACCGGATTTTCCGCGGTGCCCGCCAggtggtgggaa     <  1:49359/150‑4 (MQ=255)
                                                                                                                                              tccagctGGAAACCACTGACCAGCTCAACGCCGAAGTGGCGACCGTTCTGAACATCAGCGAAGACCACATGGACCGCTATAGCGGCCTTCCGGCCTACCACCTGGCCAAACACCGGATTTTCCGCGGTGCCCGCCAGGTGGTGGTCAAcc  <  1:1018215/150‑1 (MQ=255)
                                                                                                                                                |                                                                                                                                                   
GCCAAGAGCACCGTCACGACGCTGGTCGGTGAAATGGCCGCCAAGGCCGGCAAACGCGTGGCCGTGGGCGGTAACCTCGGTACCCCGGCACTCGACCTGCTGGCCGATGACATCGAGCTGTACGTGCTGGAGCTTTCCAGCTTCCAGCTGGAAACCACTGACCAGCTCAACGCCGAAGTGGCGACCGTTCTGAACATCAGCGAAGACCACATGGACCGCTATAGCGGCCTTCCGGCCTACCACCTGGCCAAACACCGGATTTTCCGCGGTGCCCGCCAGGTGGTGGTCAACC  >  NC_002947/1521146‑1521437

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑