Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,155,364 C→G 6.9% G300A (GGC→GCC)  PP_2771 ← sigma‑54 dependent transcriptional regulator

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029473,155,3640CG 6.9% 215.3 / 4.4 69G300A (GGC→GCC) PP_2771sigma‑54 dependent transcriptional regulator
Reads supporting (aligned to +/- strand):  ref base C (34/30);  new base G (2/3);  total (36/33)
Fisher's exact test for biased strand distribution p-value = 6.65e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.36e-01

ATTGACCACCAGTTCGCCGATGGCGATCAGTCGGGCGCGGGTGATGTTGCGGCTCAGGCCCAGCAGGCTGGCCGTGTGCACCTGGTTGCAGTGGCAGAAACGATAGGCACTGCGTAGCAGGCTGTCCTCGACCTTCTCGTGCAGGCCGCCGGCCTGCTCCTCGAACAGGCGCTGAAAGGCTCGCAGCAACTGTTCTTCGGCGCTTTCGGCGGTGCTCACCGGCGACTTGTCCTGGCGCTCG‑ATGCGCAGGTGCGACAGGCGCAGGTCGTCGTCCTGCACCAGGCCATTGCGGCACA  >  NC_002947/3155216‑3155511
                                                                                                                                                    |                                                                                                                                                    
