Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,744,607 G→T 6.8% G73G (GGC→GGA PP_4198 ← OmpA family protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,744,6070GT 6.8% 212.7 / 3.4 68G73G (GGC→GGAPP_4198OmpA family protein
Reads supporting (aligned to +/- strand):  ref base G (33/30);  new base T (2/3);  total (35/33)
Fisher's exact test for biased strand distribution p-value = 6.68e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GGCCTCACGCAAGGCCAGCTGCAGGCGCTGGCGCTCCAGGTCGAGCTTGGCCTGGCGCTCCTGGTTCAGCGCCAGCTTGGCATGCTCGCGGGCAATCTCGCTGTAGCGCTGGCTGAGGTAAGCATAATGCCGCACATCGGAACCGGTGCCGATGTAACTGGACAGGCGCTCGGCGCGAGCCAGCGACTCACCCGCCCGAATCACGTCACGCGGCGCGCTGCGCAGCACATCGGAATCATCTTTGACCCCCTGGAAGGCAACGGTGGCCTCATCCAGCGCGGCACTGCTGCGCTGA  >  NC_002947/4744460‑4744754
                                                                                                                                                    |                                                                                                                                                    
ggCCTCACGCAAGGCCAGCTGCAGGCGCTGGCGCTCCAGGTCGAGCTTGGCCTGGCGCTCCTGGTTCAGCGCCAGCTTGGCATGCTCGCGGGCAATCTCGCTGTAGCGCTGGCTGAGGTAAGCATAATGCCGCACATCGGAACCGGTGc                                                                                                                                                     >  2:176130/1‑149 (MQ=255)
 gCCTCACGCAAGGCCAGCTGCAGGCGCTGGCGCTCCAGGTCGAGCTTGGCCTGGCGCTCCTGGTTCAGCGCCAGCTTGGCATGCTCGCGGGCAATCTCGCTGTAGCGCTGGCTGAGGTAAGCATAATGCCGCACATCGGCACCGGTTcc                                                                                                                                                    <  4:145629/149‑1 (MQ=255)
      aCGCAAGGCCAGCTGCAGGCGCTGGCGCTCCAGGTCGAGCTTGGCCTGGCGCTCCTGGTTCAGCGCCAGCTTGGCATGCTCGCGGGCAATCTCGCTGTAGCGCTGGCTGAGGTAAGCATAATGCCGCACATCGGAACCGGTGCCGAt                                                                                                                                                 <  8:356324/147‑1 (MQ=255)
      aCGCAAGGCCAGCTGCAGGCGCTGGCGCTCCAGGTCGAGCTTGGCCTGGCGCTCCTGGTTCAGCGCCAGCTTGGCATGCTCGCGGGCAATCTCGCTGTAGCGCTGGCTGAGGTAAGCATAATGCCGCACATCGGAACCGGTGCCGATg                                                                                                                                                <  2:1993/148‑1 (MQ=255)
        gCAAGGCCAGCTGCAGGCGCTGGCGCTCCAGGTCGAGCTTGGCCTGGCGCTCCTGGTTCAGCGCCAGCTTGGCATGCTCGCGGGCAATCTCGCTGTAGCGCTGGCGGAGGTAAGCCTAATGCCGCACATCGGCAGCGGTGGCGAGGtcg                                                                                                                                             >  4:79612/1‑147 (MQ=255)
           aGGCCAGCTGCAGGCGCTGGCGCTCCAGGTCGAGCTTGGCCTGGCGCTCCTGGTTCAGCGCCAGCTTGGCATGCTCGCGGGCAATCTCGCTGTAGCGCTGGCTGAGGTAAGCATAATGCCGCACATCGGAACCGGTGCCGATGTAACTg                                                                                                                                          <  1:206851/149‑1 (MQ=255)
            ggCCAGCTGCAGGCGCTGGCGCTCCAGGTCGAGCTTGGCCTGGCGCTCCTGGTTCAGCGCCAGCTTGGCATGCTCGCGGGCAATCTCGCTGTAGCGCTGGCTGAGGTAAGCATAATGCCGCACATCGGAACCGGTGCCGATGTAACTg                                                                                                                                          >  2:71668/1‑148 (MQ=255)
              ccAGCTGCAGGCGCTGGCGCTCCAGGTCGAGCTTGGCCTGGCGCTCCTGGTTCAGCGCCAGCTTGGCATGCTCGCGGGCAATCTCGCTGTAGCGCTGGCTGAGGTAAGCATAATGCCGCACATCGGAACCGGTGCCGATGTAACTGGa                                                                                                                                        >  2:49249/1‑148 (MQ=255)
               cAGCTGCAGGCGCTGGCGCTCCAGGTCGAGCTTGGCCTGGCGCTCCTGGTTCAGCGCCAGCTTGGCATGCTCGCGGGCAATCTCGCTGTAGCGCTGGCTGAGGTAAGCATAATGCCGCACATCGGAACCGGTGCCGATGTAACTGGa                                                                                                                                        >  7:205294/1‑147 (MQ=255)
                aGCTGCAGGCGCTGGCGCTCCAGGTCGAGCTTGGCCTGGCGCTCCTGGTTCAGCGCCAGCTTGGCATGCTCGCGGGCAATCTCGCTGTAGCGCTGGCTGAGGTAAGCATAATGCCGCACATCGGAACCGGTGCCGATGTAACTGGACAg                                                                                                                                     <  4:173435/149‑1 (MQ=255)
                 gCTGCAGGCGCTGGCGCTCCAGGTCGAGCTTGGCCTGGCGCTCCTGGTTCAGCGCCAGCTTGGCATGCTCGCGGGCAATCTCGCTGTAGCGCTGGCTGAGGTAAGCATAATGCCGCACATCGGAACCGGTGCCGATGTAACTGGACAgg                                                                                                                                    <  2:306351/149‑1 (MQ=255)
                     cAGGCGCTGGCGCTCCAGGTCGAGCTTGGCCTGGCGCTCCTGGTTCAGCGCCAGCTTGGCATGCTCGCGGGCAATCTCGCTGTAGCGCTGGCTGAGGTAAGCATAATGCCGCACATCGGAACCGGTGCCGATGTAACTGGACAGgcg                                                                                                                                  <  6:4457/147‑1 (MQ=255)
                        gcgctggcgctCCAGGTCGAGCTTGGCCTGGCGCTCCTGGTTCAGCGCCAGCTTGGCATGCTCGCGGGCAATCTCGCTGTAGCGCTGGCTGAGGTAAGCATAATGCCGCACATCGGAACCGGTGCCGATGTAACTGGACAGGCGCTCgg                                                                                                                            <  7:14380/149‑1 (MQ=255)
                           ctggcgctCCAGGTCGAGCTTGGCCTGGCGCTCCTGGTTCAGCGCCAGCTTGGCATGCTCGCGGGCAATCTCGCTGTAGCGCTGGCTGAGGTAAGCATAATGCCGCACATCGGAACCGGTGCCGATGTAACTGGACAGGCGCTCGgcgc                                                                                                                         <  2:120297/149‑1 (MQ=255)
                              gcgctCCAGGTCGAGCTTGGCCTGGCGCTCCTGGTTCAGCGCCAGCTTGGCATGCTCGCGGGCAATCTCGCTGTAGCGCTGGCTGAGGTAAGCATAATGCCGCACATCGGAACCGGTGCCGATGTAACTGGACAGGCGCTCGGCGCGa                                                                                                                       <  8:205294/148‑1 (MQ=255)
                                    cAGGTCGAGCTTGGCCTGGCGCTCCTGGTTCAGCGCCAGCTTGGCATGCTCGCGGGCAATCTCGCTGTAGCGCTGGCTGAGGTAAGCATAATGCCGCACATCGGAACCGGTGCCGATGTAACTGGACAGGCGCTCGGCGCGAGCCAGc                                                                                                                 <  1:176130/148‑1 (MQ=255)
                                      ggTCGAGCTTGGCCTGGCGCTCCTGGTTCAGCGCCAGCTTGGCATGCTCGCGGGCAATCTCGCTGTAGCGCTGGCTGAGGTAAGCATAATGCCGCACATCGGAACCGGTGCCGATGTAACTGGACAGGCGCTCGGCGCGAGCCAGCGAc                                                                                                              >  3:297615/1‑149 (MQ=255)
                                        tCGAGCTTGGCCTGGCGCTCCTGGTTCAGCGCCAGCTTGGCATGCTCGCGGGCAATCTCGCTGTAGCGCTGGCTGAGGTAAGCATAATGCCGCACATCGGAACCGGTGCCGATGTAACTGGACAGGCGCTCGGCGCGAGCCAGCGACt                                                                                                             >  4:196840/1‑148 (MQ=255)
                                           aGCTTGGCCTGGCGCTCCTGGTTCAGCGCCAGCTTGGCATGCTCGCGGGCAATCTCGCTGTAGCGCTGGCTGAGGTAAGCATAATGCCGCACATCGGAACCGGTGCCGATGTAACTGGACAGGCGCTCGGCGCGAGCCAGCGACTCAcc                                                                                                         <  4:15590/149‑1 (MQ=255)
                                                   cTGGCGCTCCTGGTTCAGCGCCAGCTTGGCATGCTCGCGGGCAATCTCGCTGTAGCGCTGGCTGAGGTAAGCATAATGCCGCACATCGGAACCGGTGCCGATGTAACTGGACAGGCGCTCGGCGCGAGCCAGCGACTCACCCGCCCGaa                                                                                                 >  7:293977/1‑149 (MQ=255)
                                                      gcgcTCCTGGTTCAGCGCCAGCTTGGCATGCTCGCGGGCAATCTCGCTGTAGCGCTGGCTGAGGTAAGCATAATGCCGCACATCGGAACCGGTGCCGATGTAACTGGACAGGCGCTCGGCGCGAGCCAGCGACTCACCCGCCCGAAtca                                                                                              >  2:187064/1‑149 (MQ=255)
                                                        gcTCCTGGTTCAGCGCCAGCTTGGCATTCTCGCGGGCAATCTCGCTGTAGCGCTGGCTGAGGTAAGCATAATGCCGCACATCGGAACCGGTGCCGATGTAACTGGACAGGCGCTCTGGCGCGAGCCAGCGACTCACCCGCCCGAAtcac                                                                                             <  1:49249/149‑1 (MQ=255)
                                                                 tCAGCGCCAGCTTGGCATGCTCGCGGGCAATCTCGCTGTAGCGCTGGCTGAGGTAAGCATAATGCCGCACATCGGAACCGGTGCCGATGTAACTGGACAGGCGCTCGGCGCGAGCCAGCGACTCACCCGCCCGAATCACGTCACgcggc                                                                                   >  8:10239/1‑149 (MQ=255)
                                                                    gcgcCAGCTTGGCATGCTCGCGGGCAATCTCGCTGTAGCGCTGGCTGAGGTAAGCATAATGCCGCACATCGGAACCGGTGCCGATGTAACTGGACAGGCGCTCGGCGCGAGCCAGCGACTCACCCGCCCGAATCACGTCACGCGgcgcg                                                                                >  4:73197/1‑149 (MQ=255)
                                                                     cgcCAGCTTGGCATGCTCGCGGGCAATCTCGCTGTGGCGCTGGCTGAGGTAAGCATAATGCCGCACATCGGAACCGGTGCCGATGTAACTGGACAGGCGCTCGGCGCGAGCCAGCGACTCa                                                                                                           >  1:157203/1‑121 (MQ=255)
                                                                     cgcCAGCTTGGCATGCTCGCGGGCAATCTCGCTGTGGCGCTGGCTGAGGTAAGCATAATGCCGCACATCGGAACCGGTGCCGATGTAACTGGACAGGCGCTCGGCGCGAGCCAGCGACTCa                                                                                                           <  2:157203/121‑1 (MQ=255)
                                                                        cAGCTTGGCATGCTCGCGGGCAATCTCGCTGTAGCGCTGGCTGAGGTAAGCATAATGCCGCACATCGGAACCGGTGCCGATGTAACTGGACAGGCGCTCGGCGCGAGCCAGCGACTCACCCGCCCGAATCACGTCACGCGGCGCGCTgc                                                                            >  8:244468/1‑149 (MQ=255)
                                                                         aGCTTGGCATGCTCGCGGGCAATCTCGCTGTAGCGCTGGCTGAGGTAAGCATAATGCCGCACATCGGAACCGGTGCCGATGTAACTGGACAGGCGCTCGGCGCGAGCCAGCGACTCACCCGCCCGAATCACGTCACGCGGCGCGCTgcg                                                                           <  3:196840/149‑1 (MQ=255)
                                                                         aGCTTGGCATGCTCGCGGGCAATCTCGCTGTAGCGCTGGCTGAGGTAAGCATAATGCCGCACATCGGAACCGGTGCCGATGTAACTGGACAGGCGCTCGGCGCGAGCCAGCGACTCACCCGCCCGAATCACGTCACGCGGCGCGCTgcg                                                                           <  8:293977/149‑1 (MQ=255)
                                                                                 aTGCTCGCGGGCAATCTCGCTGTAGCGCTGGCTGAGGTAAGCATAATGCCGCACATCGGAACCGGTGCCGATGTAACTGGACAGGCGCTCGGCGCGAGCCAGCGCATCACCCGCCCGAATCACGTCACGCGGCGCGCTGCGCAGCACAt                                                                   >  7:104393/1‑149 (MQ=255)
                                                                                       gggggCAATCTCGCTGTATCGCTGGGTGATGGAATCATAATGCCGCCTAAAGGGAACCGGTGCCGATGTAACTGGACAGGCGCTCGGCGCGAGCCAGCGACTCACCCGCCCGAATCACGTCACGCGGCGCGCTGCGCAGCACATCg                                                                 <  1:121623/144‑1 (MQ=255)
                                                                                        cgGGCAATCTCGCTGTAGCGCTGGCTGAGGTAAGCATAATGCCGCACATCGGAACCGGTGCCGATGTAACTGGACAGGCGCTCGGCGCGAGCCAGCGACTCACCCGCCCGAATCACGTCACGCGGCGCGCTGCGCAGCACATCGGAat                                                             >  7:141230/1‑148 (MQ=255)
                                                                                           gCAATCTCGCTGTAGCGCTGGCTGAGGTAAGCATAATGCCGCACATCGGAACCGGTGCAGATGTAACTGGACAGGCGCTCGGCGCGAGCCAGCGACTCACCCGCCCGAATCACGTCACGCGGCGCGCTGCTCAGCACATCGGAatcatc                                                         >  1:103275/1‑149 (MQ=255)
                                                                                             aaTCTCGCTGTAGCGCTGGCTGAGGTAAGCATAATGCCGCACATCGGAACCGGTGCCGATGTAACTGGACAGGCGCTCGGCGCGAGCCAGCGACTCACCCGCCCGAATCACGTCACGCGGCGCGCTGCGCAGCACATCGGaa                                                              <  8:258092/142‑1 (MQ=255)
                                                                                             aaTCTCGCTGTAGCGCTGGCTGAGGTAAGCATAATGCCGCACATCGGAACCGGTGCCGATGTAACTGGACAGGCGCTCGGCGCGAGCCAGCGACTCACCCGCCCGAATCACGTCACGCGGCGCGCTGCGCAGCACATCGGaa                                                              >  7:258092/1‑142 (MQ=255)
                                                                                              aTCTCGCTGTAGCGCTGGCTGAGGTAAGCATAATGCCGCACATCGGAACCGGTGCCGATGTAACTGGACAGGCGCTCGGCGCGAGCCAGCGACTCACCCGCCCGAATCACGTCACGCGGCGCGCTGCGCAGCACATCGGAATCATCttt                                                      <  4:297615/149‑1 (MQ=255)
                                                                                              aTCTCGCTGTAGCGCTGGCTGAGGTAAGCATAATGCCGCACATCGGAACCGGTGCCGATGTAACTGGACAGGCGCTCGGCGCGAGCCAGCGACTCACCCGCCCGAATCACGTCACGCGGCGCGCTGCGCAGCACATCGGAATCATCttt                                                      <  7:10239/149‑1 (MQ=255)
                                                                                              aTCTCGCTGTAGCGCTGGCTGAGGTAAGCATAATGCCGCACATCGGAACCGGTGCCGATGTAACTGGACAGGCGCTCGGCGCGAGCCAGCGACTCACCCGCCCGAATCACGTCACGCGGCGCGCTGCGCAGCACATCGGAATCATCt                                                        >  1:131696/1‑147 (MQ=255)
                                                                                                ctcGCTGTAGCGCTGGCTGAGGTAAGCATAATGCCGCACATCGGAACCGGTGCCGATGTAACTGGACAGGCGCTCGGCGCGAGCCAGCGACTCACCCGCCCGAATCACGTCACGCGGCGCGCTGCGCAGCACATCGGAATCATCTTTGa                                                    >  1:185574/1‑149 (MQ=255)
                                                                                                     tGTAGCGCTGGCTGAGGTAAGCATAATGCCGCACATCGGAACCGGTGCCGATGTAACTGGACAGGCGCTCGGCGCGAGCCAGCGACTCACCCGCCCGAATCACGTCACGCGGCGCGCTGCGCAGCACATCGGAATCATCTTTGAccccc                                               >  5:155759/1‑149 (MQ=255)
                                                                                                       tAGCGCTGGCTGAGGTAAGCATAATGCCGCACATCGGAACCGGTGCCGATGTAACTGGACAGGCGCTCGGCGCGAGCCAGCGACTCACCCGCCCGAATCACGTCACGCGGCGCGCTGCGCAGCACATCGGAATCATCTTTGACCCCCTg                                             >  1:204472/1‑149 (MQ=255)
                                                                                                         gcgcTGGCTGAGGTAAGCATAATGCCGCACATCGGAACCGGTGCCGATGTAACTGGACAGGCGCTCGGCGCGAGCCAGCGACTCACCCGCCCGAATCACGTCACGCGGCGCGCTGCGCAGCACATCGGAATCATCTTTGACCCCCTgg                                            <  1:187064/148‑1 (MQ=255)
                                                                                                         gcgcTGGCTGAGGTAAGCATAATGCCGCACATCGGAACCGGTGCCGATGTAACTGGACAGGCGCTCGGCGCGAGCCAGCGACTCACCCGCCCGAATCACGTCACGCGGCGCGCTGCGCAGCACATCGGAATCATCTTTGACCCCCTgg                                            >  5:110095/1‑148 (MQ=255)
                                                                                                            ctggctgAGGTAAGCATAATGCCGCACATCGGAACCGGTGCCGATGTAACTGGACAGGCGCTCGGCGCGAGCCAGCGACTCACCCGCCCGAATCACGTCACGCGGCGCGCTGCGCAGCACATCGGAATCATCTTTGACCCCCTGGAAgg                                        <  5:315352/149‑1 (MQ=255)
                                                                                                                 tgAGGTAAGCATAATGCCGCACATCGGAACCGGTGCCGATGTAACTGGACAGGCGCTCGGCGCGAGCCAGCGACTCACCCGCCCGAATCACGTCACGCGGCGCGCTGCGCAGCACATCGGAATCATCt                                                        <  3:194386/128‑1 (MQ=255)
                                                                                                                 tgAGGTAAGCATAATGCCGCACATCGGAACCGGTGCCGATGTAACTGGACAGGCGCTCGGCGCGAGCCAGCGACTCACCCGCCCGAATCACGTCACGCGGCGCGCTGCGCAGCACATCGGAATCATCt                                                        >  4:194386/1‑128 (MQ=255)
                                                                                                                      tAAGCATAATGCCGCACATCGGAACCGGTGCCGATGTAACTGGACAGGCGCTCGGCGCGAGCCAGCGACTCACCCGCCCGAATCACGTCACGCGGCGCGCTGCGCTGCACATCGGAATCATCTTTGACCCCCTGGAAGGCAACGGTGGc                              <  1:71668/149‑1 (MQ=255)
                                                                                                                      tAAGCATAATGCCGCACATCGGAACCGGTGCCGATGTAACTGGACAGGCGCTCGGCGCGAGCCAGCGACTCACCCGCCCGAATCACGTCACGCGGCGCGCTGCGCAGCACATCGGAATCATCTTTGACCCCCTGGAAGGCAACGGTGGc                              >  5:235463/1‑149 (MQ=255)
                                                                                                                       aaGCATAATGCCGCACATCGGAACCGGTGCCGATGTAACTGGACAGGCGCTCGGCGCGAGCCAGCGACTCACCCGCCCGAATCACGTCACGCGGCGCGCTGCGCAGCACATCGGAATCATCTTTGACCCCCTGGAAGGCAACGGTGGcc                             >  7:9301/1‑149 (MQ=255)
                                                                                                                         gCATAATGCCGCACATCGGAACCGGTGCCGATGTAACTGGACAGGCGCTCGGCGCGAGCCAGCGACTCACCCGCCCGAATCACGTCACGCGGCGCGCTGCGCAGCACATCGGAATCATCTTTGACCCCCTGGAAGGCAACGGTGGCCt                            >  5:311139/1‑148 (MQ=255)
                                                                                                                         gCATAATGCCGCACATCGGAACCGGTGCCGATGTAACTGGACAGGCGCTCGGCGCGAGCCAGCGACTCACCCGCCCGAATCACGTCACGCGGCGCGCTGCGCAGCACATCGGAATCATCTTTGACCCCCTGGAAGGCAACGGTGGCCt                            <  6:311139/148‑1 (MQ=255)
                                                                                                                          cATAATGCCGCACATCGGAACCGGTGCCGATGTAACTGGACAGGCGCTCGGCGCGAGCCAGCGACTCCCCCGCCCGAATCACGTCACGCGGCGCGCTGCGCAGCACATCGGAATCATCTTTGACCCCCTGGAAGGCAACGGTGGCCTc                           >  7:170838/1‑148 (MQ=255)
                                                                                                                             aaTGCCGCACATCGGAACCGGTGCCGATGTAACTGGACAGGCGCTCGGCGCGAGCCAGCGACTCACCCGCCCGAATCACGTCACGCGGCGCGCTGCGCAGCACATCGGAATCATCTTTGACCCCCTGGAAGGCAACGGTGGCCTCATc                        >  4:219685/1‑148 (MQ=255)
                                                                                                                             aaTGCCGCACATCGGAACCGGTGCCGATGAAACTGGACAGGCGCTCGGCGCGAGCCAGCGACTCACCCGCCCGAATCACGTCACGCGGCGCGCTGCGCAGCACATCGGAATCATCTTTGACCCCCTGGAAGGCAACGGTGGCCTCATc                        >  5:205332/1‑148 (MQ=255)
                                                                                                                                  cGCACATCGGAACCGGTGCCGATGTAACTGGACAGGCGCTCGGCGCGAGCCAGCGACTCACCCGCCCGAATCACGTCACGCGGCGCGCTGCGCAGCACATCGGAATCATCTTTGACCCCCTGGAAGGCAACGGTGGCCTCATCCAgcgc                  <  5:20063/149‑1 (MQ=255)
                                                                                                                                   gttacaTCGGCACCGGTTCCGATGTAACTGGACAGGCGCTCGGCGCGAGCCAGCGACTCACCCGCCCGAATCACGTCACGCGGCGCGCTGCGCAGCACATCGGAATCATCTTTGACCCCCTGGAAGGCAACGGTg                                >  6:241147/4‑135 (MQ=255)
                                                                                                                                   gttacaTCGGCACCGGTTCCGATGTAACTGGACAGGCGCTCGGCGCGAGCCAGCGACTCACCCGCCCGAATCACGTCACGCGGCGCGCTGCGCAGCACATCGGAATCATCTTTGACCCCCTGGAAGGCAACGGTg                                <  5:241147/132‑1 (MQ=255)
                                                                                                                                   gttacaTCGGCACCGGTTCCGATGTAACTGGACAGGCGCTCGGCGCGAGCCAGCGACTCACCCGCCCGAATCACGTCACGCGGCGCGCTGCGCAGCACATCGGAATCATCTTTGACCCCCTGGAAGGCAACGGTGGCCTCATCCAgcgc                  >  4:266772/4‑149 (MQ=255)
                                                                                                                                   gCACATCGGAACCGGTGCCGATGTAACTGGACAGGCGCTCGGCGCGAGCCAGCGACTCACCCGCCCGAATCACGTCACGCGGCGCGCTGCGCAGCACATCGGAATCATCTTTGACCCCCTGGAAGGCAACGGTGGCCTCATCCAGcgcg                 <  2:131696/149‑1 (MQ=255)
                                                                                                                                     cacaTCGGAACCGGTGCCGATGTAACTGGACAGGCGCTCGGCGCGAGCCAGCGACTCACCCGCCCGAATCACGTCACGCGGCGCGCTGCGCAGCACATCGGAATCATCTTTGACCCCCTGGAAGGCAACGGTGGCCTCATCCAGCGCgg                >  5:222339/1‑149 (MQ=255)
                                                                                                                                        aTCGGAACCGGTGCCGATGTAACTGGACAGGCGCTCGGCGCGAGCCAGCGACTCACCCGCCCGAATCACGTCACGCGGCGCGCTGCGCAGCACATCGGAATCATCTTTGACCCCCTGGAAGGCAACGGTGGCCTCATCCAGCGCGGCAc             >  5:60972/1‑149 (MQ=255)
                                                                                                                                        aTCGGAACCGGTGCCGATGTAACTGGACAGGCGCTCGGCGCGAGCCAGCGACTCACCCGCCCGAAGCACGTCACGCGGCGCGCTGCGCAGCACATCGGAATCATCTTTGCCCCCCTGGAAGGCAACGGTGGCCTCATCCAGCGCGGCAc             >  8:288464/1‑149 (MQ=255)
                                                                                                                                           ggAACCGGTGCCGATGTAACTGGACAGGCGCTCGGCGCGAGCCAGCGACTCACCCGCCCGAATCACGTCACGCGGCGCGCTGCGCAGCACATCGGAATCATCTTTGACCCCCTGGAAGGCAACGGTGGCCTCATCCAGCGCGGCActgc          >  1:221014/1‑149 (MQ=255)
                                                                                                                                             aaCCGGTTCCGATGTAACTGGACAGGCGCGCGGCGCGAGCCAGCGACTCACCCGCCCGAATCACGTCTCGCGGCGCGCTGCGCAGCACATCGGAATCATCTTTGTCCCCCTGGAAGGCAACGGTGGCCTCATCCAGCGCGGCActgctg        <  8:104393/149‑1 (MQ=255)
                                                                                                                                             aaCCGGTGCCGATGTAACTGGACAGGCGCTCGGCGCGAGCCAGCGACTCACCCGCCCGAATCACGTCACGCGGCGCGCTGCGCAGCACATCGGAATCATCTTTGACCCCCTGGAAGGCAACGGTGGCCTCATCCAGCGCGGCActgctg        <  8:267357/149‑1 (MQ=255)
                                                                                                                                                 ggTGCCGATGTAACTGGACAGGCGCTCGGCGCGAGCCAGCGACTCACCCGCCCGAATCACGTCACGCGGCGCGCTGCGCAGCACATCGGAATCATCTTTGACCCCCTGGAAGGCAACGGTGGCCTCATCCAGCGCGGCACTGCTGCGCt    <  8:141230/149‑1 (MQ=255)
                                                                                                                                                  gTGCCGATGTAACTGGACAGGCGCTCGGCGCGAGCCAGCGACTCACCCGCCCGAATCACGTCACGCGGCGCGCTGCGCAGCACATCGGAATCATCTTTGACCCCCTGGAAGGCAACGGTGGCCTCATCCAGCGCAGCAGTGCCGCGCTg   <  2:204472/149‑1 (MQ=255)
                                                                                                                                                   tGCCGATGTAACTGGACAGGCGCTCGGCGCGAGCCAGCGACTCACCCGCCCGAATCACGTCACGCGGCGCGCTGCGCAGCACATCGGAATCATCTTTGACCCCCTGGAAGGCAACGGTGGCCTCATCCAGCGCGGCACTGCTGCGCTGa  <  6:110095/149‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                    
GGCCTCACGCAAGGCCAGCTGCAGGCGCTGGCGCTCCAGGTCGAGCTTGGCCTGGCGCTCCTGGTTCAGCGCCAGCTTGGCATGCTCGCGGGCAATCTCGCTGTAGCGCTGGCTGAGGTAAGCATAATGCCGCACATCGGAACCGGTGCCGATGTAACTGGACAGGCGCTCGGCGCGAGCCAGCGACTCACCCGCCCGAATCACGTCACGCGGCGCGCTGCGCAGCACATCGGAATCATCTTTGACCCCCTGGAAGGCAACGGTGGCCTCATCCAGCGCGGCACTGCTGCGCTGA  >  NC_002947/4744460‑4744754

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: