Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,549,606 T→A 7.5% intergenic (+19/‑153) PP_t47 → / → PP_4881 tRNA‑Leu/iron ABC transporter substrate‑binding protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,549,6060TA 7.5% 226.5 / 6.4 80intergenic (+19/‑153)PP_t47/PP_4881tRNA‑Leu/iron ABC transporter substrate‑binding protein
Reads supporting (aligned to +/- strand):  ref base T (38/36);  new base A (4/2);  total (42/38)
Fisher's exact test for biased strand distribution p-value = 6.78e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.21e-01

AAGCGAACGATAGTCGCCTTCTGCAGATGATGTAAAGCTTTGCCCAGGTGGCGGAATTGGTAGACGCGCTGGTTTCAGGTATCAGTGACCGCAAGGTCGTGGAAGTTCGAGTCTTCTCCTGGGCACCATACAAAAACCCACTAGCTTGCTAGTGGGTTTTTTCTTGCCTGCGATTTATGTGCCCCCTCCCCCCCTTGCACCATGGTCTGACCATGAATTTCTCGTCACACGGCCACTTGAGAAACAATTTCGTTTACCATTGTTTCCAAATATGTAGCCGTAAGCGAGGAAGT  >  NC_002947/5549460‑5549752
                                                                                                                                                  |                                                                                                                                                  
