Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,707,227 G→C 8.4% intergenic (+45/+97) PP_1501 → / ← PP_1502 hypothetical protein/OmpA family protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029471,707,2270GC 8.4% 180.3 / 5.6 58intergenic (+45/+97)PP_1501/PP_1502hypothetical protein/OmpA family protein
Reads supporting (aligned to +/- strand):  ref base G (27/26);  new base C (3/2);  total (30/28)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TTCGACACCCAGCAGAAAGGCCTGACCGGCGTGCTGGTGGACAACCGCAAGCGCGCGGTGCGCTGCATGGGCGAGGAAACACCGCGGCAGGTGCAGTGAGGTTTCGGCGCAGCATTTTTGGCGCCCGCGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGGCACAACATTATCG‑‑GCGAACAAAGAAAA‑‑GCCCCGCG‑TCGATCGCGGGGCTTGTCCTTTCAGCTAGCCAGCCTTAGCGCTGAACCTGGCTCACCGAGCGTGGTGCCACCGGCTGGTTGTCGTTG  >  NC_002947/1707084‑1707375
                                                                                                                                               |                                                                                                                                                         
ttCGACACCCAGCAGAAAGGCCTGACCGGCGTGCTGGTGGACAACCGCAAGCGCGCGGTGCGCTGCATGGGCGAGGAAACACCGCGGCAGGTGCAGTGAGGTTTCGGCGCAGCATTTTTGGCGCCCGCGAGATCGAGCGCCGCgcgg                                                                                                                                                        <  5:399232/147‑1 (MQ=255)
   gACACCCAGCAGAAAGGCCTGACCGGCGTGCTGGTGGACAACCGCAAGCGCGCGGTGCGCTGCATGGGCGAGGAAACACCGCGGCAGGTGCAGTGAGGTTTCGGCGCAGCATTTTTGGCGCCCGCGAGATCGAGCGCCGCGCGggcggc                                                                                                                                                   >  5:181597/1‑149 (MQ=255)
    acacCCAGCAGAAAGGCCTGACCGGCGTGCTGGTGGACAACCGCAAGCGCGCGGTGCGCTACATGGGCGAGGAAACACCGCGGCAGGTGCAGTGAGGTTTCGGCGCAGCATTTTTGGCGCCCGCGAGATCGAGCGCCGCGCGGgcggcg                                                                                                                                                  <  3:386850/149‑1 (MQ=255)
       cccAGCAGAAAGGCCTGACCGGCGTGCTGGTGGACAACCGCAAGCGCGCGGTGCGCTGCATGGGCGAGGAAACACCGCGGCAGGTGCAGTGAGGTTTCGGCGCAGCATTTTTGGCGCCCGCGAGATCGAGCGCCGCGCGGGCGGCGCTc                                                                                                                                               <  3:111644/149‑1 (MQ=255)
       cccAGCAGAAAGGCCTGACCGGCGTGCTGGTGGACAACCGCAAGCGCGCGGTGCGCTGCATGGGCGAGGAAACACCGCGGCAGGTGCAGTGAGGTTTCGGCGCAGCATTTTTGGCGCCCGCGAGATCGAGCGCCGCGCGGGCGGCGCTc                                                                                                                                               >  2:48535/1‑149 (MQ=255)
        ccAGCAGAAAGGCCTGACCGGCGTGCTGGTGGACAACCGCAAGCGCGCGGTGCGCTGCATGGGCGAGGAAACACCGCGGCAGGTGCAGTGAGGTTTCGGCGCAGCATTTTTGGCGCCCGCGAGATCGAGCGCCGCGCGGGCGGCGCTc                                                                                                                                               >  3:317808/1‑148 (MQ=255)
        ccAGCAGAAAGGCCTGACCGGCGTGCTGGTGGACAACCGCAAGCGCGCGGTGCGCTGCATGGGCGAGGAAACACCGCGGCAGGTGCAGTGAGGTTTCGGCGCAGCATTTTTGCCGCCCGCGAGATCGAGCGCCGCGCGGGCGGCGCTCg                                                                                                                                              <  4:98234/149‑1 (MQ=255)
          agcagAAAGGCCTGACCGGCGTGCTGGTGGACAACCGCAAGCGCGCGGTGCGCTGCATGGGCGAGGAAACACCGCGGCAGGTGCAGTGAGGTTTCGGCGCAGCATTTTTGGCGCCCGCGAGATCGAGCGCCGCGCGGGCGGCGCTCGAt                                                                                                                                            <  2:218830/149‑1 (MQ=255)
               aaaGGCCTGACCGGCGTGCTGGTGGACAACCGCAAGCGCGCGGTGCGCTGCATGGGCGAGGAAACACCGCGGCAGGTGCAGTGAGGTTTCGGCGCAGCATTTTTGGCGCCCGCGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTcac                                                                                                                                       <  1:238039/149‑1 (MQ=255)
                aaGGCCTGACCGGCGTGCTGGTGGACAACCGCAAGCGCGCGGTGCGCTGCATGGGCGAGGAAACACCGCGGCAGGTGCAGTGAGGTTTCGGCGCAGCATTTTTGGCGCCCGCGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTcac                                                                                                                                       <  2:129022/148‑1 (MQ=255)
                  ggCCTGACCGGCGTGCTGGTGGACAACCGCAAGCGCGCGGTGCGCTGCATGGGCGAGGAAACACCGCGGCAGGTGCAGTGAGGTTTCGGCGCAGCATTTTTGGCGCCCGCGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAgg                                                                                                                                    <  6:368553/149‑1 (MQ=255)
                  ggCCTGACCGGCGTGCTGGTGGACAACCGCAAGCGCGCGGTGCGCTGCATGGGCGAGGAAACACCGCGGCAGGTGCAGTGAGGGTTCGGCGCAGCATTTTTGGCGCCCGCGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAgg                                                                                                                                    <  3:344369/149‑1 (MQ=255)
                        aCCGGCGTGCTGGTGGACAACCGCAAGCGCGCGGTGCGCTGCATGGGCGAGGAAACACCGCGGCAGGTGCAGTGAGGTTTCGGCGCAGCATTTTTGGCGCCCGCGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGGcac                                                                                                                              <  5:119310/149‑1 (MQ=255)
                        aCCGGCGTGCTGGTGGACAACCGCAAGCGCGCGGTGCGCTGCATGGGCGAGGAAACACCGCGGCAGGTGCAGTGAGGTTTCGGCGCAGCATTTTTGGCGCCCGCGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGGca                                                                                                                               >  7:403822/1‑148 (MQ=255)
                            gCGTGCTGGTGGACAACCGCAAGCGCGCGGTGCGCTGCATGGGCGAGGAAACACCGCGGCAGGTGCAGTGAGGTTTCGGCGCAGCATTTTTGGCGCCCGCGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGGCacaaca                                                                                                                          >  3:156972/1‑149 (MQ=255)
                                 cTGGTGGACAACCGCAAGCGCGCGGTGCGCTGCATGGGCGAGGAAACACCGCGGCAGGTGCAGTGAGGTTTCGGCGCAGCATTTTTGGCGCCCGCGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGGCACAACATTATc                                                                                                                     <  8:403822/149‑1 (MQ=255)
                                  tggtggACAACCGCAAGCGCGCGGTGCGCTGCATGGGCGAGGAAACACCGCGGCAGGTGCAGTGAGGTTTCGGCGCAGCATTTTTGGCGCCCGCGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGGCACAACATTATCg                                                                                                                    >  1:64045/1‑149 (MQ=255)
                                  tggtggACAACCGCAAGCGCGCGGTGCGCTGCATGGGCGAGGAAACACCGCGGCAGGTGCAGTGAGGTTTCGGCGCAGCATTTTTGGCGCCCGCGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGGCACAACATTATCg                                                                                                                    >  3:320470/1‑149 (MQ=255)
                                    gtggACAACCGCAAGCGCGCGGTGCGCTGCATGGGCGAGGAAACACCGCGGCAGGTGCAGTGAGGTTTCGGCGCAGCATTTTTGGCGCCCGCGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGGCACAACATTATCg                                                                                                                    <  1:48535/147‑1 (MQ=255)
                                          aaCCGCAAGCGCGCGGTGCGCTGCATGGGCGAGGAAACACCGCGGCAGGTGCAGTGAGGTTTCGGCGCAGCATTTTTGGCGCCCGCGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGGCACAACATTATCG‑‑GCGAACa                                                                                                           >  6:366525/1‑148 (MQ=255)
                                                         gTGCGCTGCATGGGCGAGGAAACACCGCGGCAGGTGCAGTGAGGTTTCGGCGCAGCATTTTTGGCGCCCGCGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGGCACAACATTATCG‑‑GCGAACAAAGAAAA‑‑GCCCCGCG‑t                                                                                        <  2:407470/149‑1 (MQ=255)
                                                          tGCGCTGCATGGGCGAGGAAACACCGCGGCAGGTGCAGTGAGGTTTCGGCGCAGCATTTTTGGCGCCCGCGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGGCACAACATTATCG‑‑GCGAACAAAGAAAA‑‑GCCCCGCG‑Tc                                                                                       >  5:261877/1‑149 (MQ=255)
                                                             gcTGCATGGGCGAGGAAACACCGCGGCAGGTGCAGTGAGGTTTCGGCGCAGCATTTTTGGCGCCCGCGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGGCACAACATTATCG‑‑GCGAACAAAGAAAA‑‑GCCCCGCG‑TCGAt                                                                                    >  5:261873/1‑149 (MQ=255)
                                                               tGCATGGGCGAGGAAACACCGCGGCAGGTGCAGTGAGGTTTCGGCGCAGCATTTTTGGCGCCCGCGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGGCACAACATTATCG‑‑GCGAACAAAGAAAA‑‑GCCCcgcgatcgacg                                                                                   >  2:187490/1‑142 (MQ=255)
                                                                  aTGGGCGAGGAAACACCGCGGCAGGTGCAGTGAGGTTTCGGCGCAGCATTTTTGGCGCCCGCGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGGCACAACATTATCG‑‑GCGAACAAAGAAAA‑‑GCCCCGCG‑TCGATCGCgg                                                                               <  2:64045/149‑1 (MQ=255)
                                                                    gggCGAGGAAGCACCGCGGCAGGTGCAGTGAGGTTTCGGCGCAGCATTTTTGGCGCCCGCGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGGCACAACATTATCG‑‑GCGAACAAAGAAAA‑‑GCCCCCCG‑ACGATCCCcggg                                                                             >  4:11800/1‑149 (MQ=39)
                                                                       cGAGGAAACACCGCGGCAGGTGCAGTGAGGTTTCGGCGCAGCATTTTTGGCGCCCGCGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGGCACAACATTATCG‑‑GCGAACaa                                                                                                          >  3:144843/1‑120 (MQ=255)
                                                                       cGAGGAAACACCGCGGCAGGTGCAGTGAGGTTTCGGCGCAGCATTTTTGGCGCCCGCGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGGCACAACATTATCG‑‑GCGAACaa                                                                                                          <  4:144843/120‑1 (MQ=255)
                                                                       cGAGGAAACACCGCGGCAGGTGCAGTGAGGTTTCGGCGCAGCATTTTTGGCGCCCGCGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGGCACAACATTATCG‑‑GCGAACAAAGAAAA‑‑GCCCCGCGATCGA‑CGCGGGGCtt                                                                          <  1:187490/149‑1 (MQ=255)
                                                                       cGAGGAAACACCGCGGCAGGTGCAGTGAGGTTTCGGCGCAGCATTTTTGGCGCCCGCGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGGCACAACATTATCG‑‑GCGAACAAAGAAAA‑‑GCCCCGCG‑TCGATCGCGGGGCtt                                                                          >  8:63249/1‑149 (MQ=255)
                                                                           gAAACACCGCGGCAGGTGCAGTGAGGTTTCGGCGCAGCATTTTTGGCGCCCGCGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGGCACAACATTATCG‑‑GCGAACAAAGAAAA‑‑GCCCCGCG‑TCGATCGCGGGGCTTGTcc                                                                      <  4:156972/149‑1 (MQ=255)
                                                                           gAAACACCGCGGCAGGTGCAGTGAGGTTTCGGCGCAGCATTTTTGGCGCCCGCGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGGCACAACATTATCG‑‑GCGAACAAAGAAAA‑‑GCCCCGCG‑TCGATCGCGGGGCTTGTcc                                                                      >  6:142283/1‑149 (MQ=255)
                                                                            aaaCACCGCGGCAGGTGCAGTGAGGTTTCGGCGCAGCATTTTTGGCGCCCGCGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGGCACAACATTATCG‑‑GCGAACAAAGAAAA‑‑GCCCCGCG‑TCGATCGCGGGGCTTGTCCt                                                                     <  2:429678/149‑1 (MQ=255)
                                                                               cacCGCGGCAGGTGCAGTGAGGTTTCGGCGCAGCATTTTTGGCGCCCGCGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGGCACAACATTATCG‑‑GCGAACAAAGAAAA‑‑GCCCCGCG‑TCGATCGCGGGGCTTGTCCTTTc                                                                  <  2:383570/149‑1 (MQ=255)
                                                                                acCGCGGCAGGTGCAGTGAGGTTTCGGCGCAGCATTTTTGGCGCCCGCGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGGCACAACATTATCG‑‑GCGAACAAAGAAAA‑‑GCCCCGCG‑TCGATCGCGGGGCTTGTCCtta                                                                   >  6:136450/1‑146 (MQ=255)
                                                                                    cgGCAGGTGCAGTGAGGTTTCGGCGCAGCATTTTTGGCGCCCGCGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGGCACAACATTATCG‑‑GCGAACAAAGAAAA‑‑GCCCCGCG‑TCGATCGCGGGGCTTGTCCTTTCAGCTa                                                             <  5:52819/149‑1 (MQ=255)
                                                                                        aGGTGCAGTGAGGTTTCGGCGCAGCATTTTTGGCGCCCGCGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGGCACAACATTATCG‑‑GCGAACAAAGAAAA‑‑GCCCCGCG‑TCGATCGCGGGGCTTGTCCCTTCAGCTagcca                                                         <  6:181597/149‑1 (MQ=255)
                                                                                          gTGCAGTGAGGTTTCGGCGCAGCATTTTTGGCGCCCGCGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGGCACAACATTATCG‑‑GCGAACAAAGAAAA‑‑GCCCCGCG‑TCGATCGCGGGGCTTGTCCTTTCAGCTagccagc                                                       >  7:296712/1‑149 (MQ=255)
                                                                                               gTGAGGTTTCGGCGCAGCATTTTTGGCGCCCGCGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGGCACAACATTATCG‑‑GCGAACAAAGAAAA‑‑GCCCCGCG‑TCGATCGCGGGGCTTGTCCttt                                                                   <  7:260708/132‑1 (MQ=255)
                                                                                               gTGAGGTTTCGGCGCAGCATTTTTGGCGCCCGCGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGGCACAACATTATCG‑‑GCGAACAAAGAAAA‑‑GCCCCGCG‑TCGATCGCGGGGCTTGTCCttt                                                                   >  8:260708/1‑132 (MQ=255)
                                                                                                tGAGGTTTCGGCGCAGCATTTTTGGCGCCCGCGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGGCACAACATTATCG‑‑GCGAACAAAGAAAA‑‑GCCCCGCG‑TCGATCGCGGGGCTTGTCCTTTCAGCTAGCCAGCCTTAgc                                                 <  7:333755/149‑1 (MQ=255)
                                                                                                tGAGGTTTCGGCGCAGCATTTTTGGCGCCCGCGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGGCACAACATTATCG‑‑GCGAACAAAGAAAA‑‑GCCCCGCG‑TCGATCGCGGGGCTTGTCCTTTCAGCTAGCCAGCCTTAgc                                                 >  4:213214/1‑149 (MQ=255)
                                                                                                  aGGTTTCGGCGCAGCATTTTTGGCGCCCGCGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGGCACAACATTATCG‑‑GCGAACAAAGAAAA‑‑GCCCCGCG‑TCGATCGCGGGGCTTGTCCTTTCAGCTAGCCAGCCTTAgcg                                                <  3:11800/148‑1 (MQ=255)
                                                                                                     tttCGGCGCAGCATTTTTGGCGCCCGGGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGGCACAACATTATCG‑‑GCGAACAAAGAAAA‑‑GCCCCGCG‑TCGATCGCGGGGCTTGTCCTTTCAGCTAGCCAGCCTTAGCGCTg                                             >  1:233534/1‑148 (MQ=255)
                                                                                                         ggCGCAGCATTTTTGGCGCCCGCGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGGCACAACATTATCG‑‑GCGAACAAAGAAAA‑‑GCCCCGCG‑TCGATCGCGGGGCTTGTCCTTTCAGCTAGCCAGCCTTAGCGCTGAAcc                                         <  6:66143/148‑1 (MQ=255)
                                                                                                              agcaTTTTTGGCGCCCGCGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGGCACAACATTATGGTAGCGAACAAAGAAAAAGGCCCCGCG‑TCGATCGCGGGGCTTGTCCTTTCAGCTAGCCAGCCTTAGCGCTGAACCTg                                       >  1:353641/1‑149 (MQ=38)
                                                                                                                caTTTTTGGCGCCCGCGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGGCACAACATTATCG‑‑GCGAACAAAGAAAA‑‑GCCCCGCG‑TCGATCGCGGGGCTTGTCCTTTCAGCTAGCCAGCCTTAGCGCTGAACCTGGCTCAc                                 >  6:101157/1‑149 (MQ=255)
                                                                                                                caTTTTTGGCGCCCGCGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGGCACAACATTATCG‑‑GCGAACAAAGAAAA‑‑GCCCCGCG‑TCGATCGCGGGGCTTGTCCTTTCAGCTAGCCAGCCTTAGCGCTGAACCTGGCTCAc                                 >  6:327475/1‑149 (MQ=255)
                                                                                                                    tttGGCGCCCGCGAGATCGAGCGCCGCGCGGGCGGCGGTCGCTCTCACAGGCGGCACAACATTATCG‑‑GCGAACAAAGAAAA‑‑GCCCCGCG‑TCGATCGCGGGGCTTGTCCTTTCAGCTAGCCAGCCTTAGCGCTGAACCTGGCTCTCCGAg                             >  7:317191/1‑149 (MQ=255)
                                                                                                                        gcgcCCGCGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGGCACAACATTATCG‑‑GCGAACAAAGAAAA‑‑GCCCCGCG‑TCGATCGCGGGGCTTGTCCTTTCAGCTAGCCAGCCTTa                                                   <  8:212753/123‑1 (MQ=255)
                                                                                                                        gcgcCCGCGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGGCACAACATTATCG‑‑GCGAACAAAGAAAA‑‑GCCCCGCG‑TCGATCGCGGGGCTTGTCCTTTCAGCTAGCCAGCCTTa                                                   >  7:212753/1‑123 (MQ=255)
                                                                                                                          gcCCGCGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGGCACAACATTATCG‑‑GCGAACAAAGAAAA‑‑GCCCCGCG‑TCGATCGCGGGGCTTGTCCTTTCAGCTAGCCAGCCTTAGCGCTGAACCTGGCTCACCGAGCgtggt                       <  3:213214/149‑1 (MQ=255)
                                                                                                                                   aTCGAGCGCCGCCCGCGCGGCGCTCGATCTCACAGGCGGCACAACATTATCG‑‑GCGAACAAAGAAAA‑‑GCCCCGCG‑TCGATCGCGGGGCTTGTCCTTTCAGCTAGCCAGCCTTAGCGCTGAACCTGGCTCACCGAGCGTGGTGc                     <  3:161750/142‑1 (MQ=255)
                                                                                                                                   aTCGAGCGCCGCCCGCGCGGCGCTCGATCTCACAGGCGGCACAACATTATCG‑‑GCGAACAAAGAAAA‑‑GCCCCGCG‑TCGATCGCGGGGCTTGTCCTTTCAGCTAGCCAGCCTTAGCGCTGAACCTGGCTCACCGAGCGTGGTGc                     >  4:161750/1‑142 (MQ=255)
                                                                                                                                    tCGAGCGCCGCCCGCGCGGCGCTCGATCTCACAGGCGGCACAACATTATCG‑‑GCGAACAAAGAAAA‑‑GCCCCGCG‑TCGATGGCGGGGCTTGTCCTTTCAGCTAGCCAGCCTTAGCGCTGAACCTgg                                      >  7:57187/1‑124 (MQ=255)
                                                                                                                                    tCGAGCGCCGCCCGCGCGGCGCTCGATCTCACAGGCGGCACAACATTATCG‑‑GCGAACAAAGAAAA‑‑GCCCCGCG‑TCGATCGCGGGGCTTGTCCTTTCAGCTAGCCAGCCTTAGCGCTGAACCTgg                                      <  8:57187/124‑1 (MQ=255)
                                                                                                                                      gAGCGCCGCCTGCGCGGCGCTCGATCTCGCGGGCGGCACAACATTATCG‑‑GCGAACAAAGAAAA‑‑GCCCCGCG‑TCGATCGCGGGGCTTGTCCTTTCAGCTAGCCAGCCTTAGCGCTGAACCTGGCTCACCGAGCGTGGTGCCACCGGCTgg           >  4:367701/1‑149 (MQ=255)
                                                                                                                                       aGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGGCACAACATTATCG‑‑GCGAACAAAGAAAA‑‑GCCCCGCG‑TCGATCGCGGGGCTTGTCCTTTCAGCTAGCCAGCCTTAGCGCTGAACCTGGCTCACCGAGCGTGGTGCCACCGGCTGGt          >  5:63692/1‑149 (MQ=255)
                                                                                                                                           ccgcGCGGGCGGCGCTCGATCTCACAGGCGGCACAACATTATCG‑‑GCGAACAAAGAAAA‑‑GCCCCGCG‑TCGATCGCGGGGCTTGTCCTTTCAGCTAGCCAGCCTTAGCGCTGAACCTGGCTCACCGAGCGTGGTGCCACCGGCTGGTTGTc      >  3:5036/1‑149 (MQ=255)
                                                                                                                                            cgcgcgGGCGGCGCTCGATCTCACAGGCGGCACAACATTATCG‑‑GCGAACAAAGAAAA‑‑GCCCCGCG‑TCGATCGCGGGGCTTGTCCTTTCAGCTAGCCAGCCTTAGCGCTGAACCTGGCTCACCGAGCGTGGTGCCACCGGCTGGTTGTCg     <  4:320470/149‑1 (MQ=255)
                                                                                                                                            cgcccgCGCGGCGCTCGATCTCACAGGCGGCACAACATTATCG‑‑GCGAACAAAGAAAA‑‑GCCCCGCG‑TCGATCGCGGGGCTTGTCCTTTCAGCTAGCCAGCCTTAGCGCTGAACCTGGCt                                    >  6:289642/1‑118 (MQ=255)
                                                                                                                                            cgcccgCGCGGCGCTCGATCTCACAGGCGGCACAACATTATCG‑‑GCGAACAAAGAAAA‑‑GCCCCGCG‑TCGATCGCGGGGCTTGTCCTTTCAGCTAGCCAGCCTTAGCGCTGAACCTGGCt                                    <  5:289642/118‑1 (MQ=255)
                                                                                                                                               gcgggcggCGCTCGATCTCACAGGCGGCACAACATTATCG‑‑GCGAACAAAGAAAA‑‑GCCCCGCG‑TCGATCGCGGGGCTTGTCCTTTCAGCTAGCCAGCCTTAGCGCTGAACCTGGCTCACCGAGCGTGGTGCCACCGGCTGGTTGTCGttt  <  5:327475/149‑2 (MQ=255)
                                                                                                                                               gcgggcggCGCTCGATCTCACAGGCGGCACAACATTATCG‑‑GCGAACAAAGAAAA‑‑GCCCCGCG‑TCGATCGCGGGGCTTGTCCTTTCAGCTAGCCAGCCTTAGCGCTGAACCTGGCTCACCGAGCGTGGTGCCACCGGCTGGTTGTCGTTg  >  5:122024/1‑149 (MQ=255)
                                                                                                                                               |                                                                                                                                                         
TTCGACACCCAGCAGAAAGGCCTGACCGGCGTGCTGGTGGACAACCGCAAGCGCGCGGTGCGCTGCATGGGCGAGGAAACACCGCGGCAGGTGCAGTGAGGTTTCGGCGCAGCATTTTTGGCGCCCGCGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGGCACAACATTATCG‑‑GCGAACAAAGAAAA‑‑GCCCCGCG‑TCGATCGCGGGGCTTGTCCTTTCAGCTAGCCAGCCTTAGCGCTGAACCTGGCTCACCGAGCGTGGTGCCACCGGCTGGTTGTCGTTG  >  NC_002947/1707084‑1707375

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: ‑