aTTGACCACCAGTTCGCCGATGGCGATCAGTCGGGCGCGGGTGATGTTGCGGCTCAGGCCCAGCAGGCTGGCCGTGTGCACCTGGTTGCAGTGGCAGAAACGATAGGCACTGCGTAGCAGGCTGTCCTCGACCTTCTCGTGCAGgccgc                                                                                                                                                      >  6:100255/1‑149 (MQ=255)
    accaccAGTTCGCCGATGGCGATCAGTCGGGCGCGGGTGATGTTGCGGCTCAGGCCCAGCAGGCTGGCCGTGTGCACCTGGTTGCAGTGGCAGAAACGATAGGCACTGCGTAGCAGGCTGTCCTCGACCTTCTCGTGCAGGCCGCCGGc                                                                                                                                                  >  1:14630/1‑149 (MQ=255)
    accaccAGTTCGCCGATGGCGATCAGTCGGGCGCGGGTGATGTTGCGGCTCAGGCCCAGCAGGCTGGCCGTGTGCACCTGGTTGCAGTGGCAGAAACGATAGGCACTGCGTAGCAGGCTGTCCTCGACCTTCTCGTGCAGGCCGCCGGc                                                                                                                                                  >  8:341159/1‑149 (MQ=255)
     ccaccaGTTCGCCGATGGCGATCAGTCGGGCGCGGGTGATGTTGCGGCTCAGGCCCAGCAGGCTGGCCGTGTGCACCTGGTTGCAGTGGCAGAAACGATAGGCACTGCGTAGCAGGCTGTCCTCGACCTTCTCGTGCAGGCCGGCGGcc                                                                                                                                                 >  6:349132/1‑149 (MQ=255)
     ccaccaGTTCGCCGATGGCGATCAGTCGGGCGCGGGTGATGTTGCGGCTCAGGCCCAGCAGGCTGGCCGTGTGCACCTGGTTGCAGTGGCAGAAACGATAGGCACTGCGTAGCAGGCTGTCCTCGACCTTCTCGTGCAGGCCGCCGGcc                                                                                                                                                 >  7:349148/1‑149 (MQ=255)
      caccaGTTCGCCGATGGCGATCAGTCGGGCGCGGGTGATGTTGCGGCTCAGGCCCAGCAGGCTGGCCGTGTGCACCTGGTTGCAGTGGCAGAAACGATAGGCACTGCGTAGCAGGCTGTCCTCGACCTTCTCGTGCAGGCCGCCGGCCt                                                                                                                                                >  4:300676/1‑149 (MQ=255)
        ccaGTTCGCCGATGGCGATCAGTCGGGCGCGGGTGATGTTGCGGCTCAGGCCCAGCAGGCTGGCCGTGTGCACCTGGTTGCAGTGGCAGAAACGATAGGCACTGCGTAGCAGGCTGTCCTCGACCTTCTCGTGCAGGCCGCCGGCCTGc                                                                                                                                              >  7:272257/1‑149 (MQ=255)
           gTTCGCCGATGGCGATCAGTCGGGCGCGGGTGATGTTGCGGCTCAGGCCCAGCAGGCTGGCCGTGTGCACCTGGTTGCAGTGGCAGAAACGATAGGCACTGCGTAGCAGGCTGTCCTCGACCTTCTCGTGCAGGCCGCCGGCCTGctc                                                                                                                                            >  8:242789/1‑148 (MQ=255)
            ttCGCCGATGGCGATCAGTCGGGCGCGGGTGATGTTGCGGCTCAGGCCCAGCAGGCTGGCCGTGTGCACCTGGTTGCAGTGGCAGAAACGATAGGCACTGCGTAGCAGGCTGTCCTCGACCTTCTCGTGCAGGCCGGCGGCCTGcacga                                                                                                                                          <  5:349132/149‑5 (MQ=255)
            ttCGCCGATGGCGATCAGTCGGGCGCGGGTGATGTTGCGGCTCAGGCCCAGCAGGCTGGCCGTGTGCACCTGGTTGCAGTGGCAGAAACGATAGGCACTGCGTAGCAGGCTGTCCTCGACCTTCTCGAGCAGGCCGGCGGCCTGcacga                                                                                                                                          <  6:223510/149‑5 (MQ=255)
                ccGATGGCGATCAGTCGGGCGCGGGTGATGTTGCGGCTCAGGCCCAGCAGGCTGGCCGTGTGCACCTGGTTGCAGTGGCAGAAACGATAGGCACTGCGTAGCAGGCTGTCCTCGACCTTCTCGTGCAGGCCGCCGGCCTGCTCCTCGaa                                                                                                                                      >  5:390043/1‑149 (MQ=255)
                      gCGATCAGTCGGGCGCGGGTGATGTTGCGGCTCAGGCCCAGCAGGCTGGCCGTGTGCACCTGGTTGCAGTGGCAGAAACGATAGGCACTGCGTAGCAGGCTGTCCTCGACCTTCTCGTGCAGGCCGCCGGCCTGCTCCTCGAACAGgc                                                                                                                                 <  4:201985/148‑1 (MQ=255)
                        gATCAGTCGGGCGCGGGTGATGTTGCGGCTCAGGCCCAGCAGGCTGGCCGTGTGCACTTGGTTGCAGTGGCAGAAACGATAGGCACTGCGTAGCAGGCTGTCCTCGACCTTCTCGTGCAGGCCGCCGGCCTGCTCCTCGAACAGGCGCt                                                                                                                              <  5:387148/149‑1 (MQ=255)
                        gATCAGTCGGGCGCGGGTGATGTTGCGGCTCAGGCCCAGCAGGCTGGCCGTGTGCACCTGGTTGCAGTGGCAGAAACGATAGGCACTGCGTAGCAGGCTGTCCTCGACCTTCTCGTGCAGGCCGCCGGCCTGCTCCTCGAACAGGCGCt                                                                                                                              <  3:360380/149‑1 (MQ=255)
                            aGTCGGGCGCGGGTGATGTTGCGGCTCAGGCCCAGCAGGCTGGCCGTGTGCACCTGGTTGCAGTGGCAGAAACGATAGGCACTGCGTAGCAGGCTGTCCTCGACCTTCTCGTGCAGGCCGCCGGCCTGCTCCTCGAACAGGCGCTGaaa                                                                                                                          <  5:372642/149‑1 (MQ=255)
                            aGTCGGGCGCGGGTGATGTTGCGGCTCAGGCCCAGCAGGCTGGCCGTGTGCACCTGGTTGCAGTGGCAGAAACGATAGGCACTGCGTAGCAGGCTGTCCTCGACCTTCTCGTGCAGGCCGCCGGCCTGCTCCTCGAACAGGCGCTGaaa                                                                                                                          <  5:378265/149‑1 (MQ=255)
                                  gcgcGGGTGATGTTGCGGCTCAGGCCCAGCAGGCTGGCCGTGTGCACCTGGTTGCAGTGGCAGAAACGATAGGCACTGCGTAGCAGGCTGTCCTCGACCTTCTCGTGCAGGCCGCCGGCCTGCTCCTCGAACAGGCGCTGAAAGGCTCg                                                                                                                    <  6:140353/149‑1 (MQ=255)
                                    gcgGGTGATGTTGCGGCTCAGGCCCAGCAGGCTGGCCGTGTGCACCTGGTTGCAGTGGCAGAAACGATAGGCACTGCGTAGCAGGCTGTCCTCGACCTTCTCGTGCAGGCCGCCGGCCTGCTCCTCGAACAGGCGCTGAAAGGCTCgca                                                                                                                  <  4:179343/149‑1 (MQ=255)
                                    gcgGGTGATGTTGCGGCGCAGGCCCAGCAGGCTGGCCGTGTGCACCTGGTTGCAGTGGCAGAAACGATAGGCACTGCGTAGCAGGCTGTCCTCGACCTTCTCGTGCAGGCCGCCGGCCTGCTCCTCGAACAGGCGCTGAAAGGCTCgca                                                                                                                  <  5:313429/149‑1 (MQ=255)
                                           aTGTTGCGGCTCAGGCCCAGCAGGCTGGCCGTGTGCACCTGGTTGCAGTGGCAGAAACGATAGGCACTGCGTAGCAGGCTGTCCTCGACCTTCTCGTGCAGGCCGCCGGCCTGCTCCTCGAACAGGCGCTGAAAGGCTCGCAGCAACTg                                                                                                           <  1:53772/149‑1 (MQ=255)
                                             gTTGCGGCTCAGGCCCAGCAGGCTGGCCGTGTGCACCTGGTTGCAGTGGCAGAAACGATAGGCACTGCGTAGCAGGCTGTCCTCGACCTTCTCGTGCAGGCCGCCGGCCTGCTCCTCGAACAGGCGCTGAAAGGCTCGCAGCAACTGtt                                                                                                         >  3:265706/1‑149 (MQ=255)
                                              ttGCGGCTCAGGCCCAGCAGGCTGGCCGTGTGCACCTGGTTGCAGTGGCAGAAACGATAGGCACTGCGTAGCAGGCTGTCCTCGACCTTCTCGTGCAGGCCGCCGGCCTGCTCCTCGAACAGGCGCTGAAAgg                                                                                                                        <  4:63495/133‑1 (MQ=255)
                                              ttGCGGCTCAGGCCCAGCAGGCTGGCCGTGTGCACCTGGTTGCAGTGGCAGAAACGATAGGCACTGCGTAGCAGGCTGTCCTCGACCTTCTCGTGCAGGCCGCCGGCCTGCTCCTCGAACAGGCGCTGAAAgg                                                                                                                        >  3:63495/1‑133 (MQ=255)
                                                gCGGCTCAGGCCCAGCAGGCTGGCCGTGTGCACCTGGTTGCAGTGGCAGAAACGATAGGCACTGCGTAGCAGGCTGTCCTCGACCTTCTCGTGCAGGCCGCCGGCCTGCTCCTCGAACAGGCGCTGAAAGGCTCGCAGCAACTGttctt                                                                                                      >  6:52182/1‑149 (MQ=255)
                                                gCGGCTCAGGCCCAGCAGGCTGGCCGTGTGCACCTGGTTGCAGTGGCAGAAACGATAGGCACTGCGTAGCAGGCTGTCCTCGACCTTCTCGTGCAGGCCGCCGGCCTGCTCCTCGAACAGGCGCTGAAAGGCTCGCAGCAACTGttctt                                                                                                      >  7:362859/1‑149 (MQ=255)
                                                gCGGCTCAGGCCCAGCAGGCTGGCCGTGTGCACCTGGTTGCAGTGGCAGAAACGATAGGCACTGCGTAGCAGGCTGTCCTCGACCTTCTCGTGCAGGCCGCCGGCCTGCTCCTCGAACAGGCGCTGAAAGGCTCGCAGCAACTGttctt                                                                                                      >  6:152893/1‑149 (MQ=255)
                                                 cGGCTCAGGCCCAGCAGGCTGGCCGTGTGCACCTGGTTGCAGTGGCAGAAACGATAGGCACTGCGTAGCAGGCTGTCCTCGACCTTCTCGTGCAGGCCGCCGGCCTGCTCCTCGAACAGTTGCTGCGAGCCTTTCAGCGCCTGttcgag                                                                                                     <  3:300676/149‑4 (MQ=38)
                                                 cGGCTCAGGCCCAGCAGGCTGGCCGTGTGCACCTGGTTGCAGTGGCAGAAACGATAGGCACTGCGTAGCAGGCTGTCCTCGACCTTCTCGTGCAGGCCGCCGGCCTGCTCCTCGAACAGGCGCTGAAAGGCTCGCAGCAACTGttctt                                                                                                      <  6:390043/148‑1 (MQ=255)
                                                         gCCCAGCAGGCTGGCCGTGTGCACCTGGTTGCAGTGGCAGAAACGATAGGCACTGCGTAGCAGGCTGTCCTCGACCTTCTCGTGCAGGCCGCCGGCCTGCTCCTCGAACAGGTGCTGAAAGGCTCGCAGCAACTGTTCTTCGGCGCttt                                                                                             >  2:184129/1‑149 (MQ=255)
                                                           ccAGCAGGCTGGCCGTGTGCACCTGGTTGCAGTGGCAGAAACGATAGGCACTGCGTAGCAGGCTGTCCTCGACCTTCTCGTGCAGGCCGCCGGCCTGCTCCTCGAACAGGCGCTGAAAGGCTCGCAGCAACTGTTCTTCGGCGCTTTcg                                                                                           >  5:344530/1‑149 (MQ=255)
                                                                 ggCTGGCCGTGTGCACCTGGTTGCAGTGGCAGAAACGATAGGCACTGCGTAGCAGGCTGTCCTCGACCTTCTCGTGCAGGCCGCCGGCCTGCTCCTCGAACAGGCGCTGAAAGGCTCGCAGCAACTGTTCTTCGGCGCTTTCGGCGGTg                                                                                     <  8:272257/149‑1 (MQ=255)
                                                                  gCTGGCCGTGTGCACCTGGTTGCAGTGGCAGAAACGATAGGCACTGCGTAGCAGGCTGTCCTCGACCTTCTCGTGCAGGCCGCCGGCCTGCTCCTCGAACAGGCGCTGAAAGGCTCGCAGCAACTGTTCTTCGGCGCTTTCGGCGGTGc                                                                                    >  1:48993/1‑149 (MQ=255)
                                                                        cGTGTGCACCTGGTTGCAGTGGCAGAAACGATAGGCACTGCGTAGCAGGCTGTCCTCGACCTTCTCGTGCAGGCCGCCGGCCTGCTCCTCGAACAGGCGCTGAAAGGCTCGCAGCAACTGTTCTTCGGCGCTTTCTGCGGTGCTCACCg                                                                              >  3:321786/1‑149 (MQ=255)
                                                                           gtgCACCTGGTTGCAGTGGCACAAACGATAGGCACTGCGTAGCAAGCTGTCCTCGACCTTCTCGTGCAGGCCGCCGGCCTGCTCCTCGAACAGGCGCAGAAAGGCTCGCAGCAACTGTTCTTCGGCGCTTTCGGCGGTg                                                                                     >  1:279732/1‑139 (MQ=255)
                                                                           gagCACCTGGTTGCAGTGGCAGAAACGATAGGCACTGCGTAGCAGGCTGTCCTCGACGTTCTCGTGCAGGCCGCGGGCCTGCTCCTCGAACAGGCGCTGAAAGGCTCGCAGCAACTGTTATTCGGCGCTTTCGGCGGTg                                                                                     <  2:279732/137‑1 (MQ=255)
                                                                               aCCTGGTTGCAGTGGCAGAAACGATAGGCACTGCGTAGCAGGCTGTCCTCGACCTTCTCGTGCAGGCCGCCGGCCTGCTCCTCGAACAGGCGCTGAAAGGCTCGCAGCAACTGTTCTTCGGCGCTTTCGGCGGTGCTCACCGGCGACtt                                                                       >  7:282257/1‑149 (MQ=255)
                                                                                ccTGGTTGCAGTGGCAGAAACGATAGGCACTGCGTAGCAGGCTGTCCTCGACCTTCTCGTGCAGGCCGCCGGCCTGCTCCTCGAACAGGCGCTGAAAGGCTCGCAGCAACTGTTCTTCGGCGCTTTCGGCGGTGCTCACCGGCGACTTg                                                                      <  5:152893/149‑1 (MQ=255)
                                                                                   ggTTGCAGTGGCAGAAACGATAGGCACTGCGTAGCAGGCTGTCCTCGACCTTCTCGTGCAGGCCGCCGGCCTGCTCCTCGAACAGGCGCTGAAAGGCTCGCAGCAACTGTTCTTCGGCGCTTTCGGCGGTGCTCACCGGCGACTTGTcc                                                                   >  5:132384/1‑149 (MQ=255)
                                                                                    gTTGCAGTGGCAGAAACGATAGGCACTGCGTAGCAGGCTGTCCTCGACCTTCTCGTGCAGGCCGCCGGCCTGCTCCTCGAACAGGCGCTGAAAGGCTCGCAGCAACTGTTCTTCGGCGCTTTCGGCGGTGCTCACCGGCGACTTGTCCt                                                                  <  5:100255/149‑1 (MQ=255)
                                                                                     ttGCAGTGGCAGAAACGATAGGCACTGCGTAGCAGGCTGTCCTCGACCTTCTCGTGCAGGCCGCCGGCCTGCTCCTCGAACAGGCGCTGAAAGGCTCGCAGCAACTGTTCTTCGGCGCTTTCGGCGGTGCTCACCGGCGACTTGTCCTg                                                                 >  2:289750/1‑149 (MQ=255)
                                                                                      tgaaGTGGAAGAAACGATATGCACTGCGTAGCAGGCTGTCCTCGACCTTCTCGTGCAGGCCGCCGGCCTGCTCCTCGAACAGGCGCTGAAAGGCTCGCAGCAACTGTTCTTCGGCGCTTTCGGCGGTGCTCACCGGCGACTTGTCCTg                                                                 <  8:349148/145‑1 (MQ=255)
                                                                                               aGAAACGATAGGCACTGCGTAGCAGGCTGTCCTCGACCTTCTCGTGCAGGCCGCCGGCCTGCTCCTCGAACAGGCGCTGAAAGGCTCGCAGCAACTGTTCTTCGGCGCTTTCGGCGGTGCTCACCGGCGACTTGTCCTGGCGCTCG‑ATg                                                      >  3:74851/1‑149 (MQ=255)
                                                                                                 aaaCGATAGGCACTGCGTAGCAGGCTGTCCTCGACCTTCTCGTGCAGGCCGCCGGCCTGCTCCTCGAACAGGCGCTGAAAGGCTCGCAGCAACTGTTCTTCGGCGCTTTCGGCGGTGCTCACCGGCGACTTGTCCTGGCGCTCG‑ATgcg                                                    >  3:95766/1‑149 (MQ=255)
                                                                                                    cGATAGGCACTGCGTAGCAGGCTGTCCTCGACCTTCTCGTGCAGGCCGCCGGCCTGCTCCTCGAACAGGCGCTGAAAGGCTCGCAGCAACTGTTCTTCGGCGCTTTCGGCGGTGCTCACCGGCGACTTGTCCTGGCGCTCG‑ATGCGCAg                                                 <  7:341159/149‑1 (MQ=255)
                                                                                                     gATAGGCACTGCGTAGCAGGCTGTCCTCGACCTTCTCGTGCAGGCCGCCGGCCTGCTCCTCGAACAGGCGCTGAAAGGCTCGCAGCAACTGTTCTTCGGCGCTTTCGGCGGTGCTCACCGGCGACTTGTCCTGGCTCTAG‑ATGCGCAgg                                                >  8:353984/1‑149 (MQ=255)
                                                                                                         ggCACTGCGTAGCAGGCTGTCCTCGACCTTCTCGTGCAGGCCGCCGGCCTGCTCCTCGAACAGGCGCTGAAAGGCTCGCAGCAACTGTTCTTCGGTGCTTTCGGCGGTGCTCACCGGCGACTTGTCCTGGCGCTCG‑ATGCGCAGGTGc                                             <  8:362859/148‑1 (MQ=255)
                                                                                                          gCACTGCGTAGCAGGCTGTCCTCGACCTTCTCGTGCAGGCCGCCGGCCTGCTCCTCGAACAGGCGCTGAAAGGCTCGCAGCAACTGTTCTTCGGCGCTTTCGGCGGTGCTCACCGGCGACTTGTCCTGGCGCTCG‑At                                                       >  7:251119/1‑137 (MQ=255)
                                                                                                          gCACTGCGTAGCAGGCTGTCCTCGACCTTCTCGTGCAGGCCGCCGGCCTGCTCCTCGAACAGGCGCTGAAAGGCTCGCAGCAACTGTTCTTCGGCGCTTTCGGCGGTGCTCACCGGCGACTTGTCCTGGCGCTCG‑At                                                       <  8:251119/137‑1 (MQ=255)
                                                                                                          gCACTGCGTAGCAGGCTGTCCTCGACCTTCTCGTGCAGGCCGCCGGCCTGCTCCTCGAACAGGCGCTGAAAGGCTCGCAGCAACTGTTCTTCGGCGCTTTCGGCGGTGCTCACCGGCGACTTGTCCTGGCGCTCG‑ATGCGCAGGTGCGa                                           >  7:96477/1‑149 (MQ=255)
                                                                                                           cACTGCGTAGCAGGCTGTCCTCGACCTTCTCGTGCAGGCCGCCGGCCTGCTCCTCGAACAGGTGCTGAAAGGCTCGCAGCAACTGTTCTTCGGCGCTTTCGGCGGTGCTCACCGGCGACTTGTCCTGGCGCTCG‑ATGCGCAGGTGCGAc                                          <  1:184129/149‑1 (MQ=255)
                                                                                                                 gTAGCAGGCTGTCCTCGACCTTCTCGTGCAGGCCGCCGGCCTGCTCCTCGAACAGGCGCTGAAAGGCTCGCAGCAACTGTTCTTCGGCGCTTTCGGCGGTGCTCACCGGCGACTTGTCCTGGCGCTCG‑ATGCGCAGGTGCGACAGgcgc                                    <  8:96477/149‑1 (MQ=255)
                                                                                                                    gCAGGCTGTCCTCGACCTTCTCGTGCAGGCCGCCGGCCTGCTCCTCGAACAGGCGCTGAAAGGCTCGCAGCAACTGTTCTTCGGCGCTTTCGGCGGTGCTCACCGGCGACTTGTCCTGGCGCTCG‑ATGCGCAGGTGCGACAGGCGCAg                                  <  6:344530/148‑1 (MQ=255)
                                                                                                                    gCAGGCTGTCCTCGACCTTCTCGTGCAGGCCGCCGGCCTGCTCCTCGAACAGGCGCTGAAAGGCTCGCAGCAACTGTTCTTCGGCGCTTTCGGCGGTGCTCACCGGCGACTTGTCCTGGCACTCG‑ATGCGCAGGTGCGACAGGCGCAgg                                 <  8:282257/149‑1 (MQ=255)
                                                                                                                         cTGTCCTCGACCTTCTCGTGCAGGCCGCCGGCCTGCTCCTCGAACAGGCGCTGAAAGGCTCGCAGCAACTGTTCTTCGGCGCTTTCGGCGGTGCTCACCGGCGACTTGTCCTGGCGCTCG‑ATGCGCAGGTGCGACAGGCGCAGgtcgtc                            <  2:14630/149‑1 (MQ=255)
                                                                                                                           gTCCTCGACCTTCTCGTGCAGGCCGCCGGCCTGCTCCTCGAACAGGCGCTGAAAGGCTCGCAGCAACTGTTCTTCGGCGCTTTCGGCGGTGCTCACCGGCGACTTGTCCTGGCGCTCG‑ATGCGCAGGTGCGACAGGCGCAGgtcgtcgt                          >  4:11944/1‑149 (MQ=255)
                                                                                                                             ccTCGACCTTCTCGTGCAGGCCGCCGGCCTGCTCCTCGAACAGGCGCTGAAAGGCTCGCAGCAACTGTTCTTCGGCGCTTTCGGCGGTGCTCACCGGCGACTTGTCCTGGCGCTCG‑ATGCGCAGGTGCGACAGGCGCAGGTCGTCGTcc                        >  5:32177/1‑149 (MQ=255)
                                                                                                                              cTCGACCTTCTCGTGCAGGCCGCCGGCCTGCTCCTCGAACAGGCGCTGAAAGGCTCGCAGCAACTGTTCTTCGGCGCTTTCGGCGGTGCTCACCGGCGACTTGTCCTGGCGCTCG‑ATGCGCAGGTGCGACAGGCGCAGGTCGTCGTCCt                       >  7:194840/1‑149 (MQ=255)
                                                                                                                              cTCGACCTTCTCGTGCAGGCCGCCGGCCTGCTCCTCGAACAGGCGCTGAAAGGCTCGCAGCAACTGTTCTTCGGCGCTTTCGGCGGTGCTCACCGGCGACTTGTCCTGGCGCTCG‑ATGCGCAGGTGCGACAGGCGCAGGTCGTCGTCCt                       >  2:383699/1‑149 (MQ=255)
                                                                                                                                 gACCTTCTCGTGCAGGCCGCCGGCCTGCTCCTCGAACAGGCGCTGAAAGGCTCGCAGCAACTGTTCTTCGGCGCTTTCGGCGGTGCTCACCGGCGACTTGTCCTGGCGCTCG‑ATGCGCAGGTGCGACAGGCGCAGGTCGTCGTCCTGc                     <  8:194840/148‑1 (MQ=255)
                                                                                                                                 gACCTTCTCGTGCAGGCCGCCGGCCTGCTCCTCGAACAGGCGCTGAAAGGCTCGCAGCAACTGTTCTTCGGCGCTTTCGGCGGTGCTCACCGGCGACTTGTCCTGGCGCTCG‑ATGCGCAGGTGCGACAGGCGCAGGTCGTCGTCCTGCa                    >  3:271960/1‑149 (MQ=255)
                                                                                                                                          gagCAGGCCGGCGGCCTGCTCCTCGAACAGGCGCTGAAAGGCTCGCAGCAACTGTTCTTCGGCGCTTTCGGCGGTGCTCACCGGCGACTTGTCCTGGCGCTCG‑ATGCGCAGGTGCGACAGGCGCAGGTCGTCGTCCTg                      <  6:182071/136‑1 (MQ=255)
                                                                                                                                          gagCAGGCCGGCGGCCTGCTCCTCGAACAGGCGCTGAAAGGCTCGCAGCAACTGTTCTTCGGCGCTTTCGGCGGTGCTCACCGGCGACTTGTCCTGGCGCTCG‑ATGCGCAGGTGCGACAGGCGCAGGTCGTCGTCCTg                      >  5:182071/3‑138 (MQ=255)
                                                                                                                                          gTGCAGGCCGCCGGCCTGCTCCTCGAACAGGCGCTGAAAGGCTCGCAGCAACTGTTCTTCGGCGCTTTCGGCGGTGCTCACCGGCGACTTGTCCTGGCGCTCG‑ATGCGCAGGTGCGACAGGCGCAGGTCGTCGTCCTGCACCAGGCCAt           >  6:325874/1‑149 (MQ=255)
                                                                                                                                          gTGCAGGCCGCCGGCCTGCTCCTCGAACAGGCGCTGAAAGGCTCGCAGCAACTGTTCTTCGGCGCTTTCGGCGGTGCTCACCGGCGACTTGTCCTGGCGCTCG‑ATGCGCAGGTGCGACAGGCGCAGGTCGTCGTCCTGCACCAGGCCAt           >  1:105652/1‑149 (MQ=255)
                                                                                                                                          gTGCAGGCCGCCGGCCTGCTCCTCGAACAGGCGCTGAAAGGCTCGCAGCAACTGTTCTTCGGCGCTTTCGGCGGTGCTCACCGGCGACTTGTCCTGGCGCTCG‑ATGCGCAGGTGCGACAGGCGCAGGTCGTCGTCCTGCACCAGGCCAt           >  3:19060/1‑149 (MQ=255)
                                                                                                                                           tGCAGGCCGCCGGCCTGCTCCTCGAACAGGCGCTGAAAGGCTCGCAGCAACTGTTCTTCGGCGCTTTCGGCGGTGCTCACCGGCGACTTGTCCTGGCGCTCG‑ATGCGCAGGTGCGACAGGCGCAGGTCGTCGTCCTGCACCAGGCCAtt          <  4:265706/149‑1 (MQ=255)
                                                                                                                                             cAGGCCGCCGGCCTGCTCCTCGAACAGGCGCTGAAAGGCTCGCAGCAACTGTTCTTCGGCGCTTTCGGCGGTGCTCACCGGCGACTTGTCCTGGCGCTCG‑ATGCGCAGGTGCGACAGGCGCAGGTCGTCGTCCTGCACCAGGCCATTGc        <  3:11944/149‑1 (MQ=255)
                                                                                                                                             cAGGCCGCCGGCCTGCTCCTCGAACAGGCGCTGAAAGGCTCGCAGCAACTGTTCTTCGGCGCTTTCGGCGGTGCTCACCGGCGACTTGTCCTGGCGCTCG‑ATGCGCAGGTGCGACAGGCGCAGGTCGTCGTCCTGCACCAGGCCATTGc        >  4:359443/1‑149 (MQ=255)
                                                                                                                                              aGGCCGCCGGCCTGCTCCTCGAACAGGCGCTGAAAGGCTCGCAGCAACTGTTCTTCGGCGCTTTCGGCGGTGCTCACCGGCGACTTGTCCTGGCGCTCGCATGCGCAGGTGCGACAGGCGCAGGTCGTCGTCCTGCACCAGGCCATTGc        <  1:289750/149‑1 (MQ=255)
                                                                                                                                              aGGCCGCCGGCCTGCTCCTCGAACAGGCGCTGAAAGGCTCGCAGCAACTGTTCTTCGGCGCTTTCGGCGGTGCTCACCGGCGACTTGTCCTGGCGCTCG‑ATGCGCAGGTGCGACAGGCGCAGGTCGTCGTCCTGCACCAGGCCATTGCg       <  2:48993/149‑1 (MQ=255)
                                                                                                                                                gccgccGGCCTGCTCCTCGAACAGGCGCTGAAAGGCTCGCAGCAACTGTTCTTCGGCGCTTTCGGCGGTGCTCACCGGCGACTTGTCCTGGCGCTCG‑ATGCGCAGGTGCGACAGGCGCAGGTCGTCGTCCTGCACCAGGCCATTGCgg      >  7:139299/1‑148 (MQ=255)
                                                                                                                                                gccgccGGCCTGCTCCTCGAACAGGCGCTGAAAGGCTCGCAGCAACTGTTCTTCGGCGCTTTCGGCGGTGCTCACCGGCGACTTGTCCTGGCGCTCG‑ATGAGCAGGTGCGACAGGCGCAGGTCGTCGTCCTGCACCAGGCCATTGCgg      >  4:294611/1‑148 (MQ=255)
                                                                                                                                                  cgccgGCCTGCTCCTCGAACAGGCGCTGAAAGGCTCGCAGCAACTGTTCTTCGGCGCTTTCGGCGGTGCTCACCGGCGACTTGTCCTGGCGCTCG‑ATGCGCAGGTGCGACAGGCGCAGGTCGTCGTCCTGCACCAGGCCATTGCGGca    >  4:115309/1‑148 (MQ=255)
                                                                                                                                                   gccgGCCTGCTCCTCGAACAGGCGCTGAAAGGCTCGCAGCAACTGTTCTTCGGCGCTTTCGGCGGTGCTCACCGGCGACTTGTCCTGGCGCTCG‑ATGCGCAGGTGCGACAGGCGCAGGTCGTCGTCCTGCACCAGGCCATTGCGGcaca  >  1:125871/1‑149 (MQ=255)
                                                                                                                                                    |                                                                                                                                                    
ATTGACCACCAGTTCGCCGATGGCGATCAGTCGGGCGCGGGTGATGTTGCGGCTCAGGCCCAGCAGGCTGGCCGTGTGCACCTGGTTGCAGTGGCAGAAACGATAGGCACTGCGTAGCAGGCTGTCCTCGACCTTCTCGTGCAGGCCGCCGGCCTGCTCCTCGAACAGGCGCTGAAAGGCTCGCAGCAACTGTTCTTCGGCGCTTTCGGCGGTGCTCACCGGCGACTTGTCCTGGCGCTCG‑ATGCGCAGGTGCGACAGGCGCAGGTCGTCGTCCTGCACCAGGCCATTGCGGCACA  >  NC_002947/3155216‑3155511

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