aaGCGAACGATCGTCGCCTTCTGCAGTTGATGTAAATCTTTGCCCAGGTGGCGGAATTGGTAGTCGGGCTTGTTTAAGGTATCAGTGACCGCAAGGTCGTGGAAGTTCTAGTCTTATCCTGGTCACCATACAAAAACCCACTAGCTTGc                                                                                                                                                  <  4:38937/149‑1 (MQ=21)
aaGCGAACGATAGTCGCCTTCTGCAGATGATGTAAAGCTTTGCCCAGGCGGCGGAATTGGTAGACGCGCTGGTTTCAGGTATCAGTGACCGCAAGGTCGTGGAAGTTCGAGTCTTCTCCTGGGCACCATACAAAAACCCACTAGCTTGc                                                                                                                                                  <  6:1647/149‑1 (MQ=255)
 aGCGAACGATAGTCGCCTTCTGCAGATGATGTAAAGCTTTGCCCAGGTGGCGGAATTGGTAGACGCGCTGGTTTCAGGTATCAGTGACCGCAAGGTCGTGGAAGTTCGAGTCTTCTCCTGGGCACCATACAAAAACCCTCTAGCTTGc                                                                                                                                                  >  3:162583/1‑148 (MQ=255)
   cGAACGATAGTCGCCTTCTGCAGATGATGTAAAGCTTTGCCCAGGTGGCGGAATTGGTAGACGCGCTGGTTTCAGGTATCAGTGACCGCAAGGTCGTGGAAGTTCGAGTCTTCTCCTGGGCACCATACAAAAACCCACTAGCTTGCTAg                                                                                                                                               >  4:286452/1‑149 (MQ=255)
     aaCGATAGTCGCCTTCTGCAGATGATGTAAAGCTTTGCCCAGGTGGTGGAATTGGTAGACGCGCTGGTTTCAGGTATCAGTGACCGCAAGGTCGTGGAAGTTCGAGTCTTCTCCTGGGCACCATACAAAAACCCACTAGCTTGCTAGTg                                                                                                                                             >  6:267030/1‑149 (MQ=255)
     aaCGATAGTCGCCTTCTGCAGATGATGTAAAGCTTTGCCCAGGTGGCGGAATTGGTAGACGCGCTGGTTTCAGGTATCAGTGACCGCAAGGTCGTGGAAGTTCGAGTCTTCTCCTGGGCACCATACAAAAACCCACTAGCAAGCTAGTg                                                                                                                                             <  1:153510/149‑1 (MQ=255)
      aCGATAGTCGCCTTCTGCAGATGATGTAAAGCTTTTCCCAGGTGGCGGAATTGGTAGACGATCTGTTTTCAGTTATCAGTTACCGCAATGTCGTGGAAGTTCGAGTCTTCTCCTGGGCACCATACAAAAACCCACTAGCTTGCTAGTgg                                                                                                                                            <  2:242863/149‑1 (MQ=25)
      aCGATAGTCGCCTTCTGCAGATGATGTAAAGCTTTGCCCAGGTGGCGGAATTGGTAGACGCGCTGGTTTCAGGTATCAGTGACCGCAAGGTCGTGGAAGTTCGAGTCTTCTCCTGGGCACCATACAAAAACCCACTAGCAAGCTAGTgg                                                                                                                                            <  8:150396/149‑1 (MQ=255)
             tCGCCTTCTGCAGATGATGTAAAGCTTTGCCCAGGTGGCGGAATTGGTAGACGCGCTGGTTTCAGGTATCAGTGACCGCAAGGTCGTGGAAGTTCGAGTCTTCTCCTGGGCACCATACAAAAACCCACTAGCTTGCTAGTGGGtttttt                                                                                                                                     >  1:351120/1‑149 (MQ=255)
                 cTTCTGCAGATGATGTAAAGCTTTGCCCAGGTGGCGGAATTGGTAGACGCGCTGGTTTCAGGTATCAGTGACCGCAAGGTCGTGGAAGTTCGAGTCTTCTCCTGGGCACCATACAAAAACCCACTAGCTTGCTAGTGGGTTTTTTCTTg                                                                                                                                 <  2:355770/149‑1 (MQ=255)
                     tGCAGATGATGTAAAGCTTTGCCCAGGTGGCGGAATTGGTAGACGCGCTGGTTTCAGGTATCAGTGACCGCAAGGTCGTGGAAGTTCGAGTCTTCTCCTGGGCACCATACAAAAACCCACTAGCTTGCTAGTGGGTTTTTTCTTGCCTg                                                                                                                             >  3:294488/1‑149 (MQ=255)
                      gCAGATGATGTAAAGCTTTGCCCAGGTGGCGGAATTGGTAGACGCGCTGGTTTCAGGTATCAGTGACCGCAAGGTCGTGGAAGTTCGAGTCTTCTCCTGGGCACCATACAAAAACCCACTAGCTTGCTAGTGGGTTTTTTCTTGCCTGc                                                                                                                            >  3:157408/1‑149 (MQ=255)
                       cAGATGATGTAAAGCTTTGCCCAGGTGGCGGAATTGGTAGACGCGCTGGTTTCAGGTATCAGTGACCGCAAGGTCGTGGAAGTTCGAGTCTTCTCCTGGGCACCATACAAAAACCCACTAGCTTGCTAGTGGGTTTTTTCTTGCCTGc                                                                                                                            >  8:92501/1‑148 (MQ=255)
                         gatgatGTAAAGCTTTGCCCAGGTGGTGGAATTGGTAGACGCGCTGGTTTCAGGTATCAGTGACCGCAAGGTCGTGGAAGTTCGAGTCTTCTCCTGGGCACCATACAAAAACCCACTAGCTTGCTAGTGGGTTTTTTCTTGCCTGCGAt                                                                                                                         <  5:267030/149‑1 (MQ=255)
                         gatgatGTAAAGCTTTGCCCAGGTGGCGGAATTGGTAGACGCGCTGGTTTCAGGTATCAGTGACCGCAAGGTCGTGGAAGTTCGAGTCTTCTCCTGGGCACCATACAAAAACCCTCTAGCTTGCTAGTGGGTTTTTTCTTGCCTGCGAt                                                                                                                         <  4:162583/149‑1 (MQ=255)
                                tAAAGCTTTGCCCAGGTGGCGGAATTGGTAGACGCGCTGGTTTCAGGTATCAGTGACCGCAAGGTCGTGGAAGTTCGAGTCTTCTCCTGGGCACCATACAAAAACCCACTAGCTTGCTAGTGGGTTTTTTCTTGCCTGCGATTTAtgtg                                                                                                                  >  3:122333/1‑149 (MQ=255)
                                tAAAGCTTTGCCCAGGTGGCGGAATTGGTAGACGCGCTGGTTTCAGGTATCAGTGACCGCAAGGTCGTGGAAGTTCGAGTCTTCTCCTGGGCACCATACAAAAACCCACTAGCTTGCTAGTGGGTTTTTTCTTGCCTGCGATTTAtgt                                                                                                                   <  5:218280/148‑1 (MQ=255)
                                tAAAGCTTTGCCCAGGTGGCGGAATTGGTAGACGCGCTGGTTTCAGGTATCAGTGACCGCAAGGTCGTGGAAGTTCGAGTCTTCTCCTGGGCACCATACAAAAACCCACTAGCTTGCTAGTGGGTTTTTTCTTGCCTGCGATTTAt                                                                                                                     <  6:248545/146‑1 (MQ=255)
                                tAAAGCTTTGCCCAGGTGGCGGAATTGGTAGACGCGCTGGTTTCAGGTATCAGTGACCGCAAGGTCGTGGAAGTTCGAGTCTTCTCCTGGGCACCATACAAAAACCCACTAGCTTGCTAGTGGGTTTTTTCTTGCCTGCGATTTAt                                                                                                                     >  5:248545/1‑146 (MQ=255)
                                  aaGCTTTGCCCAGGTGGCGGAATTGGTAGACGCGCTGGTTTCAGGTATCAGTGACCGCAAGGTCGTGGAAGTTCGAGTCTTCTCCTGGGCGCCATACAAAAACCCACTAGCTTGCTAGTGGGTTTTTTCTTGCCTGCGATTTAtgtgag                                                                                                                >  7:115960/1‑147 (MQ=255)
                                  aaGCTTTGCCCAGGTGGCGGAATTGGTAGACGCGCTGGTTTCAGGTATCAGTGACCGCAAGGTCGTGGAAGTTCGAGTCTTCTCCTGGGCACCATACAAAAACCCACTAGCTTGCTAGTGGGTTTTTTCTTGCCTGCGATTTATGTGcc                                                                                                                <  4:135297/149‑1 (MQ=255)
                                     gtgtGCCCAGGTGGCGGAATTGGTAGACGCGCTGGTTTCAGGTATCAGTGACCGCAAGGTCGTGGAAGTTCGAGTCTTCTCTTGGGCACCATACAAAAACCCACTAGCTTGCTa                                                                                                                                                <  1:272532/111‑1 (MQ=255)
                                     gtgtGCCCAGGTGGCGGAATTGGTAGACGCGCTGGTTTCAGGTATCAGTGACCGCAAGGTCGTGGAAGTTCGAGTCTTCTCCTGGGCACCATACAAAAACCCACTAGCTTGCTa                                                                                                                                                >  2:272532/4‑114 (MQ=255)
                                      tttGCCCAGGTGGCGGAATTGGTAGACGCGCTGGTTTCAGGTATCAGTGACCGCAAGGTCGTGGAAGTTCGAGTCTTCTCCTGGGCACCATACAAAAACCCACTAGCTTGCTAGTGGGTTTTTTCTTGCCTGCGATTTATGTGCCCCCt                                                                                                            >  3:356405/1‑149 (MQ=255)
                                       ttGCCCAGGTGGCGGAATTGGTAGACGCGCTGGTTTCAGGTATCAGTGACCGCAAGGTCGTGGAAGTTCGAGTCTTCTCCTGGGCACCATACAAAAACCCACTAGCTTGCTAGTGGGTTTTTTCTTg                                                                                                                                 <  4:67268/127‑1 (MQ=255)
                                       ttGCCCAGGTGGCGGAATTGGTAGACGCGCTGGTTTCAGGTATCAGTGACCGCAAGGTCGTGGAAGTTCGAGTCTTCTCCTGGGCACCATACAAAAACCCACTAGCTTGCTAGTGGGTTTTTTCTTg                                                                                                                                 >  3:67268/1‑127 (MQ=255)
                                          cccAGGTGGCGGAATTGGTAGACGCGCTGGTTTCAGGTATCAGTGACCGCAAGGTCGTGGAAGTTCGAGTCTTCTCCTGGGCACCATACAAAAACCCACTAGCTTGCTAGTGGGTTTTTTCTTGCCTGCGATTTATGTGCGCGCTcccc                                                                                                        >  5:42570/1‑149 (MQ=255)
                                            cAGGTGGCGGAATTGGTAGACGCGCTGGTTTCAGGTATCAGTGACCGCAAGGTCGTGGAAGTTCGAGTCTTCTCCTGGGCACCATACAAAAACCCACTAGCTTGCTAGTGGGTTTTTTCTTGCCTGCGATTTATGTGCCCCCTcccccc                                                                                                      <  2:351120/149‑1 (MQ=255)
                                                 ggCGGAATTGGTAGACGCGCTGGTTTCAGGTATCAGTGACCGCAAGGTCGTGGAAGTTCGAGTCTTCTCCTGGGCACCATACAAAAACCCACTAGCTTGCTAGTGGGTTTTATCTTGCCTGCGATTTATGTGCCCCCTCACCCCCTTg                                                                                                  >  3:95663/1‑148 (MQ=255)
                                                    ggAATTGGTAGACGCGCTGGTTTCAGGTATCAGTGACCGCAAGGTCGTGGAAGTTCGAGTCTTCTCCTGGGCACCATACAAAAACCCACTAGCTTGCTAGTGGGTTTTTTCTTGCCTGCGATTTATGTGGCCGCTCCCCCCCTGGCAcc                                                                                              >  6:316733/1‑149 (MQ=255)
                                                           gTAGACGCGCTGGTTTCAGGTATCAGTGACCGCAAGGTCGTGGAAGTTCGAGTCTTCACCTGGGCACCATACAAAAACCCACTAGCTTGCTAGTGGGTTTTTTCTTGCCTGCGATTTATGTGGCACCGAGCCCCCGTGCACCATGGTCt                                                                                       >  7:384164/1‑149 (MQ=25)
                                                              gACGCGCTGGTTTCAGGTATCAGTGACCGCAAGGTCGTGGAAGTTCGAGTCTTCTCCTGGGCACCATACAAAAACCCACTAGACAGCTAGtcg                                                                                                                                            >  8:186014/1‑91 (MQ=17)
                                                               tcgcgCTGGTTTCACGTATCAGTGAGCGCAAGGTCGTGGAAGTTCGATTCTTCTCCTGGGCACCATACAAAAACCCACTAGCTTGCTAGTGGTTTTTTTCTTGCCTGCGATTTATGTGCCCCCTCCCCCCCTTTCACCATGGTCTGAc                                                                                    <  8:305433/147‑1 (MQ=255)
                                                                 gcgcTGGTTTCAGGTATCAGTGACCGCAAGGTCGTGGAAGTTCGAGTCTTCTCCTGGGCACCATACAAAAACCCACTAGCTTGCTAGTGGGTTTTTTCTTGCCTGCGATTTATGTGGCCCCGCCCGGCCGTGCGCCCTGGGCTGACCAt                                                                                 >  7:63685/1‑149 (MQ=25)
                                                                 gcgcTGGTTTCAGGTATCAGTGACCGCAAGGTCGTGGAAGTTCGAGTCTTCTCCTGCGCACCATACAAAAACCCACTAGCTTGCTAGTGGGTTTTTTCTTGCCTGCGATTTATGTGCCCCCTCCCCCCCTTGCACCATGGTCTGACCAt                                                                                 <  7:184762/149‑1 (MQ=255)
                                                                  cgcTGGTTTCAGGTATCAGTGACCGCAAGGTCGTGGAAGTTCGAGTCTTCTCCTGGGCACCATACAAAAACCCACTAGCTTGCTAGTGGGTTTTTTCTTGCCTGCGATTTATGTGCCCCCTCCCCCCCTTGCACCATGGTCTGACCATg                                                                                <  1:122362/149‑1 (MQ=255)
                                                                   gcTGGTTTCAGGTATCAGTGACCGCAAGGTCGTGGAAGTTCGAGTCTTCTCCTGGGCACCATACAAAAACCCACTAGCTTGCTAGTGGGTTTTTTCTTGCCGGCGATTTATGTGCCCCCGCCCCCCCGTGCACCATGGTCGGACCATg                                                                                >  4:59107/1‑148 (MQ=255)
                                                                       gTTTCAGGTATCAGTGACCGCAAGGTCGTGGAAGTTCGAGTCTTCTCCTGGGCACCATACAAAAACCCACTAGCTTGCTAGTGGGTTTTTTCTTGCCTGCGATTTATGTGCCCCCTCCCCCCCTTGCACCATGGTCTGACCATGAAt                                                                             <  4:122333/147‑1 (MQ=255)
                                                                           cAGGTATCAGTGACCGCAAGGTCGTGGAAGTTCGAGTCTTCTCCTGGGCACCATACAAAAACCCACTAGCTTGCTAGTGGGTTTTTTCTTGCCTGCGATTTATGTGCCCCCTCCCCCCCTTGCACCATGGTCTGACCATGAATTtctc                                                                        <  3:286452/148‑1 (MQ=255)
                                                                            aGGTATCAGTGACCGCAAGGTCGTGGAAGTTCGAGTCTTCTCCTGGGCACCATACAAAAACCCACTAGCTTGCTAGTGGGTTTTTTCTTGCCTGCGATTTATGTGCCCCCTCCCCCCCTTGCACCATGGTCTGACCATGAATTTCTCGt                                                                      <  5:316733/149‑1 (MQ=255)
                                                                              gTATCAGTGACCGCAAGGTCGTGGAAGTTCGAGTCTTCTCCTGGGCACCATACAAAAACCCACTAGCTTGCTAGTGGGTTTTTTCTTGCCTGCGATTTATGTGGCCCCTCCCACCCTTGCCCCATGGTCTGACCATGAATTTCTCGTc                                                                     >  1:56450/1‑148 (MQ=255)
                                                                              gTATCAGTGACCGCAAGGTCGTGGAAGTTCGAGTCTTCTCCTGGGCACCATACAAAAACCCACTAGCTTGCTAGTGGGTTTTTTCTTGCCTGCGATTTATGTGCCTCCTCCCCCCCTTGCACCATGGTCTGACCATGAATTTCTCGTca                                                                    >  3:294678/1‑149 (MQ=255)
                                                                              gTATCAGTGACCGCAAGGTCGTGGAAGTTCGAGTCTTCTCCTGGGCACCATACAAAAACCCACTAGCTTGCTAGTGGGTTTTTTCTTGCCTGCGATTTATGTGCCCGCGCCCCCCCTTGCACCATGGGCGGACCATGAATGTCTCGTcc                                                                    >  4:277397/1‑148 (MQ=38)
                                                                                    gTGACCGCAAGGTCGTGGAAGTTCGAGTCTTCTCCTGGGCACCATACAAAAACCCACTAGCTTGCTAGTGGGttttt                                                                                                                                      >  7:255853/1‑77 (MQ=255)
                                                                                    gTGACCGCAAGGTCGTGGAAGTTCGAGTCTTCTCCTGGGCACCATACAAAAACCCACTAGCTTGCTAGTGGGttttt                                                                                                                                      <  8:255853/77‑1 (MQ=255)
                                                                                    gTGACCGCAAGGTCGTGGAAGTTCGAGTCTTCTCCTGGGCACCATACAAAAACCCACTAGCTTGCTAGTGGGTTTTTTCTTGCGTGCGATTTATGTGCCCCCTCCGCCCCTTGCACCATGGTCTGCCCATGAATTTCTAGTCACACGGc                                                              >  2:153252/1‑149 (MQ=255)
                                                                                    gTGACCGCAAGGTCGTGGAAGTTCGAGTCTTCTCCTGGGCACCATACAAAAACCCACTAGCTTGCTAGTGGGTTTTTTCTTGCCTGCGATTTATGTGGCCCCGCCCCCCCGGGCACCCTGGGCTGACCATGAATTTCTCGTCACACGGc                                                              >  2:61356/1‑149 (MQ=38)
                                                                                    gTGACCGCAAGGTCGTGGAAGTTCGAGTCTTCTCCTGGGCACCATACAAAAACCCACTAGCTTGCTAGTGGGTTTTTTCTTGCCTGCGATTTATGTGCGCCCTCCCCCCCTTGCGCCATGGTCTGCCCATGAATTTCTCGTGACACGGc                                                              >  6:29091/1‑149 (MQ=255)
                                                                                         cGCAAGGTCGTGGAAGTTCGAGTCTTCTCCTGGGCACCATACAAAAACCCACTAGCTTGCTAGTGGGTTTTTTCTTGCCTGCGATTTATGTGCCCCCTCCCCCCCTTGCACCATGGTCTGACCATGAATTTCTCGTCACACGGCCACtt                                                         <  4:157408/149‑1 (MQ=255)
                                                                                           cAAGGTCGTGGAAGTTCGAGTCTTCTCCTGGGCACCATACAAAAACCCACTAGCTTGCTAGTGGGTTTTTTCTTGCCTGCGATTTATGTGCGGCCGGCCCCCCGTGCACCATGGTCGGACCAGGAATTTCTCGTACCACGGCCACtt                                                         >  1:61374/1‑147 (MQ=38)
                                                                                                   tGGAAGTTCGAGTCTTCTCCTGGGCACCATACAAAAACCCACTAGCTTGCTAGTGGGTTTTTTCTTGCCTGCGATTTATGTGCGCGCGGCCCCCCTTGCGCCAGGCGGTGAACAAGAATTTAGCGTGACACGCCCGCTTGACaacgata                                               >  4:197733/1‑144 (MQ=255)
                                                                                                     gAAGTTCGAGTCTTCTCCTGGGCACCATACAAAAACCCACTAGCTTGCTAGTGGGTTTTTTCTTGCCTGCGATTTATATGCCCCCTCCCCCCCTTGCACCATGGTCTGACCATGAATTTCTCGTCACACGGCCACTTGAGAAACAAttt                                             <  2:81026/149‑1 (MQ=255)
                                                                                                      aaGTTCGAGTCTTCTCCTTCGCACCATACAAAAACCCACTAGCTTGCTAGTGGGTTTTTTCTTGCCTGCGATTTATGTGCCCCCTCCCCCCCTTGCACCATGGTCTGACCATGAATTTCTCGTCACACGGCCACTTGAGAAACAATTTc                                            <  8:384164/149‑1 (MQ=255)
                                                                                                      aaGTTCGAGTCTTCTCCTGGGCACCATACAAAAACCCACTAGCTTGCt                                                                                                                                                 >  1:130220/1‑48 (MQ=255)
                                                                                                      aaGTTCGAGTCTTCTCCTGGGCACCATACAAAAACCCACTAGCTTGCt                                                                                                                                                 <  2:130220/48‑1 (MQ=255)
                                                                                                          tCGAGTCTTCTCCTGGGCACCATACAAAAACCCACTAGCTTGCTAGTGGGTTTTTTCTTGCCTGCGATTTATGTGCCCCCTCCCCCCCTTGCACCATGGTCTGACCATGAATTTCTCGTCACACGGCCACTTGAGAAACAATTTCGtt                                         <  5:326894/148‑1 (MQ=255)
                                                                                                           cGAGTCTTCTCCTGGGCACCATACAAAAACCCACTAGCTTGCTAGTGGGTTTTTTCTTGCCTGCGATTTATGTGCCCCCTCCCCCCCTTGCACCATGGTCTGACCATGAATTTCTCGTCACACGGCCACTTGAGAAACAATTTCGttt                                        <  4:356405/148‑1 (MQ=255)
                                                                                                            gAGTCTTCTCCTGGGCACCATACAAAAACCCACTAGCTTGCTAGTGGGTTTTTTCTTGCCTGCGATTTATGTGCCCCCTCCCCCCCTTGCACCATGGTCTGACCATGAATTTCTCGTCAAACGGCCACTTGAGAAACAATTTCGTTTa                                       >  2:429345/1‑148 (MQ=255)
                                                                                                            gAGTCTTCTCCTGGGCACCATACAAAAACCCACTAGCTTGCTAGTGGGTTTTTTCTTGCCTGCGATTTATGTGCCCCCTCCCCCCCGTGCACCATGGTCTGACCATGAATTTCTCGTCACACGGCCACTTGAGAAACAATTTCGTTTa                                       >  2:429350/1‑148 (MQ=255)
                                                                                                             aGTCTTCTCCTGGTCACCATACAAAAACCCTCTAGCTTGCTAGTGGGTTTTTTCTTGCCTGCGATTTATGTGCCCCCTCCCCCGCTTGCACCATGGTCTGACCATGAATTTCTCGTCACACGGCCACTTGAGAAACAATTTCGTTTAcc                                     <  8:115960/149‑1 (MQ=255)
                                                                                                                   ctcCTGGGCACCATACAAAAACCCACTAGCTTGCTAGTGGGTTTTTTCTTGCCTGCGATTTATGTGCCTCCTCCCCCCCTTGCACCATGGTCTGACCATGAATTTCTCGTCACACGGCCACTTGAGAAACAATTTCGTTTACCAtt                                  <  4:294678/146‑1 (MQ=255)
                                                                                                                    tcCTGGGCACCATACAAAAACCCACTAGCTTGCTAGTGGGTTTTTTCTTGCCTGCGATTTATGTGCCCCCTCCCCCCCTTGCACCATGGTCTGACCATGAATTTCTCGTCACACGGCCACTTGAGAAACAATTTc                                            <  3:71513/135‑1 (MQ=255)
                                                                                                                    tcCTGGGCACCATACAAAAACCCACTAGCTTGCTAGTGGGTTTTTTCTTGCCTGCGATTTATGTGCCCCCTCCCCCCCATGCACCCTGGTCTGACCATGAATTTCTCGTCACCAGGCGACTTGAGAAACAAattg                                            >  4:71513/1‑134 (MQ=255)
                                                                                                                      cTGGGCACCATACAAAAACCCACTAGCTTGCTAGTGGGTTTTTTCTTGCCTGCGATTTATGTGGCCCCTCCCCCCCTTGCACCATGGTCTGACCATGAATTTCTAGTCAAACGGCCACTTGAGAAACAATTTCGTTTAACATTGTTTc                             >  4:144755/1‑148 (MQ=255)
                                                                                                                      cTGGGCACCATACAAAAACCCACTAGCTTGCTAGTGGGTTTTTTCTTGCCTGCGATTTATGTGCCCCCTCCCCCCCTTGCACCATGGTCTGACCATGAATTTCTCGTCACACGGCCACTTGAGaaa                                                   <  3:394098/126‑1 (MQ=255)
                                                                                                                      cTGGGCACCATACAAAAACCCACTAGCTTGCTAGTGGGTTTTTTCTTGCCTGCGATTTATGTGCCCCCTCCCCCCCTTGCACCATGGTCTGACCATGAATTTCTCGTCACACGGCCACTTGAGAAACAATTTCGTTTACCATTGTTTcc                            <  5:248933/149‑1 (MQ=255)
                                                                                                                        gggCACCATACAAAAACCCACTAGCTTGCTAGTGGGTTTTTTCTTGCCTGCGATTTATGTGCCCCCTCCCCCCCTTGCACCATGGTCTGACCATGAATTTCTCGTCACGCTGCCACTTGAGAAACAATTTCGTTTACCATTGTTTCCa                           >  6:435512/1‑148 (MQ=255)
                                                                                                                          tgacgaTACAAAAACCTAGTAGCTTGCTAGTGTGTTTTTTCTTGCCTGCGATTTATGTTCGCCCTCCCCCCGTAGCACGATGGTCTGACCATGAATTTCTCGTCACACGGCCACTTTAGAAACAATTTCGTTTACCATTGTTTCCaaa                         <  3:309132/143‑1 (MQ=255)
                                                                                                                            aCCATACAAAAACCCACTAGCTTGCTAGTGGGTTTTTTCTTGCCTGCGATTTATGTGCCCCCTCCCCCCCTTGCACCATGGTCTGACCATGAATTTCTCGTCACACGGCCACTTGAGAAACAATTTCGTTTACCATTGTTTCCAAatat                      >  3:340061/1‑149 (MQ=255)
                                                                                                                               aTACAAAAACCCACTAGCTTGCTAGTGGGTTTTTTCTTGCCTGCGATTTATGTGCCCCCTCCCCCCCTTGCACCATGGTCTGACCATGAATTTCTCGTCACACGGCCACTTGAGAAACAATTTCGTTTACCATTGTTTCCAAATATGTa                   <  5:29091/149‑1 (MQ=255)
                                                                                                                                tACAAAAACCCACTAGCTTGCTAGTGGGTTTTTTCTTGCCTGCGATTTATGTGCCCCCTCCCCCCCTTGCACCATGGTCTGACCATGAATTTCTCGTCACACGGCCACTTGAGAAACAATTTCGTTTACCATTGTTTCCAAATATGTAg                  <  4:294488/149‑1 (MQ=255)
                                                                                                                                 gaaaaaaCCCACTAGCAAGCTAGTGGGTTTTTTCTTGCCTGCGATTTATGTGCCCCCTCCCCCCCTTGCACCATGGTCTGAGCATGAATTTCTCGTCACACGGCCACTTGAGAAACAATTTCGTTTACCATTGTTTCCAAATATGTAGc                 >  4:294898/3‑149 (MQ=255)
                                                                                                                                     aaaCCCACTAGCTTGCTAGTGGGTTTTTTCTTGCCTGCGATTTATGTGCCCCCTCCCCCCCTGGCACCATGGTCTGACCATGAATTTCTCGTCACACGGCCACTTGAGAAACAATTTCGTTTACCATTGTTTCCAAATATGTAGCCGt              >  2:106605/1‑148 (MQ=255)
                                                                                                                                       aCCCACTAGCTTGCTAGTGGGTTTTTTCTTGCCTGCGATTTATGTGCCCCCTCCCCCCCTTGCACCATGGTGTGACCATGAATTTCTCGTCACACGGCCACTTGAGAAACAATTTCGTTTACCATTGTTTCCAAATATGTAGCCGTAAg           >  5:65316/1‑149 (MQ=255)
                                                                                                                                        cccACTAGCAAGCTAGTGGGTTTTTTCTTGCCTGCGATTTATGTGCCCCCTCCCCCCCTGGCACCAAGGTATGACCATGAATTTGTCGTCACACGGCCACTTGAGAAACAATTTCGTTTAAAATTGTTTCCAAATATATAGCCGTaa            >  8:252135/1‑147 (MQ=255)
                                                                                                                                         ccACTAGCAAGCTAGTGGGTTTTTTCTTGCCTGCGATTTATGTGCACCCTCCCCCCCTTGCACAATGGAATCACAATGAATTTCTCATCACACCTCCACTTAATAAACAAATTCGTTTACAATTGTTTACAAATATGTACCAGTAAGc          >  2:271970/1‑148 (MQ=255)
                                                                                                                                          cACTAGCTTGCTAGTGGGTTTTTTCTTGCCTGCGATTTATGTGCCCCCTCCCCCCCTTGCACCATGGTCTGACCATGAATTTCTCGTCACACGGCCACTTGAGAAACAATTTCGTTTACCATTGTTTCCAAATATGTAGCCGTAAGCGa        <  3:277397/149‑1 (MQ=255)
                                                                                                                                          cACTAGCTTGCTAGTGGGTTTTTTCTTGCCTGCGATTTATGTGCCCCCTCCCCCCCTTGCACCATGGTCTGACCATGAATTTCTCGTCACACGGCCACTTGAGAAACAATTTCGTTTACCATTGTTTCCAAATATGTAGCCGTAAGCGa        <  7:92501/149‑1 (MQ=255)
                                                                                                                                          cACTAGCTTGCTAGAGGCTTTTTTCTTGCCTGCGATTTATGTGCCCCCTCCCCCCCATGCACCAAGGACTCACCATAAATTTCTCGTAACAAGGACACTACAGAAACAATATCCTATACCATTATTACCAAATATGAAGCAGTAAga          >  4:246241/1‑146 (MQ=255)
                                                                                                                                                cTTGCTAGTGGGTTTTTTCTTGCCTGCGATTTATGTGCCCCCTCCCCCCCTTGCACCATGGTCTGACCATGAATTTCTCGTCACACGGCCACTTGAGAAACAATTTCGTTTACCATTGTTTCCAAATATGTAGCCGTAAGCGAGGAAGt  <  1:153252/149‑1 (MQ=255)
                                                                                                                                                  |                                                                                                                                                  
AAGCGAACGATAGTCGCCTTCTGCAGATGATGTAAAGCTTTGCCCAGGTGGCGGAATTGGTAGACGCGCTGGTTTCAGGTATCAGTGACCGCAAGGTCGTGGAAGTTCGAGTCTTCTCCTGGGCACCATACAAAAACCCACTAGCTTGCTAGTGGGTTTTTTCTTGCCTGCGATTTATGTGCCCCCTCCCCCCCTTGCACCATGGTCTGACCATGAATTTCTCGTCACACGGCCACTTGAGAAACAATTTCGTTTACCATTGTTTCCAAATATGTAGCCGTAAGCGAGGAAGT  >  NC_002947/5549460‑5549752

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