Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,835,637 A→G 7.3% intergenic (‑79/+53) 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,6370AG 7.3% 182.9 / 4.7 69intergenic (‑79/+53)PP_1639/yohCSprT family protein/inner membrane protein
Reads supporting (aligned to +/- strand):  ref base A (31/32);  new base G (3/2);  total (34/35)
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.82e-01

AGCGGCGTTTGAAAAAGGCTTCGGCTTGCTGGTAACAGGTTTCGACGCGTTGTTTGAGCTGCTCGGGCATGGCGGGGCGGTTCTCCTGACCAGGCATTATGCCGCAAGCGCAGGTGGCTTGCCGAGTCACCGGTCAGCGCACAGCCAGAGGCCGCCTTGCGGCGGCCCCTGGTTGATGCCCCAGTGTTGGTTCGGTCTGTTCTAGTTGGTGTAGACGGGCCCCACGCCCAGTCCCCAGATAATCACGGTAGCGGCCATGATTGCCACCAGTACTACCAAGCCCACTGCCAGCACC  >  NC_002947/1835489‑1835783
                                                                                                                                                    |                                                                                                                                                  
aGCGGCGTTTGAAAAAGGCTTCGGCTTGCTGGTAACAGGTTTCGACGCGTTGTTTGAGCTGCTCGGGCATGGCGGGGCGGTTCTCCTGACCAGGCATTATGCCGCAAGCGCAGGTGGCTTGCCGAGTCACCGGTCAGCGCACAGCCaga                                                                                                                                                    <  7:8954/149‑1 (MQ=255)
 gcggcgTTTGAAAAAGGCTTCGGCTTGCTGGTAACAGGTTTCGACGCGTTGTTTGAGCTGCTCGGGCATGGCGGGGCGGTTCTCCTGACCAGGCATTATGCCGCAAGCGCAGGTGGCTTGCCGAGTCACCGGTCAGCGCACAGCCagag                                                                                                                                                   <  5:242038/149‑1 (MQ=255)
  cggcgTTTGAAAAAGGCTTCGGCTTGCTGGTAACAGGTTTCGACGCGTTGTTTGAGCTGCTCGTGCATGGCGGGGCGGTTCTCCTGACCAGGCATTATGCCTCAAGCGCAGGTGGCTTGCCGAGTCACCGGTCAGCGCACAGCCAGAgg                                                                                                                                                  <  6:286612/149‑1 (MQ=255)
  cggcgTTTGAAAAAGGCTTCGGCTTGCTGGTAACAGGTTTCGACGCGTTGTTTGAGCTGCTCGGGCATGGCGGGGCGGTTCTCCTGACCAGGCATTATGCCGCAAGCGCAGGTGGCTTGCCGAGTCACCGGTCAGCGCACAGCCAGAgg                                                                                                                                                  <  8:261340/149‑1 (MQ=255)
  cggcgTTTGAAAAAGGCTTCGGCTTGCTGGTAACAGGTTTCGACGCGTTGTTTGAGCTGCTCGGGCATGGCGGGGCGGTTCTCCTGACCAGGCATTATGCCGCAAGCGCAGGTGGCTTGCCGAGTCACCGGTCAGCGCACAGCCAGAgg                                                                                                                                                  >  4:267951/1‑149 (MQ=255)
  cggcgTTTGAAAAAGGCTTCGGCTTGCTGGTAACAGGTTTCGACGCGTTGTTTGAGCTGCTCGGGCATGGCGGGGCGGTTCTCCTGACCAGGCATTATGCCGCAAGCGCAGGTGGCTTGCCGAGTCACCGGTCAGCGCACAGCCAGAgg                                                                                                                                                  <  1:49155/149‑1 (MQ=255)
          gAAAAAGGCTTCGGCTTGCTGGTAACAGGTTTCGACGCGTTGTTTGAGCTGCTCGGGCATGGCGGGGCGGTTCTCCTGACCAGGCATTATGCCGCAAGCGCAGGTGGCTTGCCGAGTCACCGGTCAGCGCACAGCCAGAGGCCGCCTTg                                                                                                                                          >  7:386275/1‑149 (MQ=255)
           aaaaaGGCTTCGGCTTGCTGGTAACAGGTTTCGACGCGTTGTTTGAGCTGCTCGGGCATGGCGGGGCGGTTCTCCTGACCAGGCATTATGCCGCAAGCGCAGGTGGCTTGCCGAGTCACCGGTCAGCGCACAGCCAGAGGCCGCCTTgc                                                                                                                                         >  1:173034/1‑149 (MQ=255)
           aaaaaGGCTTCGGCTTGCTGGTAACAGGTTTCGACGCGTTGTTTGAGCTGCTCGGGCATGGCGGGGCGGTTCTCCTGACCAGGCATTATGCCGCAAGCGCAGGTGGCTTGCCGAGTCACCGGTCAGCGCACAGCCAGAGGCCGCCTTgc                                                                                                                                         <  3:19828/149‑1 (MQ=255)
             aaaGGCTTCGGCTTGCTGGTAACAGGTTTCGACGCGTTCTTTGAGCTGCTCGGGCATGGCGGGGCGGTTCTCCTGACCAGGCATTATGCCGCAAGCGCAGGTGGCTTGCCGAGTCACCGGTCAGCGCACAGCCAGGGgccgccgcaagg                                                                                                                                       >  5:246598/1‑143 (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                                                                                                       >  1:312057/1‑130 (MQ=255)
                                                                gggCATGGCGGGGCGGTTCTCCTGACCAGGCATTATGCCGCAAGCGCAGGTGGCTTGCCGAGTCACCGGTCAGCGCACAGCCAGAGGCCGCCTTGCGGCGGCCCCTGGTTGATGCCCCAGTGTTGGTTCg                                                                                                       <  2:312057/130‑1 (MQ=255)
                                                                gggCATGGCGGGGCGGTTCTCCTGACCAGGCATTATGCCGCAAGCGCAGGTGGCTTGCCGAGTCACCGGTCAGCGCACAGCCAGAGGCCGCCTTGCGGCGGCCCCTGGTTGATGCCCCAGTGTTGGTTCGGTCTGTTCTAGTTGgtg                                                                                      <  1:277612/147‑1 (MQ=255)
                                                                gggCATGGCGGGGCGGGTCTCCTTCCCAGGCATTATGCCGCAAGCGCAGGTGGCTTGCCGAGTCACCGGTCAGCGCACAGCCAGAGGCCGCCTTGCGGCGGCCCCTGGTTGATGCCCCAGTGTTGGTTCGGTCTGTTCTAGTTGgtg                                                                                      <  8:272952/147‑1 (MQ=255)
                                                                  gCATGGCGGGGCGGTTCTCCTGACCAGGCATTATGCCGCAAGCGCAGGTGGCTTGCCGAGTCACCGGTCAGCGCACAGCCAGAGGCCGCCTTGCGGCGGCCCCTGGTTGATGCCCCAGTGTTGGTTCGGTCTGTTCTAGTTGGTGTAGa                                                                                  <  8:386275/149‑1 (MQ=255)
                                                                    aTGGCGGGGCGGTTCTCCTGACCAGGCATTATGCCGCAAGCGCAGGTGGCTTGCCGAGTCACCGGTCAGCGCACAGCCAGAGGCCGCCTTGCGGCGGCCCCTGGTTGATGCCCCAGTGTTGGTTCGGTCTGTTCTAGtt                                                                                          >  1:240734/1‑139 (MQ=255)
                                                                    aTGGCGGGGCGGTTCTCCTGACCAGGCATTATGCCGCAAGCGCAGGTGGCTTGCCGAGTCACCGGTCAGCGCACAGCCAGAGGCCGCCTTGCGGCGGCCCCTGGTTGATGCCCCAGTGTTGGTTCGGTCTGTTCTAGtt                                                                                          <  2:240734/139‑1 (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:373658/1‑149 (MQ=255)
                                                                                                      cGCAAGCGCAGGTGGCTTGCCGAGTCACCGGTCAGCGCACAGCCAGAGGCCGCCTTGCGGCGGCCCCTGGTTGATGCCCCAGTGTTGGTTCGGTCTGTTCTAGTTGGTGTAGACGGGCCCCACGCCCAGTCCCCAGATAATCACGGTAg                                              >  1:88725/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)
                                                                                                                                                    |                                                                                                                                                  
AGCGGCGTTTGAAAAAGGCTTCGGCTTGCTGGTAACAGGTTTCGACGCGTTGTTTGAGCTGCTCGGGCATGGCGGGGCGGTTCTCCTGACCAGGCATTATGCCGCAAGCGCAGGTGGCTTGCCGAGTCACCGGTCAGCGCACAGCCAGAGGCCGCCTTGCGGCGGCCCCTGGTTGATGCCCCAGTGTTGGTTCGGTCTGTTCTAGTTGGTGTAGACGGGCCCCACGCCCAGTCCCCAGATAATCACGGTAGCGGCCATGATTGCCACCAGTACTACCAAGCCCACTGCCAGCACC  >  NC_002947/1835489‑1835783

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