Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,563,782 C→T 12.1% D27N (GAT→AAT)  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,7820CT12.1% 167.0 / 11.8 66D27N (GAT→AAT) vgrG‑IIItype VI secretion system protein
Reads supporting (aligned to +/- strand):  ref base C (32/26);  new base T (4/4);  total (36/30)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.77e-02

CTAGGCCATGCACATGACGCAGTGGCCGGTCTGCGCTGTAAAGGGTGAATAGCGCAGGCTTGTCGAGGACGTGGCCGAAATCGATATCACTTTCGTAGCTCACAAGTTCCAGGTTCAGCACAAATGGCTGGTTGAGGGTTTCATCCAGCTTGAATGAAACCACTTCAAATTCCGTGCGCCCGGCGAGCGCCTGAAAACTGAAGCGAAAGTCGGACTGTTTCATGAGTTGCTCTACCCGCCATCATGTTTGCCGGCGAAAAACAACCACACCCTGGCTTTCTTCAAGATGAG  >  NC_002947/4563638‑4563928
                                                                                                                                                |                                                                                                                                                  
cTAGGCCATGCACATGACGCAGTGGCCGGTCTGCGCTGTAAAGGGTGAATAGCGCAGGCTTGTCGAGGACGTGGCCGAAATCGATATCACTTTCGTAGCTCACAAGTTCCAGGTTCAGCACAAATGGCTGGTTGAGGGTTTCATCCAGc                                                                                                                                                <  4:16294/149‑1 (MQ=255)
 tAGGCCATGCACATGACGCAGTGGCCGGTCTGCGCTGTAAAGGGTGAATAGCGCAGGCTTGTCGAGGACGTGGCCGAAATCGATATCACTTTCGTAGCTCACAAGTTCCAGGTTCAGCACAAATGGCTGGTTGAGGGTTTTATTCaagc                                                                                                                                               >  5:69379/1‑146 (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                                                                                                                      >  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                                                                                                                 <  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                                                             >  1:59221/1‑149 (MQ=255)
                                                                                   atatCACTTTCGTAGCTCACAAGTTCCAGGTTCAGCACAAATGGCTGGTTGAGGGTTTCATCCAGCTTGAATGAAACCACTTCAAATTCCGTGCGCCCGGCGAGCGCCTGAAAACTGAAGCGAAAGTCGGACTGTTTCATGAGTTGctc                                                             >  3:331183/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                            <  1:185791/149‑1 (MQ=255)
                                                                                                                    aGCACAAATGGCTGGTTGAGGGTTTCATCCAGCTTGAATGAAACCACTTCAAATTCCGTGCGCCCGGCGAGCGCCTGAAAACTGAAGCGAAAGTCGGACTGTTTCATGAGTTGCTCTACCCGCCATCATGTTTGCCGGCGAAAaacaac                            <  4:331183/149‑1 (MQ=255)
                                                                                                                       acaAATGGCTGGTTGAGGGTTTCATCCAGCTTGAATGAAACCACTTCAAATTCCGTGCGCCCGGCGAGCGCCTGAAAACTGAAGCGAAAGTCGGACTGTTTCATGAGTTGCTCTACCCGCCATCATGTTTGCCGGCGAAAAACAACca                          >  6:194047/1‑148 (MQ=255)
                                                                                                                           aTGGCTGGTTGAGGGTTTCATCCAGCTTGAATGAAACCACTTCAAATTCCGTGCGCCCGGCGAGCGCCTGAAAACTGAAGCGAAAGTCGGACTGTTTCATGAGTTGCTCTACCCGCCATCATGTTTGCCGGCGAAAAACAACCACAccc                     <  6:255621/149‑1 (MQ=255)
                                                                                                                           aTGGCTGGTTGAGGGTTTCATCCAGCTTGAATGAAACCACTTCAAATTCCGTGCGCCCGGCGAGCGCCTGAAAACTGAAGCGAAAGTCGGACTGTTTCATGAGTTGCTCTACCCGCCATCATGTTTGCCGGCGAAAAACAACCACAccc                     <  7:97322/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           >  3:103243/1‑148 (MQ=255)
                                                                                                                                      aGGGTTTCATCCAGCTTGAATGAAACCACTTCAAATTCCGTGCGCCCGGCGAGCGCCTGAAAACTGAAGCGAAAGTCGGACTGTTTCATGAGTTGCTCTACCCGCCATCATGTTTGCCGGCGAAAAACAACCACACCCTGGCTttctt           <  5:194047/148‑1 (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)
                                                                                                                                                |                                                                                                                                                  
CTAGGCCATGCACATGACGCAGTGGCCGGTCTGCGCTGTAAAGGGTGAATAGCGCAGGCTTGTCGAGGACGTGGCCGAAATCGATATCACTTTCGTAGCTCACAAGTTCCAGGTTCAGCACAAATGGCTGGTTGAGGGTTTCATCCAGCTTGAATGAAACCACTTCAAATTCCGTGCGCCCGGCGAGCGCCTGAAAACTGAAGCGAAAGTCGGACTGTTTCATGAGTTGCTCTACCCGCCATCATGTTTGCCGGCGAAAAACAACCACACCCTGGCTTTCTTCAAGATGAG  >  NC_002947/4563638‑4563928

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: