Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,995,018 T→A 7.1% intergenic (+146/+31) PP_5251 → / ← PP_5252 transporter/hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,995,0180TA 7.1% 258.3 / 7.9 85intergenic (+146/+31)PP_5251/PP_5252transporter/hypothetical protein
Reads supporting (aligned to +/- strand):  ref base T (38/41);  new base A (3/3);  total (41/44)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.33e-01

CCGGCCCTTACGGCCCATTCGCGGGTGAACCCGCTCCTACACAGATCTGCGGTAGACACACACCGGACGTTGGAAACAGGCCGTTGTAGGAGCGGGTTTACCCGCGAAAGGGCCGGCACTGAAGACACAAGAAAGGGGCCAGTTATTGGCCCCCGTTCGTCACATCAGCTCAGCTTAGAAAGCGAAGCACGAACAGCCCAGCGCACCCCAGAAGCCCTGGAAGTCATTTACCGGCACGCTCGAATGCCGCGCCTTTTCGTGGCTGTGCGCATGCACACCAC  >  NC_002947/5994875‑5995155
                                                                                                                                               |                                                                                                                                         
gcGGCCCTTACGGCCCATTCGCGGGTCAACCGGCTCCTACTCAGTTCTGCGGTAGACACACACCGGACGTTGGAAACAGCCCGTTGTAGGAGCGGGTTTACCCGCGAAAGGGCCGGCACTGAAGACACAAGAAAGGGGCCAGTTATTg                                                                                                                                       <  3:370399/147‑1 (MQ=37)
 cGGCCCTTACGGCCCATTCGCGGGTGAACCCGCTCCTACACAGATCTGCGGTAGACACACACCGGACGTTGGAAACAGGCCGTTGTAGGAGCGGGTTTACCCGCGAAAGGGCCGGCACTGAAGACACAAGAAAGGGGCCAGTTATTGGc                                                                                                                                     >  1:263832/1‑149 (MQ=255)
   gCCCTTACGGCCCATTCGCGGGTGAACCCGCTCCTACACAGATCTGCGGTAGACACACACCGGACGTTGGAAACAGGCCGTTGTAGGAGCGGGTTTACCCGCGAAAGGGCCGGCACTGAAGACACAAGAAAGGGGCCAGTTATTGGccc                                                                                                                                   >  6:255221/1‑149 (MQ=255)
   gCCCTTACGGCCCATTCGCGGGTGAACCCGCTCCTACACAGATCTGCGGTAGACACACACCGGACGTTGGAAACAGGCCGTTGTAGGAGCGGGTTTACCCGCGAAAGGGCCGGCACTGAAGACACAAGAAAGGGGCCAATAACTGGccc                                                                                                                                   >  4:332632/1‑149 (MQ=255)
     ccTTACGGCCCATTCGCGGGTGAACCCGCTCCTACACAGATCTGCGGTAGACACACACCGGACGTTGGAAACAGGCCGTTGTAGGAGCGGGTTTACCCGCGAAAGGGCCGGCACTGAAGACACAAGAAAGGGGCCAGTTATTGGccccc                                                                                                                                 >  6:69538/1‑149 (MQ=255)
     ccTTACGGCCCATTCGCGGGTGAACCCGCTCCTACACAGATCTGCGGTAGACACACACCGGACGTTGGAAACAGGCCGTTGTAGGAGCGGGTTTACCCGCGAAAGGGCCGGCACTGAAGACACAAGAAAGGGGCCAATAACTGGcccct                                                                                                                                 <  3:332632/149‑2 (MQ=255)
       ttACGGCCCATTCGCGGGTGAACCCGCTCCTACACAGATCTGCGGTAGACACACACCGGACGTTGGAAACAGGCCGTTGTAGGAGCGGGTTTACCCGCGAAAGGGCCGGCACTGAAGACACAAGAAAGGGGCCAGTTATTGGCCCCCGt                                                                                                                               >  6:227420/1‑149 (MQ=255)
          cGGCCCATTCGCGGGTGAACCCTCTCCTACACAGATCTGCGGTAGACACACACCGGACGTTGGAAACAGTCCGTTGTAGGAGCGGGTTTACCCGCGAAAGGGCCGGCACTGAAGACACAAGAAAGGGGCCAATAACTGGcccctttctt                                                                                                                            <  4:181656/149‑7 (MQ=37)
           ggCCCATTCGCGGGTGAACCCGCTCCTACACAGATCTGCGGTAGACACACACCGGACGTTGGAAACAGGCCGTTGTAGGAGCGGGTTTACCCGCGAAAGGGCCGGCACTGAAGACACAAGAAAGGGGCCAGTTATTGGCCCCCGTTCGt                                                                                                                           <  4:119871/149‑1 (MQ=255)
             cccATTCGCGGGTGAACCCGCTCCTACACAGATCTGCGGTAGACACACACCGGACGTTGGAAACAGGCCGTTGTAGGAGCGGGTTTACCCGCGAAAGGGCCGGCACTGAAGACACAAGAAAGGGGCCAGTTATTGGCCCCCGTTCGTca                                                                                                                         >  6:284336/1‑149 (MQ=255)
              ccATTCGCGGGTGAACCCGCTCCTACACAGATCTGCGGTAGACACACACCGGACGTTGGAAACAGGCCGTTGTAGGAGCGGGTTTACCCGCGAAAGGGCCGGCACTGAAGACACAAGAAAGGGGCCAGTTATTGGc                                                                                                                                     <  2:346172/136‑1 (MQ=255)
               cATTCGCGGGTGAACCCGCTCCTACACAGATCTGCGGTAGACACACACCGGACGTTGGAAACAGGCCGTTGTAGGAGCGGGTTTACCCGCGAAAGGGCCGGCACTGAAGACACAAGAAAGGGGCCAGTTATTGGCCCCCGTTCGTcaca                                                                                                                       >  5:403115/1‑149 (MQ=255)
                    gcgGGTGAACCCGCTCCTACACAGATCTGCGGTAGACACACACCGGACGTTGGAAACAGGCCGTTGTAGGAGCGGGTTTACCCGCGAAAGGGCCGGCACTGAAGACACAAGAAAGGGGCCAGTTATTGGCCCCCGTTCGTCACAtcag                                                                                                                   <  5:255221/148‑1 (MQ=255)
                             cccGCTCCTACACAGATCTGCGGTAGACACACACCGGACGTTGGAAACAGGCCGTTGTAGGAGCGGGTTTACCCGCGAAAGGGCCGGCACTGAAGACACAAGAAAGGGGCCAGTTATTGGCCCCCGTTCGTCACATcagctcagct                                                                                                            >  6:9769/1‑146 (MQ=255)
                              ccGCTCCTACACAGATCTGCGGTAGACACACACCGGACGTTGGAAACAGGCCGTTGTAGGAGCGGGTTTACCCGCGAAAGGGCCGGCACTGAAGACACAAGAAAGGGGCCAGTTATTGGCCCCCGTTCGTCACATCAGCTCAGCTTAg                                                                                                         >  1:137639/1‑148 (MQ=255)
                                    cTACACAGATCTGCGGTAGACACACACCGGACGTTGGAAACAGGCCGTTGTAGGAGGGGGTTTACCCGCGAAAGGGCCGGCACTGAAGACACAAGAAAGGGGCCAGTTATTGGCCCCCGTTCGTCACATCAGCTCAGCTTAGAAAGCg                                                                                                   <  5:41226/148‑1 (MQ=255)
                                    cTACACAGATCTGCGGTAGACACACACCGGACGTTGGAAACAGGCCGTTGTAGGAGCGGGTTTACCCGCGAAAGGGCCGGCACTGAAGACACAAGAAAGGGGCCAGTTATTGGCCCCCGTTCGTcaca                                                                                                                       <  6:283106/128‑1 (MQ=255)
                                    cTACACAGATCTGCGGTAGACACACACCGGACGTTGGAAACAGGCCGTTGTAGGAGCGGGTTTACCCGCGAAAGGGCCGGCACTGAAGACACAAGAAAGGGGCCAGTTATTGGCCCCCGTTCGTcaca                                                                                                                       >  5:283106/1‑128 (MQ=255)
                                      acacAGATCTGCGGTAGACACACACCGGACGTTGGAAACAGGCCGTTGTAGGAGCGGGTTTACCCGCGAAAGGGCCGGCACTGAAGACACAAGAAAGGGGCCAGTTATTGGCCCCCGTTCGTca                                                                                                                         >  3:78806/1‑124 (MQ=255)
                                      acacAGATCTGCGGTAGACACACACCGGACGTTGGAAACAGGCCGTTGTAGGAGCGGGTTTACCCGCGAAAGGGCCGGCACTGAAGACACAAGAAAGGGGCCAGTTATTGGCCCCCGTTCGTca                                                                                                                         <  4:78806/124‑1 (MQ=255)
                                       cacaGATCTGCGGTAGACACACACCGGAGGGTGGAAACAGGCCGTTGTAGGAGCGGGTTTACCCGCGAAAGGGCCGGCACTGAAGACACAAGAAAGGGGCCAGTTATTGGCCCCCGTTCGTCACATCAGCTCAGCTTAGAAAGCGaa                                                                                                 <  1:359820/147‑1 (MQ=255)
                                       cacaGATCTGCGGTAGACACACACCGGACGTTGGAAACAGGCCGTTGTAGGAGCGGGTTTACCCGCGAAAGGGCCGGCACTGAAGACACAAGAAAGGGGCCAGTTATTGGCCCCCGTTCGTCACATCAGCTCAGCTTAGAAAGCGAAg                                                                                                >  5:170900/1‑148 (MQ=255)
                                        acaGATCTGCGGTAGACACACACCGGACGTTGGAAACAGGCCGTTGTAGGAGCGGGTTTACCCGCGAAAGGGCCGGCACTGAAGACACAAGAAAGGGGCCAGTTATTGGCCCCCGTTCGTCACATCAGCTCAGCTTAGAAAGCGAAGCa                                                                                              <  5:284336/149‑1 (MQ=255)
                                        acaGATCTGCGGTAGACACACACCGGACGTTGGAAACAGGCCGTTGTAGGAGCGGGTTTACCCGCGAAAGGGCCGGCACTGAAGACACAAGAAAGGGGCCAGTTATTGGCCCCCGTTCGTCACATCAGCTCAGCTTAGAAAGCGAAGCa                                                                                              >  7:207876/1‑149 (MQ=255)
                                         caGATCTGCGGTAGACACACACCGGACGTTGGAAACAGGCCGTTGTAGGAGCGGGTTTACCCGCGAAAGGGCCGGCACTGAAGACACAAGAAAGGGGCCAGTTATTGGCCCCCGTTCGTCACATCAGCTCAGCTTAGAAAGCGAAGCAc                                                                                             <  5:227420/149‑1 (MQ=255)
                                                 cGGTAGACACACACCGGACGTTGGAAACAGGCCGTTGTAGGAGCGGGTTTACCCGCGAAAGGGCCGGCACTGAAGACACAAGAAAGGGGCCAGTTATTGGCCCCCGTTCGTCACATCAGCTCAGCTTAGAAAGCGAAGCACGAACAGc                                                                                      <  3:324143/148‑1 (MQ=255)
                                                  ggTAGACACACACCGGACGTTGGAAACAGGCCGTTGTAGGAGCGGGTTTACCCGCGAAAGGGCCGGCACTGAAGACACAAGAAAGGGGCCAGTTATTGGCCCCCGTTCGTCACATCAGCTCAGCTTAGAAAGCGAAGCACGAACAGccc                                                                                    <  2:263832/149‑1 (MQ=255)
                                                   gTAGACACACACCGGACGTTGGAAACAGGCCGTTGTAGGAGCGGGTTTACCCGCGAAAGGGCCGGCACTGAAGACACAAGAAAGGGGCCAGTTATTGGCCCCCGTTCGTCACATCAGCTCAGCTTAGAAAGCGAAGCACGAACAGcc                                                                                     >  8:20831/1‑147 (MQ=255)
                                                   gTAGACACACACCGGACGTTGGAAACAGGCCGTTGTAGGAGCGGGTTTACCCGCGAAAGGGCCGGCACTGAAGACACAAGAAAGGGGCCAGTTATTGGCCCCCGTTCGTCACATCAGCTCAGCTTAGAAAGCGAAGCACGAACAGCCCa                                                                                   <  7:121769/149‑1 (MQ=255)
                                                   gTAGACACACACCGGACGGTGTACACAGGCCGTTGTAGGAGCGGGTTTACCCGCGAAAGGGCCGGCACTGAAGACACAAGAAAGGGGCCAGTTATTGGCCCCCGTTCGTCACATCAGCTCAGCTTAGAAAGCGAAGCACGAACAGCCCa                                                                                   <  6:239802/149‑1 (MQ=255)
                                                       acacacacCGGACGTTGGAAACAGGCCGTTGTAGGAGCGGGTTTACCCGCGAAAGGGCCGGCACTGAAGACACAAGAAAGGGGCCAGTTATTGGCCCCCGTTCGCCACATCAGCTCAGCTTAGAAAGCGAAGCAc                                                                                             <  7:237351/135‑1 (MQ=255)
                                                       acacacacCGGACGTTGGAAACAGGCCGTTGTAGGAGCGGGTTTACCCGCGAAAGGGCCGGCACTGAAGACACAAGAAAGGGGCCAGTTATTGGCCCCCGTTCGCCACATCAGCTCAGCTTAGAAAGCGAAGCAc                                                                                             >  8:237351/1‑135 (MQ=255)
                                                        cacacacCGGACGTTGGAAACAGGCCGTTGTAGGAGCGGGTTTACCCGCGAAAGGGCCGGCACTGAAGACACAAGAAAGGGGCCAGTTATTGGTCCCCGTTCGTCACATCAGCTCAGCTTAGAAAGCGAAGCACGAACAGCCCAGCGCa                                                                              >  8:48538/1‑149 (MQ=255)
                                                         acacacCGGACGTTGGAAACAGGCCGTTGTAGGAGCGGGTTTACCCGCGAAAGGGCCGGCACTGAAGACACAAGAAAGGGGCCAGTTATTGGCCCCCGTTCGTCACATCAGCTCCGCTTAGAAAGCGAAGc                                                                                               >  4:301634/1‑131 (MQ=255)
                                                         acacacCGGACGTTGGAAACAGGCCGTTGTAGGAGCGGGTTTACCCGCGAAAGGGCCGGCACTGAAGACACAAGAAAGGGGCCAGTTATTGGCCCCCGTTCGTCACATCAGCTCCGCTTAGAAAGCGAAGc                                                                                               <  3:301634/131‑1 (MQ=255)
                                                           acacCGGACGTTGGAAACAGGCCGTTGTAGGAGCGGGTTTACCCGCGAAAGGGCCGGCACTGAAGACACAAGAAAGGGGCCAGTTATTGGCCCCCGTTCGTCACATCAGCTCAGCTTAGAAAGCGAAGCACGAACAGCCCAGCGCAccc                                                                           >  8:233424/1‑149 (MQ=255)
                                                            cacCGGACGTTGGAAACAGGACGTTGTAGGAGCGGTTTTACCCGCGAAAGGGCCGGCACTGAAGACACAAGAAAGGGGCCAGTTATTGGCCCCCGTTCGTCACATCAGCTCAGCTTAGAAAGCGAAGCACGAACAGCCCAGCGCAcccc                                                                          <  3:364400/149‑1 (MQ=255)
                                                                   cGTTGGAAACAGGCCGTTGTAGGAGCGGGTTTACCCGCGAAAGGGCCGGCACTGAAGACACAAGAAAGGGGCCAGTTATTGGCCCCCGTTCGTCACATCAGCTCAGCTTAGAAAGCGAAGCACGAACAGCCCAGCGCACCCCAGAAGcc                                                                   <  2:137639/149‑1 (MQ=255)
                                                                    gTTGGAAACAGGCCGTTGTAGGAGCGGGTTTACCCGCGAAAGGGCCGGCACTGAAGACACAAGAAAGGGGCCAGTTATTGGCCCCCGTTCGTCACATCAGCTCAGCTTAGAAAGCGAAGCACGAACAGCCCAGCGCACCCCAGAAGccc                                                                  <  6:51943/149‑1 (MQ=255)
                                                                      tGGAAACAGGCCGTTGTAGGAGCGGGTTTACCCGCGAAAGGGCCGGTACTGAAGACACAAGAAAGGGGCCAGTTATTGGCCCCCGTTCGTCACATCAGCTCAGCTTAGAAAGCGAAGCACGAACAGCCCAGCGCAccc                                                                           >  6:36054/1‑138 (MQ=255)
                                                                      tGGAAACAGGCCGTTGTAGGAGCGGGTTTACCCGCGAAAGGGCCGGTACTGAAGACACAAGAAAGGGGCCAGTTATTGGCCCCCGTTCGTCACATCAGCTCAGCTTAGAAAGCGAAGCACGAACAGCCCAGCGCAccc                                                                           <  5:36054/138‑1 (MQ=255)
                                                                        gAAACAGGCCGTTGTAGGAGCGGGTTTACCCGCGAAAGGGCCGGCACTGAAGACACAAGAAAGGGGCCAGTTATTGGCCCCCGTACGTCACATCAGCTCAGCTTAGAAAGCGAAGCACGAACAGCCCAGCGCACCCCAGAAGCCCTgg                                                               <  5:69538/148‑1 (MQ=255)
                                                                         aaaCAGGCCGTTGTAGGAGCGGGTTTACCCGCGAAAGGGCCGGCACTGAAGACACAAGAAAGGGGCCAGTTATTGGCCCCCGTTCGTCACATCAGCTCAGCTTAGAAAGCGAAGCACGAACAGCCCAGCGCACCCCAGAAGCCCTGGaa                                                             <  7:233424/149‑1 (MQ=255)
                                                                         aaaCAGGCAGTTGTAGGAGCGGGTTTACCCGCGAAAGGGCCGGCACTGAAGACACAAGAAAGGGGCCAGTTATTGGCCCCCGTTCGTCACATCAGCTCAGCTTAGAAAGCGAAGCACGAACAGCCCAGCGCACCCCAGAAGCCCTGGaa                                                             <  1:411956/149‑1 (MQ=255)
                                                                             aGGCCGTTGTAGGAGCGGGTTTACCCGCGAAAGGGCCGGCACTGAAGACACAAGAAAGGGGCCAGTTATTGGCCCCCGTTCGTCACATCAGCTCAGCTTAGAAAGCGAAGCACGAACAGCCCAGCGCATCCCAGAAGCCCTGGAAGTCa                                                         <  6:403115/149‑1 (MQ=255)
                                                                               gCCGTTGTAGGAGCGGGTTTACCCGCGAAAGGGCCGGCACTGAAGACACAAGAAAGGGGCCAGTTATTGGCCCCCGTTCGTCACATCAGCTCAGCTTAGAAAGCGAAGCACGAACAGCCCAGCGCAccc                                                                           <  2:200422/129‑1 (MQ=255)
                                                                               gCCGTTGTAGGAGCGGGTTTACCCGCGAAAGGGCCGGCACTGAAGACACAAGAAAGGGGCCAGTTATTGGCCCCCGTTCGTCACATCAGCTCAGCTTAGAAAGCGAAGCACGAACAGCCCAGCGCAccc                                                                           >  1:200422/1‑129 (MQ=255)
                                                                               gCCGTTGTAGGAGCGGGTTTACCCGCGAAAGGGCCGGCACTGAAGACACAAGAAAGGGGCCAGTTATTGGCCCCCGTTCGTCACATCAGCTCAGCTTAGAAAGCGAAGCACGAACAGCCCAGCGCACCCCAGAAGCCCTGGAAGTCAtt                                                       <  7:341271/149‑1 (MQ=255)
                                                                               gCCGTTGTAGGAGCGGGTTTACCCGCGAAAGGGCCGGCACTGAAGACACAAGAAAGGGGCCAGTTATTGGCCCCCGTTCGTCACATCAGCTCAGCTTAGAAAGCGAAGCACGAACAGCCCAGCGCACCCCAGAAGCCCTGGAAGTCAtt                                                       >  3:263622/1‑149 (MQ=255)
                                                                                     gTAGGAGCGGGTTTACCCGCGAAAGGGCCGGCACTGAAGACACAAGAAAGGGGCCAGTTATTGGCCCCCGTTCGTCACATCAGCTCAGCTTAGAAAGCGAAGCACGAACAGCCCAGCGCACCCCAGAAGCCCTGGAAGTCATTTACCgg                                                 >  3:221935/1‑149 (MQ=255)
                                                                                      tAGGAGCGGGTTTACCCGCGAAAGGGCCGGCACTGAAGACACAAGAAAGGGGCCAGTTATTGGCCCCCGTTCGTCACATCAGCTCAGCTTAGAAAGCGAAGCACGAACAGCCCAGCGCACCCCAGAAGCCCTGGAAGTCATTTACCgg                                                 >  4:144118/1‑148 (MQ=255)
                                                                                       aGGAGCGGGTTTACCCGCGAAAGGGCCGGCACTGAAGACACAAGAAAGGGGCCAGTTATTGGCCCCCGTTCGTCACATCAGCTCAGCTTAGAAAGCGAAGCACGAACAGCCCAGCGCACCCCAGAAGCCCTGGAAGTCATTTACCGGCa                                               >  8:254120/1‑149 (MQ=255)
                                                                                         gAGCGGGTTTACCCGCGAAAGGGCCGGCACTGAAGACACAAGAAAGGGGCCAGTTATTGGCCCCCGTTCGTCACATCAGCTCAGCTTAGAAAGCGAAGCACGAACAGCCCAGCGCACCCCAGAAGCCCTGGAAGTCATTTACCGGCACg                                             >  3:306488/1‑149 (MQ=255)
                                                                                           gCGGGTTTACCCGCGAAAGGGCCGGCACTGAAGACACAAGAAAGGGGCCAGTTATTGGCCCCCGTTCGTCACATCAGCTCAGCTTa                                                                                                          <  6:169239/86‑1 (MQ=255)
                                                                                           gCGGGTTTACCCGCGAAAGGGCCGGCACTGAAGACACAAGAAAGGGGCCAGTTATTGGCCCCCGTTCGTCACATCAGCTCAGCTTa                                                                                                          >  5:169239/1‑86 (MQ=255)
                                                                                               gTTTACCCGCGAAAGGGCCGGCACTGAAGACACAAGAAAGGGGCCAGTTATTGGCCCCCGTTCGTCACATCAGCTCAGCTTAGAAAGCGAAGCACGAACAGCCCAGCGCACCCCAGAAGCCCTGGAAGTCATTTACCGGCACGCTCGaa                                       <  4:350759/149‑1 (MQ=255)
                                                                                               gTTTACCCGCGAAAGGGCCGGCACTGAAGACACAAGAAAGGGGCCAGTTATTGGCCCCCGTTCGTCACATCAGCTCAGCTTAGAAAGCGAAGCACGAACAGCCCAGCGCACCCCAGAAGCCCTGGAAGTCATTTACCGGCACGCTCGaa                                       >  6:76294/1‑149 (MQ=255)
                                                                                                tttACCCGCGAAAGGGCCGGCACTGAAGACACAAGAAAGGGGCCAGTTATTGGTCCCCGTTCGTCACATCAGCTCAGCTTAGAAAGCGAAGCACGAACAGCCCAGCGCACCCCAGAAGCCCTGGAAGTCATTTACCGGCACGCTCGAAt                                      <  7:48538/149‑1 (MQ=255)
                                                                                                tttACCCGCGAAAGGGCCGGCACTGAAGACACAAGAAAGGGGCCAGTTATTGGCCCCCGTTCGTCACATCAGCTCAGCTTAGAAAGCGAAGCACGAACAGCCCAGCGCACCCCAGAAGCCCTGGAAGTCATTTACCGGTACGCTCGAAt                                      >  6:102833/1‑149 (MQ=255)
                                                                                                   aCCCGCGAAAGGGCCGGCACTGAAGACACAAGAAAGGGGCCAGTTATTGGCCCCCGTTCGTCACATCAGCTCAGCTTAGAAAGCGAAGCACGAACAGCCCAGCGCACCCCAGAAGCCCTGGAAGTCATTTACCGGCACGCTCGAATGcc                                   >  2:147646/1‑149 (MQ=255)
                                                                                                    cccGCGAAAGGGCCGGCACTGAAGACACAAGAAAGGGGCCAGTTATTGGCCCCCGTTCGTCACATCAGCTCAGCTTAGAAAGCGAAGCACGAACAGCCCAGCGCACCCCAGAAGCCCTGGAAGTCATTTACCGGCACGc                                            >  2:155950/1‑139 (MQ=255)
                                                                                                    cccGCGAAAGGGCCGGCACTGAAGACACAAGAAAGGGGCCAGTTATTGGCCCCCGTTCGTCACATCAGCTCAGCTTAGAAAGCGAAGCACGAACAGCCCAGCGCACCCCAGAAGCCCTGGAAGTCATTTACCGGCACGc                                            <  1:155950/139‑1 (MQ=255)
                                                                                                     ccGCGAAAGGGCCGGCACTGAAGACACAAGAAAGGGGCCAGTTATTGGCCCCCGTTCGTCACATCAGCTCAGCTTAGAAAGCGAAGCACGAACAGCCCAGCGCACCCCAGAAGCCCTGGAACTCATTTACCGGCACGCTCGAATGCcgc                                 >  6:54695/1‑149 (MQ=255)
                                                                                                       gcgAAAGGGCCGGCACTGAAGACACAAGAAAGGGGCCAGTTATTGGCCCCCGTTCGTCACATCAGCTCAGCTTAGAAAGCGAAGCACGAACAGCCCAGCGCACCCCAGAAGCCCTGGaa                                                             <  2:96804/119‑1 (MQ=255)
                                                                                                       gcgAAAGGGCCGGCACTGAAGACACAAGAAAGGGGCCAGTTATTGGCCCCCGTTCGTCACATCAGCTCAGCTTAGAAAGCGAAGCACGAACAGCCCAGCGCACCCCAGAAGCCCTGGaa                                                             >  1:96804/1‑119 (MQ=255)
                                                                                                          aaaGGGCCGGCACTGAAGACACAAGAAAGGGGCCAGTTATTGGCCCCCGTTCGTCACATCAGCTCAGCTTAGAAAGCGAAGCACGAACAGCCCAGCGCACCCCAGAAGCCCTGGa                                                              <  4:290033/115‑1 (MQ=255)
                                                                                                          aaaGGGCCGGCACTGAAGACACAAGAAAGGGGCCAGTTATTGGCCCCCGTTCGTCACATCAGCTCAGCTTAGAAAGCGAAGCACGAACAGCCCAGCGCACCCCAGAAGCCCTGGa                                                              >  3:290033/1‑115 (MQ=255)
                                                                                                                 cGGCACTGAAGACACAAGAAAGGGGCCAGTTATTGGCCCCCGTTCGTCACATCAGCTCAGCTTAGAAAGCGAAGCACGAACAGCCCAGCGCACCCCAg                                                                        >  2:179485/1‑98 (MQ=255)
                                                                                                                 cGGCACTGAAGACACAAGAAAGGGGCCAGTTATTGGCCCCCGTTCGTCACATCAGCTCAGCTTAGAAAGCGAAGCACGAACAGCCCAGCGCACCCCAg                                                                        <  1:179485/98‑1 (MQ=255)
                                                                                                                     aCTGAAGACACAAGAAAGGGGCCAGTTATTGGCCCCCGTTCGTCACATCAGCTCAGCTTAGAAAGCGAAGCACGAACAGCCCAGCGCACCCCAGAAGCCCTGGAAGTCATTTACCGGCACGCTCGAATGCCGCGCCTTTTCGTGGCtgt                 <  7:20831/149‑1 (MQ=255)
                                                                                                                      cTGAAGACACAAGAAAGGGGCCAGTTATTGGCCCCCGTTCGTCACATCAGCTCAGCTTAGAAAGCGAAGCACGAACAGCCCAGCGCACCCCAGAAGCCCTGGAAGTCATTTACCGGCACGCTCGAATGCCGCGCCTTTTCGTGGCtgtg                <  3:393818/149‑1 (MQ=255)
                                                                                                                         aaGACACAAGAAAGGGGCCAGTTATTGGCCCCCGTTCGTCACATCAGCTCAGCTTAGAAAGCGAAGCACGAACAGCCCAGCGCACCCCAGAAGCCCTGGAAGTCATTTACCGGCACGCTCGAATGCCGCGCCTTTTCGTGGCTGTgcg              >  8:14739/1‑148 (MQ=255)
                                                                                                                           gACACAAGAAAGGGGCCAGTTATTGGCCCCCGTTCGTCACATCAGCTCAGCTTAGAAAGCGAAGCACGAACAGCCCAGCGCACCCCAGAAGCCCTGGAAGTCATTTACCGGCACGCTCGAATGCCGCGCCTTTTCGTGGCTGTGCGCAt           <  8:207876/149‑1 (MQ=255)
                                                                                                                             cacaAGAAAGGGGCCAGTTATTGGCCCCCGTTCGTCACATCAGCTCAGCTTAGAAAGCGAAGCACGAACAGCCCAGCGCACCCCAGAAGCCCTGGAAGTCATTTACCGGCACGCTCGAATGCCGCGCCTTTTCGTGGCTGTGCGCATg          <  6:170900/148‑1 (MQ=255)
                                                                                                                              acaAGAAAGGGGCCAGTTATTGGCCCCCGTTCGTCACATCAGCTCAGCTTAGAAAGCGAAGCACGAACAGCCCAGCGCACCCCAGAAGCCCTGGAAGTCATTTACCGGCACGCTCGAATGCCGCGCCTTTTCGTgg                     >  4:56705/1‑136 (MQ=255)
                                                                                                                              acaAGAAAGGGGCCAGTTATTGGCCCCCGTTCGTCACATCAGCTCAGCTTAGAAAGCGAAGCACGAACAGCCCAGCGCACCCCAGAAGCCCTGGAAGTCATTTACCGGCACGCTCGAATGCCGCGCCTTTTCGTgg                     <  3:56705/136‑1 (MQ=255)
                                                                                                                               acgaacgggggCCAATAACTGGCCCCTGTTCGTCACATCAGCTCAGCTTAGAAAGCGAAGCACGAACAGCCCAGCGCACCCCAGAAGCCCTGGAAGTCATTTACCGGCACGCTCGAATGCCGCGCCTTTTCGTGGCTGTGCGCATGca        >  7:193241/8‑148 (MQ=255)
                                                                                                                                aaGAAAGGGGCCAGTTATTGGCCCCCGTTCGTCACATCAGCTCAGCTTAGAAAGCGAAGCACGAACAGCCCAgcg                                                                                <  5:375012/75‑1 (MQ=255)
                                                                                                                                aaGAAAGGGGCCAGTTATTGGCCCCCGTTCGTCACATCAGCTCAGCTTAGAAAGCGAAGCACGAACAGCCCAgcg                                                                                >  6:375012/1‑75 (MQ=255)
                                                                                                                                aaGAAAGGGGCCAGTTATTGGCCCCCGTTCGTCACATCAGCTCAGCTTAGAAAGCGAAGCACGAACAGCCCAGCGCACCCCAGAAGCCCTGGAAGTCATTTACCGGCACGCTCGAATGCCGCGCCTTTTCGTGGCTGTGCGCATGcaca      >  4:396415/1‑149 (MQ=255)
                                                                                                                                 aGAAAGGGGCCAGTTATTGGCCCCCGTTCGTCACATCAGCTCAGCTTAGAAAGCGAAGCACGAACAGCCCAGCGCACCCCAGAAGCCCTGGAAGTCATTTACCGGCACGCTCGAATGCCGCGCCTTTTCGTGGCTGTGCGCATGCacac     <  4:306488/149‑1 (MQ=255)
                                                                                                                                  gAAAGGGGCCAGTTATTGGCCCCCGTTCGTCACATCAGCTCAGCTTAGAAAGCGAAGCAc                                                                                             >  5:253159/1‑60 (MQ=255)
                                                                                                                                  gAAAGGGGCCAGTTATTGGCCCCCGTTCGTCACATCAGCTCAGCTTAGAAAGCGAAGCAc                                                                                             <  6:253159/60‑1 (MQ=255)
                                                                                                                                    cgggggCCAATAACTGGCCCCTGTTCGTCACATCAGCTCAGCTTAGAAAGCGAAGCACGAACAGCCCAGCGCACCCCAGAAGCCCTGGAAGTCATTTACCGGCACGCTCGAATGCCGCGCCTTTTCGTGGCTGTGCGCATGCAcaccac  <  8:193241/147‑1 (MQ=255)
                                                                                                                                    cgggggCCAATAACTGGCCCCCGTTCGTCACATCAGCTCAGCTTAGAAAGCGAAGCACGAACAGCCCAGCGCACCCCAGAAGCCCTGGAAGACATTTACCGGCACGCTCGAATGCCCCGCCTTTTCTTGGATGTGCGAATGCAcacca   >  2:308864/3‑148 (MQ=255)
                                                                                                                                               |                                                                                                                                         
CCGGCCCTTACGGCCCATTCGCGGGTGAACCCGCTCCTACACAGATCTGCGGTAGACACACACCGGACGTTGGAAACAGGCCGTTGTAGGAGCGGGTTTACCCGCGAAAGGGCCGGCACTGAAGACACAAGAAAGGGGCCAGTTATTGGCCCCCGTTCGTCACATCAGCTCAGCTTAGAAAGCGAAGCACGAACAGCCCAGCGCACCCCAGAAGCCCTGGAAGTCATTTACCGGCACGCTCGAATGCCGCGCCTTTTCGTGGCTGTGCGCATGCACACCAC  >  NC_002947/5994875‑5995155

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: