Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,563,787 T→C 10.6% K25R (AAG→AGG) ‡ vgrG‑III ← type VI secretion system protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,563,7870TC10.6% 172.7 / 10.4 66K25R (AAG→AGG) ‡vgrG‑IIItype VI secretion system protein
Reads supporting (aligned to +/- strand):  ref base T (33/26);  new base C (3/4);  total (36/30)
Fisher's exact test for biased strand distribution p-value = 6.93e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.14e-01

GCCATGCACATGACGCAGTGGCCGGTCTGCGCTGTAAAGGGTGAATAGCGCAGGCTTGTCGAGGACGTGGCCGAAATCGATATCACTTTCGTAGCTCACAAGTTCCAGGTTCAGCACAAATGGCTGGTTGAGGGTTTCATCCAGCTTGAATGAAACCACTTCAAATTCCGTGCGCCCGGCGAGCGCCTGAAAACTGAAGCGAAAGTCGGACTGTTTCATGAGTTGCTCTACCCGCCATCATGTTTGCCGGCGAAAAACAACCACACCCTGGCTTTCTTCAAGATGAGAATATC  >  NC_002947/4563642‑4563934
                                                                                                                                                 |                                                                                                                                                   
gCCATGCACATGACGCAGTGGCCGGTCTGCGCTGTAAAGGGTGAATAGCGCAGGCTTGTCGAGGACGTGGCCGAAATCGATATCACTTTCGTAGCTCACAAGTTCCAGGTTCAGCACAAATGGCTGGTTGAGGGTTTCATCCAGCTtga                                                                                                                                                  <  3:130361/149‑1 (MQ=255)
gCCATGCACATGACGCAGTGGCCGGTCTGCGCTGTAAAGGGTGAATAGCGCAGGCTTGTCGAGGACGTGGCCGAAATCGATATCACTTTCGTAGCTCACAAGTTCCAGGTTCAGCACAAATGGCTGGTTGAGGGTTTCATCCAGCTtga                                                                                                                                                  <  8:129105/149‑1 (MQ=255)
  cATGCACATGACGCAGTGGCCGGTCTGCGCTGTAAAGGGTGAATAGCGCAGGCTTGTCGAGGACGTGGCCGAAATCGATATCACTTTCGTAGCTCACAAGTTCCAGGTTCAGCACAAATGGCTGGTTGAGGGTTTCATCCAGCTtgaat                                                                                                                                                <  6:138073/149‑1 (MQ=255)
       acaTGACGCAGTGGCCGGTCTGCGCTGTAAAGGGTGAATAGCGCAGGCTTGTCGAGGACGTGGCCGAAATCGATATCACTTTCGTAGCTCACAAGTTCCAGGTTCAGCACAAATGGCTGGTTGAGGGTTTCATTCAAGCTGGATGAAAc                                                                                                                                           <  2:38004/149‑1 (MQ=255)
         aTGACGCAGTGGCCGGTCTGCGCTGTAAAGGGTGAATAGCGCAGGCTTGTCGAGGACGTGGCCGAAATCGATATCACTTTCGTAGCTCACAAGTTCCAGGTTCAGCACAAATGGCTGGTTGAGGGTTTCATCCAGCTTGAATGAAACCa                                                                                                                                         >  7:105379/1‑149 (MQ=255)
          tGACGCAGTGGCCGGTCTGCGCTGTAAAGGGGGAATAGCGCAGGCTTGTCGAGGACGTGGCCGAAATCGATATCACTTTCGTAGCTCACAAGTGCCAGGTTCAGCACAAATGGCTGGTTGAGGGTTTCATTCAAGCTGGATGAAAccct                                                                                                                                        <  6:153695/149‑3 (MQ=255)
           gACGCAGTGGCCGGTCTGCGCTGTAAAGGGTGAATAGCGCAGGCTTGTCGAGGACGTGGCCGAAATCGATATCACTTTCGTAGCTCACAAGTTCCAGGTTCAGCACAAATGGCTGGTTGCGGGTTTCATCCAGCTTGAATGAAACCACt                                                                                                                                       >  5:233115/1‑149 (MQ=255)
            aCGCAGTGGCCGGTCTGCGCTGTAAAGGGTGAATAGCGCAGGCTTGTCGAGGACGTGGCCGAAATCGATATCACTTTCGTAGCTCACAAGTTCCAGGTTCAGCACAAATGGCTGGTTGAGGGTTTCATCCAGCTTGAATGAAACCACtt                                                                                                                                      >  3:139359/1‑149 (MQ=255)
             cGCAGTGGCCGGTCTGCGCTGTAAAGGGTGAATAGCGCAGGCTTGTCGAGGACGTGGCCGAAATCGATATCACTTTCGTAGCTCACAAGTTCCAGGTTCAGCACAAATGGCTGGTTGAGGGTTTTATTCAAGCTGGATGAAAccctcaa                                                                                                                                     <  6:69379/149‑6 (MQ=255)
                     ccGGTCTGCGCTGTAAAGGGTGAATAGCGCAGGCTTGTCGAGGACGTGGCCGAAATCGATATCACTTTCGTAGCTCACAAGTTCCAGGTTCAGCACAAATGGCTGGTTGAGGGTTTCATCCAGCTTGAATGAAACCACTTCAAATTCCg                                                                                                                             <  4:188486/149‑1 (MQ=255)
                      cGGTCTGCGCTGTAAAGGGTGAATAGCGCAGGCTTGTCGAGGACGTGGCCGAAATCGATATCACTTTCGTAGCTCACAAGTTCCAGGTTCAGCACAAATGGCTGGTTGAGGGTTTCATCCAGCTTGAATGAAACCACTTCAAATTCCGt                                                                                                                            <  5:202616/149‑1 (MQ=255)
                      cGGTCTGCGCTGTAAAGGGTGAATAGCGCAGGCTTGTCGAGGACGTGGCCGAAATCGATATCACTTTCGTAGCTCACAAGTTCCAGGTTCAGCACAAATGGCTGGTTGAGGGTTTCATCCAGCTTGAATGAAACCACTTCAAATTCCGt                                                                                                                            >  2:185791/1‑149 (MQ=255)
                      cGGTCTGCGCTGTAAAGCGTGAATAGCGCAGGCTTGTCGAGGACGTGGCCGAAATCGATATCACTTTCGTAGCTCACAAGTTCCAGGTTCAGCACAAATGGCTGCTTGAGGGTTTCATCCAGCTTGAATGAAACCACTTCAAATTCCGt                                                                                                                            >  6:165646/1‑149 (MQ=255)
                               cTGTAAAGGGTGAATAGCGCAGGCTTGTCGAGGACGTGGCCGAAATCGATATCACTTTCGTAGCTCACAAGTTCCAGGTTCAGCACAAATGGCTGGTTGAGGGTTTCATCCAGCTTGAATGAAACCACTTCAAATTCCGTGCGCCCGGc                                                                                                                   <  6:233115/149‑1 (MQ=255)
                                       ggTGAATAGCGCAGGCTTGTCGAGGACGTGGCCGAAATCGATATCACTTTCGTAGCTCACAAGTTCCAGGTTCAGCACAAATGGCTGGTTGAGGGTTTCATCCAGCTTGAATGAAACCACTTCAAATTCCGTGCGcc                                                                                                                       <  2:387290/137‑1 (MQ=255)
                                       ggTGAATAGCGCAGGCTTGTCGAGGACGTGGCCGAAATCGATATCACTTTCGTAGCTCACAAGTTCCAGGTTCAGCACAAATGGCTGGTTGAGGGTTTCATCCAGCTTGAATGAAACCACTTCAAATTCCGTGCGcc                                                                                                                       >  1:387290/1‑137 (MQ=255)
                                           aaTAGCGCATGCTTGTCGAGGACGTGGCCGAAATCGATATAACTTTCGTAGCTCAATAGTTCCAGGTTCAGCACAAATGGCAGGTTGAGGGTTTCATCCATCTTGAATGAAAACACTTCAAATTACGAGAGACCGGCGAGCGCCTCaa                                                                                                        >  8:179656/1‑148 (MQ=255)
                                           aaTAGCGCAGGCTTGTCGAGGACGTGGCCGAAATCGATATCACTTTCGTAGCTCACAAGTTCCAGGTTCAGCACAAATGGCTGGTTGAGGGTTTCATCCAGCTTGAATGAAACCACTTCAAATTCCGTGCGACCGGCGAGCTCCTGaa                                                                                                        >  8:139274/1‑148 (MQ=255)
                                               gcgcAGGCTTGTCGAGGACGTGGCCGAAATCGATATCACTTTCGTAGCTCACAAGTTCCAGGTTCAGCACAAATGGCTGGTTGAGGGTTTCATCCAGCTTGAATGAAACCACTTCAAATTCCGTGCGCCCGGCGAGCGCCTGAAAACt                                                                                                    >  2:403453/1‑148 (MQ=255)
                                                 gcAGGCTTGTCGAGGACGTGGCCGAAAGCGATATCACTTGCGTAGCTCACAAGTTCCAGGTTCAGCACAAATGGCTGGTTGAGGGTTTCATCCAGCTTGAATGAAACCACTTCAAATTCCGTGCGCCCGGCGAGCGCCTGAAAACTGaa                                                                                                 <  4:240744/149‑1 (MQ=255)
                                                 gcAGGCTTGTCGAGGACGGTGCCGAAATCGATATCACTTTCGTAGCTCACAAGTTCCAGGTTCAGCACAAATGGCTGGTTGAGGGTTTCATCCAGCTTGAATGAAACCACTTCAAATTCCGTGCGCCCGGCGAGCGCCTGAAAACTGaa                                                                                                 <  4:250117/149‑1 (MQ=255)
                                                  cAGGCTTGTCGAGGACGTGGCCGAAATCGATATCACTTTCGTAGCTCACAAGTTCCAGGTTCAGCACAAATGGCTGGTTGAGGGTTTCATCCAGCTTGAATGAAACCACTTCAAATTCCGTGCGCCCGGCGAGCGCCTGAAAACTGAAg                                                                                                >  1:403326/1‑149 (MQ=255)
                                                     gCTTGTCGAGGACGTGGCCGAAATCGATATCACTTTCGTAGCTCACAAGTTCCAGGTTCAGCACAAATGGCTGGTTGAGGGTTTCATCCAGCTTGAATGAAACCACTTCAAATTCCGTGCGCCCGGCGAGCGCCTGAAAACTGAAGCg                                                                                              >  5:345441/1‑148 (MQ=255)
                                                       ttGTCGAGGACGTGGCCGAAATCGATATCACTTTCGTAGCTCACAAGTTCCAGGTTCAGCACAAATGGCTGGTTGAGGGTTTCATCCAGCTTGAATGAAACCACTTCAAATTCCGTGCGCCCGGCGAGCGCCTGAAAACTGAAGCGa                                                                                             >  3:228050/1‑147 (MQ=255)
                                                            gAGGACGTGGCCGAAATCGATATCACTTTCGTAGCTCACAAGTTCCAGGTTCAGCACAAATGGCTGGTTGAGGGTTTCATCCGGCTTGAATGAAACCACTTCAAATTCCGTGCGCCCGGCGAGCGCCTGAAAACTGAAGCGAAAGTCg                                                                                       >  3:159257/1‑148 (MQ=255)
                                                             aGGACGTGGACGAAAACGATATAAATTTCTTAGCTCACAAGATCGAGGTTCAGCAGAAATGGCTGGTTGAGGGTTTCATCAAGATAGAAAGAAACCACTTCAAATTCCGTGCGACCGCCGAGCGCCTGAAAACTGAAGCTAAAGTCgg                                                                                      >  2:133208/1‑148 (MQ=255)
                                                                     gCCGAAATCGATATCACTTTCGTAGCTCACAAGTTCCAGGTTCAGCACAAATGGCTGGTTGAGGGTTTCATCCGGCTTGAATGAAACCACTTCAAATTCCGTGCGCCCGGCGAGCGCCTGAAAACTGAAGCGAAAGTCGGACTGTTTCa                                                                             <  4:159257/149‑1 (MQ=255)
                                                                     gCCGAAATCGATATCAATTTCTTAGCTCACAAGTTCCATGTGCAGCACAAATGGCTGGTTGAGTGTTTCATCCAGCTTGAATGAAACCACTTCAAATTCCGTGCGCCCGGCGAGCGCCTGAAAACTGAAGCGAAAGTCGGACTGTTTCa                                                                             <  7:179656/149‑1 (MQ=255)
                                                                        gAAATCGATATCACTTTCGTAGCTCACAAGTTCCAGGTTCAGCACAAATGGCTGGTTGAGGGTTTCATCCAGCTTGAATGAAACCACTTCAAATTCCGTGCGCCCGGCGAGCGCCTGAAAACTGAAGCGAAAGTCGGACTGTTTCATGa                                                                          <  8:105379/149‑1 (MQ=255)
                                                                        gAAATCGATATCACTTTCGTAGCTCACAAGTTCCAGGTTCAGCACAAATGGCTGGTTGAGGGTTTCATCCAGCTTGAATGAAACCACTTCAAATTCCGTGCGCCCGGCGAGCGCCTGAAAACTGAAGCGAAAGTCGGACTGTTTCATGa                                                                          <  2:337452/149‑1 (MQ=255)
                                                                               atatCACTTTCGTAGCTCACAAGTTCCAGGTTCAGCACAAATGGCTGGTTGAGGGTTTCATCCAGCTTGAATGAAACCACTTCAAATTCCGTGCGCCCGGCGAGCGCCTGAAAACTGAAGCGAAAGTCGGACTGTTTCATGAGTTGctc                                                                   >  3:331183/1‑149 (MQ=255)
                                                                               atatCACTTTCGTAGCTCACAAGTTCCAGGTTCAGCACAAATGGCTGGTTGAGGGTTTCATCCAGCTTGAATGAAACCACTTCAAATTCCGTGCGCCCGGCGAGCGCCTGAAAACTGAAGCGAAAGTCGGACTGTTTCATGAGTTGctc                                                                   >  1:59221/1‑149 (MQ=255)
                                                                                tatCACTTTCGTAGCTCACAAGTTCCAGGTTCAGCACAAATGGCTGGTTGAGGGTTTCATCCAGCTTGAATGAAACCACTTCAAATTCCGTGCGCCCGGCGAGCGCCTGAAAACTGAAGCGAAAGTCGGACTGTTTCATGAGTTGctc                                                                   >  6:159367/1‑148 (MQ=255)
                                                                                 atCACTTTCGTAGCTCACAAGTTCCAGGTTCAGCACAAATGGCTGGTTGAGGGTTTCATCCAGCTTGAATGAAACCACTTCAAATTCCGTGCGCCCGGCGAGCGCCTGAAAACTGAAGCGAAAGTCg                                                                                       <  6:154107/127‑1 (MQ=255)
                                                                                 atCACTTTCGTAGCTCACAAGTTCCAGGTTCAGCACAAATGGCTGGTTGAGGGTTTCATCCAGCTTGAATGAAACCACTTCAAATTCCGTGCGCCCGGCGAGCGCCTGAAAACTGAAGCGAAAGTCg                                                                                       >  5:154107/1‑127 (MQ=255)
                                                                                    aCTTTCGTAGCTCACAAGTTCCAGGTTCAGCACAAATGGCTGGTTGAGGGTTTCATCCAGCTTGAATGAAACCACTTCAAATTCCGTGCGCCCGGCGAGCGCCTGAAAACTGAAGCGAAAGTCGGACTGTTTCATGAGTTGCTCTAccc                                                              >  6:250073/1‑149 (MQ=255)
                                                                                           tAGCTCACAAGTTCCAGGTTCAGCACAAATGGCTGGTTGAGGGTTTCATCCAGCTTGAATGAAACCACTTCAAATTCAGTGCGCCCGGCGAGCGCCTGAAAACTGAAGCGAAAGTCGGACTGTTTCATGAGTTGCTCTACCCGCcatca                                                       >  5:228255/1‑149 (MQ=255)
                                                                                              cTCACAAGTTCCAGGTTCAGCACAAATGGCTGGTTGAGGGTTTCATCCAGCTTGAATGAAACCACTTCAAATTCCGTGCGCCCGGCGAGCGCCTGAAAACTGAAGCGAAAGTCGGACTGTTTCATGAGTTGCTCTACCCGCCATCATGt                                                    >  5:255621/1‑149 (MQ=255)
                                                                                               tCACAAGTTCCAGGTTCAGCACAAATGGCTGGTTGAGGGTTTCATCCAGCTTGAATGAAACCACTTCAAATTCCGTGCGCCCGGCGAGCGCCTGAAAACTGAAGCGAAAGTCGGACTGTTTCATGAGTTGCTCTACCCGCCATCATGtt                                                   >  7:17762/1‑149 (MQ=255)
                                                                                               tCACAAGTTCCAGGTTCAGCACAAATGGCTGGTTGAGGGTTTCATCCAGCTTGAATGAAACCACTTCAAATTCCGTGCGCCCGGCGAGCGCCTGAAAACTGAAGCGAAAGTCGGACTGTTTCATGAGTTGCTCTACCCGCCATCATGtt                                                   >  8:97322/1‑149 (MQ=255)
                                                                                               tCACAAGTTCCAGGTTCAGCACAAATCGATGGTTGAGGGTTTCATCCAGCTTGAATGAAACCACTTCAAATTCCGTGCGCCAGGCGAGCGCCTGAAAACTGAAGCGAAAGTCCGACTGTTTCATGAGTTGCTCTACCCGCCATCATGtt                                                   >  3:366200/1‑149 (MQ=255)
                                                                                                 acaAGTTCCAGGTTCAGCACAAATGGCTGGTTGAGGGTTTCATCCAGCTTGAATGAAACCACTTCAAATTCCGTGCGCCCGGCGAGCGCCTGAAAACTGAAGCGAAAGTCGGACTGTTTCATGAGTTGCTCTACCCGCCATCATGttt                                                  <  6:345441/148‑1 (MQ=255)
                                                                                                 acaAGTTCCAGGTTCAGCACAAATGGCTGGTTGAGGGTTTCATCCAGCTTGAATGAAACCACTTCAAATTCCGTGCGCCCGGCGAGCGCCTGAAAACTGAAGCGAAAGTCGGACTGTTTCATGAGTTGCTCTACCCGCCATCATGTTTg                                                 >  7:43546/1‑149 (MQ=255)
                                                                                                          aGGTTCAGCACAAATGGCTGGTTGAGGGTTTCATCCAGCTTGAATGAAACCACTTCAAATTCCGTGCGCCCGGCGAGCGCCTGAAAACTGAAGCGAAAGTCGGACTGTTTCATGAGTTGCTCTACCCGCCATCATGTTTGCCGGCGaa                                         >  1:382146/1‑148 (MQ=255)
                                                                                                           ggTTCAGCACAAATGGCTGGTTGAGGGTTTCATCCAGCTTGAATGAAACCACTTCAAATTCCGTGCGCCCGGCGAGCGCCTGAAAACTGAAGCGAAAGTCGGACTGTTTCATGAGTTGCTCTACCCGCCATCATGTTTGCCGGCGaaaa                                       <  1:133208/149‑1 (MQ=255)
                                                                                                                aGCACAAATGGCTGGTTGAGGGTTTCATCCAGCTTGAATGAAACCACTTCAAATTCCGTGCGCCCGGCGAGCGCCTGAAAACTGAAGCGAAAGTCGGACTGTTTCATGAGTTGCTCTACCCGCCATCATGTTTGCCGGCGAAAaacaac                                  <  4:331183/149‑1 (MQ=255)
                                                                                                                aGCACAAATGGCTGGTTGAGGGTTTCATCCAGCTTGAATGAAACCACTTCAAATTCCGTGCGCCCGGCGAGCGCCTGAAAACTGAAGCGAAAGTCGGACTGTTTCATGAGTTGCTCTACCCGCCATCATGTTTGCCGGCGAAAaacaac                                  <  1:185791/149‑1 (MQ=255)
                                                                                                                   acaAATGGCTGGTTGAGGGTTTCATCCAGCTTGAATGAAACCACTTCAAATTCCGTGCGCCCGGCGAGCGCCTGAAAACTGAAGCGAAAGTCGGACTGTTTCATGAGTTGCTCTACCCGCCATCATGTTTGCCGGCGAAAAACAACca                                >  6:194047/1‑148 (MQ=255)
                                                                                                                       aTGGCTGGTTGAGGGTTTCATCCAGCTTGAATGAAACCACTTCAAATTCCGTGCGCCCGGCGAGCGCCTGAAAACTGAAGCGAAAGTCGGACTGTTTCATGAGTTGCTCTACCCGCCATCATGTTTGCCGGCGAAAAACAACCACAccc                           <  7:97322/149‑1 (MQ=255)
                                                                                                                       aTGGCTGGTTGAGGGTTTCATCCAGCTTGAATGAAACCACTTCAAATTCCGTGCGCCCGGCGAGCGCCTGAAAACTGAAGCGAAAGTCGGACTGTTTCATGAGTTGCTCTACCCGCCATCATGTTTGCCGGCGAAAAACAACCACAccc                           <  6:255621/149‑1 (MQ=255)
                                                                                                                        tGGCTGGTTGAGGGTTTCATCCAGCTTGAATGTAACCACTTCTAATTCCGTGCGCCCGGCGAGCGCCTGAAAACTGAAGAGAAAGTCGGACTGTTTCATGAGTTGCTCTTCCCGCCATCATGTTTGCCGGCGAAAAACAACCACACCCt                          <  4:139359/149‑1 (MQ=255)
                                                                                                                          gCTGGTTGAGGGTTTCATCCAGCTTGAATGAAACCACTTCAAATTCCGTGCGCCCGGCGAGCGCCTGAAAACTGAAGCGAAAGTCGGACTGTTTCATGAGTTGCTCTACCCGCCATCATGTTTGCCGGCGAAAAACAACCACACCCTgg                        <  8:43546/149‑1 (MQ=255)
                                                                                                                             ggTTGAGGGTTTCATCCAGCTTGAATGAAACCACTTCAAATTCCGTGCGCCCGGCGAGCGCCTGAAAACTGAAGCGAAAGTCGGACTGTTTCATGAGTTGCTCTACCCGCCATCATGTTTGCCGGCGAAAAACAACCACACCCTGGCtt                     >  7:66681/1‑149 (MQ=255)
                                                                                                                              tgaagtggTTTCATTCAAGCTGGATGAAACCACTTCAAATTCCGTGCGCCCGGCGAGCGCCTGAAAACTGAAGCGAAAGTCGGACTGTTTCATGAGTTGCTCTACCCGCAATCATGTTTGCCGGCGAAAAACAACCACACCCTGGCtt                     >  4:343951/7‑148 (MQ=255)
                                                                                                                               ttGAGGGTTTCATCCAGCTTGAATGAAACCACTTCAAATTCCGTGCGCCCGGCGAGCGCCTGAAAACTGAAGCGAAAGTCGGACTGTTTCATGAGTTGCTCTACCCGCCATCATGTTTGCCGGCGAAAAACAACCACACCCTGGCTttc                   <  5:250073/149‑1 (MQ=255)
                                                                                                                               ttGAGGGTTTCATCCAGCTTGAATGAAACCACTTCAAATTCCGTGCGCCCGGCGAGCGCCTGAAAACTGAAGCGAAAGTCGGACTGTTTCATGAGTTGCTCTACCCGCCATCATGTTTGCCGGCGAAAAACAACCACACCCTGGCTttc                   >  4:318417/1‑149 (MQ=255)
                                                                                                                                  aGGGTTTCATCCAGCTTGAATGAAACCACTTCAAATTCCGTGCGCCCGGCGAGCGCCTGAAAACTGAAGCGAAAGTCGGACTGTTTCATGAGTTGCTCTACCCGCCATCATGTTTGCCGGCGAAAAACAACCACACCCTGGCTttctt                 <  5:194047/148‑1 (MQ=255)
                                                                                                                                  aGGGTTTCATCCAGCTTGAATGAAACCACTTCAAATTCCGTGCGCCCGGCGAGCGCCTGAAAACTGAAGCGAAAGTCGGACTGTTTCATGAGTTGCTCTACCCGCCATCATGTTTGCCGGCGAAAAACAACCACACCCTGGCTttctt                 >  3:103243/1‑148 (MQ=255)
                                                                                                                                   tggTTTCATTCAAGCTGGATGAAACCACTTCAAATTCCGTGCGCCCGGCGAGCGCCTGAAAACTGAAGCGAAAGTCGGACTGTTTCATGAGTTGCTCTACCCGCCATCATGTTTGCCGGCGAAAaaca                                    <  2:239988/127‑1 (MQ=255)
                                                                                                                                   tggTTTCATTCAAGCTGGATGAAACCACTTCAAATTCCGTGCGCCCGGCGAGCGCCTGAAAACTGAAGCGAAAGTCGGACTGTTTCATGAGTTGCTCTACCCGCCATCATGTTTGCCGGCGAAAaaca                                    >  1:239988/2‑128 (MQ=255)
                                                                                                                                     gTTTCATTCAAGCTGGATGAAACCACTTCAAATTCCGTGCGCCCGGCGAGCGCCTGAAAACTGAAGCGAAAGTCGGACTGTTTCATGAGTTGCTCTACCCGCCATCATGTTTGCCGGCGAAAAACAACCACACCCTGGCTttcttc                >  4:342948/1‑146 (MQ=255)
                                                                                                                                      tttCATCCAGCTTGAATGAAACCACTTCAAATTCCGTGCGCCCGGCGAGCGCCTGAAAACTGAAGCGAAAGTCGGACTGTTTCATGAGTTGCTCTACCCGCCATCATGTTTGCCGGCGAAAAACAACCACACCCTGGCTttcttc                >  2:77910/1‑145 (MQ=255)
                                                                                                                                          aTCCAGCTTGAATGAAACCACTTCAAATTCCGTGCGCCCGGCGAGCGCCTGAAAACTGAAGCGAAAGTCGGACTGTTTCATGAGTTGCTCTACCCGCCATCATGTTTGCCGGCGAAAAACAACCACACCCTGGCTTTCTTCAAGATga         <  2:59221/148‑1 (MQ=255)
                                                                                                                                           tCCAGCTTGAATGAAACCACTTCAAATTCCGTGCGCCCGGCGAGCGCCTGAAAACTGAAGCGAAAGTCGGACTGTTTCATGAGTTGCTCTACCCGCCATCATGTTTGCCGGCGAAAAACAACCACACCCTGGCTTTCTTCAAGATgag        >  8:97958/1‑148 (MQ=255)
                                                                                                                                                cTTGAATGAAACCACTTCAAATTCCGTGCGCCCGGCGAGCGCCTGAAAACTGAAGCGAAAGTCGGACTGTTTCATGAGTTGCTCTACCCGCCATCATGTTTGCCGGCGAAAAACAACCACACCCTGGCTTTCTTCAAGATGAGAATATc  >  4:68634/1‑149 (MQ=255)
                                                                                                                                                cTTGAATGAAACCACTTCAAATTCCGTGCGCCCGGCGAGCGCCTGAAAACTGAAGCGAAAGTCGGACTGTTTCATGAGTTGCTCTACCCGCCATCATGTTTGCCGGCGAAAAACAACCACACCCTGGCTTTCTTCAAGATGAGAATATc  <  2:403326/149‑1 (MQ=255)
                                                                                                                                                 ttGAATGAAACCACTTCAAATTCCGTGCGCCCGGCGAGCGCCTGAAAACTGAAGCGAAAGTCGGACTGTTTCATGAGTTGCTCTACCCGCCATCATGTTTGCCGGCGaaa                                        <  2:379187/110‑1 (MQ=255)
                                                                                                                                                 ttGAATGAAACCACTTCAAATTCCGTGCGCCCGGCGAGCGCCTGAAAACTGAAGCGAAAGTCGGACTGTTTCATGAGTTGCTCTACCCGCCATCATGTTTGCCGGCGaaa                                        >  1:379187/1‑110 (MQ=255)
                                                                                                                                                 |                                                                                                                                                   
GCCATGCACATGACGCAGTGGCCGGTCTGCGCTGTAAAGGGTGAATAGCGCAGGCTTGTCGAGGACGTGGCCGAAATCGATATCACTTTCGTAGCTCACAAGTTCCAGGTTCAGCACAAATGGCTGGTTGAGGGTTTCATCCAGCTTGAATGAAACCACTTCAAATTCCGTGCGCCCGGCGAGCGCCTGAAAACTGAAGCGAAAGTCGGACTGTTTCATGAGTTGCTCTACCCGCCATCATGTTTGCCGGCGAAAAACAACCACACCCTGGCTTTCTTCAAGATGAGAATATC  >  NC_002947/4563642‑4563934

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: