Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,507,607 G→C 10.3% A290P (GCG→CCG)  cycA → D‑alanine/beta‑alanine/D‑serine/glycine permease

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,507,6070GC10.3% 193.2 / 10.1 66A290P (GCG→CCG) cycAD‑alanine/beta‑alanine/D‑serine/glycine permease
Reads supporting (aligned to +/- strand):  ref base G (27/32);  new base C (5/2);  total (32/34)
Fisher's exact test for biased strand distribution p-value = 2.51e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.32e-01

GGCGATCAACTCGATCCCGATCCGTATCATCGTGTTCTATGTGCTGGCGCTGATCGCGATCATGGCCGTGACCCCGTGGCGCGATGTGGTGCCCGGAAAAAGCCCGTTCGTTGAGTTGTTCGTGCTGGCAGGCTTGCCGGCGGCGGCAAGCATCATCAACTTCGTGGTGCTGACCTCGGCGGCTTCGTCGGCCAACAGCGGCGTGTTCTCCACCAGCCGCATGCTCTACGGCCTGTCGCAGGAAGGTGACGCGCCCAAGGCTTTCGAGAAGCTGTCGAGCCGTTCGGTGC  >  NC_002947/5507465‑5507754
                                                                                                                                              |                                                                                                                                                    
ggCGATCAACTCGATCCCGATCCGGATCATCGTGGTCTATGTGCTGGCGCTGATCGCGATCATGGCCGTGACCCCGTGGCGCGATGTGGTGCCCGGAAAAAGCCCGTTCGTTGAGTTGTTCGTGCTGGCAGGCTTGCCGGCGGCGGCaa                                                                                                                                                <  7:197588/149‑1 (MQ=255)
  cGATCAACTCGATCCCGATCCGTATCATCGTGTTCTATGTGCTGGCGCTGATCGCGATCATGGCCGTGACCCCGTGGCGCGATGTGGTGCCCGGAAAAAGCCCGTTCGTTGAGTTGTTCGTGCTGGCAGGCTTGCCGGCGGCGGCAAGc                                                                                                                                              <  7:197120/149‑1 (MQ=255)
  cGATCAACTCGATCCCGATCCGTATCATCGTGTTCTATGTGCTGGCGCTGATCGCGATCATGGCCGTGACCCCGTGGCGCGATGTGGTGCCCGGAAAAAGCCCGTTCGTTGAGTTGTTCGTGCTGGCAGGCTTGCCGCCGCCGGCAAGc                                                                                                                                              <  3:96318/149‑1 (MQ=255)
     tCAACTCGATCCCGATCCGTATCATCGTGTTCTATGTGCTGGCGCTGATCGCGATCATGGCCGTGACCCCGTGGCGCGATGTGGTGCCCGGAAAAAGCCCGTTCGTTGAGTTGTTCGTGCTGGCAGGCTTGCCGGCGGCGGCAAGcat                                                                                                                                            >  3:36052/1‑148 (MQ=255)
      cAACTCGATCCCGATCCGTATCATCGTGTTCTATGTGCTGGCGCTGATCGCGATCATGGCCGTGACCCCGTGGCGCGATGTGGTGCCCGGAAAAAGCCCGTTCGTTGAGTTGTTCGTGCTGGCAGGCTTGCCGGCGGCGGCAAGcatca                                                                                                                                          >  5:96183/1‑149 (MQ=255)
       aaCTCGATCCCGATCCGTATCATCGTGTTCTATGTGCTGGCGCTGATCGCGATCATGGCCGTGACCCCGTGGCGCGATGTGGTGCCCGGAAAAAGCCCGTTCGTTGAGTTGTTCGTGCTGGCAGGCTTGCCGGCGGCGGCAAGcatcat                                                                                                                                         <  7:147864/149‑1 (MQ=255)
       aaCTCGATCCCGATCCGTATCATCGTGTTCTATGTGCTGGCGCTGATCGCGATCATGGCCGTGACCCCGTGGCGCGATGTGGTGCCCGGAAAAAGCCCGTTCGTTGAGTTGTTCGTGCTGGCAGGCTTGCCGGCGGCGGCAAGcatcat                                                                                                                                         <  6:53232/149‑1 (MQ=255)
       aaCTCGATCCCGATCCGTATCATCGTGTTCTATGTGCTGGCGCTGATCGCGATCATGGCCGTGACCCCGTGGCGCGATGTGGTGCCCGGAAAAAGCCCGTTCGTTGAGTTGTTCGTGCTGGCAGGCTTGCCGGCGGCGGCAAGcatca                                                                                                                                          >  3:69902/1‑148 (MQ=255)
        aCTCGATCCCGATCCGTATCATCGTGGTCTATGTGCTGGCGCTGATCGCGATCATGGCCGTGACCCCGTGGCGCGATGTGGTGCCCGGAAAAAGCCCGTGCGTTGAGTTGTGCGTGCTGGCAGGCTTGCCGCCGCCGGCAAGcctgcca                                                                                                                                        <  1:126076/149‑7 (MQ=255)
            gatcccgatccGTATCATCGTGTTCTATGTGCTGGCGCTGATCGCGATCATGGCCGTGACCCCGTGGCGCGATGTGGTGCCCGGAAAAAGCCCGTTCGTTGAGTTGTTCGTGCTGGCAGGCTTGCCGGCGGCGGCaa                                                                                                                                                >  7:117999/1‑137 (MQ=255)
            gatcccgatccGTATCATCGTGTTCTATGTGCTGGCGCTGATCGCGATCATGGCCGTGACCCCGTGGCGCGATGTGGTGCCCGGAAAAAGCCCGTTCGTTGAGTTGTTCGTGCTGGCAGGCTTGCCGGCGGCGGCaa                                                                                                                                                <  8:117999/137‑1 (MQ=255)
            gatcccgatccGTATCATCGTGTTCTATGTGCTGGCGCTGATCGCGATCATGGCCGTGACCCCGTGGCGCGATGTGGTGCCCGGAAAAAGCCCGTTCGTTGAGTTGTTCGTGCTGGCAGGCTTGCCGGCGGCGGCAAGCATCATCaa                                                                                                                                      <  2:66630/147‑1 (MQ=255)
            gatcccgatccGTATCATCGTGTTCTATGTGCTGGCGCTGATCGCGATCATGGCCGTGACCCCGTGGCGCGATGTGGTGCCCGGAAAAAGCCCGTTCGTTGAGTTGTTCGTGCTGGCAGGCTTGCCGGCGGCGGCAAGCATCATCAACt                                                                                                                                    <  4:58423/149‑1 (MQ=255)
                    tccGTATCATCGTGTTCTATGTGCTGGCGCTGATCGCGATCATGGCCGTGACCCCGTGGCGCGATGTGGTGCCCGGAAAAAGCCCGTTCGTTGAGTTGTTCGTGCTGGCAGGCTTGCCGGCGGCGGCAAGCATCATCAACTTCgtggtg                                                                                                                            >  5:69347/1‑149 (MQ=255)
                        tATCATCGTGTTCTATGTGCTGGCGCTGATCGCGATCATGGCCGTGACCCCGTGGCGCGATGTGGTGCCCGGAAAAAGCCCGTTCGTTGAGTTGTTCGTGCTGGCAGGCTTGCCGGCGGCGGCAAGCATCATCAACTTCGTGGTGCTGa                                                                                                                        <  7:130477/149‑1 (MQ=255)
                                 gtTCTATGTGCTGGCGCTGATCGCGATCATGGCCGTGACCCCGTGGCGCGATGTGGTGCCCGGAAAAAGCCCGTTCGTTGAGTTGTTCGTGCTGGCAGGCTTGCCGGCGGCGGCAAGCATCATCAACTTCGTGGTGCTGACCTcggcgg                                                                                                              >  4:39346/1‑149 (MQ=255)
                                    cTATGTGCTGGCGCTGATCGCGATCATGGCCGTGACCCCGTGGCGCGATGTGGTGCCCGGAAAAAGCCCGTTCGTTGAGTTGTTCGTGCTGGCAGGCTTGCCGGCGGCGGCAAGCATCATCAACTTCGTGGTGCTGACCTCGGCGGCtt                                                                                                           >  6:177554/1‑149 (MQ=255)
                                    cTATGTGCTGGCGCTGATCGCGATCATGGCCGTGACCCCGTGGCGCGATGTGGTGCCCGGAAAAAGCCCGTTCGTTGAGTTGTTCGTGCTGGCAGGCTTGCCGGCGGCGGCAAGCATCATCAACTTCGTGGTGCTGACCTCGGCGGCt                                                                                                            <  7:39993/148‑1 (MQ=255)
                                         tgCTGGCGCTGATCGCGATCATGGCCGTGACCCCGTGGCGCGATGTGGTGCCCGGAAAAAGCCCGTTCGTTGAGTTGTTCGTGCTGGCAGGCTTGCCGGCGGCGGCAAGCATCATCAACTTCGTGGTGCTGACCTCGGCGGCTtcgtcg                                                                                                      >  2:178693/1‑149 (MQ=255)
                                          gCTGGCGCTGATCGCGATCATGGCCGTGACCCCGTGGCGCGATGTGGTGCCCGGAAAAAGCCCGTTCGTTGAGTTGTTCGTGCTGGCAGGCTTGCCGGCGGCGGCAAGCATCATCAACTTCGTGGTGCTGACCTCGGCGGCTTCGTCgg                                                                                                     <  4:69902/149‑1 (MQ=255)
                                                    aTCGCGATCATGGCCGTGACCCCGTGGCGCGATGTGGTGCCCGGAAAAAGCCCGTTCGTTGAGTTGTTCGTGCTGGCAGGCTTGCCGGCGGCGGCAAGCATCATCAACTTCGTGGTGCTGACCTCGGCGGCTTCGTCGGCCAACAgcg                                                                                            <  5:177554/148‑1 (MQ=255)
                                                    aTCGCGATCATGGCCGTGACCCCGAGGCGCGATGTGGTGCCCGGAAAAAGCCCGTTCGTTGAGTTGTTCGTGCTGGCAGGCTTGCCGGCGGCGGCAAGCATCATCAACTTCGTGGTGCTGACCTCGGCGGCTTCGTCGGCCAACAgcgg                                                                                           >  2:92855/1‑149 (MQ=255)
                                                     tCGCGATCATGGCCGTGACCCCGTGGCGCGATGTGGTGCCCGGAAAAAGCCCGTTCGTTGAGTTGTTCGTGCTGGCAGGCTTGCCGGCGGCGGCAAGCATCATCAACTTCGTGGTGCTGAc                                                                                                                       <  6:104973/121‑1 (MQ=255)
                                                     tCGCGATCATGGCCGTGACCCCGTGGCGCGATGTGGTGCCCGGAAAAAGCCCGTTCGTTGAGTTGTTCGTGCTGGCAGGCTTGCCGGCGGCGGCAAGCATCATCAACTTCGTGGTGCTGAc                                                                                                                       >  5:104973/1‑121 (MQ=255)
                                                     tCGCGATCATGGCCGTGACCCCGTGGCGCGATGTGGTGCCCGGAAAAAGCCCGTTCGTTGAGTTGTTCGTGCTGGCAGGCTTGCCGGCGGCGGCAAGCATCATCAACTTCGTGGTGCTGACCTCGGCGGCTTCGTCGGCCAACAgcggc                                                                                          >  7:84023/1‑149 (MQ=255)
                                                     tCGCGATCATGGCCGTGACCCCGTGGCGCGATGTGGTGCCCGGAAAAAGCCCGTTCGTTGAGTTGTTCGTGCTGGCAGGCTTGCCGGCGGCGGCAAGCATCATCAACTTCGTGGTGCTGACCTCGGCGGCTTCGTCGGCCAACAgcggc                                                                                          >  3:299360/1‑149 (MQ=255)
                                                     gcgcgcgcATGGCCGTGCCCCCGGGGCGCGATGTGGGGCCCGGAAAAAGCCCGTTCGTTGCGTTGTTCGTGCTGGCAGGCTTGCCGGCGGCGGCAAGCATCATCAACTTCGTGGTGCTGACCTCGGCGGCTTCGTCGGCCAACAgcggc                                                                                          <  8:55796/142‑1 (MQ=255)
                                                         gATCATGGCCGTGACCCCGTGGCGCGATGTGGTGCCCGGAAAAAGCCCGTTCGTTGAGTTGTTCGTGCTGGCAGGCTTGCCGGCGGCGGCAAGCATCATCAACTTCGAGGTGCTGACCTCGGCGGCTTCGTCGGCCAACAGCGGCgtgt                                                                                      >  5:313241/1‑149 (MQ=255)
                                                                gCCGTGACCCCGTGGCGCGATGTGGTGCCCGGAAAAAGCCCGTTCGTTGAGTTGTTCGTGCTGGCAGGCTTGCCGGCGGCGGCAAGCATCATCAACTTCGTGGTGCTGACCTCGGCGGCTTCGTCGGCCAACAGCGGCGTGTTCTccac                                                                               <  6:69347/149‑1 (MQ=255)
                                                                      aCCCCGTGGCGCGATGTGGTGCCCGGAAAAAGCCCGTTCGTTGAGTTGTTCGTGCTGGCAGGCTTGCCGGCGGCGGCAAGCATCATCAACTTCGTGGTGCTGACCTCGGCGGCTTCGTCGGCCAACAGCGGCGTGTTCTCCACCAGCCg                                                                         >  8:327584/1‑149 (MQ=255)
                                                                       ccccGTGGCGCGATGTGGTGCCCGGAAAAAGCCCGTTCGTTGAGTTGTTCGTGCTGGCAGGCTTGCCGGCGGCGGCAAGCATCATCAACTTCGTGGTGCTGACCTCGGCGGCTTCGTCGGCCAACAGCGGCGTGTTCTCCACCAGCCGc                                                                        <  1:178693/149‑1 (MQ=255)
                                                                          cGTGGCGCGATGTGGTGCCCGGAAAAAGCCCGTTCGTTGAGTTGTTCGTGCTGGCAGGCTTGCCGGCGGCGGCAAGCATCATCAACTTCGTGGTGCTGACCTCGGCGGCTTCGTCGGCCAACAGCGGCGTGTTCTCCACCAGCCGCATg                                                                     <  7:31701/149‑1 (MQ=255)
                                                                          cGTGGCGCGATGTGGTGCCCGGAAAAAGCCCGTTCGTTGAGTTGTTCGTGCTGGCAGGCTTGCCGGCGGCGGCAAGCATCATCAACTTCGTGGTGCTGACCTCGGCGGCTTCGTCGGCCAACAGCGGCGTGTTCTCCACCAGCCGCATg                                                                     <  3:263950/149‑1 (MQ=255)
                                                                                gcgATGTGGTGCCCGGAAAAAGCCCGTTCGTTGAGTTGTTCGTGCTGGCAGGCTTGCCGGCGGCGGCAAGCATCATCAACTTCGTGGTGCTGACCTCGGCGGCTTCGTCGGCCAACAGCGGCGTGTTCTccac                                                                               >  2:190142/1‑133 (MQ=255)
                                                                                gcgATGTGGTGCCCGGAAAAAGCCCGTTCGTTGAGTTGTTCGTGCTGGCAGGCTTGCCGGCGGCGGCAAGCATCATCAACTTCGTGGTGCTGACCTCGGCGGCTTCGTCGGCCAACAGCGGCGTGTTCTccac                                                                               <  1:190142/133‑1 (MQ=255)
                                                                                    tgtgGTGCCCGGAAAAAGCCCGTTCGTTGAGTTGTTCGTGCTGGCAGGCTTGCCGGCGGCGGCAAGCATCATCAACTTCGTGGTGCTGACCTCGGCGGCTTCGTCGGCCAACAGCGGCGTGTTCTCCACCAGCCGCATGCTCTACGGcc                                                           <  7:327584/149‑1 (MQ=255)
                                                                                      tggtgCCCGGAAAAAGCCCGTTCGTTGAGTTGTTCGTGCTGGCAGGCTTGCCGGCGGCGGCAAGCATCATCAACTTCGTGGTGCTGACCTCGGCGGCTTCGTCGGCCAACAGCGGCGTGTTCTCCACCAGCCGCATGCTCTACGGCCt                                                          <  6:96183/148‑1 (MQ=255)
                                                                                       ggtgCCCGGAAAAAGCCCGTTCGTTGAGTTGTTCGTGCTGGCAGGCTTGCCGGCGGCGGCAAGCATCATCAACTTCGTGGTGCTGACCTCGGCGGCTTCGTCGGCCAACAGCGGCGTGTTCTCCACCAGCCGCATGCTCTACGGCCTGt                                                        >  3:196135/1‑149 (MQ=255)
                                                                                          gCCCGGAAAAAGCCCGTTCGTTGAGTTGTTCGTGCTGGCAGGCTTGCCGGCGGCGGCAAGCATCATCAACTTCGTGGTGCTGACCTCGGCGGCTTCGTCGGCCAACAGCGGCGTGTTCTCCACCAGCCGCATGCTCTACGGCCTGTCg                                                      >  7:56156/1‑148 (MQ=255)
                                                                                           cccGGAAAAAGCCCGTTCGTTGAGTTGTTCGTGCTGGCAGGCTTGCCGGCGGCGGCAAGCATCATCAACTTCGTGGTGCTGACCTCGGCGGCTTCGTCGGCCAACAGCGGCGTGTTCTCCACCAGCCGCATGCTCTACGGCCTGTCGc                                                     <  8:84023/148‑1 (MQ=255)
                                                                                             cGGAAAAAGCCCGTTCGTTGAGTTGTTCGTGCTGGCAGGCTTGCCGGCGGCGGCAAGCATCATCAACTTCGTGGTGCTGACCTCGGCGGCTTCGTCGGCCAACAGCGGCGTGTTCTCCACCAGCCGCATGCTCTACGGCCTGTCGCAg                                                   >  5:130040/1‑148 (MQ=255)
                                                                                              ggAAAAAGCCCGTTCGTTGAGTTGTTCGTGCTGGCAGGCTTGCCGGCGGCGGCAAGCATCATCAACTTCGTGGTGCTGACCTCGGCGGCTTCGTCGGCCAACAGCGGCGTGTTCTCCACCAGCCGCATGCTCTACGGCCTGTCGCAGGa                                                 <  4:196135/149‑1 (MQ=255)
                                                                                                aaaaaGCCCGTTCGTTGAGTTGTTCGTGCTGGCAGGCTTGCCGGCGGCGGCAAGCATCATCAACTTCGTGGTGCTGACCTCGGCGGCTTCGTCGGCCAACAGCGGCGTGTTCTCCACCAGCCGCATGCTCTACGGCCTGTCGCAGGAAg                                               <  7:92069/149‑1 (MQ=255)
                                                                                                 aaaaGCCCGTTCGTTGAGTTGTTCGTGCTGGCAGGCTTGCCGGCTGCGGCAAGCATCATCAACTTCGTGGTGCTGACCTCGGCGGCTTCGTCGGCCAACAGCGGCGTGTTCTCCACCAGCCGCATGCTCTACGGCCTGTCGCa                                                    <  5:255341/143‑1 (MQ=255)
                                                                                                 aaaaGCCCGTTCGTTGAGTTGTTCGTGCTGGCAGGCTTGCCGGCGCCGGCAAGCATCATCAACTTCGTGGTGCTGACCTCGGCGGCTTCGTCGGCCAACAGCGGAGTGTTCTCCACCAGACGCATGCTCTACGGCCTGGCGCa                                                    >  6:255341/1‑143 (MQ=255)
                                                                                                      cccGTTCGTTGAGTTGTTCGTGCTGGCAGGCTTGCCGGCGGCGGCAAGCATCATCAACTTCGTGGTGCTGACCTCGGCGGCTTCGTCGGCCAACAGCGGCGTGTTCTCCACCAGCCGCATGCTCTACCGCCTGTCGCAGGAAGGTGAcg                                         >  6:110012/1‑149 (MQ=255)
                                                                                                      cccGTTCGTTGAGTTGTTCGTGCTGGCAGGCTTGCCGGCGGCGGCAAGCATCATCAACTTCGTGGTGCTGACCTCGGCGGCTTCGTCGGCCAACAGCGGCGTGTTCTACACCAGCCGCATGCTCTACGGCCTGTCGCAGGAAGGTGAcg                                         >  5:239535/1‑149 (MQ=255)
                                                                                                       ccGTTCGTTGAGTTGTTCGTGCTGGCAGGCTTGCCGGCGGCGGCAAGCATCATCAACTTCGTGGTGCTGACCTCGGCGGCTTCGTCGGCCAACAGCGGCgtg                                                                                       >  7:4678/1‑102 (MQ=255)
                                                                                                       ccGTTCGTTGAGTTGTTCGTGCTGGCAGGCTTGCCGGCGGCGGCAAGCATCATCAACTTCGTGGTGCTGACCTCGGCGGCTTCGTCGGCCAACAGCGGCGTGTTCTCCACCAGCCGCATGCTCTACGGCCTGTCGCa                                                    <  4:308179/137‑1 (MQ=255)
                                                                                                       ccGTTCGTTGAGTTGTTCGTGCTGGCAGGCTTGCCGGCGGCGGCAAGCATCATCAACTTCGTGGTGCTGACCTCGGCGGCTTCGTCGGCCAACAGCGGCGTGTTCTCCACCAGCCGCATGCTCTACGGCCTGTCGCa                                                    >  3:308179/1‑137 (MQ=255)
                                                                                                       ccGTTCGTTGAGTTGTTCGTGCTGGCAGGCTTGCCGGCGGCGGCAAGCATCATCAACTTCGTGGTGCTGACCTCGGCGGCTTCGTCGGCCAACAGCGGCGTGTTCTCCACCAGCCGCATGCTCTACGGCCTGTCGCAGGAAGGTGAcgc                                        <  3:65992/149‑1 (MQ=255)
                                                                                                          ttcgttGAGTTGTTCGTGCTGGCAGGCTTGCCGGCGGCGGCAAGCATCATCAACTTCGTGGTGCTGACCTCGGCGGCTTCGTCGGCCAACAGCGGCGTGTTCTCCACCAGCCGCATGCTCTACGGCCTGTCGCAGGAAGGTGACGCGcc                                     <  4:299360/149‑1 (MQ=255)
                                                                                                                  gttgttCGTGCTGGCAGGCTTGCCGGCGGCGGCAAGCATCATCAACTTCGTGGTGCTGACCTCGGCGGCTTCGTCGGCCAACAGCGGCGTGTTCTCCACCAGCCGCATGCTCTACGGCCTGTCGCAGGAAGGTGACGCGCCCAAGGCtt                             >  6:135395/1‑149 (MQ=255)
                                                                                                                      ttCGTGCTGGCAGGCTTGCCGGCGGCGGCAAGCATCATCAACTTCGTGGTGCTGACCTCGGCGGCTTCGTCGGCCAACAGCGGCGTGTTCTCCACCAGCCGCATGCTCTACGGCCTGTCGCAGGAAGGTGACGCGCCCAAGGCTTTc                           >  2:64212/1‑147 (MQ=255)
                                                                                                                         gTGCTGGCCGGCTTGCCGGCGGCGGCAAGCATCATCAACTTCGTGGTGCTGACCTCGGCGGCTTCGTCGGCCAACAGCGGCGTGTTCTCCACCAGCCGCATGCTCTACGGCCTGTCGCAGGAAGGTGACGCGCCCAAGGCTTTCGAGaa                      <  5:70125/149‑1 (MQ=255)
                                                                                                                         gTGCTGGCAGGCTTGCCGGCGGCGGCAAGCATCATCAACTTCGTTGTGCTGACCTCGGCGGCTTCGTCGGCCAACAGCGTCGTGTTCTCCACCAGCCGCATGCTCTACGGCCTGTCGCAGGAAGGTGACGCGCCCAAGGCTTTCGAGaa                      <  7:208577/149‑1 (MQ=255)
                                                                                                                         gTGCTGGCAGGCTTGCCGGCGGCGGCAAGCATCATCAACTTCGTGGTGCTGACCTCGGCGGCTTCGTCGGCCAACAGCGGCGTGTTCTCCACGAGCCGCATGCTCTACGGCCTGTCGCAGGAAGGTGACGCGCCCAAGGCTTTCGAGaa                      <  3:39346/149‑1 (MQ=255)
                                                                                                                         gTGCTGGCAGGCTTGCCGGCGGCGGCAAGCATCATCAACTTCGTGGTGCTGACCTCGGCGGCTTCGTCGGCCAACAGCGGCGTGTTCTCCACCAGCCGCATGCTCTACGGCCTGTCGCAGGAAGGTGACGCGCCCAAGGCTTTCgag                        >  5:137498/1‑147 (MQ=255)
                                                                                                                             gatgatgCTTGCCGCCGCCGGCAAGCCTCATCAACTTCGTGGTGCTGACCTCGGCGGCTTCGTCGGCCAACAGCGGCGTGTTCTCCACCAGCCGCATGCTCTACGGCCTGTCGCAGGAAGGTGACGCGCCCAAGGCTTTCGAGAAGCTg                  >  2:306577/7‑149 (MQ=255)
                                                                                                                               gCAGGCTTGCCGGCGGCGGCAAGCATCATCAACTTCGTGGTGCTGACCTCGGCGGCTTCGTCGGCCAACAGCGGCGTGTTCTCCACCAGCCGCATGCTCTACGGCCTGTCGCAGGAAGGTGACGCGCCCAAGGCTTTCGAGAAGCTGTc                >  3:22044/1‑149 (MQ=255)
                                                                                                                                 aGGCTTGCCGGCGGCGGCAAGCATCATCAACTTCGAGGTGCTGACATAGGGCGGCTTCGTCGGCCAACAGCGGCGTGTTCTCCACCAGCCGCATGCTCTACGGCCTGTCGCAGGAAGGTGACGCGCCCAAGGCTTTCGAGAAGCTGTc                >  1:266913/1‑148 (MQ=255)
                                                                                                                                  tgCTTGCCGCCGCCGGCAAGCATCATCAACTTCGTGGTGCTGACCTCGGCGGCTTCGTCGGCCAACAGCGGCGTGTTCTCCACCAGCCGCATGCTCTACGGCCTGTCGCAGGAAGGTGACGCGCCCAAGGCTTTCGAGAAGCTGTCGAg             >  7:232439/2‑149 (MQ=255)
                                                                                                                                  tgCTTGCCGCCGCCGGCAAGCATCATCAACTTCGTGGTGCTGACCTCGGCGGCTTCGTCGGCCAACAGCGGCGTGTTCTCCACCAGCCGCATGCTCTACGGCCTGTCGCAGGAAGGTGACGCGCCCAAGGCTTTCGAGAAGCTGTCGAg             >  4:301253/2‑149 (MQ=255)
                                                                                                                                  ggCTTGCCGGCGGCGGCAAGCATCATCAACTTCGTGGTGCTGACCTCGGCGGCTTCGTCGGCCAACAGCGGCGTGTTCTCCACCAGCCGCATGCTCTACGGCCTGTCGCAGGAAGGTGACGCGCCCAAGGCTTTCGAGAAGCTGTCGAg             <  5:110012/149‑1 (MQ=255)
                                                                                                                                   gCTTGCCGCCGCCGGCAAGCATCATCAACTTCGTGGTGCTGACCTCGGCGGCTTCGTCGGCCAACAGCGGCGTGTTCTCCACCAGCCGCATGCTCTACGGCCTGTCGCAGGAAGGTGACGCGCCCAAGGCTTTCGAGTAGCTGTCGAGc            >  7:337219/1‑149 (MQ=255)
                                                                                                                                        ccGGCGGCGGCAAGCATCATCAACTTCGTGGTGCTGACCTCGGCGGCTTCGTCGGCCAACAGCGGCGTGTTCTCCACCAGCCGCATGCTCTACGGCCTGTCGCAGGAAGGTGACGCGCCCAAGGCTTTCTAGAAGCTGTCGAGCCGtt        <  8:56156/148‑1 (MQ=255)
                                                                                                                                           gcggcggcAAGCATCATCAACTTCGTGGTGCTGACCTCGGCGGCTTCGTCGGCCAACAGCGGCGTGTTCTCCACCAGCCGCATGCTCTACGGCCTGTCGCAGGAAGGTGACGCGCCCAAGGCTTTCGAGAAGCTGTCGAGCCGTTCGGt    <  6:130040/149‑1 (MQ=255)
                                                                                                                                             ggcggcAAGCATCATCAACTTCGTGGTGCTGACCTCGGCGGCTTCGTCGGCCAACAGCGGCGTGTTCTCCACCAGCCGCATGCTCTACGGCCTGTCGCAGGAAGGTGACGCGCCCAAGGCTTTCGAGAAGCTGTCGAGCCGTTCGGTGc  >  5:103720/1‑149 (MQ=255)
                                                                                                                                              |                                                                                                                                                    
GGCGATCAACTCGATCCCGATCCGTATCATCGTGTTCTATGTGCTGGCGCTGATCGCGATCATGGCCGTGACCCCGTGGCGCGATGTGGTGCCCGGAAAAAGCCCGTTCGTTGAGTTGTTCGTGCTGGCAGGCTTGCCGGCGGCGGCAAGCATCATCAACTTCGTGGTGCTGACCTCGGCGGCTTCGTCGGCCAACAGCGGCGTGTTCTCCACCAGCCGCATGCTCTACGGCCTGTCGCAGGAAGGTGACGCGCCCAAGGCTTTCGAGAAGCTGTCGAGCCGTTCGGTGC  >  NC_002947/5507465‑5507754

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: