Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,835,648 C→A 5.6% intergenic (‑90/+42) PP_1639 ← / ← yohC SprT family protein/inner membrane protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029471,835,6480CA 5.6% 231.6 / 3.4 72intergenic (‑90/+42)PP_1639/yohCSprT family protein/inner membrane protein
Reads supporting (aligned to +/- strand):  ref base C (34/34);  new base A (2/2);  total (36/36)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.08e-01

AAAAAGGCTTCGGCTTGCTGGTAACAGGTTTCGACGCGTTGTTTGAGCTGCTCGGGCATGGCGGGGCGGTTCTCCTGACCAGGCATTATGCCGCAAGCGCAGGTGGCTTGCCGAGTCACCGGTCAGCGCACAGCCAGAGGCCGCCTTGCGGCGGCCCCTGGTTGATGCCCCAGTGTTGGTTCGGTCTGTTCTAGTTGGTGTAGACGGGCCCCACGCCCAGTCCCCAGATAATCACGGTAGCGGCCATGATTGCCACCAGTACTACCAAGCCCACTGCCAGCACCGAGCTGGAGAAC  >  NC_002947/1835500‑1835795
                                                                                                                                                    |                                                                                                                                                   
aaaaaGGCTTCGGCTTGCTGGTAACAGGTTTCGACGCGTTGTTTGAGCTGCTCGGGCATGGCGGGGCGGTTCTCCTGACCAGGCATTATGCCGCAAGCGCAGGTGGCTTGCCGAGTCACCGGTCAGCGCACAGCCAGAGGCCGCCTTgc                                                                                                                                                     >  1:173034/1‑149 (MQ=255)
aaaaaGGCTTCGGCTTGCTGGTAACAGGTTTCGACGCGTTGTTTGAGCTGCTCGGGCATGGCGGGGCGGTTCTCCTGACCAGGCATTATGCCGCAAGCGCAGGTGGCTTGCCGAGTCACCGGTCAGCGCACAGCCAGAGGCCGCCTTgc                                                                                                                                                     <  3:19828/149‑1 (MQ=255)
        ttCGGCTTGCTGGTAACAGGTTTCGACGCGTTGTTTGAGCTGCTCGGGCATGGCGGGGCGGTTCTCCTGACCAGGCATTATGCCGCAAGCGCAGGTGGCTTGCCGAGTCACCGGTCAGCGCACAGCCAGAGGCCGCCTTGCGGCGGccc                                                                                                                                             >  3:113012/1‑149 (MQ=255)
           ggCTTGCTGGTAACAGGTTTCGACGCGTTGTTTGAGCTGCTCGGGCATGGCGGGGCGGTTCTCCTGACCAGGCATTATGCCGCAAGCGCAGGTGGCTTGCCGAGTCACCGGTCAGCGCACAGCCAGAGGCCGCCTTGCGGCGGCCCCt                                                                                                                                           <  8:40823/148‑1 (MQ=255)
           ggCTTGCTGGTAACAGGTTTCGACGCGTTGTTTGAGCTGCTCGGGCATGGCGGGGCGGTTCTCCTGACCAGGCATTATGCCGCAAGCGCAGGTGGCTTGCCGAGTCACCGGTCAGCGCACAGCCAGAGGCCGCCTTGCGGCGGCCCCt                                                                                                                                           <  3:267951/148‑1 (MQ=255)
            gCTTGCTGGTAACAGGTTTCGACGCGTTGTTTGAGCTGCTCGGGCATGGCGGGGCGGTTCTCCTGACCAGGCATTATGCCGCAAGCGCAGGTGGCTTGCCGAGTCACCGGTCAGCGCACAGCCAGGGGCCGCCGCAAGGCGGCCTCTgg                                                                                                                                         <  6:246598/149‑1 (MQ=255)
            gCTTGCTGGTAACAGGTTTCGACGCGTTGTTTGAGCTGCTCGGGCATGGCGGGGCGGTTCTCCTGACCAGGCATTATGCCGCAAGCGCAGGTGGCTTGCCGAGTCACCGGTCAGCGCACAGCCAGAGGCCGCCTTGCGGCGGCCCCTgg                                                                                                                                         <  1:18291/149‑1 (MQ=255)
             cTTGCTGGTAACAGGTTTCGACGCGTTGTTTGAGCTGCTCGGGCATGGCGGGGCGGTTCTCCTGACCAGGCATTATGCCGCAAGCGCAGGTGGCTTGCCGAGTCACCGGTCAGCGCACAGCCAGAGGCCGCCTTGCGGCGGCCCCTGGt                                                                                                                                        <  5:424414/149‑1 (MQ=255)
              gtGCTGGTCACAGGTCTCGACTCGGTGTTTGAGCTGCTCGGGCATGGCGGGGCGGTTCTCCTGACCAGGCATTATGCCGCAAGCACAGGTGGCTTGCCGAGTCACCGGTCAGCGCACAGCCAGAGGCCGCCTTGCGGCGGCCCCTGGtt                                                                                                                                       <  8:180716/148‑1 (MQ=255)
                gCTGGTAACAGGGTTCGACGCGTTGTTTGAGCTGCTCGGGCAGGGCGGGGCGGTTCTCCTGACCAGGCATTATGCCGCAAGCGCAGGTGGCTTGCCGAGTCACCGGTCAGCGCACAGCCAGGGGCCGCCGCAAGGCGGCCTCTggct                                                                                                                                       <  7:200703/147‑3 (MQ=255)
                    gTAACAGGTTTCGACGCGTTGTTTGAGCTGCTCGGGCATGGCGGGGCGGTTCTCCTGACCAGGCATTATGCCGCAAGCGCAGGTGGCTTGCCGAGTCACCGGTCAGCGCACAGCCAGAGGCCGCCTTGCGGCGGCCCCTGGTTGATGcc                                                                                                                                 <  1:236890/149‑1 (MQ=255)
                    gTAACAGGTTTCGACGCGTTGTTTGAGCTGCTCGGGCATGGCGGGGCGGTTCTCCTGACCAGGCATTATGCCGCAAGCGCAGGTGGCTTGCCGAGTCACCGGTCAGCGCACAGCCAGAGGCCGCCTTGCGGCGGCCCCTGGTTGATGc                                                                                                                                  <  4:173448/148‑1 (MQ=255)
                      aaCAGGTTTCGACGCGTTGTTTGAGCTGCTCGGGCATGGCGGGGCGGTTCTCCTGACCAGGCATTATGCCGCAAGCGCAGGTGGCTTGCCGAGTCCCCGGTCAGCGCACAGCCAGAGGCCGCCTTGCGGCGGCCCCTGGTTGATGccc                                                                                                                                >  6:210165/1‑148 (MQ=255)
                         aGGTTTCGACGCGTTGTTTGAGCTGCTCGGGCATGGCGGGGCGGTTCTCCTGACCAGGCATTATGCCGCAAGCGCAGGTGGCTTGCCGAGTCACCGGTCAGCGCACAGCCAGAGGCCGCCTTGCGGCGGCCCCTGGTTGATGCCCCAgt                                                                                                                            >  5:175428/1‑149 (MQ=255)
                                gACGCGTTGTTTGAGCTGCTCGGGCATGGCGGGGCGGTTCTCCTGACCAGGCATTATGCCGCAAGCGCAGGTGGCTTGCCGAGTCACCGGTCAGCGCACAGCCAGAGGCCGCCTTGCGGCGGCCCCTGGTTGATGCCCCAGTGTTGGtt                                                                                                                     >  1:261899/1‑149 (MQ=255)
                                gACGCGTTGTTTGAGCTGCTCGGGCATGGCGGGGCGGTTCTCCTGACCAGGCATTATGCCGCAAGCGCAGGTGGCTTGCCGAGTCACCGGTCAGCGCACAGCCAGAGGCCGCCTTGCGGCGGCCCCTGGTTGATGCCCCAGTGTTGGtt                                                                                                                     >  1:188436/1‑149 (MQ=255)
                                 aCGCGTTGTTTGAGCTGCTCGGGCATGGCGGGGCGGTTCTCCTGACCAGGCATTATGCCGCAAGCGCAGGTGGCTTGCCGAGTCACCGGTCAGCGCACAGCCAGAGGCCGCCTTGCGGCGGCCCCTGGTTGATCCCCCAGTGTTggcgc                                                                                                                    >  8:371739/1‑146 (MQ=255)
                                  cgcgTTGTTTGAGCTGCTCGGGCATGGCGGGGCGGTTCTCCTGACCAGGCATTATGCCGCAAGCGCAGGTGGCTTGCCGAGTCACCGGTCAGCGCACAGCCAGAGGCCGCCTTGCGGCGGCCCCTGGTTGATGCCCCAGTGTTGGTTCg                                                                                                                   <  3:56941/149‑1 (MQ=255)
                                   gcgTTGTTTGAGCTGCTCGGGCATGGCGGGGCGGTTCTCCTGACCAGGCATTATGCCGCAAGCGCAGGTGGCTTGCCGAGTCACCGGTCAGCGCACAGCCAGAGGCCGCCTTGCGGCGGCCCCTGGTTGa                                                                                                                                     <  6:243532/130‑1 (MQ=255)
                                   gcgTTGTTTGAGCTGCTCGGGCATGGCGGGGCGGTTCTCCTGACCAGGCATTATGCCGCAAGCGCAGGTGGCTTGCCGAGTCACCGGTCAGCGCACAGCCAGAGGCCGCCTTGCGGCGGCCCCTGGTTGa                                                                                                                                     >  5:243532/1‑130 (MQ=255)
                                     gttgttTGAGCTGCTCGGGCATGGCGGGGCGGTTCTCCTGACCAGGCATTATGCCGCAAGCGCAGGTGGCTTGCCGAGTCACCGGTCAGCGCACAGCCAGAGGCCGCCTTGCGGCGGCCCCTGGTTGATGCCCCAGTGTTGGTTCGGTc                                                                                                                >  3:257428/1‑149 (MQ=255)
                                     gttgttTGAGCTGCTCGGGCATGGCGGGGCGGTTCTCCTGACCAGGCATTATGCCGCAAGCGCAGGTGGCTTGCCGAGTCACCGGTCAGCGCACAGCCAGAGGCCGCCTTGCGGCGGCCCCTGGTTGATGCCCCAGTGTTGGTTCGGTc                                                                                                                >  7:223365/1‑149 (MQ=255)
                                             aGCTGCTCGGGCATGGCGGGGCGGTTCTCCTGACCAGGCATTATGCCGCAAGCGCAGGTGGCTTGCCGAGTCACCGGTCAGCGCACAGCCAGAGGCCGCCTTGCGGCGGCCCCTGGTTGATGCCCCAGTGTTGGTTCGGTCTGTTCTa                                                                                                         <  6:276994/148‑1 (MQ=255)
                                                  ctCGGGCATGGCGGGGCGGTTCTCCTGACCAGGCATTATGCCGCAAGCGCAGGTGGCTTGCCGAGTCACCGGTCAGCGCACAGCCAGAGGCCGCCTTGCGGCGGCCCCTGGTTGATGCCCCAGTGTTGGTTCGGTCTGTTCTAGTTgg                                                                                                    <  4:113012/148‑1 (MQ=255)
                                                     gggCATGGCGGGGCGGTTCTCCTGACCAGGCATTATGCCGCAAGCGCAGGTGGCTTGCCGAGTCACCGGTCAGCGCACAGCCAGAGGCCGCCTTGCGGCGGCCCCTGGTTGATGCCCCAGTGTTGGTTCg                                                                                                                   <  2:312057/130‑1 (MQ=255)
                                                     gggCATGGCGGGGCGGTTCTCCTGACCAGGCATTATGCCGCAAGCGCAGGTGGCTTGCCGAGTCACCGGTCAGCGCACAGCCAGAGGCCGCCTTGCGGCGGCCCCTGGTTGATGCCCCAGTGTTGGTTCg                                                                                                                   >  1:312057/1‑130 (MQ=255)
                                                     gggCATGGCGGGGCGGTTCTCCTGACCAGGCATTATGCCGCAAGCGCAGGTGGCTTGCCGAGTCACCGGTCAGCGCACAGCCAGAGGCCGCCTTGCGGCGGCCCCTGGTTGATGCCCCAGTGTTGGTTCGGTCTGTTCTAGTTGgtg                                                                                                  <  1:277612/147‑1 (MQ=255)
                                                     gggCATGGCGGGGCGGGTCTCCTTCCCAGGCATTATGCCGCAAGCGCAGGTGGCTTGCCGAGTCACCGGTCAGCGCACAGCCAGAGGCCGCCTTGCGGCGGCCCCTGGTTGATGCCCCAGTGTTGGTTCGGTCTGTTCTAGTTGgtg                                                                                                  <  8:272952/147‑1 (MQ=255)
                                                       gCATGGCGGGGCGGTTCTCCTGACCAGGCATTATGCCGCAAGCGCAGGTGGCTTGCCGAGTCACCGGTCAGCGCACAGCCAGAGGCCGCCTTGCGGCGGCCCCTGGTTGATGCCCCAGTGTTGGTTCGGTCTGTTCTAGTTGGTGTAGa                                                                                              <  8:386275/149‑1 (MQ=255)
                                                         aTGGCGGGGCGGTTCTCCTGACCAGGCATTATGCCGCAAGCGCAGGTGGCTTGCCGAGTCACCGGTCAGCGCACAGCCAGAGGCCGCCTTGCGGCGGCCCCTGGTTGATGCCCCAGTGTTGGTTCGGTCTGTTCTAGtt                                                                                                      <  2:240734/139‑1 (MQ=255)
                                                         aTGGCGGGGCGGTTCTCCTGACCAGGCATTATGCCGCAAGCGCAGGTGGCTTGCCGAGTCACCGGTCAGCGCACAGCCAGAGGCCGCCTTGCGGCGGCCCCTGGTTGATGCCCCAGTGTTGGTTCGGTCTGTTCTAGtt                                                                                                      >  1:240734/1‑139 (MQ=255)
                                                         aTGGCGGGGCGGTTCTCCTGACCAGGCATTATGCCGCAAGCGCAGGTGGCTTGCCGAGTCACCGGTCAGCGCACAGCCAGAGGCCGCCTTGCGGCGGCCCCTGGTTGATGCCCCAGTGTTGGTTCGGTCTGTTCTAGTTGGTGTAGACg                                                                                            >  1:223089/1‑149 (MQ=255)
                                                              ggggcggTTCTCCTGACCAGGCATTATGCCGCAAGCGCAGGTGGCTTGCCGAGTCACCGGTCAGCGCACAGCCAGAGGCCGCCTTGCGGCGGCCCCTGGTTGATGCCCCAGTGTTGGTTCGGTCTGTTCTAGTTGGTGTAGACGGGccc                                                                                       <  8:26134/149‑1 (MQ=255)
                                                              ggggcggTTCTCCTGACCAGGCATTATGCCGCAAGCGCAGGTGGCTTGCCGAGTCACCGGTCAGCGCACAGCCAGAGGCCGCCTTGCGGCGGCCCCTGGTTGATGCCCCAGTGTTGGTTCGGTCTGTTCTAGTTGGTGTAGACGGGccc                                                                                       <  2:223716/149‑1 (MQ=255)
                                                                    gTTCTCCTGACCAGGCATTATGCCGCAAGCGCAGGTGGCTTGCCGAGTCACCGGTCAGCGCACAGCCAGAGGCCGCCTTGCGGCGGCCCCTGGTTGATGCCCCAGTGTTGGTTCGGTCTGTTCTAGTTGGTGTAGACGGGCCCCACGcc                                                                                 <  2:261899/149‑1 (MQ=255)
                                                                      tctcCTGACCAGGCATTATGCCGCAAGCGCAGGTGGCTTGCCGAGTCACCGGTCAGCGCACAGCCAGAGGCCGCCTTGCGGCGGCCCCTGGTTGATGCCCCAGTGTTGGTTCGGTCTGTTCTAGTTGGTGTAGACGGGCCCCACGcc                                                                                 >  3:258837/1‑147 (MQ=255)
                                                                             aCCAGGCATTATGCCGCAAGCGTAGGTGGATTGCCGAGTCACCGGTCAGCGCACAGCCAGAGGCCGCATTGCGGCGGCCCATGGTTGATGCCCCAGTGTTGGTTCGGTGTGTTGTAGTTGGTGTAGACGGGCCCCACGCCCAGTccc                                                                          >  7:65145/1‑147 (MQ=255)
                                                                                 ggCATTATGCCGCAAGCGCAGGTGGCTTGCCGAGTCACCGGTCAGCGCACAGCCAGAGGCCGCCTTGCGGCGGCCCCTGGTTGATGCCCCAGTGTTGGTTCGGTCTGTTCTAGTTGGTGTAGACGGGCCCCACGCCCAGTCCCCAGATa                                                                    <  1:289394/149‑1 (MQ=255)
                                                                                    aTTATGCCGCAAGCGCAGGTGGCTTGCCGAGTCACCGGTCAGCGCACAGCCAGAGGCCGCCTTGCGGCGGCCCCTGGTTGATGCCCCAGTGTTGGTTCGGTCTGTTCTAGTTGGTGTAGACGGGCCCCACGCCCAGTCCCCAGATAAt                                                                  >  5:330516/1‑148 (MQ=255)
                                                                                       aTGCCGCAAGCGCAGGTGGCTTGCCGAGTCACCGGTCAGCGCACAGCCAGAGGCCGCCTTGCGGCGGCCCCTGGTTGATGCCCCAGTGTTGGTTCGGTCTGTTCTAGTTGGTGTAGACGGGCCCCCCGCCCAGTCCCCAGATAATCACg                                                              >  7:307596/1‑149 (MQ=255)
                                                                                        tGCCGCAAGCGCAGGTGGCTTGCCGAGTCACCGGTCAGCGCACAGCCAGAGGCCGCCTTGCGGCGGCCCCTGGTTGATGCCCCAGTGTTGGTTCGGTCTGTTCTAGTTGGTGTAGACGGGCCCCACGCCCAGTCCCCAGATAATCACgg                                                             >  8:281097/1‑149 (MQ=255)
                                                                                           cGCAAGCGCAGGTGGCTTGCCGAGTCACCGGTCAGCGCACAGCCAGAGGCCGCCTTGCGGCGGCCCCTGGTTGATGCCCCAGTGTTGGTTCGGTCTGTTCTAGTTGGTGTAGACGGGCCCCACGCCCAGTCCCCAGATAATCACGGTAg                                                          >  1:88725/1‑149 (MQ=255)
                                                                                           cGCAAGCGCAGGTGGCTTGCCGAGTCACCGGTCAGCGCACAGCCAGAGGCCGCCTTGCGGCGGCCCCTGGTTGATGCCCCAGTGTTGGTTCGGTCTGTTCTAGTTGGTGTAGACGGGCCCCACGCCCAGTCCCCAGATAATCACGGTAg                                                          >  1:373658/1‑149 (MQ=255)
                                                                                             cAAGCGCAGGTGGCTTGCCGAGTCACCGGTCAGCGCACAGCCAGAGGCCGCCTTGCGGCGGCCCCTGGTTGATGCCCCAGTGTTGGTTCGGTCTGTTCTAGTTGGTGTAGACGGGCCCCACGCCCAGTCCCCAGATAATCACGGTAGCg                                                        <  2:88725/149‑1 (MQ=255)
                                                                                              aaGCGCAGGTGGCTTGCCGAGTCACCGGTCAGCGCACAGCCAGAGGCCGCCTTGCGGCGGCCCCTGGTTGATGCCCCAGTGTTGGTTCGGTCTGTTCTAGTTGGTGTAGACGGGCCCCACGCCCAGTCCCCAGATAATCACGGTAGCgg                                                       >  5:316220/1‑149 (MQ=255)
                                                                                                gcgcAGGTGGCTTGCCGAGTCACCGGTCAGCGCACAGCCAGAGGCCGCCTTGCGGCGGCCCCTGGTTGATGCCCCAGTGTTGGTTCGGTCTGTTCTAGTTGGTGTAGACGGGCCCCACGCCCAGTCCCCAGATAATCACGGTAGCGGc                                                      >  3:303568/1‑148 (MQ=255)
                                                                                                  gcAGGTGGCTTGCCGAGTCACCGGTCAGCGCACAGCCAGAGGCCGCCATGCCGCGGCCCCTGGTTGATCCCCCAGTGGTGGTTCGGTCAGATCAAGTTGGTGTAGACGGGCCCCACCCCCAGTCCCCAGATAATCACCCTAGCGGcc                                                     >  6:37211/1‑147 (MQ=255)
                                                                                                    aGGTGGCTTGCCGAGTCACCGGTCAGCGCACAGCCAGAGGCCGCCTTGCGGCGGCCCCTGGTTGATGCCCCAGTGTTGGTTCGGTCTGTTCTAGTTGGTGTAGACGGGCCCCACGCCCAGTCCCCAGATAATCACGGTAGCGGCCATGa                                                 <  2:223089/149‑1 (MQ=255)
                                                                                                      gTGGCTTGCCGAGTCACCGGTCAGCGCACAGCCAGAGGCCGCCTTGCGGCGGCCCCTGGTTGATGCCCCAGTGTTGGTTCGGTCTGTTCTAGTTGGTGTAGACGGGCCCCACGCCCAGTCCCCAGATAATCACGGTAGCGGCCATGAtt                                               <  6:330516/149‑1 (MQ=255)
                                                                                                      gTGGCTTGCCGAGTCACCGGTCAGCGCACAGCCAGAGGCCGCCTTGCGGCGGCCCCTGGTTGATGACCCAGTGTTGGTTCGGTCTGTTCTAGTTGGTGTAGACGGGCCCCACGCCCAGTCCCCAGATAATCACGGTAGCGGCCATGAtt                                               <  7:281097/149‑1 (MQ=255)
                                                                                                            tGCCGAGTCACCGGTCAGCGCACAGCCAGAGGCCGCCTTGCGGCGGCCCCTGGTTGATGCCCCAGTGTTGGTTCGGTCTGTTCTAGTTGGTGTAGACGGGCCCCACGCCCAGTCCCCAGATAATCACGGTAGCGGCCATGATTGccac                                          >  7:160781/1‑148 (MQ=255)
                                                                                                                gAGTCACCGGTCAGCGCACAGCCAGAGGCCGCCTTGCGGCGGCCCCTGGTTGATGCCCCAGTGTTGGTTCGGTCTGTTCTAGTTGGTGTAGACGGGCCCCACGCCCAGTCCCCAGATAATCACGGTAGCGGCCATGATTGCCACCAGta                                     >  8:117095/1‑149 (MQ=255)
                                                                                                                   tCACCGGTCAGCGCACAGCCAGAGGCCGCCTTGCGGCGGCCCCTGGTTGATGCCCCAGTGTTGGTTCGGTCTGTTCTAGTTGGTGTAGACGGGCCCCACGCCCAGTCCCCAGATAATCACGGTAGCGGCCATGATTGCCACCAGtacta                                  <  6:175428/149‑1 (MQ=255)
                                                                                                                     aCCGGTCAGCGCACAGCCAGAGGCCGCCTTGCGGCGGCCCCTGGTTGATGCCCCAGTGTTGGTTCGGTCTGTTCTAGTTGGTGTAGACGGGCCCCACGCCCAGTcccc                                                                         <  2:308372/108‑1 (MQ=255)
                                                                                                                     aCCGGTCAGCGCACAGCCAGAGGCCGCCTTGCGGCGGCCCCTGGGTGATGCCCCAGTGTTGGTTCGGTCTGTTCTAGTTGGTGTAGACGGGCCCCACGCCCAGTcccc                                                                         >  1:308372/1‑108 (MQ=255)
                                                                                                                      ccGGTCAGCGCACAGCCAGAGGCCGCCTTGCGGCGGCCCCTGGTTGATGCCCCAGTGTTGGTTCGGTCTGTTCTAGTTGGTGTAGACGGGCCCCACGCCCAGTCCCCAGATAATCACGGTAGCGGCCATGATTGCCACCAGTACTAcc                                >  4:233504/1‑148 (MQ=255)
                                                                                                                       cGGTCAGCGCACAGCCAGAGGCCGCCTTGCGGCGGCCCCTGGTTGATGCCCGAGTGTTGGTTCGGTCTGTTCTAGTTGGTGTAGACGGGCCCCACGCCCAGTCCCCAGATAATCACGGTAGCGGCCATGATTGCCACCAGTACTACCaa                              <  2:188436/149‑1 (MQ=255)
                                                                                                                       cGGTCAGCGCACAGCCAGAGGCCGCCTTGCGGCGGCCCCTGGTTGATGCCCCAGTGTTGGTTCGGTCTGTTCTAGTTGGTGTAGACGGGCCCCACGCCCAGTCCCCAGATAATCACGGTAGCGGCCATGATTGCCACCAGTACTAcc                                >  1:388784/1‑147 (MQ=255)
                                                                                                                       cGGTCAGCGCACAGCCAGAGGCCGCCTTGCGGCGGCCCCTGGTTGATGCCCCAGTGTTGGTTCGGTCTGTTCTAGTTGGTGTAGACGGGCCCCACGCCCAGTCCCCAGATAATCACGGTAGCGGCCATGATTGCCACCAGTACTACCaa                              <  2:373658/149‑1 (MQ=255)
                                                                                                                            aGCGCACAGCCAGAGGCCGCCTTGCGGCGGCCCCTGGTTGATGCCCCAGTGTTGGTTCGGTCTGTTCTAGTTGGTGTAGACGGGCCCCACGCCCAGTCCCCAGATAATCACGGTAGCGGCCATGATTGCCACCAGTACTACCAAGCCCa                         >  5:114518/1‑149 (MQ=255)
                                                                                                                              cgcACAGCCAGCGGCCGCCTTGCGGCGGCCCCTGGTTGATGCCCCAGTGTTGGTTCGGTCTGTTCTAGTTGGTGTAGACGGGCCCCACGCCCAGTCCCCAGATAATCACGGTAGCGGCCATGATTGCCACCAGTACTACCAAGCCCACt                       <  4:258837/149‑1 (MQ=255)
                                                                                                                                    accAGGGGCCGCCGCAAGGCGGCCTCTGGTTGATGCCCCAGTGTTGGTTCGGTCTGTTCTAGTTGGTGTAGACGGGCCCCACGCCCAGTCCCCAGATAATCACGGTAGCGGCCATGATTGCCACCAGTACTACCAAGCCCACTGCCAg                  >  3:60573/2‑148 (MQ=255)
                                                                                                                                     ccAGGGGCCGCCGCAAGGCGGCCTCTGGTTGATGCCCCAGTGTTGGTTCGGTCTGTTCTAGTAGATGTAGACGGGCCCCACGCCCAGTCCCCAGATAATCACGGTAGCGGCCATGATTTCCACCAGTACTACCAAGCCCACTGCCAGCa                >  1:131184/1‑149 (MQ=255)
                                                                                                                                       agagGCCGCCTTGCGGCGGCCCCTGGTTGATGCCCCAGTGTTGGTTCGGTCTGTTCTAGTTGGTGTAGACGGGCCCCACGCCCAGTCCCCAGATAATCACGGTAGCGGCCATGATTGCCACCAGTACTACCAAGCCCACTGCCAGCAcc              >  4:166117/1‑149 (MQ=255)
                                                                                                                                          ggCCGCCTTGCGGCGGCCCCTGGTTGATGCCCCAGTGTTGGTTCGGTCTGTTCTAGTTGGTGTAGACGGGCCCCACGCCCAGTCCCCAGATAATCACGGTAGCGGCCATGATTGccac                                          <  2:200591/118‑1 (MQ=255)
                                                                                                                                          ggCCGCCTTGCGGCGGCCCCTGGTTGATGCCCCAGTGTTGGTTCGGTCTGTTCTAGTTGGTGTAGACGGGCCCCACGCCCAGTCCCCAGATAATCACGGTAGCGGCCATGATTGccac                                          >  1:200591/1‑118 (MQ=255)
                                                                                                                                           gccgccTTGCGGCGGCCCCTGGTTGATGCCCCAGTGTTGGTTCGGTCTGTTCTAGTTGGTGTAGACGGGCCCCACGCCCAGTCCCCAGATAATCACGGTAGCGGCCATGATTGCCACCAGTACTACCAAGCCCACTGcc                    <  4:145627/139‑1 (MQ=255)
                                                                                                                                           gccgccTTGCGGCGGCCCCTGGTTGATGCCCCAGTGTTGGTTCGGTCTGTTCTAGTTGGTGTAGACGGGCCCCACGCCCAGTCCCCAGATAATCACGGTAGCGGCCATGATTGCCACCAGTACTACCAAGCCCACTGcc                    >  3:145627/1‑139 (MQ=255)
                                                                                                                                           gccgccTTGCGGCGGCCCCTGGTTGATGCCCCAGTGTTGGTTCGGTCTGTTCTAGTTGGTGTAGACGGGCCCCACGCCCAGTCCCCAGATAATCACGGTAGCGGCCATGATTGCCACCAGTACTACCAAGCCCACTGCCAGCACCGAGc          >  8:38737/1‑149 (MQ=255)
                                                                                                                                             cgccTTGCGGCGGCCCCTGGTTGATGCCCCAGTGTTGGTTCGGTCTGTTCTAGTTGGTGTAGACGGGCCCCACGCCCAGTCCCCAGATAATCACGGTAGCGGCCATGATTGCCACCAGTACTACCAAGCCCACTGCCAg                  >  7:236433/1‑139 (MQ=255)
                                                                                                                                             cgccTTGCGGCGGCCCCTGGTTGATGCCCCAGTGTTGGTTCGGTCTGTTCTAGTTGGTGTAGACGGGCCCCACGCCCAGTCCCCAGATAATCACGGTAGCGGCCATGATTGCCACCAGTACTACCAAGCCCACTGCCAg                  <  8:236433/139‑1 (MQ=255)
                                                                                                                                              gccTTGCGGCGGCCCCTGGTTGATGCCCCAGTGTTGGTTCGGTCTGTTCTAGTTGGTGTAGACGGGCCCCACGCCCAGTCCCCAGATAATCACGGTAGCGGCCATGATTGCCACCAGTACTACCAAGCCCACTGCCAGCACCGAGCTgg       <  8:65145/149‑1 (MQ=255)
                                                                                                                                              gccTTGCGGCGGCCCCTGGTTGATGCCCCAGTGTTGGTTCGGTCTGTTCTAGTTGGTGTAGACGGGCCCCACGCCCAGTCCCCAGATAATCACGGTAGCGGCCATGATTGCCACCAGTACTACCAAGCCCACTGCCAGCACCGAGCTg        >  1:259450/1‑148 (MQ=255)
                                                                                                                                                cTTGCGGCGGCCCCTGGTTGATGCCCCAGTGTTGGTTCGGTCTGTTCTAGTTGGTGTAGACGGGCCCCACGCCCAGTCCCCAGATAATCACGGTAGCGGCCATGATTGCCACCAGTACTACCAAGCCCACTGCCAGCACCGAGCTGgag     <  3:166117/149‑1 (MQ=255)
                                                                                                                                                    cggcggCCCCTGGTTGATGCCCCAGTGTTGGTTCGGTCTGTTCTAGTTGGTGTAGACGGGCCCCACGCCCAGTCCCCAGATAATCACGGTAGCGGCCATGATTGCCACCAGTACTACCAAGCCCACTGCCAGCTCGGTGCTggcagtg  <  2:388784/148‑6 (MQ=255)
                                                                                                                                                    |                                                                                                                                                   
AAAAAGGCTTCGGCTTGCTGGTAACAGGTTTCGACGCGTTGTTTGAGCTGCTCGGGCATGGCGGGGCGGTTCTCCTGACCAGGCATTATGCCGCAAGCGCAGGTGGCTTGCCGAGTCACCGGTCAGCGCACAGCCAGAGGCCGCCTTGCGGCGGCCCCTGGTTGATGCCCCAGTGTTGGTTCGGTCTGTTCTAGTTGGTGTAGACGGGCCCCACGCCCAGTCCCCAGATAATCACGGTAGCGGCCATGATTGCCACCAGTACTACCAAGCCCACTGCCAGCACCGAGCTGGAGAAC  >  NC_002947/1835500‑1835795

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: