Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,909,295 A→T 9.4% intergenic (‑75/+94) spuD ← / ← spuC‑II spermidine/putrescine ABC transporter substrate‑binding protein/polyamine:pyruvate transaminase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,909,2950AT 9.4% 176.9 / 8.0 64intergenic (‑75/+94)spuD/spuC‑IIspermidine/putrescine ABC transporter substrate‑binding protein/polyamine:pyruvate transaminase
Reads supporting (aligned to +/- strand):  ref base A (27/31);  new base T (4/2);  total (31/33)
Fisher's exact test for biased strand distribution p-value = 4.19e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.28e-01

CTTCAGCCTGAACAGCGGAGGCCATGGCACCCATCAGGGCGGCGGCCAGCAACGTCTTGCCCATTTTATTCATGCGTAATGCTCCAGAATTTTTATTTAGCAACTGTCCAGCAGCCAGAACCCACGGTGCAGAGCGCGCGGCGACTGAAACAGTCGCTAGTCTGGCAAGTTACAAGGCCCTGTATCAAGGAAAGCAGGGCCATGTAACGACTTAATGTCACAACGCTAGCCTAGCAAACGCTTATTGAATCGCCTCAAGGGTCAGGTCCAGGCACTTGCGCGCCTTTTCCACCAGT  >  NC_002947/5909148‑5909443
                                                                                                                                                   |                                                                                                                                                    
