Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 2,699,833 G→C 8.6% P226P (CCC→CCG) ‡ fsr‑II ← fosmidomycin efflux system

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029472,699,8330GC 8.6% 179.8 / 7.0 58P226P (CCC→CCG) ‡fsr‑IIfosmidomycin efflux system
Reads supporting (aligned to +/- strand):  ref base G (27/26);  new base C (3/2);  total (30/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 = 2.72e-01

GTTGGCTGTGCTGCACACTCAGGCCGAAATGCTCGATCAGGTAGAAGGTGAAGAAGTTGGTGAACGAAGCGATGTAGACGAACTTGGCGAACATCAGCACGGCGATCACCGCCACAGCGCGCCACATGGCACCCTGCGACAGGCCGGGGGCCTGTTGGCCCGCGAAGGTCTTGAGCCGGGCCTGGCCGTGGCGCACGCTCCAGCCCGTCACCCGCAGTAGCACCAGAACAGCCAGCGCAGCGGCCAGCATGAACCAGGCAATGGCCGGTTGGCCGTGCGGAATGATAA  >  NC_002947/2699688‑2699975
                                                                                                                                                 |                                                                                                                                              
gTTGGCTGTGCTGCACACTCAGGCCGAAATGCTCGATCAGGTAGAAGGTGAAGAAGTTGGTGAACGAAGCGATGTAGACGAACTTGGCGAACATCAGCACGGCGATCACCGCCACAGCGCGCCACATGGCACCCTGCGACAGGCCgggg                                                                                                                                             <  1:109342/149‑1 (MQ=255)
gTTGGCTGTGCTGCACACTCAGGCCGAAATGCTCGATCAGGTAGAAGGTGAAGAAGTTGGTGAACGAAGCGATGTAGACGAACTTGGCGAACATCAGCACGGCGATCACCGCCACAGCGCGCCACATGGCACCCTGCGACAGGCCgggg                                                                                                                                             <  8:179587/149‑1 (MQ=255)
gTTGGCTGTGCTGCACACTCAGGCCGAAATGCTCGATCAGGTAGAAGGTGAAGAAGTTGGTGAACGAAGCGATGTAGACGAACTTGGCGAACATCAGCACGGCGATCACCGCCACAGCGCGCCACATGGCACCCTGCGACAGGCCgggg                                                                                                                                             <  6:341384/149‑1 (MQ=255)
gTTGGCTGTGCTGCACACTCAGGCCGAAATGCTCGATCAGGTAGAAGGTGAAGAAGTTGGTGAACGAAGCGATGTAGACGAACTTGGCGAACATCAGCACGGCGATCACCGCCACAGCGCGCCACATGGCACCCTGCGACAGGCCgggg                                                                                                                                             >  7:11939/1‑149 (MQ=255)
      tgtgCTGCACACTCAGGCCGAAATGCTCGATCAGGTAGAAGGTGAAGAAGTTGGTGAACGAAGCGATGTAGACGAACTTGGCGAACATCAGCACGGCGATCACCGCCACAGCGCGCCACATGGCACCCTGCGACAGGCCGGGGGCCTGt                                                                                                                                       >  3:191039/1‑149 (MQ=255)
        tgctgcACACTCAGGCCGAAATGCTCGATCAGGTAGAAGGTGAAGAAGTTGGTGAACGAAGCGATGTAGACGAACTTGGCGAACATCAGCACGGCGATCACCGCCACAGCGCGCCACATGGCACCCTGCGACAGGCCGGGGGCCTGtt                                                                                                                                      <  8:11939/148‑1 (MQ=255)
           tgcACACTCAGGCCGAAATGCTCGATCAGGTAGAAGGTGAAGAAGTTGGTGAACGAAGCGATGTAGACGAACTTGGCGAACATCAGCACGGCGATCACCGCCACAGCGCGCCACATGGCACCCTGCGACAGGCCGGGGGCCTGTTGGcc                                                                                                                                  >  1:234628/1‑149 (MQ=255)
           tgcACACTCAGGCCGAAATGCTCGATCAGGTAGAAGGTGAAGAAGTTGGTGAACGAAGCGATGTAGACGAACTTGGCGAACATCAGCACGGCGATCACCGCCACAGCGCGCCACATGGCACCCTGCAACAGGCCCCCGGCCTGtcgcag                                                                                                                                  <  2:137027/149‑6 (MQ=255)
                    aGGCCGAAATGCTCGATCAGGTAGAAGGTGAAGAAGTTGGTGAACGAAGCGATGTAGACGAACTTGGCGAACATCAGCACGGCGATCACCGCCACAGCGCGCCACATGGCACCCTGCGACAGGCCGGGGGCCTGTTg                                                                                                                                     <  1:216428/137‑1 (MQ=255)
                    aGGCCGAAATGCTCGATCAGGTAGAAGGTGAAGAAGTTGGTGAACGAAGCGATGTAGACGAACTTGGCGAACATCAGCACGGCGATCACCGCCACAGCGCGCCACATGGCACCCTGCGACAGGCCGGGGGCCTGTTg                                                                                                                                     >  2:216428/1‑137 (MQ=255)
                     ggCCGAAATGCTCGATCAGGTAGAAGGTGAAGAAGTTGGTGAACGAAGCGATGTAGACGAACTTGGCGAACATCAGCACGGCGATCACCGCCACAGCGCGCCACATGGCACCCTGCGACAGGCCGGGGGCCTGTTGGCCCGCGAAGGTc                                                                                                                        <  3:206320/149‑1 (MQ=255)
                     ggCCGAAATGCTCGATCAGGTAGAAGGTGAAGAAGTTGGTGAACGAAGCGATGTAGACGAACTTGGCGAACATCAGCACGGCGATCACCGCCACAGCGCGCCACATGGCACCCTGCGACAGGCCGGGGGCCTGTTGGCCCGCGAAGGTc                                                                                                                        <  2:234877/149‑1 (MQ=255)
                      gCCGAAATGCTCGATCAGGTAGCAGGGGCAGAAGTTGGTGAACGAAGCGATGTAGACGAACTTGGCGAACATCAGCACGGCGATCACCGCCACAGCGCGCCACATGGCACCCTGCGACAGGCCGGGGGCCTGTTGGCCCGCGAAGGTCt                                                                                                                       <  6:184867/149‑1 (MQ=255)
                        cGAAATGCTCGATCAGGTAGAAGGTGAAGAAGTTGGTGAACGAAGCGATGTAGACGAACTTGGCGAACATCAGCACGGCGATCACCGCCACAGCGCGCCACATGGCACCCTGCGACAGGCCGGGGGCCTGTTGGCCCGCGAAGGTCtt                                                                                                                      >  1:102505/1‑148 (MQ=255)
                         gAAATGCTCGATCAGGTAGAAGGTGAAGAAGTTGGTGAACGAAGCGATGTAGACGAACTTGGCGAACATCAGCACGGCGATCACCGCCACAGCGCGCCACATGGCACCCTGCGACAGGCCGGGGGCCTGTTGGCCCGCGAAGGTCtt                                                                                                                      >  5:45750/1‑147 (MQ=255)
                                  gATCAGGTAGAAGGTGAAGAAGTTGGTGAACGAAGCGATGTAGACGAACTTGGCGAACATCAGCACGGCGATCACCGCCACAGCGCGCCACATGGCACCCTGCGACAGGCCGGGGGCCTGTTGGCCCGCGAAGGTCTTGAGCCGGGCCt                                                                                                           >  7:45252/1‑149 (MQ=255)
                                    tCAGGTAGAAGGTGAAGAAGTTGGTGAACGAAGCGATGTAGACGAACTTGGCGAACATCAGCACGGCGATCACCGCCACAGCGCGCCACATGGCACCCTGCGACAGGCCGGGGGCCTGTTGGCCCGCGAAGGTCTTg                                                                                                                     >  1:182884/1‑137 (MQ=255)
                                    tCAGGTAGAAGGTGAAGAAGTTGGTGAACGAAGCGATGTAGACGAACTTGGCGAACATCAGCACGGCGATCACCGCCACAGCGCGCCACATGGCACCCTGCGACAGGCCGGGGGCCTGTTGGCCCGCGAAGGTCTTg                                                                                                                     <  2:182884/137‑1 (MQ=255)
                                      atgTAGAAGGTGAAGAAGTTGGTGAACGAAGCGATGTAGACGAACTTGGCGAACATCAGCACGGCGATCACCGCCACAGCGCGCCACATGGCACCCTGCGACAGGCCGGGGGCCTGTTGGCCCGCGAAGGTCTTGAGCCGGGCCTGGc                                                                                                        <  2:102505/146‑1 (MQ=255)
                                      aGGTAGAAGGTGAAGAAGTTGGTGAACGAAGCGATGTAGACGAACTTGGCGAACATCAGCACGGCGATCACCGCCACAGCGCGCCACATGGCACCCTGCGACAGGCCGGGGGCCTGTTGGCCCGCGAAGGTCTTGAGCCGGGCCTGGc                                                                                                        <  6:45750/148‑1 (MQ=255)
                                        gTAGAAGGTGAAGAAGTTGGTGAACGAAGCGATGTAGACGAACTTGGCGAACATCAGCACGGCGATCACCGCCACAGCGCGCCACATGGCACCCTGCGACAGGCCGGGGGCCTGTTGGCCCGCGAAGGTCTTGAGCCGGGCCTGGCCGt                                                                                                     >  3:209442/1‑149 (MQ=255)
                                        gTAGAAGGTGAAGAAGTTGGTGAACGAAGCGATGTAGACGAACTTGGCGAACATCAGCACGGCGATCACCGCCACAGCGCGCCACATGGCACCCTGCGACAGGCCGGGGGCCTGTTGGCCCGCGAAGGTCTTGAGCCGGGCCTGGCCGt                                                                                                     <  3:287521/149‑1 (MQ=255)
                                                 gaagaaGTTGGTGAACGAAGCGATGTAGACGAACTTGGCGAACATCAGCACGGCGATCACCGCCACAGCGCGCCACATGGCACCCTGCGACAGGCCGGGGGCCTGTTGGCCCGCGAAGGTCTTGAGCCGGGCCTGGCCGTg                                                                                                    <  5:117158/141‑1 (MQ=255)
                                                 gaagaaGTTGGTGAACGAAGCGATGTAGACGAACTTGGCGAACATCAGCACGGCGATCACCGCCACAGCGCGCCACATGGCACCCTGCGACAGGCCGGGGGCCTGTTGGCCCGCGAAGGTCTTGAGCCGGGCCTGGCCGTg                                                                                                    >  6:117158/1‑141 (MQ=255)
                                                  aagaagTTGGTGAACGAAGCGATGTAGACGAACTTGGCGAACATCAGCACGGCGATCACCGCCACAGCGCGCCACATGGCACCCTGCGACAGGCCGGGGGCCTGTTGGCCCGCGAAGGTCTTGAGCCGGGCCTGGCCGTGGCGCACGCt                                                                                           >  7:311958/1‑149 (MQ=255)
                                                   agaagTTGGTGAACGAAGCGATGTAGACGAACTTGGCGAACATCAGCACGGCGATCACCGCCACAGCGCGCCACATGGCACCCTGCGACAGGCCGGGGGCCTGTTGGCCCGCGAAGGTCTTGAGCCGGGCCTGGCCGTGGCGCACGCTc                                                                                          <  8:45252/149‑1 (MQ=255)
                                                    gaagTTGGTGAACGAAGCGATGTAGACGAACTTGGCGAACATCAGCACGGCGATCACCGCCACAGCGCGCCACATGGCACCCTGCGACAGGCCGGGGGCCTGTTGGCCCGCGAAGGTCTTGAGCCGGGCCTGGCCGTGGCGCACGCTc                                                                                          >  3:29165/1‑148 (MQ=255)
                                                     aagTTGGTGAACGAAGCGATGTAGACGAACTTGGCGAACATCAGCACGGCGATCACCGCCACAGCGCGCCACATGGCACCCTGCGACAGGCCGGGGGCCTGTTGGCCCGCGAAGGTCTTGAGCCGGGCCTGGCCGTGGCGCACGCTCCa                                                                                        >  3:287903/1‑149 (MQ=255)
                                                             gAACGAAGCGATGTAGACGAACTTGGCGAACATCAGCACGGCGATCACCGCCACAGCGCGCCACATGGCACCCTGCGACAGGCCGGGGGCCTGTTGGCCCGCGAAGGTCTTGAGCCGGGCCTGGCCGTGGCGCACg                                                                                             >  2:68755/1‑136 (MQ=255)
                                                             gAACGAAGCGATGTAGACGAACTTGGCGAACATCAGCACGGCGATCACCGCCACAGCGCGCCACATGGCACCCTGCGACAGGCCGGGGGCCTGTTGGCCCGCGAAGGTCTTGAGCCGGGCCTGGCCGTGGCGCACg                                                                                             <  1:68755/136‑1 (MQ=255)
                                                                  aaGCGATGTAGACGAACTTGGCGAACATCAGCACGGCGATCACCGCCACAGCGCGCCACATGGCACCCTGCGACAGGCCGGGGGCCTGTTGGCCCGCGAAGGTCTTGAGCCGGGCCTGGCCGTGGCGCACGCTCCAGCCCGTCACCCGc                                                                           <  5:330826/149‑1 (MQ=255)
                                                                    gCGATGTAGACGAACTTGGCGAACATCAGCACGGCGATCACCGCCACAGCGCGCCACATGGCACCCTGCGACAGGCCGGGGGCCTGTTGGCCCGCGAAGGTCTTGAGCCGGGCCTGGCCGTGGCGCACGCTCCAGCCCGTCACCCGCAg                                                                         <  4:29165/149‑1 (MQ=255)
                                                                    gCGATGTAGACGAACTTGGCGAACATCAGCACGGCGATCACCGCCACAGCGCGCCACATGGCACCCTGCGACAGGCCGGGGGCCTGTTGGCCCGCGAAGGTCTTGAGCCGGGCCTGGCCGTGGCGCACGCTCCAGCCCGTCACCCGCAg                                                                         <  8:311958/149‑1 (MQ=255)
                                                                       aTGTAGACGAACTTGGCGAACATCAGCACGGCGATCACCGCCACAGCGCGCCACATGGCACCCTGCGACAGGCCGGGGGCCTGTTGGCCCGCGAAGGTCTTGAGCCGGGCCTGGCCGTGGCGCACGCTCCAGCCCGTCACCCGCAGTAg                                                                      >  2:122104/1‑149 (MQ=255)
                                                                               gAACTTGGCGAACATCAGCACGGCGATCACCGCCACAGCGCGCCAGATGGCACCCTGCGACAGGCCGGGGGCCTGTTGGCCCGCGAAGGTCTTGAGCCGGGCCTGGCCGTGGCGCACGCTCCAGCCCGTCACCCGCAGTAGCACCAGaa                                                              >  1:262253/1‑149 (MQ=255)
                                                                               gAACTTGGCGAACATCAGCACGGCGATCACCGCCACAGCGCGCCACATGGCACCCTGCGACAGGCCGGGGGCCTGTTGGCCCGCGAAGGTCTTGAGCCGGGCCTGGCCGTGGCGCACGCTCCAGCCCGTCACCCGCAGTAGCACCAGaa                                                              <  4:209442/149‑1 (MQ=255)
                                                                                       cGAACATCAGCACGGCGATCACCGCCACAGCGCGCCACATGGCACCCTGCGACAGGCCGGGGGCCTGTTGGCCCGCGAAGGTCTTGAGCCGGGCCTGGCCGTGGCGCACGCTCCAGCCCGTCACCCGCAGTAGCACCAGAAcagccagc                                                      >  1:17834/1‑149 (MQ=255)
                                                                                        gAACATCAGCACGGCGATCACCGCCACAGCGCGCCACATGGCACCCTGCGACAGGCCGGGGGCCTGTTGGCCCGCGAAGGTCTTGAGCCGGGCCTGGCCGTGGCGCACGCTCCAGCCCGTCACCCGCAGTAGc                                                                     <  6:377796/133‑1 (MQ=255)
                                                                                        gAACATCAGCACGGCGATCACCGCCACAGCGCGCCACATGGCACCCTGCGACAGGCCGGGGGCCTGTTGGCCCGCGAAGGTCTTGAGCCGGGCCTGGCCGTGGCGCACGCTCCAGCCCGTCACCCGCAGTAGCACCAGAACAGCcagcg                                                     >  2:175021/1‑149 (MQ=255)
                                                                                        gAACATCAGCACGGCGATCACCGCCACAGCGCGCCACATGGCACCCTGCGACAGGCCGCGGGCCTGTTGGCCCGCGAAGGTCTTGAGCCGGGCCTGGCCGTGGCGCACGCTCCAGCCCGTCACCCGCAGTAGc                                                                     >  5:377796/1‑133 (MQ=255)
                                                                                           cATCAGCACGGCGATCACCGCCACAGCGCGCCACATGGCACCCTGCGACAGGCCGGGGGCCTGTTGGCCCGCGAAGGTCTTGAGCCGGGCCTGGCCGTGGCGCACGCTCCAGCCCGTCACCCGCAGTAGCACCAGAACAGCcagcgcag                                                  <  2:234628/149‑1 (MQ=255)
                                                                                              cAGCACGGCGATCACCGCCACAGCGCGCCACATGGCACCCTGCGACAGGCCGGGGGCCTGTTGGCCCGCGAAGGTCTTAAGCCGGGCCTGGCCGTGGCGCACGCTCCAGCCCGTCACCCGCAGTAGCACCAGAACAGCCAGCGCAGCgg                                               >  6:370846/1‑149 (MQ=255)
                                                                                               aGCACGGCGATCACCGCCACAGCGCGCCACATGGCACCCTGCGACAGGCCGGGGGCCTGTTGGCCCGCGAAGGTCTTAAGCCGGGCCTGGCCGTGGCGCACGCTCCAGCCCGTCACCCGCAGTAGCACCAGAACAGCCAGCGCAGCgg                                               <  5:370846/148‑1 (MQ=255)
                                                                                                            ccGCCACAGCGCGCCACATGGCACCCTGCGACAGGCCGGGGGCCTGTTGGCCCGCGAAGGTCTTGAGCCGGGCCTGGCCGTGGCGCACGCTCCAGCCCGTCACCCGCAGTAGCACCAGAACAGCCAGCGCAGCGGCCAGCATGAACCAg                                 >  8:36122/1‑149 (MQ=255)
                                                                                                             cGCCACAGCGCGCCACATGGCACCCTGCGACAGGCCGGGGGCCTGTTGGCCCGCGAAGGTCTTGAGCCGGGCCTGGCCGTGGCGCACGCTCCAGCCCGTCACCCGCAGTAGCACCAGAACAGCCAGCGCAGCGGCCAGCATGAACCAgg                                <  1:122104/149‑1 (MQ=255)
                                                                                                               ccACAGCGCGCCACATGGCACCCTGCGACAGGCCGGGGGCCTGTTGGCCCGCGAAGGTCTTGAGCCGGGCCTGGCCGTGGCGCACGCTCCAGCCCGTCACCCGCAGTAGCACCAGAACAGCCAGCGCAGCGGCCAGCATGAACCAGGc                               >  3:241459/1‑148 (MQ=255)
                                                                                                                cacaGCGCGCCACATGGCACCCTGCGACAGGCCGGGGGCCTGTTGGCCCGCGAAGGTCTTGAGCCGGGCCTGGCCGTGGCGCACGCTCCAGCCCGTCACCCGCAGTAGCACCAGAACAGCCAGCGCAGCGGCCAGCATGAACCAGGCaa                             >  5:79012/1‑149 (MQ=255)
                                                                                                                      gcgcCACATGGCACCCTGCGACAGGCCGGGGGCCTGTTGGCCCGCGAAGGTCTTGAGCCGGGCCTGGCCGTGGCGCACGCTCCAGCCCGTCACCCGCAGTAGCACCAGAACAGCCAGCGCAGCGGCCAGCATGAACCAGGCAATGGCCg                       <  4:241459/149‑1 (MQ=255)
                                                                                                                       cgcCACATGGCACCCTGCGACAGGCCGGGGGCCTGTTGGCCCGCGAAGGTCTTGAGCCGGGCCTGGCCGTGGCGCACGCTCCAGCCCGTCACCCGCCGTAGCACCAGAACAGCCAGCGCAGCGGCCAGCATGAACCAGGCAATGGCCgg                      <  6:162207/149‑1 (MQ=255)
                                                                                                                           acaTGGCACCCTGCGACAGGGCGGGGGCCTGTTGGCCCGCGAAGGTCTTGAGCCGGACCTGGCCGTGGCGCACGCTCCAGCCCGTCACCCGCAGTAGCACCAGAACAGCCAGCGCAGCGGCCAGCATGAACCAGGCAATGGCCGGTTgg                  <  4:191039/149‑1 (MQ=255)
                                                                                                                                 cACCCTGCGACAGGCCGGGGGCCTGTTGGCCCGCGAAGGTCTTGAGCCGGGCCTGGCCGTGGCGCACGCTCCAGCCCGTCACCCGCAGTAGCACCAGAACAGCCAGCGCAGCGGCCAGCATGAACCAGGCAATGGCCGGTTGGCCGTGc            >  7:80821/1‑149 (MQ=255)
                                                                                                                                    ccTGCGACAGGCCGGGGGCCTGTTGGCCCGCGAAGGTCTTGAGCCGGGCCTGGCCGTGGCGCACGCTCCAGCCCGTCACCCGCAGTAGCACCAGAACAGCCAGCGCAGCGGCCAGCATGAACCAGGCAATGGCCGGTTGGCCGTGCgg          >  7:311727/1‑148 (MQ=255)
                                                                                                                                    ccTGCGACAGGCCGGGGGCCTGTTGGCCCGCGAAGGTCTTGAGCCGGGCCTGGCCGTGGCGCACGCTCCAGCCCGTCACCCGCAGTAGCACCAGAACAGCCAGCGCAGCGGCCAGCATGAACCAGGCAATGGCCGGTTGGCCGTGCgg          >  3:313321/1‑148 (MQ=255)
                                                                                                                                    ccTGCGACAGGCCGGGGGCCTGTTGGCCCGCGAAGGTCTTGAGCCGGGCCTGGCCGTGGCGCACGCTCCAGCCCGTCACCCGCAGTAGCACCAGAACAGCCAGCGCAGCGGCCAGCATGAACCAGGCAATGGCCGGTTGGCCGTGCgg          >  7:107025/1‑148 (MQ=255)
                                                                                                                                      gccaaCAGGCCCCCGGCCTGTCGGCCCGCGAAGGTCTTGAGCCGGGCCTGGCCGTGGCGCACGCTCCAGCCCGTCACCCGCAGTAGCACCAGAACAGCCAGCGCAGCGGCCAGCATGAACCAGGCAATGGCCGGTTGGCCGTGCGGAAt       >  5:122033/5‑149 (MQ=255)
                                                                                                                                         gACAGGCCGGGGGCCTGTTGGCCCGCGAAGGTCTTGAGCCGGGCCTGGCCGTGGCGCACGCTCCAGCCCGTCACCCGCAGTAGCACCAGAACAGCCAGCGCAGCGGCCAGCATGAACCAGGCAATGGCCGGTTGGCCGTGCGGAATGAt    <  2:262253/149‑1 (MQ=255)
                                                                                                                                          aCAGGCCCCCGGCCTGTTGGCCCGCGAAGGTCTTGAGCCGGGCCTGGCCGTGGCGCACGCTCCAGCCCGTCACCCGCAGTAGCACCAGAACAGCCAGCGCAGCGGCCAGCATGAAcc                                   >  2:53748/1‑117 (MQ=255)
                                                                                                                                          aCAGGCCCCCGGCCTGTTGGCCCGCGAAGGTCTTGAGCCGGGCCTGGCCGTGGCGCACGCTCCAGCCCGTCACCCGCAGTAGCACCAGAACAGCCAGCGCAGCGGCCAGCATGAAcc                                   <  1:53748/117‑1 (MQ=255)
                                                                                                                                          aCAGGCCCCCGGCCTGTTGGCCCGCGAAGGTCTTGAGCCGGGCCTGGCCGTGGCGCACGCTCCAGCCCGTCACCCGCAGTAGCACCAGAACAGCCAGCGCAGCGGCCAGCATGAACCAGGCAATGGCCGGTTGGCCGTGCGGAATGAt    >  3:291633/1‑148 (MQ=255)
                                                                                                                                           cAGGCCGGGGGCCTGTTGGCCCGCGAAGGTCTTGAGCCGGGCCTGGCCGTGGCGCACGCTCCAGCCCGTCACCCGCAGTAGCACCAGAACAGCCAGCGCAGCGGCCAGCATGAACCAGGCAATGGCCGGTTGGCCGTGCGGAATGATag  <  6:79012/149‑2 (MQ=255)
                                                                                                                                              gCCGGGGGCCTGTTGGCCCGCGAAGGTCTTGAGCCGGGCCTGGCCGTGGCGCACGCTCCAGCCCGTCACCCGCAGTAGCACCAGAACAGCCAGCGCAGCGGCCAGCATGAACCAGGCAATGGCCGGTTGGCCGTGc            <  1:372427/136‑1 (MQ=255)
                                                                                                                                              gCCGGGGGCCTGTTGGCCCGCGAAGGTCTTGAGCCGGGCCTGGCCGTGGCGCACGCTCCAGCCCGTCACCCGCAGTAGCACCAGAACAGCCAGCGCAGCGGCCAGCATGAACCAGGCAATGGCCGGTTGGCCGTGc            >  2:372427/1‑136 (MQ=255)
                                                                                                                                                 |                                                                                                                                              
GTTGGCTGTGCTGCACACTCAGGCCGAAATGCTCGATCAGGTAGAAGGTGAAGAAGTTGGTGAACGAAGCGATGTAGACGAACTTGGCGAACATCAGCACGGCGATCACCGCCACAGCGCGCCACATGGCACCCTGCGACAGGCCGGGGGCCTGTTGGCCCGCGAAGGTCTTGAGCCGGGCCTGGCCGTGGCGCACGCTCCAGCCCGTCACCCGCAGTAGCACCAGAACAGCCAGCGCAGCGGCCAGCATGAACCAGGCAATGGCCGGTTGGCCGTGCGGAATGATAA  >  NC_002947/2699688‑2699975

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