Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,497,981 G→A 5.9% G376S (GGC→AGC)  opgD → OPG biosynthetic periplasmic protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,497,9810GA 5.9% 208.9 / 2.4 67G376S (GGC→AGC) opgDOPG biosynthetic periplasmic protein
Reads supporting (aligned to +/- strand):  ref base G (33/30);  new base A (2/2);  total (35/32)
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.12e-01

ATGCTTACACCGACAACAACCCGAAAGGGTTTGGTTTATTGCAACTGGATCGTGACTTCTCCCATTATCAGGACATTATGGGCTGGTATAACAAACGCCCAAGTCTGTGGGTGGAACCGCGTAACAAGTGGGGTAAGGGCACCATCGGCCTGATGGAAATCCCAACAACGGGCGAAACGCTGGATAACATTGTCTGCTTCTGGCAGCCAGAAAAAGCTGTAAAAGCAGGTGATGAGTTTGCATTCCAGTATCGTCTGTACTGGAGTGCGCAACCGCCTGTTCATTGCCCA  >  NC_000913/1497835‑1498124
                                                                                                                                                  |                                                                                                                                               
aTGCTTACACCGACAACAACCCGAAAGGGTTTGGTTTATTGCAACTGGATCGTGACTTCTCCCATTATCAGGACATTATGGGCTGGTATAACAAACGCCCAAGTCTGTGGGTGGAACCGCGTAACAAGTGGGGTAAGGGCACCATCGGc                                                                                                                                               <  8:276895/149‑1 (MQ=255)
 tGCTTACACCGACAACAACCCGAAAGGGTTTGGTTTATTGCAACTGGATCGTGACTTCTCCCATTATCAGGACATTATGGGCTGGTATAACAAACGCCCAAGTCTGTGGGTGGAACCGCGTAACAAGTGGGGTAAGGGCACCATCGGcc                                                                                                                                              >  5:325768/1‑149 (MQ=255)
      acacCGACAACAACCCGAAAGGGTTTGGTTTATTGCAACTGGATCGTGACTTCTCCCATTATCAGGACATTATGGGCTGGTATAACAAACGCCCAAGTCTGTGGGTGGAACCGCGTAACAAGTGGGGTAAGGGCACCATCGGCCTGATg                                                                                                                                         >  8:117110/1‑149 (MQ=255)
         ccGACAACCACCCGCAAGGGTTTGGTTTATTGCAACTGGAGCGGGACGTCGCCCATTAGCAGGACATGATGGGCTCGTATAACACACGCCCAAGTCTGTGGGTGGAACCCCGTAACAATTGGGGTAAGGGCACCAGCGGACTGAGGGa                                                                                                                                       >  6:160629/1‑148 (MQ=255)
          cGACAACAACCCGAAAGGGTTTGGTTTATTGCAACTGGATCGTGACTTCTCCCATTATCAGGACATTATGGGCTGGTATAACAAACGCCCAAGTCTGTGGGTGGAACCGCGTAACAAGTGGGGTAAGGGCACCATCGGCCTGATgg                                                                                                                                        <  8:186594/146‑1 (MQ=255)
          cGACAACAACCCGAAAGGGTTTGGTTTATTGCAACTGGATCGTGACTTCTCCCATTATCAGGACATTATGGGCTGGTATAACAAACGCCCAAGTCTGTGGGTGGAACCGCGTAACAAGTGGGGTAAGGGCACCATCGGCCTGATgg                                                                                                                                        >  7:186594/1‑146 (MQ=255)
            acaacaACCCGAAAGGGTTTGGTTTATTGCAACTGGATCGTGACTTCTCCCATTATCAGGACATTATGGGCTGGTATAACAAACGCCCAAGTCTGTGGGTGGAACCGCGTAACAAGTGGGGTAAGGGCACCATCAGCCTGATGGAAATc                                                                                                                                   >  1:20406/1‑149 (MQ=255)
            acaacaACCCGAAAGGGTTTGGTTTATTGCAACTGGATCGTGACTTCTCCCATTATCAGGACATTATGGGCTGGTATAACAAACGCCCAAGTCTGTGGGTGGAACCGCGTAACAAGTGGGGTAAGGGCACCATCAGCCTGATGGAAATc                                                                                                                                   >  1:14818/1‑149 (MQ=255)
              aacaacCCGAAAGGGTTTGGTTTATTGCAACTGGATCGTTACTTCTCCCAGTATCAGGACATTATGGGCTGGTATAACAAACGCCCAAGTCTGTGGGGGGAACCGCGTAACAAGTGGGGTAAGGGCACCATCGGCCTGATGGAAATccc                                                                                                                                 <  8:240129/149‑1 (MQ=255)
              aacaacCCGAAAGGGTTTGGTTTATTGCAACTGGATCGTGACTTCTCCCATTATCAGGACATTATGGGCTGGTATAACAAACGCCCAAGTCTGTGGGTGGAACCGCGTAACAAGTGGGGTAAGGGCACCATCGGCCTGATGGAAATccc                                                                                                                                 >  1:88927/1‑149 (MQ=255)
                      gAAAGGGTTTGGTTTATTGCAACTGGATCGTGACTTCTCCCATTATCAGGACATTATGGGCTGGTATAACAAACGCCCAAGTCTGTGGGTGGAACCGCGTAACAAGTGGGGTAAGGGCACCATCGGCCTGATGGAAATCCCAACAACgg                                                                                                                         >  5:141060/1‑149 (MQ=255)
                      gAAAGGGTTTGGTTTATTGCAACTGGATCGTGACTTCTCCCATTATCAGGACATTATGGGCTGGTATAACAAACGCCCAAGTCTGTGGGTGGAACCGCGTAACAAGTGGGGTAAGGGCACCATCGGCCTGATGGAAATCCCAACAACgg                                                                                                                         >  1:31839/1‑149 (MQ=255)
                      gAAAGGGTTTGGTTTATTGCAACTGGATCGTGACTTCTCCCATTATCAGGACATTATGGGCTGGTATAACAAACGCCCAAGTCTGTGGGTGGAACCGCGTAACAAGTGGGGTAAGGGCACCATCGGCCTGATGGAAATCCCAACAACgg                                                                                                                         >  6:161504/1‑149 (MQ=255)
                                 gtttATTGCAACTGGATCGTGACTTCTCCCATTATCAGGACATTATGGGCTGGTATAACAAACGCCCAAGTCTGTGGGTGGAACCGCGTAACAAGTGGGGTAAGGGCACCATCGGCCTGATGGAAATCCCAACAACGGGCGAAACGCt                                                                                                               <  7:304707/148‑1 (MQ=255)
                                      ttGCAACTGGATCGTGACTTCTCCCATTATCAGGACATTATGGGCTGGTATAACAAACGCCCAAGTCTGTGGGTGGAACCGCGTAACAAGTGGGGTAAGGGCACCATCGGCCTGATGGAAATCCCAACAACGGGCGAAACGCTGGATaa                                                                                                         >  2:276644/1‑149 (MQ=255)
                                          aaCTGGATCGTGACTTCTCCCATTATCAGGACATTATGGGCTGGTATAACAAACGCCCAAGTCTGTGGGTGGAACCGCGTAACAAGTGGGGTAAGGGCACCATCGGCCTGATGGAAATCCCAACAACGGGCGAAACGCTGGATAACAtt                                                                                                     >  1:356942/1‑149 (MQ=255)
                                          aaCTGGATCGTGACTTCTCCCATTATCAGGACATTATGGGCTGGTATAACAAACGCCCAAGTCTGTGGGTGGAACAGCGTAACAAGTGGGGTAAGGGCACCATCGGCCTGATGGAAATCCCAACAACGGGCGAAACGCTGGATAACAtt                                                                                                     >  5:137503/1‑149 (MQ=255)
                                               gttCGTGACTTCTCCCATTATCAGGACATTATGGGCTGGTATAACAAACGCCCAAGTCTGTGGGTGGAACCGCGTAACAAGTGGGGTAAGGGCACCATCGGCCTGATGGAAATCCCAACAACGGGCGAAACGCTGGATAACATTGTc                                                                                                  <  4:117135/145‑1 (MQ=255)
                                                    tGACTTCTCCCATTATCAGGACATTATGGGCTGGTATAACAAACGCCCAAGTCTGTGGGTGGAACCGCGTAACAAGTGGGGTAAGGGCCCCAGCGGCCTGATGGAAATCCCAACACCGGGCGCAACGCGGGATAACATTGTCTGCTTCt                                                                                           >  2:285642/1‑149 (MQ=255)
                                                    tGAATTCTCCCATTATCAGGACATTATTGGCTGGTATAACAAACGCCCAAGTCTGGGGGTGGAACCGCGTAACAAGTGGGGTAAGGGCACCATCGGCCTGAGGGAAATCCCAACAACGGGCGAAACGCTGGATAACATTGTCTGCTTCt                                                                                           <  5:23192/149‑1 (MQ=255)
                                                      aCTTCTCCCATTATCAGGACATTATGGGCTGGTATAACAAACGCCCAAGTCTGTGGGTGGAACCGCGTAACAAGTGGGGTAAGGGCACCATCGGCCTGATGGAAATCCCAACAACGGGCGAAACGCTGGATAACATTGTCTGCTTCTg                                                                                          <  2:23546/148‑1 (MQ=255)
                                                        ttCTCCCATTATCAGGACATTATGGGCTGGTATAACAAACGCCCAAGTCTGTGGGTGGAACCGCGTAACAAGTGGGGTAAGGGCACCATCGGCCTGATGGAAATCCCAACAACGGGCGAAACGCTGGATAACATTGTCTGCTTCTgg                                                                                         >  4:60540/1‑147 (MQ=255)
                                                        ttCTCCCATTATCAGGACATTATGGGCTGGTATAACAAACGCCCAAGTCTGTGGGTGGAACCGCGTAACAAGTGGGGTAAGGGCACCATCGGCCTGATGGAAATCCCAACAACGGGCGAAACGCTGGATAACATTGTCTGCTTCTgg                                                                                         >  2:11860/1‑147 (MQ=255)
                                                        ttCTCCCATTATCAGGACATTATGGGCTGGTATAACAAACGCCCAAGTCTGTGGGTGGAACCGCGTAACAAGTGGGGTAAGGGCACCATCGGCCTGATGGAAATCCCAACAACGGGCGAAACGCTGGATAACATTGTCTGCTTCTGGCa                                                                                       >  2:333391/1‑149 (MQ=255)
                                                            cccATTATCAGGACATTATGGGCTGGTATAACAAACGCCCAAGTCTGTGGGTGGAACCGCGTAACAAGTGGGGTAAGGGCACCATCGGCCTGATGGAAATCCCAACAACGGGCGAAACGCTGGATAACATTGTCTGCTTCTg                                                                                          <  5:49641/142‑1 (MQ=255)
                                                            cccATTATCAGGACATTATGGGCTGGTATAACAAACGCCCAAGTCTGTGGGTGGAACCGCGTAACAAGTGGGGTAAGGGCACCATCGGCCTGATGGAAATCCCAACAACGGGCGAAACGCTGGAAAACATTGTCTGCTTCTg                                                                                          >  6:49641/1‑142 (MQ=255)
                                                             ccATTATCAGGACATTATGGGCTGGTATAACAAACGCCCAAGTCTGTGGGTGGAACCGCGTAACAAGTGGGGTAAGGGCCCCATCGGCCGGATGGAAATCCCAACAACGGGCGAAACGCTGGATAACATTGTCTGCTTCTGGCAGc                                                                                     >  6:245022/1‑146 (MQ=255)
                                                             ccATTATCAGGACATTATGGGCTGGTATAACAAACGCCCAAGTCTGTGGGTGGAACCGCGTAACAAGTGGGGTAAGGGCACCATCGGCCTGATGGAAATCCCAACAACGGGCGAAACGCTGGATAACATTGTCTGCTTCTGGCAGc                                                                                     >  5:285137/1‑146 (MQ=255)
                                                              cATTATCAGGACATTATGGGCTGGTATAACAAACGCCCAAGTCTGTGGGTGGAACCGCGTAACAAGTGGTGTAAGGGCACCATCGGCCTGATGGAAATCCCAACAACGGGCGAAACGCTGGATAACATTGTCTGCTTCTGGCAGCCAg                                                                                  >  8:90766/1‑148 (MQ=255)
                                                                     aGGACATTATGGGCTGGTATAACAAACGCCCAAGTCTGTGGGTGGAACCGCGTAACAAGTGGGGTAAGGGCACCATCGGCCTGATGGAAATCCCAACAACGGGCGAAACGCTGGATAACATTGTCTGCTTCTGGCAGCCAGAAAAAGCt                                                                          >  7:66143/1‑149 (MQ=255)
                                                                     aGGACATTATGGGCTGGTATAACAAACGCCCAAGTCTGTGGGTGGAACCGCGTAACAAGTGGGGTAAGGGCACCATCGGCCTGATGGAAATCCCAACAACGGGCGAAACGCTGGATAACATTGTCTGCTTCTGGCAGCCAGAAAAAGCt                                                                          >  6:148417/1‑149 (MQ=255)
                                                                      ggACATTATGGGCTGGTATAACAAACGCCCAAGTCTGTGGGTGGAACCGCGTAACAAGTGGGGTAAGGGCACCATCGGCCTGATGGAAATCCCAACAACGGGCGAAACGCTGGATAACATTGTCTGCTTCTgg                                                                                         <  5:156101/133‑1 (MQ=255)
                                                                      ggACATTATGGGCTGGTATAACAAACGCCCAAGTCTGTGGGTGGAACCGCGTAACAAGTGGGGTAAGGGCACCATCGGCCTGATGGAAATCCCAACAACGGGCGAAACGCTGGATAACATTGTCTGCTTCTgg                                                                                         >  6:156101/1‑133 (MQ=255)
                                                                      ggACATTATGGGCTGGTATAACAAACGCCCAAGTCTGTGGGTGGAACCGCGTAACAAGTGGGGTAAGGGCACCATCGGCCTGATGGAAATCCCAACAACGGGCGAAACGCTGGATAACATTGTCTGCTTCTGGCAGCCAGAAAAAGCTg                                                                         >  6:125977/1‑149 (MQ=255)
                                                                           ttATGGGCTGGTATAACAAACGCCCAAGTCTGTGGGTGGAACCGCGTAACAAGTGGGGTAAGGGCACCATCGGCCTGATGGAAATCCCAACAACGGGCGAAACGCTGGATAACATTGTCTGCTTCTGGCAGCCAGAAAAAGCTGTaaaa                                                                    <  5:161504/149‑1 (MQ=255)
                                                                           ttATGGGCTGGTATAACAAACGCCCAAGTCTGTGGGTGGAACCGCGTAACAAGTGGGGTAAGGGCACCATCGGCCTGATGGAAATCCCAACAACGGGCGAAACGCTGGATAACATTGTCTGCTTCTGGCAGCCAGAAAAAGCTGTaaaa                                                                    <  8:328311/149‑1 (MQ=255)
                                                                           ttATGGGCTGGTATAACAAACGCCCAAGTCTGTGGGTGGAACCGCGTAACAAGTGGGGTAAGGGCACCATCGGCCTGATGGAAATCCCAACAACGGGCGAAACGCTGGATAACATTGTCTGCTTCTGGCAGCCAGAAAAAGCTGTaaaa                                                                    <  7:279257/149‑1 (MQ=255)
                                                                           gtctGGGCTGGGATAACAAACGCCCAAGTCTGTGGGTGTAACCACGTAACAAGTGGGGTAAGGGCACCATCGGCCTGCTGGAAATCCCAACAACGGGCGAAACGAGGGATAACATTGTCTGCTTCTGGCAGCCAGAAAAAGCTGTaaaa                                                                    <  5:160629/146‑1 (MQ=255)
                                                                                ggCTGGTATAACAAACGCCCAAGTCTGTGGGTGGAACCGCGTAACAAGTGGGGTAAGGGCACCATCGGCCTGATGGAAATCCCAACAACGGGCGAAACGCTGGATAACATTGTCTGCTTCTGGCAGCCAGAAAAAGCTGTAAAAGCAg                                                                <  7:117110/148‑1 (MQ=255)
                                                                                ggCTGGTATAACAAACGCCCAAGTCTGTGGGTGGAACCGCGTAACAAGTGGGGTAAGGGCACCATCAGCCTGATGGAAATCCCAACAACGGGCGTAACGCTGGATAACATTGTCTGCTTCTGGCAGCCAGAAAAAGCTGTAAAAGCAgg                                                               <  2:20406/149‑1 (MQ=255)
                                                                                ggCTGGTATAACAAACGCCCAAGTCTGTGGGTGGAACCGCGTAACAAGTGGGGTAAGGGCACCATCAGCCTGATGGAAATCCCAACAACGGGCGTAACGCTGGATAACATTGTCTGCTTCTGGCAGCCAGAAAAAGCTGTAAAAGCAgg                                                               <  2:14818/149‑1 (MQ=255)
                                                                                     gTATAACAAACGCCCAAGTCTGTGGGTGGAACCGCGTAACAAGTGGGGTAAGGGCACCATCGGCCTGATGGAAATCCCAACAACGGGCGAAACGCTGGATAACATTGTCTGCTTCTGGCAGCCAGAAAAAGCTGTAAAAGCAGGtgat                                                           <  6:325768/148‑1 (MQ=255)
                                                                                        taACAAACGCCCAAGTCTGTGGGTGGAACCGCGTAACAAGTGGGGTAAGGGCACCATCGGCCTGATGGAAATCCCAACAACGGGCGAAACGCTGGATAACATTGTCTGCTTCTGGCAGCCAGAAAAAGCTGTAAAAGCAGGTGATGAGt                                                       >  8:129838/1‑149 (MQ=255)
                                                                                          aCAAACGCCCAAGTCTGTGGGTGGAACCGCGTAACAAGTGGGGTAAGGGCACCATCGGCCTGATGGAAATCCCAACAACGGGCGAAACGCTGGATAACATTGTCTGCTTCTGGCAGCCAGAAAAAGCTGTAAAAGc                                                                  >  6:118330/1‑136 (MQ=255)
                                                                                          aCAAACGCCCAAGTCTGTGGGTGGAACCGCGTAACAAGTGGGGTAAGGGCACCATCGGCCTGATGGAAATCCCAACAACGGGCGAAACGCTGGATAACATTGTCTGCTTCTGGCAGCCAGAAAAAGCTGTAAAAGc                                                                  <  5:118330/136‑1 (MQ=255)
                                                                                          aCAAACGCCCAAGTCTGTGGGTGGAACCGCGTAACAAGTGGGGTAAGGGCACCATCGGCCTGATGGAAATCCCAACAACGGGCGAAACGCTGGATAACATTGTCTGCTTCTGGCAGCCAGAAAAAGCTGTAAAAGCAGGTGATGAGttt                                                     <  3:114400/149‑1 (MQ=255)
                                                                                             aaCGCCCAAGTCTGTGGGTGGAACCGCGTAACAAGTGGGGTAAGGGCACCATCGGCCTGATGGAAATCCCAACAACGGGCGAAACGCTGGATAACATTGTCTGCTTCTGGCAGCCAGAAAAAGCTGTAAAAGCAGGTGATGAGTTTGCa                                                  <  8:49836/149‑1 (MQ=255)
                                                                                                  ccAAGTCTGTGGGTGGAACCGCGTAACAAGTGGGGTAAGGGCACCATCGGCCTGATGGAAATCCCAACAACGGGCGAAACGCTGGATAACa                                                                                                       <  2:124331/91‑1 (MQ=255)
                                                                                                  ccAAGTCTGTGGGTGGAACCGCGTAACAAGTGGGGTAAGGGCACCATCGGCCTGATGGAAATCCCAACAACGGGCGAAACGCTGGATAACATTGTCTGCTTCTGGCAGCCAGAAAAAGCTGTAAAAGCAGGTGATGAGTTTGCATTCCa                                             <  5:125977/149‑1 (MQ=255)
                                                                                                     aGTCTGTGGGTGGAACCGCGTAACAAGTGGGGTAAGGGCACCATCGGCCTGATGGAAATCCCAACAACGGGCGAAACGCTGGATAACATTGTCTGCTTCTGGCAGCCAGAAAAAGCTGTAAAAGCAGGTGATGAGTTTGCATTCCAGTa                                          <  6:285137/149‑1 (MQ=255)
                                                                                                              gtggAACCGCGTAACAAGTGGGGTAAGGGCACCATCGGCCTGATGGAAATCCCAACAACGGGCGAAACGCTGGATAACATTGTCTGCTTCTGGCAGCCAGAAAAAGCTGTAAAAGCAGGTGATGAGTTTGCATTCCAGTATCGTCTGTa                                 >  3:61333/1‑149 (MQ=255)
                                                                                                               tggAACCGCGTAACAAGTGGGGTAAGGGCACCATCGGCCTGATGGAAATCCCAACAACGGGCGAAACGCTGGATAACATTGTCTGCTTCTGGCAGCCAGAAAAAGCTGTAAAAGCAGGTGATGAGTTTGCATTCCAGTATCGTCTGTAc                                <  3:60540/149‑1 (MQ=255)
                                                                                                                 gAACCGCGTAACAAGTGGGGTAAGGGCACCATCGGCCTGATGGAAATCCCAACAACGGGCGAAACGCTGGATAACATTGTCTGCTTCTGGCAGCCAGAAAAAGCTGTAAAAGCAGGTGATGAGTTTGCATTCCAGTATCGTCTGTACTg                              <  5:148417/149‑1 (MQ=255)
                                                                                                                   aCCGCGTAACTAGTGGGGTAAGGGCACCATCGGCCTGATGGAAATCGCAACAACGGGCGAAACGCTGGATAACATTGTCTGCTTCTGGCAGCCAGAAAAAGCTGTAAAAGCAGGTGATGAGTTTGCATTCCAGTATCGTCTGTACTGGa                            <  5:93806/149‑1 (MQ=255)
                                                                                                                   aCCGCGTAACAAGTGGGGTAAGGGCACCATCGGCCTGATGGAAATCCCAACAACGGGCGAAACGCTGGATAACATTGTCTGCTTCTGGCAGCCAGAAAAAGCTGTAAAAGCAGGTGATGAGTTTGCATTCCAGTATCGTCTGTACTGGa                            <  4:184182/149‑1 (MQ=255)
                                                                                                                     cgcgTAACAAGTGGGGTAAGGGCACCATCGGCCTGATGGAAATCCCAACAACGGGCGAAACGCTGGATAACATTGTCTGCTTCTGGCAGCCAGAAAAAGCTGTAAAAGCAGGTGATGAGTTTGCATTCCAGTATCGTCTGTACTGGAGt                          <  1:333391/149‑1 (MQ=255)
                                                                                                                       cgTAACAAGTGGGGTAAGGGCACCATCGGCCTGATGGAAATCCCAACAACGGGCGAAACGCTGGATAACATTGTCTGCTTCTGGCAGCCAGAAAAAGCTGTAAAAGCAGGTGATGAGTTTGCATTCCAGTATCGTCTGTACTGGAGTgc                        <  8:321536/149‑1 (MQ=255)
                                                                                                                          aaCAAGTGGGGTAAGGGCACCATCGGCCTGATGGAAATCCCAACAACGGGCGAAACGCTGGATAACATTGTCTGCTTCTGGCAGCCAGAAAAAGCTGTAAAAGCAGGTGATGAGTTTGCATTCCAGTATCGTCTGTACTGGAGTgcgc                      >  1:304529/1‑148 (MQ=255)
                                                                                                                          aaCAAGTGGGGTAAGGGCACAATCGGCCTGATGGAAATCCCAACAACGGGCGAAACGCTGGATAACATTGTCTGCTTCTGGCAGCCAGAAAAAGCTGTAAAAGCAGGTGATGAGTTTGCATTCCAGTATCGTCTGTACTGGAGTgcgc                      >  2:234877/1‑148 (MQ=255)
                                                                                                                           aCAAGTGGGGTAAGGGCACCATCGGCCTGATGGAAATCCCAACAACGGGCGAAACGCTGGATAACATTGTCTGCTTCTGGCAGCCAGAAAAAGCTGTAAAAGCAGGTGATGAGTTTGCATTCCAGTATCGTCTGTACTGGAGTGCGCa                     <  2:356942/148‑1 (MQ=255)
                                                                                                                                  gggTAAGGGCACCATCGGCCTGATGGAAATCCCAACAACGGGCGAAACGCTGGATAACATTGTCTGCTTCTGGCAGCCAGAAAAAGCTGTAAAAGCAGGTGATGAGTTTGCATTCCAGTATCGTCTGTACTGGAGTGCGCAACCGCCTg             >  6:322497/1‑149 (MQ=255)
                                                                                                                                  gggTAAGGGCACCATCGGCCTGATGGAAATCCCAACAACGGGCGAAACGCTGGATAACATTGTCTGCTTCTGGCAGCCAGAAAAAGCTGTAAAAGCAGGTGATCAGTTTGCATTCCAGTATCGTCTGTACTGGAGTGCGCAACCGCCTg             <  1:285642/149‑1 (MQ=255)
                                                                                                                                    gTAAGGGCACCATCGGCCTGATGGAAATCCCAACAACGGGCGAAACGCTGGATAACATTGTCTGCTTCTGGCAGCCAGAAAAAGCTGTAAAAGCAGGTGATGAGTTTGCATTCCAGTATCGTCTGTACTGGAGTGCGCAACCGCCTGt            >  5:172967/1‑148 (MQ=255)
                                                                                                                                          gCACCATCGGCCTGATGGAAATCCCAACAACGGGCGAAACGCTGGATAACATTGTCTGCTTCTGGCAGCCAGAAAAAGCTGTAAAAGCAGGTGATGAGTTTGCATTCCAGTATCGTCTGTACTGGAGTGCGCAACCGCCTGTTCATTGc     <  1:276644/149‑1 (MQ=255)
                                                                                                                                             ccATCGGCCTGATGGAAATCCCAACAACGGGCGAAACGCTGGATAACATTGTCTGCTTCTGGCAGCCAGAAAAAGCTGTAAAAGCAGGTGATGAGTTTGCATTCCAGTATCGTCTGTACTGGAGTGCGCAAcc                  >  8:743/1‑133 (MQ=255)
                                                                                                                                             ccATCGGCCTGATGGAAATCCCAACAACGGGCGAAACGCTGGATAACATTGTCTGCTTCTGGCAGCCAGAAAAAGCTGTAAAAGCAGGTGATGAGTTTGCATTCCAGTATCGTCTGTACTGGAGTGCGCAACCGTCTGTTCATTGCCCa  >  4:163299/1‑149 (MQ=255)
                                                                                                                                             ccATCGGCCTGATGGAAATCCCAACAACGGGCGAAACGCTGGATAACATTGTCTGCTTCTGGCAGCCAGAAAAAGCTGTAAAAGCAGGTGATGAGTTTGCATTCCAGTATCGTCTGTACTGGAGTGCGCAACCGCCTGTTCATTGCCCa  >  5:159359/1‑149 (MQ=255)
                                                                                                                                                  |                                                                                                                                               
ATGCTTACACCGACAACAACCCGAAAGGGTTTGGTTTATTGCAACTGGATCGTGACTTCTCCCATTATCAGGACATTATGGGCTGGTATAACAAACGCCCAAGTCTGTGGGTGGAACCGCGTAACAAGTGGGGTAAGGGCACCATCGGCCTGATGGAAATCCCAACAACGGGCGAAACGCTGGATAACATTGTCTGCTTCTGGCAGCCAGAAAAAGCTGTAAAAGCAGGTGATGAGTTTGCATTCCAGTATCGTCTGTACTGGAGTGCGCAACCGCCTGTTCATTGCCCA  >  NC_000913/1497835‑1498124

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