cTTCAGCCTGAACAGCGGAGGCCATGGCACCCATCAGGGCGGCGGCCAGCAACGTCTTGCCCATTTTATTCATGCGTAATGCTCCAGAATTTTTATTTAGCAACTGTCCAGCAGCCAGAACCCACGGTGCAGAGCGCGCGGCGACTGa                                                                                                                                                      >  7:121795/1‑148 (MQ=255)
   cAGCCTGAACAGCGGAGGCCATGGCACCCATCAGGGCGGCGGCCAGCAACGTCTTGCCCATTTTATTCATGCGTAATGCTCCAGAATTTTTATTTAGCAACTGTCCAGCAGCCAGAACCCACGGTGCAGAGCGCGCGGCGACTGaaa                                                                                                                                                    <  5:92777/147‑1 (MQ=255)
    aGCCTGAACAGCGGAGGCCATGGCACCCATCAGGGCGGCGGCCAGCAACGTCTTGCCCATTTTATTCATGCGTAATGCTCCAGAATTTTTATTTAGCAACTGTCCAGCAACCAGAACCCACGGTGCAGAGCGCGCGGCGACTGAAACa                                                                                                                                                  <  6:328080/148‑1 (MQ=255)
       cTGAACAGCGGAGGCCATGGCACCCATCAGGGCGGCGGCCAGCAACGTCTTGCCCATTTTATTCATGCGTAATGCTCCAGAATTTTTATTTAGCAACTGTCCAGCAGCCAGAACCCACGGTGCAGAGCGCGCGGCGACTGTTTCAGTCg                                                                                                                                              >  8:238767/1‑149 (MQ=255)
          caCAGCGGAGGCCATGGCACCCATCAGGGCGGCGGCCAGCAACGTCTTGCCCATTTTATTCATGCGTAATGCTCCAGAATTTTTATTTAGCAACTGTCCAGCAGCCAGAACCCACGGTGCAGAGCGCGCGGCGACTGAAACAGTCGc                                                                                                                                             <  8:121795/146‑1 (MQ=255)
           aCAGCGGAGGCCATGGCACCCATCAGGGCGGCGGCCAGCAACGTCTTGCCCATTTTATTCATGCGTAATGCTCCAGAATTTTTATTTAGCAACTGTCCAGCAGCCAGAACCCACGGTGCAGAGCGCGCGGCGACTGTTTCAGTCGccgc                                                                                                                                          <  7:238767/149‑4 (MQ=255)
              gCGGAGGCCATGGCACCCATCAGGGCGGCGGCCAGCAACGTCTTGCCCATTTTATTCATGCGTAATGCTCCAGAATTTTTATTTAGCAACTGTCCAGCAGCCAGAACCCACGGTGCAGAGCGCGCGGCGACTGAAACAGTCGCTAGTCt                                                                                                                                       <  1:333944/149‑1 (MQ=255)
                   ggCCATGGCATGCATCAGGGCGGCGGTCAGCAAGGTCTTGGCCATTTTATTCATGCGTAATGCTCCAGTATTTTTATTTAGCTTCTGTCGATCAGGGTGATCCCAATGGGATGAGCGGGTGTCGACTGATTGAGTCGCGAGTCTgtgaa                                                                                                                                  <  6:69901/149‑5 (MQ=255)
                    gCCATGGCACCCATCAGGGCGGCGGCCAGCAACGTCTTGCCCATTTTATTCATGCGTAATGCTCCAGAATTTTTATTTAGCAACTGTCCAGCAGCCAGAACCCACGGTGCAGAGCGCGCGGCGACTGAAACAGTCGCTAGTCTGGCAAg                                                                                                                                 >  3:226548/1‑149 (MQ=255)
                     ccATGGCACCCATCAGGGCGGCGGCCAGCAACGTCTTGCCCATTTTATTCATGCGTAATGCTCCAGAATTTTTATTTAGCAACTGTCCAGCAGCCAGAACCCACGGTGCAGAGCGCGCGGCGACTGAAACAGTCGCTAGTCTGGCAAGt                                                                                                                                <  4:226548/149‑1 (MQ=255)
                      cATGGCACCCATCAGGGCGGCGGCCAGCAACGTCTTGCCCATTTTATTCATGCGTAATGCTCCAGAATTTTTATTTAGCAACTGTCCAGCAGCCAGAACCCACGGTGCAGAGCGCGCGGCGACTGAAACAGTCGCTAGTCTGGCAAGt                                                                                                                                >  5:111893/1‑148 (MQ=255)
                       aTGGCACCCATCAGGGCGGCGGCCAGCAACGTCTTGCCCATTTTATTCATGCGTAATGCTCAAGAATTTTTATTTAGCAACTGTCCAGCAGCCAGAACCCACGGTGCAGAGCGCGCGGCGACTGAAACAGTCGCTAGTCTGGCAAGTTa                                                                                                                              >  1:285822/1‑149 (MQ=255)
                            aCCCATCAGGGCGGCGGCCAGCAACGTCTTGCCCATTTTATTCATGCGTAATGCTCCAGAATTTTTATTTAGCAACTGTCCAGCAGCCAGAACCCACGGTGCAGAGCGCGCGGCGACTGa                                                                                                                                                      <  5:167863/120‑1 (MQ=255)
                            aCCCATCAGGGCGGCGGCCAGCAACGTCTTGCCCATTTTATTCATGCGTAATGCTCCAGAATTTTTATTTAGCAACTGTCCAGCAGCCAGAACCCACGGTGCAGAGCGCGCGGCGACTGa                                                                                                                                                      >  6:167863/1‑120 (MQ=255)
                              ccATCAGGGCGGCGGCCAGCAACGTCTTGCCCATTTTATTCATGCGTAATGCTCCAGAATTTTTATTTAGCAACTGTCCAGCAGCCAGAACCCACGGTGCAGAGCGCGCGGCGACTGAAACAGTCGCTAGTCTGGCAAGTTACAAGGcc                                                                                                                       >  2:128381/1‑149 (MQ=255)
                                    gggCGGCGGCCAGCAACGTCTTGCCCATTTTATTCATGCGTAATGCTCCAGAATTTTTATTTAGCAACTGTCCAGCAGCCAGAACCCACGGTGCAGAGCGCGCGGCGACTGAAACAGTCGCTAGTCTGGCAAGTTACAAGGCCTTGTAt                                                                                                                 >  4:63466/1‑149 (MQ=255)
                                        ggcggcCAGCAACGTCTTGCCCATTTTATTCATGCGTAATGCTCCAGAATTTTTATTTAGCAACTGTCCAGCAGCCAGAACCCACGGTGCAGAGCGCGCGGCGACTGAAACAGTCGCTAGTCTGGCAAGTTACAAGGCCCTGTATCAAg                                                                                                             <  8:114308/149‑1 (MQ=255)
                                                  aaCGTCTTGCCCATTTTATTCATGCGTAATGCTCCAGAATTTTTATTTAGCAACTGTCCAGCAGCCAGAACCCACGGTGCAGAGCGCGCGGCGACTGAAACAGTCGCTAGTCTGGCAAGTTACAAGGCCCTGTATCAAGGAAAGCAggg                                                                                                   <  6:181990/149‑1 (MQ=255)
                                                     gTCTTGCCCATTTTATTCATGCGTAATGCTCCAGAATTTTTATTTAGCAACTGTCCAGCAGCCAGAACCCACGGTGCAGAGCGCGCGGCGACTGAAACAGTCGCTAGTCTGGCAAGTTACAAGGCCCTGTATCAAGGAAAGCAGGGCCa                                                                                                >  5:41883/1‑149 (MQ=255)
                                                      tCTTGCCCATTTTATTCATGCGTAATGCTCCAGAATTTTTATTTAGCAACTGTCCAGCAGCCAGAACCCACGGTGCAGAGCGCGCGGCGACTGAAACAGTCGCTAGTCTGGCAAGTTACAAGGCCTTGTATCAAGGAAAGCAGGGCCAt                                                                                               <  3:63466/149‑1 (MQ=255)
                                                      tCTTGCCCATTTTATGCATGCGTACTGCTGCTGAATTTTTATTTAGCAACTGTCCAGCAGCCAGAACCCACGGTGCAGAGCGCGCGGCGACTGAAACAGTCGTTAGTCTGGCAAGTTACAAGGCCCTGTATCAAGGAAAGCAGGGCCAt                                                                                               <  6:307934/149‑1 (MQ=255)
                                                        ttGCCCATTTTATTCATGCGTAATGCTCCAGAATTTTTATTTAGCAACTGTCCAGCAGCCAGAACCCACGGTGCAGAGCGCGCGGCGACTGAAACAGTCGCTAGTCTGGCaa                                                                                                                                  >  3:200091/1‑112 (MQ=255)
                                                        ttGCCCATTTTATTCATGCGTAATGCTCCAGAATTTTTATTTAGCAACTGTCCAGCAGCCAGAACCCACGGTGCAGAGCGCGCGGCGACTGAAACAGTCGCTAGTCTGGCaa                                                                                                                                  <  4:200091/112‑1 (MQ=255)
                                                        ttGCCCATTTTATTCATGCGTAATGCTCCAGAATTTTTATTTAGCAACTGTCCAGCAGCCAGAACCCACGGTGCAGAGCGCGCGGCGACTGAAACAGTCGCTAGTCTGGCAAGTTACAAGGCCCTGTATCaa                                                                                                              <  7:322309/132‑1 (MQ=255)
                                                        ttGCCCATTTTATTCATGCGTAATGCTCCAGAATTTTTATTTAGCAACTGTCCAGCAGCCAGAACCCACGGTGCAGAGCGCGCGGCGACTGAAACAGTCGCTAGTCTGGCAAGTTACAAGGCCCTGTATCaa                                                                                                              >  8:322309/1‑132 (MQ=255)
                                                          gCCCATTTTATTCATGCGTAATGCTCCAGAATTTTTATTTAGCAACTGTCCAGCAGCCAGAACCCACGGTGCAGAGCGCGCGGCGACTGAAACAGTCGCTAGTCTGGCAAGTTACAAGGCCCTGTATCAAGGAAAGCAGGGCCATGTaa                                                                                           >  2:268618/1‑149 (MQ=255)
                                                          gCCCATTTTATTCATGCGTAATGCTCCAGAATTTTTATTTAGCAACTGTCCAGCAGCCAGAACCCACGGTGCAGAGCGCGCGGCGACTGAAACAGTCGCTAGTCTGGCAAGTTACAAGGCCCTGTATCAAGGAAAGCAGGGCCATGTa                                                                                            <  6:41883/148‑1 (MQ=255)
                                                             cATTTTATTCATGCGTAATGCTCCAGAATTTTTATTTAGCAACTGTCCAGCAGCCAGAACCCACGGTGCAGAGCGCGCGGCGACTGAAACAGt                                                                                                                                                <  6:232430/93‑1 (MQ=255)
                                                             cATTTTATTCATGCGTAATGCTCCAGAATTTTTATTTAGCAACTGTCCAGCAGCCAGAACCCACGGTGCAGAGCGCGCGGCGACTGAAACAGt                                                                                                                                                >  5:232430/1‑93 (MQ=255)
                                                             cATTTTATTCATGCGTAATGCTCCAGAATTTTTATTTAGCAACTGTCCAGCAGCCAGAACCCACGGTGCAGAGCGCGCGGCGACTGAAACAGTCGCTAGTCTGGCAAGTTACAAGGCCCTGTATCAAGGAAAGCAGGGCCATGTAACGa                                                                                        >  4:146201/1‑149 (MQ=255)
                                                              aTTTTATTCATGCGTAATGCTCCAGAATTTTTATTTAGCAACTGTCCAGCAGCCAGAACCCACGGTGCAGAGCGCGCGGCGACTGAAACAGTCGCTAGTCTGGCAAGTTACAAGGCCCTGTATCAAGGAAAGCAGGGCCATGTAAc                                                                                          >  7:177001/1‑146 (MQ=255)
                                                              aTTTTATTCATGCGTAATGCTCCAGAATTTTTATTTAGCAACTGTCCAGCAGCCAGAACCCACGGTGCAGAGCGCGCGGCGACTGAAACAGTCGCTAGTCTGGCAAGTTACAAGGCCCTGTATCAAGGAAAGCAGGGCCATGTAAc                                                                                          <  8:177001/146‑1 (MQ=255)
                                                                   aTTCATGCGTAATGCTCCAGAATTTTTATTTAGCAACTGTCCAGCAGCCAGAACCCACGGTGCAGAGCGCGCGGCGACTGAAACAGTCGCTAGTCTGGCAAGTTACAAGGCCCTGTATCAAGGAAAGCAGGGCCATGTAACGACTTaa                                                                                   <  1:268618/148‑1 (MQ=255)
                                                                       aTGCGTAATGCTCCAGAATTTTTATTTAGCAACTGTCCAGCAGCCAGAACCCACGGTGCAGAGCGCGCGGCGACTGAAACAGTCGCTAGTCTGGCAAGt                                                                                                                                >  8:80900/1‑99 (MQ=255)
                                                                       aTGCGTAATGCTCCAGAATTTTTATTTAGCAACTGTCCAGCAGCCAGAACCCACGGTGCAGAGCGCGCGGCGACTGAAACAGTCGCTAGTCTGGCAAGt                                                                                                                                <  7:80900/99‑1 (MQ=255)
                                                                       aTGCGTAATGCTCCAGAATTTTTATTTAGCAACTGTCCAGCAGCCAGAACCCACGGTGCAGAGCGCGCGGCGACTGAAACAGTCGCTAGTCTGGCAAGTTACAAGGCCCTGTATCAAGGAAAGCAGGGCCATGTAACGACTTaa                                                                                   <  1:138246/144‑1 (MQ=255)
                                                                       aTGCGTAATGCTCCAGAATTTTTATTTAGCAACTGTCCAGCAGCCAGAACCCACGGTGCAGAGCGCGCGGCGACTGAAACAGTCGCTAGTCTGGCAAGTTACAAGGCCCTGTATCAAGGAAAGCAGGGCCATGTAACGACTTaa                                                                                   >  2:138246/1‑144 (MQ=255)
                                                                               tGCTCCAGAATTTTTATTTAGCAACTGTCCAGCAGCCAGAACCCACGGTGCAGAGCGCGCGGCGACTGAAACAGTCGCTAGTCTGGCAAGTTACAAGGCCCTGTATCAAGGAAAGCAGGGCCATGTAACGACTTAATGTCACAACGc                                                                        >  3:64393/1‑147 (MQ=255)
                                                                                  tCCAGAATTTTTATTTAGCAACTGTCCAGCAGCCAGAACCCACGGTGCAGAGCGCGCGGCGACTGAAACAGTCGCTAGTCTGGCAAGTTACAAGGCCCTGTATCAAGGAAAGCAGGGCCATGTAAc                                                                                          >  5:104148/1‑126 (MQ=255)
                                                                                  tCCAGAATTTTTATTTAGCAACTGTCCAGCAGCCAGAACCCACGGTGCAGAGCGCGCGGCGACTGAAACAGTCGCTAGTCTGGCAAGTTACAAGGCCCTGTATCAAGGAAAGCAGGGCCATGTAAc                                                                                          <  6:104148/126‑1 (MQ=255)
                                                                                  tCCAGAATTTTTATTTAGCAACTGTCCAGCAGCCAGAACCCACGGTGCAGAGCGCGCGGCGACTGAAACAGTCGCTAGTCTGGCAAGTTACAAGGCCCTGTATCAAGGAAAGCAGGGCCATGTAACGACTTAATGTCACAACGctagcc                                                                   >  5:129451/1‑149 (MQ=255)
                                                                                  tCCAGAATTTTTATTTAGCAACTGTCCAGCAGCCAGAACCCACGGTGCAGAGCGCGCGGCGACTGAAACAGTCGCTAGTCTGGCAAGTTACAAGGCCCTGTATCAAGGAAAGCAGGGCCATGTAACGACTTAATGTCACAACGctagcc                                                                   >  2:204213/1‑149 (MQ=255)
                                                                                    cAGAATTTTTATTTAGCAACTGTCCAGCAGCCAGAACCCACGGTGCAGAGCGCGCGGCGACTGAAACAGTCGCTAGTCTGGCAAGTTACAAGGCCCTGTATCAAGGAAAGCAGGGCCATGTAACGACTTAATGTCACAACGctagcc                                                                   >  4:224769/1‑147 (MQ=255)
                                                                                          ttttATTTAGCAACTGTCCAGCAGCCAGAACCCACGGTGCAGAGCGCGCGGCGACTGAAACAGTCGCTAGTCTGGCAAGTTACAAGGCCCTGTATCAAGGAAAGCAGGGCCATGTAACGACTTAATGTCACAACGCTAGCCTAGCAAAc                                                           >  6:282352/1‑149 (MQ=255)
                                                                                                  aGCAACTGTCCAGCAGCCAGAACCCACGGTGCAGAGCGCGCGGCGACTGAAACAGTCGCTAGTCTGGCAAGTTACAAGGCCCGGTATCAAGGAAAGCAGGGCCATGTAACGACTTAATGTCACAACGCTAGCCTAGCAAACGCTTAtt                                                    >  6:260033/1‑148 (MQ=255)
                                                                                                     aaCTGTCCAGCAGCCAGAACCCACGGTGCAGAGCGCGCGGCGACTGAAACAGTCGCTAGTCTGGCAAGTTACAAGGCCCTGTATCAAGGAAAGCAGGGCCATGTAACGACTTAATGTCACAACGCTAGCCTAGCAAACGCTTATTGAAt                                                <  4:64393/149‑1 (MQ=255)
                                                                                                               cagcCAGAACCCACGGTGCAGAGCGCGCGGCGACTGAAACAGTCGCTAGTCTGGCAAGTTACAAGGCCCTGTATCAAGGAAAGCAGGGCCATGTAACGACTTAATGTCACAACGCTAGCCTAGCAAACGCTTATTGAATCGCCTCAAgg                                      <  1:204213/149‑1 (MQ=255)
                                                                                                                agcCAGAACCCACGGTGCAGAGCGCGCGGCGACTGAAACAGTCGCTAGTCTGGCAAGTTACAAGGCCCTGTATCAAGGAAAGCAGGGCCATGTAACGACTTAATGTCACAACGCTAGCCTAGCAAACGCTTATTGAATCGCCTCAAggg                                     >  1:131685/1‑149 (MQ=255)
                                                                                                                   cAGAACCCACGGTGCAGAGCGCGCGGCGACTGAAACAGTCGCTAGTCTGGCAAGTTACAAGGCCCTGTATCAAGGAAAGCAGGGCCATGTAACGACTTAATGTCACAACGCTAGCCTAGCAAACGCTTAtt                                                    >  3:69642/1‑131 (MQ=255)
                                                                                                                   cAGAACCCACGGTGCAGAGCGCGCGGCGACTGAAACAGTCGCTAGTCTGGCAAGTTACAAGGCCCTGTATCAAGGAAAGCAGGGCCATGTAACGACTTAATGTCACAACGCTAGCCTAGCAAACGCTTAtt                                                    <  4:69642/131‑1 (MQ=255)
                                                                                                                   cAGAACCCACGGTGCAGAGCGCGCGGCGACTGAAACAGTCGCTAGTCTGGCAAGTTACAAGGCCCTGTATCAAGGAAAGCAGGGCCATGTAACGACTTAATGTCACAACGCTAGCCTAGCAAACGCTTATTGAATCGCCTCAAGGGTc                                   <  1:128381/148‑1 (MQ=255)
                                                                                                                     gAACCCACGGTGCAGAGCGCGCGGCGACTGAAACAGTCGCTAGTCTGGCAAGTTACAAGGCCCTGTATCAAGGAAAGCAGGGCCATGTAACGACTTAATGTCACAACGCTAGCCTAGCAAACGCTTATTGAATc                                               <  1:323068/134‑1 (MQ=255)
                                                                                                                     gAACCCACGGTGCAGAGCGCGCGGCGACTGAAACAGTCGCTAGTCTGGCAAGTTACAAGGCCCTGTATCAAGGAAAGCAGGGCCATGTAACGACTTAATGTCACAACGCTAGCCTAGCAAACGCTTATTGAATc                                               >  2:323068/1‑134 (MQ=255)
                                                                                                                       aCCCACGGTGCAGAGCGCGCGGCGACTGAAACAGTCGCTAGTCTGGCAAGTTACAAGGCCCTGTATCAAGGAAAGCAGGGCCATGTAACGACTTAATGTCACAACGCTAGCCTAGCAAACGCTTATTGAATCGCCTCAAGGGTCAGGTc                              <  6:129451/149‑1 (MQ=255)
                                                                                                                          cACGGTGCAGAGCGCGCGGCGACTGAAACAGTCGCTAGTCTGGCAAGTTACAAGGCCCTGTATCAAGGAAAGCAGGGCCATGTAACGACTTAATGTCACAACGCTAGCCTAGCAAACGCTTATTGAATCGCCTCAAGGGTCAGGTCCAg                           <  2:285822/149‑1 (MQ=255)
                                                                                                                                  agagCGCGCGGCGACTGAAACAGTCGCTAGTCTGGCAAGTTACAAGGCCCTGTATCAAGGAAAGCAGGGCCATGTAACGACTTAATGTCACAACGCTAGCCTAGCAAACGCTTATTGAATCGCCTCAAGGGTCAGGTCCAGGCACTTgc                   <  3:146201/149‑1 (MQ=255)
                                                                                                                                   gagCGCGCGGCGACTGAAACAGTCGCTAGTCTGGCAAGTTACAAGGCCCTGTATCAAGGAAAGCAGGGCCATGTAACGACTTAATGTCACAACGCTAGCCTAGCAAACGCTTATTGAATCGCCTCAAGGGTCAGGt                               <  8:21853/136‑1 (MQ=255)
                                                                                                                                   gagCGCGCGGCGACTGAAACAGTCGCTAGTCTGGCAAGTTACAAGGCCCTGTATCAAGGAAAGCAGGGCCATGTAACGACTTAATGTCACAACGCTAGCCTAGCAAACGCTTATTGAATCGCCTCAAGGGTCAGGt                               >  7:21853/1‑136 (MQ=255)
                                                                                                                                    agCGCGCGGCGACTGAAACAGTCGCTAGTCTGGCAAGTTACAAGGCCCTGTATCAAGGAAAGCAGGGCCATGTAACGACTTAATGTCACAACGCTAGCCTAGCAAACGCTTATTGAATCGCCTCAAGGGTCAGGTCCAGGCACTTgcgc                 <  3:224769/149‑1 (MQ=255)
                                                                                                                                      cgcgcgGCGACTGAAACAGTCGCTAGTCTGGCAAGTTACAAGGCCCTGTATCAAGGAAAGCAGGGCCATGTAACGACTTAATGTCACAACGCTAGCCTAGCAAACGCTTATTGAATCGCCTCAAGGGTCAGGTCCAGGCACTTgcgcgc               >  4:168162/1‑149 (MQ=255)
                                                                                                                                        cgcgGCGACTGAAACAGTCGCTAGTCTGGCAAGTTACAAGGCCCTGTATCAAGGAAAGCAGGGCCATGTAACGACTTAATGTCACAACGCTAGCCTAGCAAACGCTTATTGAATCGCCTCAAGGGTCAGGTCCAGGCACTTGCGCGCCt             <  2:131685/149‑1 (MQ=255)
                                                                                                                                          tagcgACTGTTTCAGTCGCTAGTCTGGCAAGTTACAAGGCCCTGTATCAAGGAAAGCAGGGCCATGTAACGACTTAATGTCACAACGCTAGCCTAGCAAACGCTTATTGa                                                  <  6:124040/108‑1 (MQ=255)
                                                                                                                                          tagcgACTGTTTCAGTCGCTAGTCTGGCAAGTTACAAGGCCCTGTATCAAGGAAAGCAGGGCCATGTAACGACTTAATGTCACAACGCTAGCCTAGCAAACGCTTATTGa                                                  >  5:124040/3‑110 (MQ=255)
                                                                                                                                          tagcgACTGTTTCAGTCGCTAGTCTGGCAAGTTACAAGGCCCTGTATCAAGGAAAGCAGGGCCATGTAACGACTTAATGTCACAACGCTAGCCTAGCAAACGCTTATTGAATCGCCTCAAGGGTCAGGTCCAGGCACTTGCGCGCCttt           >  1:126716/3‑149 (MQ=255)
                                                                                                                                          cggcgACTGAAACAGTCGCTAGTCTGGCAAGTTACAAGGCCCTGTATCAAGGAAAGCAGGGCCATGTAACGACTTAATGTCACAACGCTAGCCTAGCAAACGCTTATTGAATCGCCTCAAGGGTCAGGTCCAGGCACTTGCGCGCCttt           <  6:111893/149‑1 (MQ=255)
                                                                                                                                            gcgACTGTTTCAGTCGCTAGTCTGGCAAGTTACAAGGCCCTGTATCAAGGAAAGCAGGGCCATGTAACGACTTAATGTCACAACGCTAGCCTAGCAAACGCTTATTGAATCGCCTCAAGGGTCAGGTCCAGGCACTTGCGCGCCTTTTc         >  4:174578/1‑149 (MQ=255)
                                                                                                                                            gcgACTGAAACAGTCGCTAGTCTGGCAAGTTACAAGGCCCTGTATCAAGGAAAGCAGGGCCATGTAACGACTTAATGTCACAACGCTAGCCTAGCAAACGCTTATTGAATCGCCTCAAGGGTCAGGTCCAGGCACTTGCGCGCCTTTTc         <  5:282352/149‑1 (MQ=255)
                                                                                                                                             cgACTGAAACAGTCGCTAGTCTGGCAAGTTACAAGGCCCTGTATCAAGGAAAGCAGGGCCATGTAACGACTTAATGTCACAACGCTAGCCTAGCAAACGCTTATTGAATCGCCTCAAGGGTCAGGTCCAGGCACTTGCGCGCCTTTTc         >  4:8921/1‑148 (MQ=255)
                                                                                                                                                   aaaCAGTCGCTAGTCTGGCAAGTTACAAGGCCCTGTATCAAGGAAAGCAGGGCCATGTAACGACTTAATGTCACAACGCTAGCCTAGCAAACGCTTATTGAATCGCCTCAAGGGTCAGGTCCAGGCACTTGCGCGCCTTTTCCACCAGt  >  8:161005/1‑149 (MQ=255)
                                                                                                                                                   |                                                                                                                                                    
CTTCAGCCTGAACAGCGGAGGCCATGGCACCCATCAGGGCGGCGGCCAGCAACGTCTTGCCCATTTTATTCATGCGTAATGCTCCAGAATTTTTATTTAGCAACTGTCCAGCAGCCAGAACCCACGGTGCAGAGCGCGCGGCGACTGAAACAGTCGCTAGTCTGGCAAGTTACAAGGCCCTGTATCAAGGAAAGCAGGGCCATGTAACGACTTAATGTCACAACGCTAGCCTAGCAAACGCTTATTGAATCGCCTCAAGGGTCAGGTCCAGGCACTTGCGCGCCTTTTCCACCAGT  >  NC_002947/5909148‑5909443

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: ‑