Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,907,906:1 +G 5.5% coding (61/1095 nt) spuE ← spermidine/putrescine ABC transporter substrate‑binding protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,907,9061.G 5.5% 373.5 / 13.5 91coding (61/1095 nt)spuEspermidine/putrescine ABC transporter substrate‑binding protein
Reads supporting (aligned to +/- strand):  ref base . (42/44);  new base G (2/3);  total (44/47)
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.06e-01

CAGCAGCTTGCCTTCCAGGGTTTCGTTGGAGTCGAACACATCTTGCACAGGCTTGATGCCCGTGGCTTTCTCGAAGTCCGCGAGTGTGGTGGGGCCGATGTAGTCGGACCAGTTATAGAAGTGCACCGTCGGCGCCGCTTGGACGCTGCATGCCAGCGTCAGGCCCGCGCCAGCCATCAAGGCCTTGCGGAATACAGAAATAGACAAGTGGGGGTCCTCACAATCAGTGCCTGGGTAGCGACAATGCTGCCGCACACGCAAAACCGGCGCGCAACTTAACTTCGCAGCTTCGAT  >  NC_002947/5907760‑5908053
                                                                                                                                                   |                                                                                                                                                   
cagcagCTTGCCTTCCAGGGTTTCGTTGGAGTCGAACACATCTTGCACAGGCTTGATGCCCGTGGCTTTCTCGAAGTCCGCGAGTGTGGTGGGGCCGATGTAGTCGGACCAGTTATAGAAGTGCACCGTCGGCGCCGCTTGGACGCTg                                                                                                                                                    >  5:43933/1‑148 (MQ=255)
 agcagcTTGCCTTCCAGGGTTTCGTTGGAGTCGAACACATCTTGCACAGGCTTGATGCCCGTGGCTTTCTCGAAGTCCGCGAGTGTGGTGGGGCCGATGTAGTCGGACCAGTTATAGAAGTGCACCGTCGGCGCCGCTTGGACGCTGCa                                                                                                                                                  >  5:14209/1‑149 (MQ=255)
      cTTGCCTTCCAGGGTTTCGTTGGAGTCGAACACATCTTGCACAGGCTTGATGCCCGTGGCTTTCTCGAAGTCCGCGAGTGTGGTGGGGCCGATGTAGTCGGACCAGTTATAGAAGTGCACCGTCGGCGCCGCTTGGACGCTGCATGc                                                                                                                                               >  7:170426/1‑147 (MQ=255)
      cTTGCCTTCCAGGGTTTCGTTGGAGTCGAACACATCTTGCACAGGCTTGATGCCCGTGGCTTTCTCGAAGTCCGCGAGTGTGGTGGGGCCGATGTAGTCGGACCAGTTATAGAAGTGCACCGTCGGCGCCGCTTGGACGCTGCATGCCa                                                                                                                                             <  7:44997/149‑1 (MQ=255)
         gCCTTCCAGGGTTTCGTTGGAGTCGAACACATCTTGCACAGGCTTGATGCCCGTGGCTTTCTCGAAGTCCGCGAGTGTGGTGGGGCCGATGTAGTCGGACCAGTTATAGAAGTGCACCGTCGGCGCCGCTTGGACGCTGCATGCCAGCg                                                                                                                                          >  5:294889/1‑149 (MQ=255)
           cTTCCAGGGTTTCGTTGGAGTCGAACACATCTTGCACAGGCTTGATGCCCGTGGCTTTCTCGAAGTCCGCGAGTGTGGTGGGGCCGATGTAGTCGGACCAGTTATAGAAGTGCACCGTCGGCGCCGCTTGGACGCTGGCATGCAGCGTc                                                                                                                                        <  2:54197/149‑1 (MQ=255)
              ccAGGGTTTCGTTGGAGTCGAACACATCTTGCACAGGCTTGATGCCCGTGGCTTTCTCGAAGTCCGCGAGTGTGGTGGGGCCGATGTAGTCGGACCAGTTATAGAAGTGCACCGTCGGCGCCGCTTGGACGCTGCATGCCAGCGTCAg                                                                                                                                      <  8:13014/148‑1 (MQ=255)
               cAGGGTTTCGTTGGAGTCGAACACATCTTGCACAGGCTTGATGCCCGTGGCTTTCTCGAAGTCCGCGAGTGTGGTGGGGCCGATGTAGTCGGACCAGTTATAGAAGTGCACCGTCGGCGCCGCTTGGACGCTGGCATGCAGCGTccaa                                                                                                                                     <  8:55814/148‑4 (MQ=255)
               cAGGGTTTCGTTGGAGTCGAACACATCTTGCACAGGCTTGATGCCCGTGGCTTTCTCGAAGTCCGCGAGTGTGGTGGGGCCGATGTAGTCGGACCAGTTATAGAAGTGCACCGTCGGCGCCGCTTGGACGCTGGCATGCAGCGTccaa                                                                                                                                     <  2:111266/148‑4 (MQ=255)
               cAGGGTTTCGTTGGAGTCGAACACATCTTGCACAGGCTTGATGCCCGTGGCTTTCTCGAAGTCCGCGAGTGTGGTGGGGCCGATGTAGTCGGACCAGTTATAGAAGTGCACCGTCGGCGCCGCTTGGACGCTCCATGCGCGCGGCAGGc                                                                                                                                    >  8:109515/1‑149 (MQ=255)
                 gggTTTCGTTGGAGTCGAACACATCTTGCACAGGCTTCATGCCCGTGGCTTTCTCGAAGTCCGCGAGTGTGGTGGGGCCGATGTAGTCGCACCACTTATAGAAATGCACCGTCGGCGCCCCTTGGACGCTGCATGCCAGCGTCAGGccc                                                                                                                                  >  5:288195/1‑149 (MQ=255)
                        gTTGGAGTCGAACACATCTTGCACAGGCTTGATGCCCGTGGCTTTCTCGAAGTCCGCGAGTGTGGTGGGGCCGATGTAGTCGGACCAGTTATAGAAGTGCACCGTCGGCGCCGCTTGGACGCTGCATGCCAGCGTCAGGCCCGCgccag                                                                                                                           >  4:178632/1‑149 (MQ=255)
                        gTTGGAGTCGAACACATCTTGCACAGGCTTGATGCCCGTGGCTTTCTCGAAGTCCGCGAGTGTGGTGGGGCCGATGTAGTCGGACCAGTTATAGAAGTGCACCGTCGGCGCCGCTTGGACGCTGCATGCCAGCGTCAGGCCCGCgccag                                                                                                                           >  2:393459/1‑149 (MQ=255)
                        gTTGGAGTCGAACACATCTTGCACAGGCTTGATGCCCGTGGCTTTCTCGAAGTCCGCGAGTGTGGTGGGGCCGATGTAGTCGGACCAGTTATAGAAGTGCACCGTCGGCGCCGCTTGGACGCTGCATGCCAGCGTCAGGCCCGCgcca                                                                                                                            <  6:72034/148‑1 (MQ=255)
                            gAGTCGAACACATCTTGCACAGGCTTGATGCCCGTGGCTTTCTCGAAGTCCGCGAGTGTGGTGGGGCCGATGTAGTCGGACCAGTTATAGAAGTGCACCGTCGGCGCCGCTTGGACGCTGCATGCCAGCGTCAGGCCCGCgccagcc                                                                                                                         >  5:354261/1‑147 (MQ=255)
                            gAGTCGAACACATCTTGCACAGGCTTGATGCCCGTGGCTTTCTCGAAGTCCGCGAGTGTGGTGGGGCCGATGTAGTCGGACCAGTTATAGAAGTGCACCGTCGGCGCCGCTTGGACGCTGCATGCCAGCGTCAGGCCCGCgccagcc                                                                                                                         <  6:354261/147‑1 (MQ=255)
                                cGAACACATCTTGCACAGGCTTGATGCCCGTGGCTTTCTCGAAGTCCGCGAGTGTGGTGGGGCCGATGTAGTCGGACCAGTTATAGAAGTGCACCGTCGGCGCCGCTTGGACGCTGCATGCCAGCGTCAGGCCCGCGCCAGCCATCAAg                                                                                                                   <  8:10467/149‑1 (MQ=255)
                                cGAACACATCTTGCACAGGCTTGATGCCCGTGGCTTTCTCGAAGTCCGCGAGTGTGGTGGGGCCGATGTAGTCGGACCAGTTATAGAAGTGCACCGTCGGCGCCGCTTGGACGCTGCATGCCAGCGTCAGGCCCGCGCCAGCCATCAAg                                                                                                                   <  2:275726/149‑1 (MQ=255)
                                   acacATCTTGCACAGGGTTGATGCCCGTGGCTTTCTCGAAGTCCGGGAGTGTGGTGGGGCCGATGTAGTCGGACCAGTTATAGAAGTGCACCGTCGGCGCCGCTTGGACGCTGCATGCCAGCGTCAGGCCCGCGCCAGCCATCAAGGc                                                                                                                 <  4:148305/148‑1 (MQ=255)
                                       aTCTTGCACAGGCTTGATGCCCGTGGCTTTCTCGAAGTCCGCGAGTGTGGTGGGGCCGATGTAGTCGGACCAGTTATAGAAGTGCACCGTCGGCGCCGCTTGGACGCTGCATGCCAGCGTCAGGCCCGCGCCAGCCATCaa                                                                                                                    >  2:395726/1‑141 (MQ=255)
                                       aTCTTGCACAGGCTTGATGCCCGTGGCTTTCTCGAAGTCCGCGAGTGTGGTGGGGCCGATGTAGTCGGACCAGTTATAGAAGTGCACCGTCGGCGCCGCTTGGACGCTGCATGCCAGCGTCAGGCCCGCGCCAGCCATCaa                                                                                                                    <  1:395726/141‑1 (MQ=255)
                                       aTCTTGCACAGGCTTGATGCCCGTGGCTTTCTCGAAGTCCGCGAGTGTGGTGGGGCCGATGTAGTCGGACCAGTTATAGAAGTGCACCGTCGGCGCCGCTTGGACGCTGCATGCCAGCGTCAGGCCCGCGCCAGCCATCAAGGCCTTGc                                                                                                            <  2:19224/149‑1 (MQ=255)
                                       aTCTTGCACAGGCTTGATGCCCGTGGCTTTCTCGAAGTCCGCGAGTGTGGTGGGGCCGATGTAGTCGGACCAGTTATAGAAGTGCACCGTCGGCGCCGCTTGGACGCTGCATGCCAGCGTCAGGCCCGCGCCAGCCATCAAGGCCTTGc                                                                                                            >  4:86555/1‑149 (MQ=255)
                                       aTCTTGCACAGGCTTGATGCCCGTGGCTTTCTCGAAGTCCGCGAGTGTGGTGGGGCCGATGTAGTCGGACCAGTTATAGAAGTGCACCGTCGGCGCCGCTTGGACGCTGCATGCCAGCGTCAGGCCCGCGCCAGCCATCAAGGCCTTGc                                                                                                            >  1:311391/1‑149 (MQ=255)
                                        tCTTGCACAGGCTTGATGCCCGTGGCTTTCTCGAAGTCCGCGAGTGTGGTGGGGCCGATGTAGTCGGACCAGTTATAGAAGTGCACCGTCGGCGCCGCTTGGACGCTGCATGCCAGCGTCAGGCCCGCGCCAGCCATCAAGGCCTTGCg                                                                                                           >  3:24421/1‑149 (MQ=255)
                                            gCACAGGCTTGATGCCCGTGGCTTTCTCGAAGTCCGCGAGTGTGGTGGGGCCGATGTAGTCGGACCAGTTATAGAAGTGCACCGTCGGCGCCGCTTGGACGCTGCATGCCAGCGTCAGGCCCGCGCCAGCCATCAAGGCCTTGCg                                                                                                           <  6:246724/145‑1 (MQ=255)
                                            gCACAGGCTTGATGCCCGTGGCTTTCTCGAAGTCCGCGAGTGTGGTGGGGCCGATGTAGTCGGACCAGTTATAGAAGTGCACCGTCGGCGCCGCTTGGACGCTGCATGCCAGCGTCAGGCCCGCGCCAGCCATCAAGGCCTTGCGGaa                                                                                                        <  7:176035/148‑1 (MQ=255)
                                            gCACAGGCTTGATGCCCGTGGCTTTCTCGAAGTCCGCGAGTGTGGTGGGGCCGATGTAGTCGGACCAGTTATAGAAGTGCACCGTCGGCGCCGCTTGGACGCTGCATGCCAGCGTCAGGCCCGCGCCAGCCATCAAGGCCTTGCGGaa                                                                                                        >  7:338880/1‑148 (MQ=255)
                                             cacaGGCTTGATGCCCGTGGCTTTCTCGAAGTCCGCGAGTGTGGTGGGGCCGATGTAGTCGGACCAGTTATAGAAGTGCACCGTCGGCGCCGCTTGGACGCTGCATGCCAGCGTCAGGCCCGCGCCAGCCATCAAGGCCTTGCGGAAt                                                                                                       <  6:356251/148‑1 (MQ=255)
                                             cacaGGCTTGATGCCCGTGGCTTTCTCGAAGTCCGCGAGTGTGGTGGGGCCGATGTAGTCGGACCAGTTATAGAAGTGCACCGTCGGCGCCGCTTGGACGCTGCATGCCAGCGTCAGGCCCGCGCCAGCCATCAAGGCCTTGCGGAAt                                                                                                       <  1:393459/148‑1 (MQ=255)
                                                 ggCTTGATGCCCGTGGCTTTCTCGAAGTCCGCGAGTGTGGTGGGGCCGATGTAGTCGGACCAGTTATAGAAGTGCACCGTCGGCGCCGCTTGGACGCTGCATGCCAGCGTCAGGCCCGCGCCAGCCATCAAGGCCTTGCGGAATACAGa                                                                                                  <  3:86555/149‑1 (MQ=255)
                                                  gCTTGATGCCCGTGGCTTTCTCGAAGTCCGCGAGTGTGGTGGGGCCGATGTAGTCGGACCAGTTATAGAAGTGCACCGTCGGCGCCGCTTGGACGCTGCATGCCAGCGTCAGGCCCGCGCCAGCCATCaa                                                                                                                    >  7:363036/1‑130 (MQ=255)
                                                  gCTTGATGCCCGTGGCTTTCTCGAAGTCCGCGAGTGTGGTGGGGCCGATGTAGTCGGACCAGTTATAGAAGTGCACCGTCGGCGCCGCTTGGACGCTGCATGCCAGCGTCAGGCCCGCGCCAGCCATCaa                                                                                                                    <  8:363036/130‑1 (MQ=255)
                                                    ttGATGCCCGTGGCTTTCTCGAAGTCCGCGAGTGTGGTGGGGCCGATGTAGTCGGACCAGTTATAGAAGTGCACCGTCGGCGCCGCTTGGACGCTGCATGCCAGCGTCAGGCCCGCGCCAGCCATCAAGGCCTTGCGGAATACAGAAAt                                                                                               <  6:14209/149‑1 (MQ=255)
                                                      gATGCCCGTGGCTTTCTCGAAGTCCGCGAGTGTGGTGGGGCCGATGTAGTCGGACCAGTTATAGAAGTGCACCGTCGGCGCCGCTTGGACGCTGCATGCCAGCGTCAGGCCCGCGCCAGCCATCAAGGCCTTGCGGAATACAGAAAt                                                                                               <  5:41121/147‑1 (MQ=255)
                                                      gATGCCCGTGGCTTTCTCGAAGTCCGCGAGTGTGGTGGGGCCGATGTAGTCGGACCAGTTATAGAAGTGCACCGTCGGCGCCGCTTGGACGCTGCATGCCAGCGTCAGGCCCGCGCCAGCCATCAAGGCCTTGCGGAATACAGAAATa                                                                                              <  6:294889/148‑1 (MQ=255)
                                                        tGCCCGTGGCTTTCTCGAAGTCCGCGAGTGTGGTGGGGCCGATGTAGTCGGACCAGTTATAGAAGTGCACCGTCGGCGCCGCTTGGACGCTGCATGCCAGCGTCAGGCCCGCGCCAGCCATCAAGGCCTTGCGGAATACAGAAATAGAc                                                                                           <  1:356700/149‑1 (MQ=255)
                                                        tGCCCGTGGCTTTCTCGAAGTCCGCGAGTGTGGTGGGGCCGATGTAGTCGGACCAGTTATAGAAGTGCACCGTCGGCGCCGCTTGGACGCTGCATGCCAGCGTCAGGCCCGCGCCAGCCATCAAGGCCTTGCGGAATACAGAAATAGAc                                                                                           <  6:43933/149‑1 (MQ=255)
                                                           ccGTGGCTTTCTCGAAGTCCGCGAGTGTGGTGGGGCCGATGTAGTCGGACCAGTTATAGAAGTGCACCGTCGGCGCCCCTTGGACCCTGCATGCCAGCGTCAGGCCCGCGCCAGCCATCAAGGCCTTGCGGAATACAGAAATAGACAAg                                                                                        >  2:326401/1‑149 (MQ=255)
                                                              tGGCTTTCTCGAAGTCCGCGAGTGTGGTGGGGCCGATGTAGTCGGACCAGTTATAGAAGTGCACCGTCGGCGCCGCTTGGACGCTGCATGCCAGCGTCAGGCCCGCGCCAGCCATCAAGGCCTTGCGGAATACAGAAATAGACAAGTgg                                                                                     >  4:211777/1‑149 (MQ=255)
                                                              tGGCTTTCTCGAAGTCCGCGAGTGTGGTGGGGCCGATGTAGTCGGACCAGTTATAGAAGTGCACCGTCGGCGCCGCTTGGACGCTGCATGCCAGCGTCAGGCCCGCGCCAGCCATCAAGGCCTTGCGGAATACAGAAATAGACAAGTgg                                                                                     <  4:24421/149‑1 (MQ=255)
                                                              tGGCTTTCTCGAAGTCCGCGAGTGTGGTGGGGCCGATGTAGTCGGACCAGTTATAGAAGTGCACCGTCGGCGCCGCTTGGACGCTGCATGCCAGCGTCAGGCCCGCGCCAGCCATCAAGGCCTTGCGGAATACAGAAATAGACAAGTgg                                                                                     <  8:36488/149‑1 (MQ=255)
                                                              tGGCTTTCTCGAAGTCCGCGAGTGTGGTGGGGCCGATGTAGTCGGACCAGTTATAGAAGTGCACCGTCGGCGCCGCTTGGACGCTGCATGCCAGCGTCAGGCCCGCGCCAGCCATCAAGGCCTTGCGGAATACAGAAATAGACAAGTgg                                                                                     <  1:49333/149‑1 (MQ=255)
                                                                gCTTTCTCGAAGTCCGCGAGTGTGGTGGGGCCGATGTAGTCGGACCAGTTATAGAAGTGCACCGTCGGCGCCGCTTGGACGCTGCATGCCAGCGTCAGGCCCGCGCCAGCCATCAAGGCCTTGCGGAATACAGAAATAGACAAGTgggg                                                                                   <  1:395478/149‑1 (MQ=255)
                                                                gCTTTCTCGAAGTCCGCGAGTGTGGTGGGGCCGATGTAGTCGGACCAGTTATAGAAGTGCACCGTCGGCGCCGCTTGGACGCTGCATGCCAGCGTCAGGCCCGCGCCAGCCATCAAGGCCTTGCGGAATACAGAAATAGACAAGTgggg                                                                                   <  3:178632/149‑1 (MQ=255)
                                                                gCTTTCTCGAAGTACGCGAGTGTGGGGGGGCCGATGTAGTCGGACCAGTTATAGAAGTGCACCGTCGGCGCCGCTTGGACGCTGCATGCCAGCGTCAGGCCCGCGCCAGCCATCAAGGCCTTGCGGAATACAGAAATAGACAAGTgggg                                                                                   <  7:347485/149‑1 (MQ=255)
                                                                  tttCTCGAAGTCCGCGAGTGTGGTGGGGCCGATGTAGTCGGACCAGTTATAGAAGTGCACCGTCGGCGCCGCTTGGACGCTGCATGCCAGCGTCAGGCCCGCGCCAGCCATCAAGGCCTTGCGGAATACAGAAATAGACAAGTGGGGGt                                                                                 >  7:13237/1‑149 (MQ=255)
                                                                  tttCTCGAAGTCCGCGAGTGTGGTGGGGCCGATGTAGTCGGACCAGTTATAGAAGTGCACCGTCGGCGCCGCTTGGACGCTGCATGCCAGCGTCAAGCCCGCGCCAGCCATCAAGGCCTTGCGGAATACAGAAATAGACAAGTGGGGGt                                                                                 >  6:173710/1‑149 (MQ=255)
                                                                    tctcGAAGTCCGCGAGTGTGGTGGGGCCGATGTAGTCGGACCAGTTATAGAAGTGCACCGTCGGCGCCGCTTGGACGCTGCATGCCAGCGTCAGGCCCGCGCCAGCCATCAAGGCCTTGCGGAATACAGAAATAGACAAGTGGGGGTcc                                                                               >  2:145458/1‑149 (MQ=255)
                                                                      tcGAAGTCCGCGAGTGTGGTGGGGCCGATGTAGTCGGACCAGTTATAGAAGTGCACCGTCGGCGCCGCTTGGACGCTGCATGCCAGCGTCAGGCCCGCGCCAGCCATCAAGGCCTTGCGGAATACAGAAATAGACAAGTGGGGGTCCt                                                                              >  4:159876/1‑148 (MQ=255)
                                                                      tcGAAGTCCGCGAGTGTGGTGGGGCCGATGTAGTCGGACCAGTTATAGAAGTGCACCGTCGGCGCCGCTTGGACGCTGCATGCCAGCGTCAGGCCCGCGCCAGCCATCAAGGCCTTGCGGAATACAGAAATAGACAAGTGGGGGTCCt                                                                              >  3:203888/1‑148 (MQ=255)
                                                                        gAAGTCCGCGAGTGTGGTGGGGCCGATGTAGTCGGACCAGTTATAGAAGTGCACCGTCGGCGCCGCTTGGACGCTGCATGCCAGCGTCAGGCCCGCGCCAGCCATCAAGGCCTTGCGGAATACAGAAATAGACAAGTGGGGGTCCTca                                                                            >  5:18373/1‑148 (MQ=255)
                                                                        gAAGTCCGCGAGTGTGGTGGGGCCGATGTAGTCGGACCAGTTATAGAAGTGCACCGTCGGCGCCGCTTGGACGCTGCATGCCAGCGTCAGGCCCGCGCCAGCCATCAAGGCCTTGCGGAATACAGAAATAGACAAGTGGGGGTCCTc                                                                             <  3:159876/147‑1 (MQ=255)
                                                                         aaGTCCGCGAGTGTGGTGGGGCCGATGTAGTCGGACCAGTTATAGAAGTGCACCGTCGGCGCCGCTTGGACGCTGCATGCCAGCGTCAGGCCCGCGCCAGCCATc                                                                                                                      <  5:326533/105‑1 (MQ=255)
                                                                         aaGTCCGCGAGTGTGGTGGGGCCGATGTAGTCGGACCAGTTATAGAAGTGCACCGTCGGCGCCGCTTGGACGCTGCATGCCAGCGTCAGGCCCGCGCCAGCCATc                                                                                                                      >  6:326533/1‑105 (MQ=255)
                                                                              cgcgAGTGTGGTGGGGCCGATGTAGTCGGACCAGTTATAGAAGTGCACCGTCGGCGCGGCTTGGACGCTGCATGCCAGCGTCAGGCCCGCGCCAGCCATCATGGCCTTGCGGAATACAGAAATAGACAAGTGGGGGTCCTCACAAtga                                                                      <  6:288195/148‑3 (MQ=255)
                                                                              cgcgAGTGTGGTGGGGCCGATGTAGTCGGACCAGTTATAGAAGTGCACCGTCGGCGCCGCTTGGACGCTGCATGCCAGCGTCAGGCCCGCGCCAGCCATCAAGGCCTTGCGGAATACAGAAATAGAc                                                                                           >  6:362197/1‑127 (MQ=255)
                                                                              cgcgAGTGTGGTGGGGCCGATGTAGTCGGACCAGTTATAGAAGTGCACCGTCGGCGCCGCTTGGACGCTGCATGCCAGCGTCAGGCCCGCGCCAGCCATCAAGGCCTTGCGGAATACAGAAATAGAc                                                                                           <  5:362197/127‑1 (MQ=255)
                                                                               gcgAGTGTGGTGGGGCCGATGTAGTCGGACCAGTTATAGAAGTGCACCGTCGGCGCCGCTTGGACGCTGCATGCCAGCGTCAGGCCCGCGCCAGCCATCAAGGCCTTGCGGAATACAGAAATAGACAAGTGGGGGTCCTCACAATCAGt                                                                    <  3:211777/149‑1 (MQ=255)
                                                                                 gAGTGTGGTGGGGCCGATGTAGTCGGACCAGTTATAGAAGTGCACCGTCGGCGCCGCTTGGACGCTGCATGCCAGCGTCAGGCCCGCGCCAGCCATCAAGGCCTTGCGGAATACAGAAATAGACAAGTGGGGGTCCTCACAATCAGTGc                                                                  >  1:143087/1‑149 (MQ=255)
                                                                                          ggggCCGATGTAGTCGGACCAGTTATAGAAGTGCACCGTCGGCGCCGCTTGGACGCTGCATGCCAGCGTCAGGCCCGCGCCAGCCATCAAGGCCTTGCGGAATACAGAAATAGACAAGTGGGGGTCCTCACAATCAGTGCCTGGGTAGc                                                         <  1:326401/149‑1 (MQ=255)
                                                                                           gggCCGATGTAGTCGGACCAGTTATAGAAGTGCACCGTCGGCGCCGCTTGGACGCTGCATGCCAGCGTCAGGCCCGCGCCAGCCATCAAGGCCTTGCGGAATAGTGAAATAGACAAGTGGGGGTCCTCACAATCAGTGCCTGGGTAGCg                                                        <  1:278054/149‑1 (MQ=255)
                                                                                               cGATGTAGTCGGACCAGTTATAGAAGTGCACCGTCGGCGCCGCTTGGACGCTGCATGCCAGCGTCAGGCCCGCGCCAGCCATCAAGGCCTTGCGGAATACAGAAATAGACAAGTGGGGGTCCTCACAATCAGTGCCTGGGTAGCGACaa                                                    >  7:20848/1‑149 (MQ=255)
                                                                                                gATGTAGTCGGACCAGTTATAGAAGTGCACCGTCGGCGCCGCTTGGACGCTGCATGCCAGCGTCAGGCCCGCGCCAGCCATCAAGGCCTTGCGGAATACAGAAATAGACAAGTGGGGGTCCTCACAATCAGTGCCTGGGTAGCGACAAt                                                   <  8:170426/149‑1 (MQ=255)
                                                                                                      gTCGGACCAGTTATAGAAGTGCACCGTCGGCGCCGCTTGGACGCTGCATGCCAGCGTCAGGCCCGCGCCATCCATCAAGGCCTTGCGGAATACAGAAATAGACAAGTGGGGGTCCTCACAATCAGTGCCTGGGTAGCGACAAtgctgc                                              <  5:200766/148‑1 (MQ=255)
                                                                                                      gTCGGACCAGTTATAGAAGTGCACCGTCGGCGCCGCTTGGACGCTGCATGCCAGCGTCAGGCCCGCGCCAGCCATCAAGGCCTTGCGGAATACAGAAATAGACAAGTGGGGGTCCTCACAATCAGTGCCTGGGTAGCGACAAtgctgc                                              >  6:200766/1‑148 (MQ=255)
                                                                                                      gTCGGACCAGTTATAGAAGTGCACCGTCGGCGCCGCTTGGACGCTGCATGCCAGCGTCAGGCCCGCGCCAGCCATCAAGGCCTTGCGGAATACAGAAATAGACAAGTGGGGGTCCTCACAATCAGTGCCTGGGTAGCGACAATGCTGcc                                             >  4:19572/1‑149 (MQ=255)
                                                                                                      gTCGGACCAGTTATAGAAGTGCACCGTCGGCGCCGCTTGGACGCTGCATGCCAGCGTCAGGCCCGCGCCAGCCATCAAGGCCTTGCGGAATACAGAAATAGACAAGTGGGGGTCCTCACAATCAGTGCCTGGGTAGCGACAATGCTGcc                                             >  5:83865/1‑149 (MQ=255)
                                                                                                       tCGGACCAGTTATAGAAGTGCACCGTCGGCGCCGCTTGGACGCTGCATGCCAGCGTCAGGCCCGCGCCAGCCATCAAGGCCTTGCGGAATACAGAAATAGACAAGTGGGGGTCCTCACAATCAGTGCCTGGGTAGCGACAATGCTGCCg                                            >  5:219655/1‑149 (MQ=255)
                                                                                                          gACCAGTTATAGAAGTGCACCGTCGGCGCCGCTTGGACGCTGCATGCCAGCGTCAGGCCCGCGCCAGCCATCAACGCCTTGCGGAATACAGAAATAGACAAGTGGGGGTCCTCACAATCAGTGCCTCGGTAGCCACAATGCTGCCGca                                          >  5:31453/1‑148 (MQ=255)
                                                                                                              aGTTATAGAAGTGCACCGTCGGCGCCGCTTGGACGCTGCATGCCAGCGTCAGGCCCGCGCCAGCCATCAAGGCCTTGCGGAATACAGAAATAGACAAGTGGGGGTCCTCACAATCAGTGCCTGGGTAGCGACAATGCTGCCGCACACg                                      >  5:197993/1‑148 (MQ=255)
                                                                                                              aGTTATAGAAGTGCACCGTCGGCGCCGCTTGGACGCTGCATGCCAGCGTCAGGCCCGCGCCAGCCATCAAGGCCTTGCGGAATACAGAAATAGACAAGTGGGGGTCCTCACAATCAGTGCCTGGGTAGCGACAATGCTGCCGCACACGc                                     <  7:109515/149‑1 (MQ=255)
                                                                                                               gTTATAGAAGTGCACCGTCGGCGCCGCTTGGACGCTGCATGCCAGCGTCAGGCCCGCGCCAGCCATCAAGGCCTTGCGGAATACAGAAATAGACAAGTGGGGGTCCTCACAATCAGTGCCTGGGTAGCGACAATGCTGCCGCACACGc                                     >  5:157615/1‑148 (MQ=255)
                                                                                                               gTTATAGAAGTGCACCGTCGGCGCCGCTTGGACGCTGCATGCCAGCGTCAGGCCCGCGCCAGCCATCAAGGCCTTGCGGAATACAGAAATAGACAAGTGGGGGTCCTCACAATCAGTGCCTGGGTAGCGACAATGCTGCCGCACACGc                                     >  6:391732/1‑148 (MQ=255)
                                                                                                                   taGAAGTGCACCGTCGGCGCCGCTTGGACGCTGCATGCCAGCGTCAGGCCCGCGCCAGCCATCAAGGCCTTGCGGAATACAGAAATAGACAAGTGGGGGTCCTCACAATCAGTGCCTGGGTAGCGACAATGCTGCCGCACACGCaaaa                                 <  8:338880/148‑1 (MQ=255)
                                                                                                                   taGAAGTGCACCGTCGGCGCCGCTTGGACGCTGCATGCCAGCGTCAGGCCCGCGCCAGCCATCAAGGCCTTGCGGAATACAGAAATAGACAAGTGGGGGTCCTCACAATCAGTGCCTGGGTAGCGACAATGCTGCCGCACACGCAAAAc                                >  7:258488/1‑149 (MQ=255)
                                                                                                                   taGAAGTGCACCGTCGGCGCCGCTTGGACGCTGCATGCCAGCGTCAGGCCCGCGCCAGCCATCAAGGCCTTGCGGAATACAGAAATAGACAAGTGGGGGTCCTCACAATCAGTGCCTGGGTAGCGACAATGCTGCCGCACACGCAAAAc                                >  6:359416/1‑149 (MQ=255)
                                                                                                                     gAAGTGCACCGTCGGCGCCGCTTGGACGCTGCATGCCAGCGTCAGGCCCGCGCCAGCCATCAAGGCCTTGCGGAATACAGAAATAGACAAGTGGGGGTCCTCACAATCAGTGCCTGGGTAGCGACAATGCTGCCGCACACGCAAAACCg                              <  5:250630/149‑1 (MQ=255)
                                                                                                                           cACCGTCGGCGCCGCTTGGACGCTGCATGCCAGCGTCAGGCCCGCGCCAGCCATCAAGGCCTTGCGGAATACAGAAATAGACAAGTGGGGGTCCTCACAATCAGTGCCTGGGTAGCGACAATGCTGCCGCACACGCAAAACCGgcgcg                         >  1:341247/1‑148 (MQ=255)
                                                                                                                                tCGGCGCCGCTTGGACGCTGCATGCCAGCGTCAGGCCCGCGCCAGCCATCAAGGCCTTGCGGAATACAGAAATAGACAAGTGGGGGTCCTCACAATCAGTGCCTGGGTAGCGACAATGCTGCCGCACACGCAAAACCGGCGCGCaactt                   <  1:160193/149‑1 (MQ=255)
                                                                                                                                  ggCGCCGCTTGGACGCTGCATGCCAGCGTCAGGCCCGCGCCAGCCATCAAGGCCTTGCGGAATACAGAAATAGACAAGTGGGGGTCCTCACAATCAGTGCCTGGGTAGCGACAATGCTGCCGCACACGCAAAACCGGCGCGCaacttaa                 <  5:391732/149‑1 (MQ=255)
                                                                                                                                  ggCGCCGCTTGGACGCTGCATGCCAGCGTCAGGCCCGCGCCAGCCATCAAGGCCTTGCGGAATACAGAAATAGACAAGTGGGGGTCCTCACAATCAGTGCCTGGGTAGCGACAATGCTGCCGCACACGCAAAACCGGCGCGCaactta                  >  6:122484/1‑148 (MQ=255)
                                                                                                                                        gcTTGGACGCTGCATGCCAGCGTCAGGCCCGCGCCAGCCATCAAGGCCTTGCGGAATACAGAAATAGACAAGTGGGGGTCCTCACAATCAGTGCCTGGGTAGCGACAATGCTGCCGCACACGCAAAACCGGCGCGCAACTTAACTTCGc           <  1:145458/149‑1 (MQ=255)
                                                                                                                                           ctgACGCTGGCATGCAGCGTCAGGCCCGCGCCAGCCATCAAGGCCTTGCGGAATACAGAAATAGACAAGTGGGGGTCCTCACAATCAGTGCCTGGGTAGCGACAATGCTGCCGCACACGCAAAACCGGCGCGCAACTTAACTTCGCAGc        >  4:298044/3‑149 (MQ=255)
                                                                                                                                             gACGCTGGCATGCAGCGTCAGGCCCGCGCCAGCCATCAAGGCCTTGCGGAATACAGAAATAGACAAGTGGGGGTCCTCACAATCAGTGCCTGGGTAGTGACAATGCTGCCGCACACGCAAAACCGGCGCGCAACTTAACTTCGCAGCtt      >  1:8109/1‑149 (MQ=255)
                                                                                                                                              aCGCTGCATGCCAGCGTCAGGCCCGCGCCAGCCATCAAGGCCTTGCGGAATACAGAAATAGACAAGTGGGGGTCCTCACAATCAGTGCCTGGGTAGCGACAATGCTGCCGCACACGCAAAACCGgc                            >  5:98992/1‑126 (MQ=255)
                                                                                                                                              aCGCTGCATGCCAGCGTCAGGCCCGCGCCAGCCATCAAGGCCTTGCGGAATACAGAAATAGACAAGTGGGGGTCCTCACAATCAGTGCCTGGGTAGCGACAATGCTGCCGCACACGCAAAACCGgc                            <  6:98992/126‑1 (MQ=255)
                                                                                                                                              aCGCTGCATGCCAGCGTCAGGCCCGCGCCAGCCATCAAGGCCTTGCGGAATACAGAAATAGACAAGTGGGGGTCCTCACAATCAGTGCCTGGGTAGCGACAATGCTGCCGCACACGCAAAACCGGCGCGCAACTTAACTTCGCAGCTTc     >  7:58204/1‑149 (MQ=255)
                                                                                                                                                gCTGCATGCCAGCGTCAGGCCCGCGCCAGCCATCAAGGCCTTGCGGAATACAGAAATAGACAAGTGGGGGTCCTCACAATCAGTGCCTGGGTAGCGACAATGCTGCCGCACACGCAAAACCGGCGCGCAACTTAACTTCGCAGCTTCGa   >  2:200180/1‑149 (MQ=255)
                                                                                                                                                gCTGCATGCCAGCGTCAGGCCCGCGCCAGCCATCAAGGCCTTGCGGAATACAGAAATAGACAAGTGGGGGTCCTCACAATCAGTGCCTGGGTAGCGACAATGCTGCCGCACACGCAAAACCGGCGCGCAACTTAACTTCGCAGCTTCGa   <  8:258488/149‑1 (MQ=255)
                                                                                                                                                gCTGCATGCCAGCGTCAGGCCCGCGCCAGCCATCAAGGCCTTGCGGAATACAGAAATAGACAAGTGGGGGTCCTCACAATCAGTGCCTGGGTAGCGACAATGCTGCCGCACACGCAAAACCGGCGCGCAACTTAACTTCGCAGCTTCGa   >  1:84576/1‑149 (MQ=255)
                                                                                                                                                 cTGCATGCCAGCGTCAGGCCCGCGCCAGCCATCAAGGCCTTGCGGAATACAGAAATAGACAAGTGGGGGTCCTCACAATCAGTGCCTGGGTAGCGACAATGCTGCCGCACACGCAAAACCGGCGCGCAACTTAACTTCGCAGCTTCGAt  <  1:200180/149‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                   
CAGCAGCTTGCCTTCCAGGGTTTCGTTGGAGTCGAACACATCTTGCACAGGCTTGATGCCCGTGGCTTTCTCGAAGTCCGCGAGTGTGGTGGGGCCGATGTAGTCGGACCAGTTATAGAAGTGCACCGTCGGCGCCGCTTGGACGCTGCATGCCAGCGTCAGGCCCGCGCCAGCCATCAAGGCCTTGCGGAATACAGAAATAGACAAGTGGGGGTCCTCACAATCAGTGCCTGGGTAGCGACAATGCTGCCGCACACGCAAAACCGGCGCGCAACTTAACTTCGCAGCTTCGAT  >  NC_002947/5907760‑5908053

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: