Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,905,973 C→G 5.6% P276P (CCG→CCC spuF ← spermidine/putrescine ABC transporter ATP binding protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,905,9730CG 5.6% 227.1 / 3.7 71P276P (CCG→CCCspuFspermidine/putrescine ABC transporter ATP binding protein
Reads supporting (aligned to +/- strand):  ref base C (39/28);  new base G (2/2);  total (41/30)
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.50e-01

TGGATCTTGCCGCGCGACCAGTTGTGCTCGCAGGTTGGCTGCTGGGTGGTGACCAGCATCTTTTCCGGGCGCAGGGCGTAGGTGATGTGCTTGTCTTCCACCGAGGTGGTGATGCCGTGGCCCACGTAGATCTTGCGCTCCAGCTCCGGGCTGGCAATGATCGCGTAGCCTTCGGCGTCGTCGACCACGTCACCTTCGAACAGGTTGACGTTACCGATGAACTCGCACACCAGGCGGCTGGTAGGCGTCTCGTAGATGTCCACTGGCGAGCCGATCTGGGCGATCCAGCCCAGGT  >  NC_002947/5905827‑5906121
                                                                                                                                                  |                                                                                                                                                    
tGGATCTTGCCGCGCGACCAGTTGTGCTCGCAGGTTGGCTGCTGGGTGGTGACCAGCATCTTTTCCGGGCGCAGGGCGTAGGTGATGTGCTTGTCTTCCACCGAGGTGGTGATGCCGTGGCCCACGTAGATCTTGCGCTCCAGCTCCgg                                                                                                                                                    <  1:49316/149‑1 (MQ=255)
 ggATCTTGCCGCGCGACCAGTTGTGCTCGCAGGTTGGCTGCTGGGTGGTGACCAGCATCTTTTCCGGGCGCAGGGCGTAGGTGATGTGCTTGTCTTCCACCGAGGTGGTGATGCCGTGGCCCACGTAGATCTTGCGCTCCAGCTCCggg                                                                                                                                                   >  1:169224/1‑149 (MQ=255)
    tCTTGCCGCGCGACCAGTTGTGCTCGCAGGTTGGCTGCTGGGTGGTGACCAGCATCTTTTCCGGGCGCAGGGCGTAGGTGATGTGCTTGTCTTCCACCGAGGTGGTGATGCCGTGGCCCACGTAGATCTTGCGCTCCAGCTCCGGGCTg                                                                                                                                                >  8:298287/1‑149 (MQ=255)
        gCCGCGCGACCAGTTGTGCTCGCAGGTTGGCTGCTGGGTGGTGACCAGCATCTTTTCCGGGCGCAGGGCGTAGGTGATGTGCTTGTCTTCCACCGAGGTGGTGATGCCGTGGCCCACGTAGATCTTGCGCTCCAGCTCCGGGCTgg                                                                                                                                               <  7:298287/146‑1 (MQ=255)
        gCCGCGCGACCAGTTGTGCTCGCAGGTTGGCTGCTGGGTGGTGACCAGCATCTTTTCCGGGCGCAGGGCGTAGGTGATGTGCTTGTCTTCCACCGAGGTGGTGATGCCGTGGCCCACGTAGATCTTGCGCTCCAGCTCCGGGCTgg                                                                                                                                               <  3:234417/146‑1 (MQ=255)
        gCCGCGCGACCAGTTGTGCTCGCAGGTTGGCTGCTGGGTGGTGACCAGCATCTTTTCCGGGCGCAGGGCGTAGGTGATGTGCTTGTCTTCCACCGAGGTGGTGATGCCGTGGCCCACGTAGATCTTGCGCTCCAGCTCCGGGCTgg                                                                                                                                               <  2:169224/146‑1 (MQ=255)
               gACCAGTTGTGCTCGCAGGTTGGCTGCTGGGTGGTGACCAGCATCTTTTCCGGGCGCAGGGCGTAGGTGATGTGCTTGTCTTCCACCGAGGTGGTGATGCCGTGGCCCACGTAGATCTTGCGCTCCAGCTCCGGGCTGGCAATGATcgc                                                                                                                                     <  1:32554/149‑1 (MQ=255)
                  cAGTTGTGCTCGCAGGTTGGCTGCTGGGTGGTGACCAGCATCTTTTCCGGGCGCAGGGCGTAGGTGATGTGCTTGTCTTCCACCGAGGTGGTGATGCCGTGGCCCACGTAGATCTTGCGCTCCAGCTCCGGGCTGGCAATGATCGCGTa                                                                                                                                  <  7:219185/149‑1 (MQ=255)
                     ttGTGCTCGCAGGTTGGCTGCTGGGTGGTGACCAGCATCTTTTCCGGGCGCAGGGCGTAGGTGATGTGCTTGTCTTCCACCGAGGTGGTGATGCCGTGGCCCACGTAGATCTTGCGCTCCAGCTCCGGGCTGGCAATGATCGCGTAGcc                                                                                                                               >  2:182212/1‑149 (MQ=255)
                          cTCGCAGGTTGGCTGCTGGGTGGTGACCAGCATCTTTTCCGGGCGCAGGGCGTAGGTGATGTGCTTGTCTTCCACCGAGGTGGTGATGCCGTGGCCCACGTAGATCTTGCGCTCCAGCTCCGGGCTGGCAATGATCGCGTAGCCTTCgg                                                                                                                          >  7:180263/1‑149 (MQ=255)
                             gCAGGTTGGCTGCTGGGTGGTGACCAGCATCTTTTCCGGGCGCAGGGCGTAGGTGATGTGCTTGTCTTCCACCGAGGTGGTGATGCCGTGGCCCACGTAGATCTTGCGCTCCAGCTCCGGGCTGGCAATGATCGCGTAGCCTTCGGcgt                                                                                                                       >  7:103557/1‑149 (MQ=255)
                                    ggCTGCTGGGTGGTGACCAGCATCTTTTCCGGGCGCAGGGCGTAGGTGATGTGCTTGTCTTCCACCGAGGTGGTGATGCCGTGGCCCACGTAGATCTTGCGCTCCAGCTCCGGGCTGGCAATGATCGCGTAGCCTTCGGCGTCGTCGa                                                                                                                 <  1:158216/148‑1 (MQ=255)
                                     gctgctGGGTGGTGACCAGCATCTTTTCCGGGCGCAGGGCGTAGGTGATGTGCTTGTCTTCCACCGAGGTGGTGATGCCGTGGCCCACGTAGATCTTGCGCTCCAGCTCCGGGCTGGCAATGATCGCGTAGCCTTCgg                                                                                                                          <  7:119892/138‑1 (MQ=255)
                                     gctgctGGGTGGTGACCAGCATCTTTTCCGGGCGCAGGGCGTAGGTGATGTGCTTGTCTTCCACCGAGGTGGTGATGCCGTGGCCCACGTAGATCTTGCGCTCCAGCTCCGGGCTGGCAATGATCGCGTAGCCTTCgg                                                                                                                          >  8:119892/1‑138 (MQ=255)
                                     gctgctGGGTGGTGACCAGCATCTTTTCCGGGCGCAGGGCGTAGGTGATGTGCTTGTCTTCCACCGAGGTGGTGATGCCGTGGCCCACGTAGATCTTGCGCTCCAGCTCCGGGCTGGCAATGATCGCGTAGCCTTCGGCGtcgtcg                                                                                                                  <  7:272223/146‑1 (MQ=255)
                                     gctgctGGGTGGTGACCAGCATCTTTTCCGGGCGCAGGGCGTAGGTGATGTGCTTGTCTTCCACCGAGGTGGTGATGCCGTGGCCCAAGTAGATCTTGCGCTCCAGCTCCGGGCTGGCAATGATCGCGTAGCCTTCGGCGtcgtcg                                                                                                                  >  8:272223/1‑146 (MQ=255)
                                       tgctgGGTGGTGACCAGCATCTTTTCCGGGCGCAGGGCGTAGGTGATGTGCTTGTCTTCCACCGAGGTGGTGATGCCGTGGCCCACGTAGATCTTGCGCTCCAGCTCCGGGCTGGCAATGATCGCGTAGCCTTCGGCGTCGTCGACCAc                                                                                                             >  1:325674/1‑149 (MQ=255)
                                           gggTGGTGACCAGCATCTTTTCCGGGCGCAGGGCGTAGGTGATGTGCTTGTCTTCCACCGAGGTGGTGATGCCGTGGCCCACGTAGATCTTGCGCTCCAGCTCCGGGCTGGCAATGATCGCGTAGCCTTCGGCGTCGTCGACCACGTCa                                                                                                         >  6:266973/1‑149 (MQ=255)
                                            ggtggtGACCAGCATCTTTTCCGGGCGCAGGGCGTAGGTGATGTGCTTGTCTTCCACCGAGGTGGTGATGCCGTGGCCCACGTAGATCTTGCGCTCCAGCTCCGGGCTGGCAATGATCGCGTAGCCTTCGGCGTCGTCGACCACGTCAc                                                                                                        >  6:28064/1‑149 (MQ=255)
                                                   aCCAGCATCTTTTCCGGGCGCAGGGCGTAGGTGATGTGCTTGTCTTCCACCGAGGTGGTGATGCCGTGGCCCACGTAGATCTTGCGCTCCAGCTCCGGGCTGGCAATGATCGCGTAGCCTTCGGCGTCGTCGACCACGTCACCTTCGaa                                                                                                 <  5:266973/149‑1 (MQ=255)
                                                       gCATCTTTTCCGGGCGCAGGGCGTAGGTGATGTGCTTGTCTTCCACCGAGGTGGTGATGCCGTGGCCCACGTAGATCTTGCGCTCCAGCTCCGGGCTGGCAATGATCGCGTAGCCTTCGGCGTCGTCGACCACGTCACCTTCGAACAgg                                                                                             >  8:59045/1‑149 (MQ=255)
                                                       gCATCTTTTCCGGGCGCAGGGCGTAGGTGATGTGCTTGTCTTCCACCGAGGTGGTGATGCCGTGGCCCACGTAGATCTTGCGCTCCAGCTCCGGGCTGGCAATGATCGCGTAGCCTTCGGCGTCGTCGACCACGTCACCTTCGAACAgg                                                                                             <  5:174748/149‑1 (MQ=255)
                                                        cATCTTTTCCGGGCGCAGGGCGTAGGTGATGTGCTTGTCTTCCACCGAGGTGGTGATGCCGTGGCCCACGTAGATCTTGCGCTCCAGCTCCGGGCTGGCAATGATCGCGTAGCCTTCGGCGTCGTCGACCACGTCACCTTCGAACAGGt                                                                                            >  1:173055/1‑149 (MQ=255)
                                                        cATCTTTTCCGGGCGCAGGGCGTAGGTGATGTGCTTGTCTTCCACCGAGGTGGTGATGCCGTGGCCCACGTAGATCTTGCGCTCCAGCTCCGGGCTGGCAATGATCGCGTAGCCTTCGGCGTCGTCGACCACGTCACCTTCGAACAGGt                                                                                            >  1:63413/1‑149 (MQ=255)
                                                        cATCTTTTCCGGGCGCAGGGCGTAGGTGATGTGCTTGTCTTCCACCGAGGTGGTGATGCCGTGGCCCACGTAGATCTTGCGCTCCAGCCCGGAGCTggagc                                                                                                                                            >  6:126843/1‑98 (MQ=255)
                                                        cATCTTTTCCGGGCGCAGGGCGTAGGTGATGTGCTTGTCTTCCACCGAGGTGGTGATGCCGTGGCCCACGTAGATCTTGCGCTCCAGCCCGGAGCTggagc                                                                                                                                            <  5:126843/101‑4 (MQ=255)
                                                          tCTTTTCCGGGCGCAGGGCGTAGGTGATGTGCTTGTCTTCCACCGAGGTGGTGATGCCGTGGCCCACGTAGATCTTGCGCTCCAGCTCCGGGCTGGCAATGATCGCGTAGCCTTCGGCGTCGTCGACCACGTCACCTTCGAACAGGtt                                                                                           >  1:295428/1‑148 (MQ=255)
                                                            ttttCCGGGCGCAGGGCGTAGGTGATGTGCTTGTCTTCCACCGAGGTGGTGATGCCGTGGCCCACGTAGATCTTGCGCTCCAGCTCCGGGCTGGCAATGATCGCGTAGCCTTCGGCGTCGTCGACCACGTCACCTTCGAACAGGTTGAc                                                                                        <  5:28064/149‑1 (MQ=255)
                                                              ttCCGGGCGCAGGGCGTAGGTGATGTGCTTGTCTTCCACCGAGGTGGTGATGCCGTGGCCCACGTAGATCTTGCGCTCCAGCTCCGGGCTGGCAATGATCGCGTAGCCTTCGGCGTCGTCGACCACGTCACCTTCGAACAGGTTGAc                                                                                        <  1:182212/147‑1 (MQ=255)
                                                                ccGGGCGCAGGGCGTAGGTGATGTGCTTGTCTTCCACAGAGGTGGTGATGCCGTGCCCCACGTAGATCTTGCGCTCCAGCTCCGGGATGGCAATGATAGCTTAGACTTCGGCGTCGTCGACCACGACACCTTc                                                                                                    >  7:86427/1‑133 (MQ=255)
                                                                      gcAGGGCGTAGGTGATGTGCTTGTCTTCCACCGAGGTGGTGATGCCGTGGCCCACGTAGATCTTGCGCTCCAGCTCCGGGCTGGCAATGATCGCGTAGCCTTCGGCGTCGTCGACCACGTCACCTTCGAACAGGTTGACGTTACCGAt                                                                               >  7:76170/1‑148 (MQ=255)
                                                                         gggCGTAGGTGATGTGCTTGTCTTCCACCGAGGTGGTGATGCCGTGGCCCACGTAGATCTTGCGCTCCAGCTCCGGGCTGGCAATGATCGCGTAGCCTTCGGCGTCGTCGACCACGTCACCTTCGAACAGGTTGACGTTACCGATGAAc                                                                           >  1:290493/1‑149 (MQ=255)
                                                                          ggCGTAGGTGATGTGCTTGTCTTCCACCGAGGTGGTGATGCCGTGGCCCACGTAGATCTTGCGCTCCAGCTCCGGGCTGGCAATGATCGCGTAGCCTTCGGCGTCGTCGACCACGTCACCTTCGAACAGGTTGACGTTACCGATGAACt                                                                          >  3:92605/1‑149 (MQ=255)
                                                                              tAGGTGATGTGCTTGTCTTCCACCGAGGTGGTGATGCCGTGGCCCACGTAGATCTTGCGCTCCAGCTCCGGGCTGGCAATGATCGCGTAGCCTTCGGCGTCGTCGACCACGTCACCTTCGAACAGGTTGACGTTACCGATGAACTCGca                                                                      >  7:231286/1‑149 (MQ=255)
                                                                              tAGGTGATGTGCTTGTCTTCCACCGAGGTGGTGATGCCGTGGCCCACGTAGATCTTGCGCTCCAGCTCCGGGCTGGCAATGATCGCGTAGCCTTCGGCGTCGTCGACCACGTCACCTTCGAACAGGTTGACGTTACCGATGAACTCGca                                                                      <  2:63413/149‑1 (MQ=255)
                                                                                  tGATGTGCTTGTCTTCCACCGAGGTGGTGATGCCGTGGCCCACGTAGATCTTGCGCTCCAGCTCCGGGCTGGCAATGATCGCGTAGCCTTCGGCGTCGTCGACCACGTCACCTTCGAACAGGTTGACGTTACCGAt                                                                               >  5:140066/1‑136 (MQ=255)
                                                                                  tGATGTGCTTGTCTTCCACCGAGGTGGTGATGCCGTGGCCCACGTAGATCTTGCGCTCCAGCTCCGGGCTGGCAATGATCGCGTAGCCTTCGGCGTCGTCGACCACGTCACCTTCGAACAGGTTGACGTTACCGAt                                                                               <  6:140066/136‑1 (MQ=255)
                                                                                  tGATGTGCTTGTCTTCCACCGAGGTGGTGATGCCGTGGCCCACGTAGATCTTGCGCTCCAGCTCCGGGCTGGCAATGATCGCGTAGCCTTCGGCGTCGTCGACCACGTCACCTTCGAACAGGTTGACGTTACCGATGAACTCGCacac                                                                   >  3:156086/1‑148 (MQ=255)
                                                                                       tgCTTGTCTTCCACCGAGGTGGTGATGCCGTGGCCCACGTAGATCTTGCGCTCCAGCTCCGGGCTGGCAATGATCGCGTAGCCTTCGGCGTCGTCGACCACGTCACCTTCGAACAGGTTGACGTTACCGATGAACTCGCGCACCAggc                                                              <  2:325674/148‑1 (MQ=255)
                                                                                        gCTTGTCTTCCACCGAGGTGGTGATGCCGTGGCCCACGTAGATCTTGCGCTCCAGCTCCGGGCTGGCAATGATCGCGTAGCCTTCGGCGTCGTCGACCACGTCACCTTCGAACAGGTTGACGTTACCGATGAACTCGCACACCAggcgg                                                            >  7:345189/1‑149 (MQ=255)
                                                                                         cTTGTCTTCCACCGAGGTGGTGATGCCGTGGCCCACGTAGATCTTGCGCTCCAGCTCCGGGCTGGCAATGATCGCGTAGCCTTCGGCGTCGTCGACCACGTCACCTTCGAACAGGTTGACGTTACCGATGAACTCGCACACCAggcggc                                                           >  6:64689/1‑149 (MQ=255)
                                                                                             tCTTCCACCGAGGTGGTGATGCCGTGGCCCACGTAGATCTTGCGCTCCAGCTCCGGGCTgg                                                                                                                                               <  1:31408/61‑1 (MQ=255)
                                                                                             tCTTCCACCGAGGTGGTGATGCCGTGGCCCACGTAGATCTTGCGCTCCAGCTCCGGGCTgg                                                                                                                                               >  2:31408/1‑61 (MQ=255)
                                                                                               ttCCACCGAGGTGGTGATGCCGTGGCCCACGTAGATCTTGCGCTCCAGCTCCGGGCTGGCAATGATCGCGTAGCCTTCGGCGTCGTCGACCACGTCACCTTCGAACAGGTTGACGTTACCGATGAACTCGCACACCAGGCGGCTGGTAg                                                     >  5:227699/1‑149 (MQ=255)
                                                                                               ttCCACCGAGGTGGTGATGCCGTGGCCCACGTAGATCTTGCGCTCCAGCTCCGGGCTGGCAATGATCGCGTAGCCTTCGGCGTCGTCGACCACGTCACCTTCGAACAGGTTGACGTTACCGATGAACTCGCACACCAGGCGGCTGGGAg                                                     >  6:304517/1‑149 (MQ=255)
                                                                                               ttCCACCGAGGTGGTGATGCCGTGGCCCACGTAGATCTTGAGCTCCAGCTCCGGGCTGGCAATGATCGCGTAGCCTTCGGCGTCGTCGACCACGTCACCTTCGAACAGGTTGa                                                                                         >  5:108434/1‑113 (MQ=255)
                                                                                                 ccACCGAGGTGGTGATGCCGTGGCCCACGTAGATCTTGCGCTCCAGCTCCGGGCTGGCAATGATCGCGTAGCCTTCGGCGTCGTCGACCACGTCACCTTCGAACAGGTTGACGTTACCGATGAACTCGCACACCAGGCGGCTGGTAg                                                     <  4:92605/147‑1 (MQ=255)
                                                                                                 ccACCGAGGTGGTGATGCCGTGGCCCACGTAGATCTTGCGCTCCAGCTCCGGGCTGGCAATGATCGCGTAGCCTTCGGCGTCGTCGACCACGTCACCTTCGAACAGGTTGACGTTACCGATGAACTCGCACACCAGGCGGCTGGTAg                                                     <  7:59045/147‑1 (MQ=255)
                                                                                                  cACCGAGGTGGTGATGCCGTGGCCCACGTAGATCTTGCGCTCCAGCTCCGGGCTGGCAATGATCGCGTAGCCTTCGGCGTCGTCGACCACGTCACCTTCGAACAGGTTGACGTTACCGATGAACTCGCACACCAGGCGGCTGGTAGGCg                                                  >  1:136708/1‑149 (MQ=255)
                                                                                                      gAGGTGGTGATGCCGTGGCCCACGTAGATCTTGCGCTCCAGCTCCGGGCTGGCAATGATCGCGTAGCCTTCGGCGTCGTCGACCACGTCACCTTCGAACAGGTTGACGTTACCGATGAACTCGCACACCAGGCGGCTGGTAGGCGtctc                                              <  8:180263/149‑1 (MQ=255)
                                                                                                          tggtgATGCCGTGGCCCACGTAGATCTTGCGCTCCAGCTCCGGGCTGGCAATGATCGCGTAGCCTTCGGCGTCGTCGACCACGTCACCTTCGAACAGGTTGACGTTACCGATGAACTCGCACACCAGGCGGCTGGTAGGCGTCTCGTAg                                          <  8:345189/149‑1 (MQ=255)
                                                                                                            gtgATGCCGTGGCCCACGTAGATCTTGCGCTCCAGCTCCGGGCTGGCAATGATCGCGTAGCCTTCGGCGTCGTCGACCACGTCACCTTCGAACAGGTTGACGTTACCGATGAACTCGCACACCAGGCGGCTGGTAGGCGTCTCGTAGAt                                        <  6:227699/149‑1 (MQ=255)
                                                                                                             tgATGCCGTGGCCCACGTAGATCTTGCGCTCCAGCTCCGGGCTGGCAATGATCGCGTAGCCTTCGGCGTCGTCGACCACGTCACCTTCGAACAGGTTGACGTTACCGATGAACTCGCACACCAGGCGGCTGGTAGGCGTCTCGTAGATg                                       >  8:294093/1‑149 (MQ=255)
                                                                                                                 gCCGTGGCCCACGTAGATCTTGCGCTCCAGCTCCGGGCTGGCAATGATCGCGTAGCCTTCGGCGTCGTCGACCACGTCACCTTCGAACAGGTTGACGTTACCGATGAACTCGCACACCAGGCGGCTGGt                                                       >  6:268100/1‑129 (MQ=255)
                                                                                                                 gCCGTGGCCCACGTAGATCTTGCGCTCCAGCTCCGGGCTGGCAATGATCGCGTAGCCTTCGGCGTCGTCGACCACGTCACCTTCGAACAGGTTGACGTTACCGATGAACTCGCACACCAGGCGGCTGGt                                                       <  5:268100/129‑1 (MQ=255)
                                                                                                                 gCCGTGGCCCACGTAGATCTTGCGCTCCAGCTCCGGGCTGGCAATGATCGCGTAGCCTTCGGCGTCGTCGACCACGTCACCTTCGAACAGGTTGACGTTACCGATGAACTCGCACACCAGGCGGCTGGTAGGCGTCTCGTAGATGTCCa                                   >  7:181910/1‑149 (MQ=255)
                                                                                                                          cACGTAGATCTTGCGCTCCAGCTCCGGGCGGGCAATGATCGCGTAGCCTTCGGCGTCGTCGACCACGTCACCTTCGAACAGGTTGACGTTACCGATGAACTCGCACACCAGGCGGCTGGTAGGCGTCTCGTAGATGTCCACTGGCGAGc                          <  6:97930/149‑1 (MQ=255)
                                                                                                                                  tCTTGCGCTCCAGCTCCGGGCTGGCAATGATCGCGTAGCCTTCGGCGTCGTCGACCACGTCACCTTc                                                                                                    >  5:332315/1‑67 (MQ=255)
                                                                                                                                  tCTTGCGCTCCAGCTCCGGGCTGGCAATGATCGCGTAGCCTTCGGCGTCGTCGACCACGTCACCTTc                                                                                                    <  6:332315/67‑1 (MQ=255)
                                                                                                                                    ttGCGCTCCAGCTCCGGGCTGGCAATGATCGCGTAGCCTTCGGCGTCGTCGACCACGTCACCTTCGAACAGGTTGACGTTACCGATGAACTCGCACACCAGGCGGCTGGTAGGCGTCTCGTAGATGTCCACTGGCGAGCCGATCTGGGc                >  6:203680/1‑149 (MQ=255)
                                                                                                                                       cgcTCCAGCTCCGGGCTGGCAATGATCGCGTAGCCTTCGGCGTCGTCGACCACGTCACCTTCGAACAGGTTGACGTTACCGATGAACTCGCACACCAGTCGGCTGGTAGGCGTCTCGTAGATGTCCACTGGCGAGCCGATCTGGGCGa              <  2:290493/148‑1 (MQ=255)
                                                                                                                                       cgcTCCAGCTCCGGGCTGGCAATGATCGCGTAGCCTTCGGCGTCGTCGACCACGTCACCTTCGAACAGGTTGACGTTACCGATGAACTCGCACACCAGGCGGCTGGTAGGCGTCTCGTAGATGTCCACTGGCGAGCCGATCTGGGCGa              <  8:324259/148‑1 (MQ=255)
                                                                                                                                        ttgccAGCCCGGAGCTGGCAATGATCGCGTAGCCTTCGGCGTCGTCGACCACGTCACCTTCGAACAGGTTGACGTTACCGATGAACTCGCACACCAGGCGGCTGGTAGGCGTCTCGTAGATGTCCACTGGc                              <  4:265821/128‑1 (MQ=255)
                                                                                                                                        ttgccAGCCCGGAGCTGGCAATGATCGCGTAGCCTTCGGCGTCGTCGACCACGTCACCTTCGAACAGGTTGACGTTACCGATGAACTCGCACACCAGGCGGCTGGTAGGCGTCTCGTAGATGTCCACTGGc                              >  3:265821/4‑131 (MQ=255)
                                                                                                                                        gcTCCAGCTCCGGGCTGGCAATGATCGCGTAGCCTTCGGCGTCGTCGACCACGTCACCTTCGAACAGGTTGACGTTACCGATGAACTCGCACACCAGGCGGCTGGTAGGCGTCTCGTAGATGTCCACTGGCGAGCCGATCTGGGCGATc            >  5:230293/1‑149 (MQ=255)
                                                                                                                                        gcTCCAGCTCCGGGCTGGCAATGATCGCGTAGCCTTCGGCGTCGTCGACCACGTCACCTTCGAACAGGTTGACGTTACCGATGAACTCGCACACCAGGCGGCTGGTAGGCGTCTCGTAGATGTCCACTGGCGAGCCGATCTGGGCGATc            <  2:295428/149‑1 (MQ=255)
                                                                                                                                           ccAGCTCCGGGCTGGCAATGATCGCGTAGCCTTCGGCGTCGTCGACCACGTCACCTTCGAACAGGTTGACGTTACCGATGAACTCGTACACCAGGCGGCTGGTAGGCGTCTCGTAGATGTCCACTGGCGAGCCGATCTGGGCGATCCAg         >  8:62545/1‑149 (MQ=255)
                                                                                                                                             aGCTCCGGGCTGGCAATGATCGCGTAGCCTTCGGCGTCGTCGACCACGTCACCTTCGAACAGGTTGACGTTACCGATGAACTCGCACACCAGGCGGCTGGTAGGCGTCTCGTAGATGTCCACTGGCGAGCCGATCTGGGCGATCCAGcc       >  7:22433/1‑149 (MQ=255)
                                                                                                                                              gCTCCGGGCTGGCAATGATCGCGTAGCCTTCGGCGTCGTCGACCACGTCACCTTCGAACAGGTTGACGTTACCGATGAACTCGCACACCAGGCGGCTGGTAGGCGTCTCGTAGATGTCCACTGGCGAGCCGATCTGGGCGATCCAGcc       >  7:75885/1‑148 (MQ=255)
                                                                                                                                               cTCCGGGCTGGCAATGATCGCGTAGCCTTCGGCGTCGTCGACCACGTCACCTTCGAACAGGTTGACGTTACCGATGAACTCGCACACCAGGCGGCTGGTAGGCGTCTCGTAGATGTCCACTGGCGAGCCGATCTGGGCGATCCAGcc       >  7:105916/1‑147 (MQ=255)
                                                                                                                                                tCCGGGCTGGCAATGATCGCGTAGCCTTCGGCGTCGTCGACCACGTCACCTTCGAACAGGTTGACGTTACCGATGAACTCGCACACCAGGCGGCTGGTAGGCGTCTCGTAGATGTCCACTGGCGAGCCGATCTGGGCGATCCAGCCCAg    >  3:355504/1‑149 (MQ=255)
                                                                                                                                                  cGGGCTGGCAATGATCGCGTAGCCTTCGGCGTCGTCGACCACGTCACCTTCGAACAGGTTGACGTTACCGATGAACTCGTACACCAGGCGGCTGGTAGGCGTCTCGTAGATGTCCACTGGCGAGCCGATCTGGGCGATCCAGCCCAg    <  7:62545/147‑1 (MQ=255)
                                                                                                                                                  cGGGCTGGCAATGATCGCGTAGCCTTCGGCGTCGTCGACCACGTCACCTTCGAACAGGTTGACGTTACCGATGAACTCGCACACCAGGCGGCTGGTAGGCGTCTCGTAGATGTCCACTGGCGAGCCGATCTGGGCGATCCAGCCCAg    <  2:173055/147‑1 (MQ=255)
                                                                                                                                                  cGGGCTGGCAATGATCGCGTAGCCTTCGGCGTCGTCGACCACGTCACCTTCGAACAGGTTGACGTTACCGATGAACTCGCACACCAGGCGGCTGGTAGGCGTCTCGTAGATGTCCACTGGCGAGCCGATCTGGGCGATCCAGCCCAGGt  <  2:136708/149‑1 (MQ=255)
                                                                                                                                                  |                                                                                                                                                    
TGGATCTTGCCGCGCGACCAGTTGTGCTCGCAGGTTGGCTGCTGGGTGGTGACCAGCATCTTTTCCGGGCGCAGGGCGTAGGTGATGTGCTTGTCTTCCACCGAGGTGGTGATGCCGTGGCCCACGTAGATCTTGCGCTCCAGCTCCGGGCTGGCAATGATCGCGTAGCCTTCGGCGTCGTCGACCACGTCACCTTCGAACAGGTTGACGTTACCGATGAACTCGCACACCAGGCGGCTGGTAGGCGTCTCGTAGATGTCCACTGGCGAGCCGATCTGGGCGATCCAGCCCAGGT  >  NC_002947/5905827‑5906121

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: