Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,526,631 A→C 9.7% E259A (GAA→GCA)  ddlB → D‑alanine‑‑D‑alanine ligase B

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029471,526,6310AC 9.7% 170.1 / 8.7 62E259A (GAA→GCA) ddlBD‑alanine‑‑D‑alanine ligase B
Reads supporting (aligned to +/- strand):  ref base A (24/32);  new base C (3/3);  total (27/35)
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.86e-01

CCTGGGTACGCCGCACGTGTTCTACGACTATGACGCCAAGTACATCGTCAACGATACCCAGTACCGCATCCCGTGCGGGCTGGATAGCGTCAAGGAACAGGAACTGATCGACCTGACCGCGCGC‑GCCTGCGATGCCATTGGCATCGAAGGCTGGGGGCGGTTGGACGTGATGCAGGACGAGCAGGGCCGGTTCTGGCTGCTCGAAGTCAACACCGCACCAGGCATGACAGATCATAGCCTGGTGCCCATGGCGGCCCGTGCGGCTGGCCTGGACTTCCAGCAACTGGTG  >  NC_002947/1526485‑1526773
                                                                                                                                                   |                                                                                                                                              
ccTGGGTACGCCGCACGTGTTCTACGACTATGACGCCAAGTACATCGTCAACGATACCCAGTACCGCATCCCGTGCGGGCTGGATAGCGTCAAGGAACAGGAACTGATCAACCTGACCGCGCCCAGCCTTCGATGCCAATGGCATCGCAg                                                                                                                                              <  2:118290/150‑1 (MQ=255)
         gccgcACGTGTTCTACGACTATGACGCCAAGTACATCGTCAACGATACCCAGTACCGCATCCCGTGCGGGCTGGATAGCGTCAAGGAACAGGAACTGATCGACCTGACCGCGCGC‑GCCTGCGATGCCATTGGCATCGAAGGCTGGGGGCg                                                                                                                                    >  2:1397231/1‑150 (MQ=255)
         gccgcACGTGTTCTACGACTATGACGCCAAGTACATCGTCAACGATACCCAGTACCGCATCCCGTGCGGGCTGGATAGCGTCAAGGAACAGGAACTGATCGACCTGACCGCGCGC‑GCCTGCGATGCCATTGGCATCGAAGGCTGGGGGCg                                                                                                                                    >  1:702590/1‑150 (MQ=255)
          ccgcACGTGTTCTACGACTATGACGCCAAGTACATCGTCAACGATACCCAGTACCGCATCCCGTGCGGGCTGGATAGCGTCAAGGAACAGGAACTGAGCGACCTGACCGCGCGC‑GCCTGCGATGCCATTGGCATCGAAGGCTGGGGGCgg                                                                                                                                   <  2:702590/150‑1 (MQ=255)
            gcACGTGTTCTACGACTATGACGCCAAGTACATCGTCAACGATACCCAGTACCGCATCCCGTGCGGGCTGGATAGCGTCAAGGAACAGGAACTGATCGACCTGACCGCGCGC‑GCCTGCGATGCCATTGGCATCGAAGGCTGGGGGCGGtt                                                                                                                                 <  2:1177650/150‑1 (MQ=255)
              aCGTGTTCTACGACTATGACGCCAAGTACATCGTCAACGATACCCAGTACCGCATCCCGTGCGGGCTGGATAGCGTCAAGGAACAGGAACTGATCGACCTGACCGCGCGC‑GCCTGCGATGCCATTGGCATCGAAGGCTGGGGGCGGTTgg                                                                                                                               >  1:1201423/1‑150 (MQ=255)
                   ttCTACGACTATGACGCCAAGTACATCGTCAACGATACCCAGTACCGCATCCCGTGCGGGCTGGATAGCGTCAAGGAACAGGAACTGATCGACCTGACCGCGCGC‑GCCTGCGATGCCATTGGCATCGAAGGCTGGGGGCGGTTGGACGTg                                                                                                                          >  2:726331/1‑150 (MQ=255)
                   ttCTACGACTATGACGCCAAGTACATCGTCAACGATACCCAGTACCGCATCCCGTGCGGGCTGGATAGCGTCAAGGAACAGGAACTGATCGACCTGACCGCGCGC‑GCCTGCGATGCCATTGGCATCGAAGGCTGGGGGCGGTTGGACGTg                                                                                                                          <  2:261927/150‑1 (MQ=255)
                    tCTACGACTATGACGCCAAGTACATCGTCAACGATACCCAGTACCGCATCCCGTGCGGGCTGGATAGCGTCAAGGAACAGGAACTGATCGACCTGACCGCGCGC‑GCCTGCGATGCCATTGGCATCGAAGGCTGGGGGCGGTTGGACGTGa                                                                                                                         <  1:559646/150‑1 (MQ=255)
                    tCTACGACTATGACGCCAAGTACATCGGCAACGATACCCAGTACCGCATCCCGTGCGGGCTGGATAGCGTCAAGGAACAGGAACTGATCGACCTGACCGCGCGC‑GCCTGCGATGCCATTGGCATCGAAGGCTGGGGGCGGTTGGACGTGa                                                                                                                         <  2:1078594/150‑1 (MQ=255)
                      tACGACTATGACGCCAAGTACATCGTCAACGATACCCAGTACCGCATCCCGTGCGGGCTGGATAGCGTCAAGGAACAGGAACTGATCGACCTGACCGCGCGC‑GCCTGCGATGCCATTGGCATCGAAGGCTGGGGGCGGTTGGACGTGATg                                                                                                                       <  2:92459/150‑1 (MQ=255)
                       aCGACTATGACGCCAAGTACATCGTCAACGATACCCAGTACCGCATCCCGTGCGGGCTGGATAGCGTCAAGGAACAGGAACTGATCGACCTGACCGCGCGC‑GCCTGCGATGCCATTGGCATCGAAGGCTGGGGGCGGTTGGACGTGATGc                                                                                                                      <  2:909459/150‑1 (MQ=255)
                        cGACTATGACGCCAAGTACATCGTCAACGATACCCAGTACCGCATCCCGTGCGGGCTTTATAGCGTCAAGGAACAGGAACTGATCGACCTGACCGCGCGC‑GCCTGCGATGCCATTGGCATCGAAGGCTGGGGGCGGTTGGACGTGATGCa                                                                                                                     <  2:1404914/150‑1 (MQ=255)
                         gACTATTACGCCAAGTACATCGTCAACGATACCCAGTACCGCATCCCGTGCGGGCTGGATAGCGTCAAGGAACAGGAACTGATCGACCTGACCGCGCGC‑GCCTGCGATGCCATTGGCATCGAAGGCTGGGGGCGGTTGGACGTGATGCAg                                                                                                                    >  2:1158806/1‑150 (MQ=255)
                         gACTATGACGCCAAGTACATCGTCAACGATACCCAGTACCGCATCCCGTGCGGGCTGGATAGCGTCAAGGAACAGGAACTGATCGACCTGACCGCGCGC‑GCCTGCGATGCCATTGGCATCGAAGGCTGGGGGCGGTTGGACGTGATGCAg                                                                                                                    <  1:1397231/150‑1 (MQ=255)
                                  gCCAAGTACATCGTCAACGATACCCAGTACCGCATCCCGTGCGGGCTGGATAGCGTCAAGGAACAGGAACTGATCGACCTGACCGCGCGC‑GCCTGCGATGCCATTGGCATCGAAGGCTGGGGGCGGTTGGACGTGATGCAGGACGAGCAg                                                                                                           <  2:1201423/150‑1 (MQ=255)
                                      aGTACATCGTCAACGATACCCAGTACCGCATCCCGTGCGGGCTGGATAGCGTCAAGGAACAGGAACTGATCGACCTGACCGCGCGC‑GCCTGCGATGCCATTGGCATCGAAGGCTGGGGGCGGTTGGACGTGATGCAGGACGAGCAGGGcc                                                                                                       <  1:745148/150‑1 (MQ=255)
                                          cATCGTCAACGATACCCAGTACCGCATCCCGTGCGGGCTGGATAGCGTCAAGGAACAGGAACTGATCGACCTGACCGCGCGC‑GCCTGCGATGCCATTGGCATCGAAGGCTGGGGGCGGTTGGACGTGATGCa                                                                                                                     >  2:783534/1‑132 (MQ=255)
                                          cATCGTCAACGATACCCAGTACCGCATCCCGTGCGGGCTGGATAGCGTCAAGGAACAGGAACTGATCGACCTGACCGCGCGC‑GCCTGCGATGCCATTGGCATCGAAGGCTGGGGGCGGTTGGACGTGATGCa                                                                                                                     <  1:783534/132‑1 (MQ=255)
                                             cGTCAACGATACCCAGTACCGCATCCCGTGCGGGCTGGATAGCGTCAAGGAACAGGAACTGATCGACCTGACCGCGCGC‑GCCTGCGATGCCATTGGCATCGAAGGCTGGGGGCGGTTGGACGTGATGCAGGACGAGCAGGGCCGGTTCTg                                                                                                <  1:726331/150‑1 (MQ=255)
                                              gTCAACGATACCCAGTACCGCATCCCGTGCGGGCTGGATAGCGTCAAGGAACAGGAACTGATCGACCTGACCGCGCGC‑TCCTGCGATGCCATTGGCATCGAAGGCTGGGGGCGGTTGGACGTGATGCAGGACGAGCAGGGCCGGTTCTgg                                                                                               <  2:1405905/150‑1 (MQ=255)
                                                  aCGATACCCAGTACCGCATCCCGTGCGGGCTGGATAGCGTCAAGGAACAGGAACTGATCGACCTGACCGCGCGC‑GCCTGCGATGCCATTGGCATCGAAGGCTGGGGGCGGTTGGACGTGATGCAGGACGAGCAGGGCCGGTTCTGgctgc                                                                                           >  2:331261/1‑150 (MQ=255)
                                                    gATACCCAGTACCGCATCCCGTGCGGGCTGGATAGCGTCAAGGAACAGGAACTGATCGACCTGACCGCGCGC‑GCCTGCGATGCCATTGGCATCGAAGGCTGGGGGCGGTTGGACGTGATGCAGGACGAGCAGGGCCGGTTCTGGCTGCTc                                                                                         <  1:874277/150‑1 (MQ=255)
                                                      tACCCAGTACCGCATCCCGTGCGGGCTGGATAGCGTCAAGGAACAGGAACTGATCGACCTGACCGCGCGC‑GCCTGCGATGCCATTGGCATCGAAGGCTggg                                                                                                                                        <  1:896726/101‑1 (MQ=255)
                                                      tACCCAGTACCGCATCCCGTGCGGGCTGGATAGCGTCAAGGAACAGGAACTGATCGACCTGACCGCGCGC‑GCCTGCGATGCCATTGGCATCGAAGGCTggg                                                                                                                                        >  2:896726/1‑101 (MQ=255)
                                                       aCCCAGTACCGCATCCCGTGCGGGCTGGATAGCGTCAAGGAACAGGAACTGATCGACCTGACCGCGCGC‑GCCTGCGATGCCATTGGCATCGAAGGCTGGGGGCGGTTGGACGTGATGCAGGACGAGCAGGGCCGGTTCTGGCTGCTCGaa                                                                                      <  1:331261/150‑1 (MQ=255)
                                                       aCCCAGTACCGCATCCCGTGCGGGCTGGATAGCGTCAAGGAACAGGAACTGATCGACCTGACCGCGCGC‑GCCTGCGATGCCATTGGCATCGAAGGCTGGGGGCGGTTGGACGTGATGCAGGACGAGCAGGGCCGGTTCTGGCTGCTCGaa                                                                                      >  2:109068/1‑150 (MQ=255)
                                                        cccAGTACCGCATCCCGTGCGGGCTGGATAGCGTCAAGGAACAGGAACTGATCGACCTGACCGCGCGC‑GCCTGCGATGCCATTGGCATCGAAGGCTGGGGGCGGTTGGACGTGATGCAGGAc                                                                                                                 <  1:696394/122‑1 (MQ=255)
                                                        cccAGTACCGCATCCCGTGCGGGCTGGATAGCGTCAAGGAACAGGAACTGATCGACCTGACCGCGCGC‑GCCTGCGATGCCATTGGCATCGAAGGCTGGGGGCGGTTGGACGTGATGCAGGAc                                                                                                                 >  2:696394/1‑122 (MQ=255)
                                                          cAGTACCGCATCCCGTGCGGGCTGGATAGCGTCAAGGAACAGGAACTGATCGACCTGACCGCGCGC‑GCCTGCGATGCCATTGGCATCGAAGGCTGGGGGCGGTTGGACGTGATGCAGGACGAGCAGGGCCGGTTCTGGCTGCTCGaa                                                                                      >  2:351147/1‑147 (MQ=255)
                                                          cAGTACCGCATCCCGTGCGGGCTGGATAGCGTCAAGGAACAGGAACTGATCGACCTGACCGCGCGC‑GCCTGCGATGCCATTGGCATCGAAGGCTGGGGGCGGTTGGACGTGATGCAGGACGAGCAGGGCCGGTTCTGGCTGCTCGaa                                                                                      <  1:351147/147‑1 (MQ=255)
                                                               ccGCATCCCGTGCGGGCTGGATAGCGTCAAGGAACAGGAACTGATCGACCTGACCGCGCGC‑GCCTGCGATGCCATTGGCATCGAAGGCTGGGGGCGGTTGGACGTGATGCAGGACGAGCAGGGCCGGTTCTGGCTGCTCGAAGTCAacac                                                                              <  1:363564/150‑1 (MQ=255)
                                                                   aTCCCGTGCGGGCTGGATAGCGTCAAGGAACAGGAACTGATCGACCTGACCGCGCGC‑GCCTGCGATGCCATTGGCATCGAAGGCTGGGGGCGGTTGGACGTGATGCAGGACGAGCAGGGCCGGTTCTGGCTGCTCGAAGTCAACACCGCa                                                                          >  2:552691/1‑150 (MQ=255)
                                                                      ccGTGCGGGCTGGATAGCGTCAAGGAACAGGAACTGATCGACCTGACCGCGCGC‑GCCTGCGATGCCATTGGCATCGAAGGCTGGGGGCGGTTGGACGTGATGCAGGACGAGCAGGGCCGGTTCTGGCTGCTCGAAGTCAACACCGCACCa                                                                       >  1:1218720/1‑150 (MQ=255)
                                                                      ccGTGCGGGCTGGATAGCGTCAAGGAACAGGAACTGATCGACCTGACCGCGCGC‑GCCTGCGATGCCAATGGCATCGCAg                                                                                                                                              <  2:1347923/79‑1 (MQ=255)
                                                                          gCGGGCTGGATAGCGTCAAGGAACAGGAACTGATCGACCTGACCGCGCGC‑GCCTGCGATGCCATTGGCATCGAAGGCTGGGGGCGGTTGGACGTGATGCAGGACGAGCAGGGCCGGTTCTGGCTGCTCGAAGTCAACACCGCACCAGGCa                                                                   <  1:270157/150‑1 (MQ=255)
                                                                             ggCTGGATAGCGTCAAGGAACAGGAACTGATCGACCTGACCGCGCGC‑GCCTGCGATGCCATTGGCATCGAAGGCTGGGGGCGGTTGGACGTGATGCAGGACGAGCAggg                                                                                                         >  2:1196747/1‑109 (MQ=255)
                                                                             ggCTGGATAGCGTCAAGGAACAGGAACTGATCGACCTGACCGCGCGC‑GCCTGCGATGCCATTGGCATCGAAGGCTGGGGGCGGTTGGACGTGATGCAGGACGAGCAggg                                                                                                         <  1:1196747/109‑1 (MQ=255)
                                                                              gCTGGATAGCGTCAAGGAACAGGAACTGATCGACCTGACCGCGCGC‑GCCTGCGATGCCATTGGCATCGAAGGCTGGGGGCGGTTGGACGTGATGCAGGACGAGCAGGGCCGGTTCTGGCTGCTCGAAGTCAACACCGCACCAGGCATGAc                                                               <  1:109068/150‑1 (MQ=255)
                                                                                tGGATAGCGTCAAGGAACAGGAACTGATCGACCTGACCGCGCGC‑GCCTGCGATGCCATTGGCATCGaa                                                                                                                                               >  2:367638/1‑68 (MQ=255)
                                                                                tGGATAGCGTCAAGGAACAGGAACTGATCGACCTGACCGCGCGC‑GCCTGCGATGCCATTGGCATCGaa                                                                                                                                               <  1:367638/68‑1 (MQ=255)
                                                                                 ggATAGCGTCAAGGAACAGGAATTGATCGACCTGACCGCGCGC‑GCCTGCGATGCCATTGGCATCGAAGGCTGGAGGCGGTTGGACGTGATGCAGGACGAGCAGGGCCGGTTCTGGCTGCTCGAAGTCAACACCGCACCAGGCATGACAGa                                                            <  1:172497/150‑1 (MQ=255)
                                                                                      gCGTCAGGGAACAGGAACTGATCGACCTGACCGCGCGC‑GCCTGCGATGCCATTGGCATCGAAGGCTGGGGGCGGTTGGACGTGATGCAGGACGAGCAGGGCCGGTTCTGGCTGCTCGAAGTCAACACCGCACCAGGCATGACAGATCATa                                                       <  2:531239/150‑1 (MQ=255)
                                                                                        gTCAAGGAACAGGAACTGATCGACCTGACCGCGCGC‑GCCTGCGATGCCATTGGCATCGAAGGCTGGGGGCGGTTGGACGTGATGCAGGACGAGCAGGGCCGGTTCTGGCTGCTCGAAGTCAACACCGCACCAGGCATGACAGATCATAGc                                                     >  1:1149190/1‑150 (MQ=255)
                                                                                               aacaggaacTGATCGACCTGACCGCGCGC‑GCCTGCGATGCCATTGGCATCGAAGGCTGGGGGCGGTTGGACGTGATGCAGGACGAGCAGGGCCGGTTCTGGCTGCTCGAAGTCAACACCGCACCAGGCATGACAGATCATAGCCTGGTGc                                              >  1:593568/1‑150 (MQ=255)
                                                                                                acaggaacTGATCGACCTGACCGCGCGC‑GCCTGCGATGCCATTGGCATCGAAGGCTGGGGGCGGTTGGACGTGATGCAGGACGAGCAGGGCCGGTTCTGGCTGCTCGAAGTCAACACCGCACCAGGCATGACAGATCATAGCCTGGTGcc                                             >  2:66283/1‑150 (MQ=255)
                                                                                                  aggaacTGATCGACCTGACCGCGCGC‑GCCTGCGATGCCATTGGCATCGAAGGCTGGGGGCGGTTGGACGTGATGCAGGACGAGCAGGGCCGGTTCTGGCTGCTCGAAGTCAACACCGCACCAGGCATGACAGATCATAGCCTGGTGCCCa                                           >  2:1225727/1‑150 (MQ=255)
                                                                                                    gaacTGATCGACCTGACCGCGCGC‑GCCTGCGATGCCATTGGCATCGAAGGCTGGGGGCGGTTGGACGTGATGCAGGACGAGCAGGGCCGGTTCTg                                                                                                >  2:517011/1‑95 (MQ=255)
                                                                                                    gaacTGATCGACCTGACCGCGCGC‑GCCTGCGATGCCATTGGCATCGAAGGCTGGGGGCGGTTGGACGTGATGCAGGACGAGCAGGGCCGGTTCTg                                                                                                <  1:517011/95‑1 (MQ=255)
                                                                                                         gATCCACCAGACCGCGCGC‑GCCTGCGATGCCATTGGCATCGAAGGCTGGGGGCGGTTGGACGTGATGCAGGACGAGCAGGGCCGGTTCTGGCTGCTCGAAGTCAACACCGCACCAGGCATGACAGATCATAGCCTGGTGCCCATggcggc                                    <  2:523136/150‑1 (MQ=255)
                                                                                                          aTCGACCTGACCGCGCGC‑GCCTGCGATGCCATTGGCATCGAAGGCTGGGGGCGGTTGGACGTGATGCAGGACGAGCAGGGCCGGTTCTGGCTGCTCGAAGTCAACACCGCACCAGGCATGACAGATCATAGCCTGGTGCCCATGGCGGcc                                   >  2:273647/1‑150 (MQ=255)
                                                                                                              aCCTGACCGCGCGC‑GCCTGCGATGCCATTGGCATCGAAGGCTGGGGGCGGTTGGACGTGATGCAGGACGAGCAGGGCCGGTTCTGGCTGCTCGAAGTCAACACCGCACCAGGCATGACAGAt                                                           <  2:1346285/122‑1 (MQ=255)
                                                                                                              aCCTGACCGCGCGC‑GCCTGCGATGCCATTGGCATCGAAGGCTGGGGGCGGTTGGACGTGATGCAGGACGAGCAGGGCCGGTTCTGGCTGCTCGAAGTCAACACCGCACCAGGCATGACAGAt                                                           >  1:1346285/1‑122 (MQ=255)
                                                                                                                  gACCGCGCGC‑GCCTGCGATGCCATTGGCATCGAAGGCTGGGGGCGGTTGGACGTGATGCAGGACGAGCAGGGCCGGTTCTGGCTGCTCGAAg                                                                                     >  1:400418/1‑92 (MQ=255)
                                                                                                                  gACCGCGCGC‑GCCTGCGATGCCATTGGCATCGAAGGCTGGGGGCGGTTGGACGTGATGCAGGACGAGCAGGGCCGGTTCTGGCTGCTCGAAg                                                                                     <  2:400418/92‑1 (MQ=255)
                                                                                                                         cgc‑gcCTCCGATGCCATTGGCATCGAAGGCTGGGGGCGGTTGGACGTGATGCAGGACGAGCAGGGCCGGTTCTGGCTGCTCGAAGTCAACACCGAACCAGGCATGACAGATCATAGCCTGGTGCCCATGGCGGCCCGTGCGGCTGGCCTg                    <  1:35785/150‑1 (MQ=255)
                                                                                                                           cagcCTTCGATGCCAATGGCATCGCAGGCTGGGGGCGGTTGGACGTGATGCAGGACGAGCAGGGCCGGTTCTGGCTGCTCGAAGTCAACACCGCACCAGGCATGACAGATCATAGCCTGGTGCCCATGGCGGCCCGTGCGGCTGGCCTgg                   >  1:209889/3‑150 (MQ=255)
                                                                                                                           c‑gcCTGCGATGCCATTGGCATCGAAGGCTGGGGGCGGTTGGACGTGATGCAGGACGAGCAGGGCCGGTTCTGGCTGCTCGAAGTCAACACCGCACCAGGCATGACAGATCATAGCCTGGTGCCCATGGCGGCCCGTGCGGCTGGCCTGGa                  <  2:593568/150‑1 (MQ=255)
                                                                                                                              ccTTCGATGCCAATGGCATCGCAGGCTGGGGGCGGTTGGACGTGATGCAGGACGAGCAGGGCCGGTTCTGGCTGCTCGAAGTCAACACCGCACCAGGCATGACAGATCATAGCCTGGTGCCCATGGCGGCCCGTGCGGCTGGCCTGGACt                >  1:292512/1‑150 (MQ=255)
                                                                                                                                ttcGATGCCAATGGCATCGCAGGCTGGGGGCGGTTGGACGTGATGCAGGACGAGCAGGGCCGGTTCTGGCTGCTCGAAGTCAACACCGCACCAGGCATGACAGATCATAGCCTGGTGCCCATGGCGGCCCGTGCGGCTGGCCTGGACTTc              >  1:388076/3‑150 (MQ=255)
                                                                                                                                 tcGATGCCAATGGCATCGCAGGCTGGGGGCGGTTGGACGTGATGCAGGACGAGCAGGGCCGGTTCTGGCTGCTCGAAGTCAACACCGCACCAGGCATGACAGATCATAGCCTGGTGCCCATGGCGGCCCGTGCGGCTGGCCTGGACTTcc             <  2:292512/149‑1 (MQ=255)
                                                                                                                                        cATTGGCATCGAAGGCTGGGGGCGGTTGGACGTGATGCAGGACGAGCAGGGCCGGTTCTGGCTGCTCGAAGTCAACACCGCACCAGGCATGACAGATCATAGCCTGGTGCCCATGGCGGCCCGTGCGGCTGGCCTGGACTTCCAGCAACt      >  2:1210173/1‑150 (MQ=255)
                                                                                                                                         aTTGGCATCGAAGGCTGGGGGCGGTTGGACGTGATGCAGGACGAGCAGGGCCGGTTCTGGCTGCTCGAAGTCACAACCCCAAACGCCAGTACACAGCATAGCCTGGTGCCCATGGCGGCCCGTGCGGCTGGCCTGGCCTTCCAGCAACTg     >  1:678266/1‑150 (MQ=255)
                                                                                                                                            ggCATCGAAGGCTGGGGGCGGTTGGACGTGATGCAGGACGAGCAGGGCCGGTTCTGGCTGCTCGAAGTCAACACCGCACCAGGCATGACAGATCATAGCCTGGTGCCCATGGCGGCCCGTGCGGCTGGCCTGGACTTCCAGCAACTGGTg  <  1:552691/150‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                              
CCTGGGTACGCCGCACGTGTTCTACGACTATGACGCCAAGTACATCGTCAACGATACCCAGTACCGCATCCCGTGCGGGCTGGATAGCGTCAAGGAACAGGAACTGATCGACCTGACCGCGCGC‑GCCTGCGATGCCATTGGCATCGAAGGCTGGGGGCGGTTGGACGTGATGCAGGACGAGCAGGGCCGGTTCTGGCTGCTCGAAGTCAACACCGCACCAGGCATGACAGATCATAGCCTGGTGCCCATGGCGGCCCGTGCGGCTGGCCTGGACTTCCAGCAACTGGTG  >  NC_002947/1526485‑1526773

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