Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,563,790 A→G 10.8% F24S (TTC→TCC)  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,7900AG10.8% 169.1 / 10.3 65F24S (TTC→TCC) vgrG‑IIItype VI secretion system protein
Reads supporting (aligned to +/- strand):  ref base A (34/24);  new base G (3/4);  total (37/28)
Fisher's exact test for biased strand distribution p-value = 4.53e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.90e-01

GCCATGCACATGACGCAGTGGCCGGTCTGCGCTGTAAAGGGTGAATAGCGCAGGCTTGTCGAGGACGTGGCCGAAATCGATATCACTTTCGTAGCTCACAAGTTCCAGGTTCAGCACAAATGGCTGGTTGAGGGTTTCATCCAGCTTGAATGAAACCACTTCAAATTCCGTGCGCCCGGCGAGCGCCTGAAAACTGAAGCGAAAGTCGGACTGTTTCATGAGTTGCTCTACCCGCCATCATGTTTGCCGGCGAAAAACAACCACACCCTGGCTTTCTTCAAGATGAGAATATCGAA  >  NC_002947/4563642‑4563937
                                                                                                                                                    |                                                                                                                                                   
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                                                                                                                               >  2:185791/1‑149 (MQ=255)
                      cGGTCTGCGCTGTAAAGGGTGAATAGCGCAGGCTTGTCGAGGACGTGGCCGAAATCGATATCACTTTCGTAGCTCACAAGTTCCAGGTTCAGCACAAATGGCTGGTTGAGGGTTTCATCCAGCTTGAATGAAACCACTTCAAATTCCGt                                                                                                                               <  5:202616/149‑1 (MQ=255)
                      cGGTCTGCGCTGTAAAGCGTGAATAGCGCAGGCTTGTCGAGGACGTGGCCGAAATCGATATCACTTTCGTAGCTCACAAGTTCCAGGTTCAGCACAAATGGCTGCTTGAGGGTTTCATCCAGCTTGAATGAAACCACTTCAAATTCCGt                                                                                                                               >  6:165646/1‑149 (MQ=255)
                               cTGTAAAGGGTGAATAGCGCAGGCTTGTCGAGGACGTGGCCGAAATCGATATCACTTTCGTAGCTCACAAGTTCCAGGTTCAGCACAAATGGCTGGTTGAGGGTTTCATCCAGCTTGAATGAAACCACTTCAAATTCCGTGCGCCCGGc                                                                                                                      <  6:233115/149‑1 (MQ=255)
                                       ggTGAATAGCGCAGGCTTGTCGAGGACGTGGCCGAAATCGATATCACTTTCGTAGCTCACAAGTTCCAGGTTCAGCACAAATGGCTGGTTGAGGGTTTCATCCAGCTTGAATGAAACCACTTCAAATTCCGTGCGcc                                                                                                                          >  1:387290/1‑137 (MQ=255)
                                       ggTGAATAGCGCAGGCTTGTCGAGGACGTGGCCGAAATCGATATCACTTTCGTAGCTCACAAGTTCCAGGTTCAGCACAAATGGCTGGTTGAGGGTTTCATCCAGCTTGAATGAAACCACTTCAAATTCCGTGCGcc                                                                                                                          <  2:387290/137‑1 (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                                                      >  8:97322/1‑149 (MQ=255)
                                                                                               tCACAAGTTCCAGGTTCAGCACAAATGGCTGGTTGAGGGTTTCATCCAGCTTGAATGAAACCACTTCAAATTCCGTGCGCCCGGCGAGCGCCTGAAAACTGAAGCGAAAGTCGGACTGTTTCATGAGTTGCTCTACCCGCCATCATGtt                                                      >  7:17762/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                                     <  1:185791/149‑1 (MQ=255)
                                                                                                                aGCACAAATGGCTGGTTGAGGGTTTCATCCAGCTTGAATGAAACCACTTCAAATTCCGTGCGCCCGGCGAGCGCCTGAAAACTGAAGCGAAAGTCGGACTGTTTCATGAGTTGCTCTACCCGCCATCATGTTTGCCGGCGAAAaacaac                                     <  4:331183/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)
                                                                                                                                                   gaatgaaACCACTTCAAATTCCGTGCGCCCGGCGAGCGCCTGAAAACTGAAGCGAAAGTCGGACTGTTTCATGAGTTGCTCTACCCGCCATCATGTTTGCCGGCGAAAAACAACCACACCCTGGCTTTCTTCAAGATGAGAATATCGa   <  1:403453/148‑1 (MQ=255)
                                                                                                                                                    aatgaaACCACTTCAAATTCCGTGCGCCCGGCGAGCGCCTGAAAACTGAAGCGAAAGTCGGACTGTTTCATGAGTTGCTCTACCCGCCATCATGTTTGCCGGCGAAAAACAACCACACCCTGGCTTTCTTCAAGATGAGAATATCGaa  <  4:228050/148‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                   
GCCATGCACATGACGCAGTGGCCGGTCTGCGCTGTAAAGGGTGAATAGCGCAGGCTTGTCGAGGACGTGGCCGAAATCGATATCACTTTCGTAGCTCACAAGTTCCAGGTTCAGCACAAATGGCTGGTTGAGGGTTTCATCCAGCTTGAATGAAACCACTTCAAATTCCGTGCGCCCGGCGAGCGCCTGAAAACTGAAGCGAAAGTCGGACTGTTTCATGAGTTGCTCTACCCGCCATCATGTTTGCCGGCGAAAAACAACCACACCCTGGCTTTCTTCAAGATGAGAATATCGAA  >  NC_002947/4563642‑4563937

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: