Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,563,786 C→G 10.9% K25N (AAG→AAC) ‡ 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,7860CG10.9% 177.2 / 11.3 64K25N (AAG→AAC) ‡vgrG‑IIItype VI secretion system protein
Reads supporting (aligned to +/- strand):  ref base C (31/25);  new base G (3/4);  total (35/29)
Fisher's exact test for biased strand distribution p-value = 6.94e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.38e-01

CTAGGCCATGCACATGACGCAGTGGCCGGTCTGCGCTGTAAAGGGTGAATAGCGCAGGCTTGTCGAGGACGTGGCCGAAATCGATATCACTTTCGTAGCTCACAAGTTCCAGGTTCAGCACAAATGGCTGGTTGAGGGTTTCATCCAGCTTGAATGAAACCACTTCAAATTCCGTGCGCCCGGCGAGCGCCTGAAAACTGAAGCGAAAGTCGGACTGTTTCATGAGTTGCTCTACCCGCCATCATGTTTGCCGGCGAAAAACAACCACACCCTGGCTTTCTTCAAGATGAGAATATC  >  NC_002947/4563638‑4563934
                                                                                                                                                    |                                                                                                                                                    
cTAGGCCATGCACATGACGCAGTGGCCGGTCTGCGCTGTAAAGGGTGAATAGCGCAGGCTTGTCGAGGACGTGGCCGAAATCGATATCACTTTCGTAGCTCACAAGTTCCAGGTTCAGCACAAATGGCTGGTTGAGGGTTTCATCCAGc                                                                                                                                                      <  4:16294/149‑1 (MQ=255)
    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                                                                                       >  5:154107/1‑127 (MQ=255)
                                                                                     atCACTTTCGTAGCTCACAAGTTCCAGGTTCAGCACAAATGGCTGGTTGAGGGTTTCATCCAGCTTGAATGAAACCACTTCAAATTCCGTGCGCCCGGCGAGCGCCTGAAAACTGAAGCGAAAGTCg                                                                                       <  6:154107/127‑1 (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)
                                                                                                                                                    |                                                                                                                                                    
CTAGGCCATGCACATGACGCAGTGGCCGGTCTGCGCTGTAAAGGGTGAATAGCGCAGGCTTGTCGAGGACGTGGCCGAAATCGATATCACTTTCGTAGCTCACAAGTTCCAGGTTCAGCACAAATGGCTGGTTGAGGGTTTCATCCAGCTTGAATGAAACCACTTCAAATTCCGTGCGCCCGGCGAGCGCCTGAAAACTGAAGCGAAAGTCGGACTGTTTCATGAGTTGCTCTACCCGCCATCATGTTTGCCGGCGAAAAACAACCACACCCTGGCTTTCTTCAAGATGAGAATATC  >  NC_002947/4563638‑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: ‑