Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 660,909 T→C 5.7% intergenic (+41/‑126) speA → / → PP_0568 biosynthetic arginine decarboxylase/hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_002947660,9090TC 5.7% 212.6 / 2.1 70intergenic (+41/‑126)speA/PP_0568biosynthetic arginine decarboxylase/hypothetical protein
Reads supporting (aligned to +/- strand):  ref base T (33/33);  new base C (2/2);  total (35/35)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.19e-01

CATTACCGCGACAAGGTTGCCAGCGCCAAGATCACTGCGCGTGAGCGGACCCAGTATCTGGATGCCTTGCGCCTGGGCTTGACCCGGTCTTCGTACCTGTCGTCGTAATTGACGGGGCCGCTCCTGCAAAAGCCTGTGCTGCCTTTTGTGGGAGCGGCCTTGCGTCGCGATATGAGGGCGAAGCCCTCACCAGCGTCACAACCGGTAAATATCCGGCCCTGGTCCATCGCACCAATGCCGCTACCTTTGAATGCAATGTTCAAAGGAGGACGATATGCAACGTCCAGGCTCC  >  NC_002947/660761‑661052
                                                                                                                                                    |                                                                                                                                               
cATTACCGCGACAATGTTGCCAGCGCCAAGATCACTGCGCGTGAGCGGACCCAGTATCTGGATGCCTTGCGCCTGGGCTTGACCCGGTCTTCGTACCTGTCGTCGTAATTGACGGGGCCGCTCCTGCAAAAGCCTGTGCTGCCTTTtgt                                                                                                                                                 <  6:277681/149‑1 (MQ=255)
    aCCGCGACAAGGTTGCCAGCGCCAAGATCACTGCGCGTGAGCGGACCCAGTATCTGGATGCCTTGCGCCTGGGCTTGACCCGGTCTTCGTACCTGTCGTCGTAATTGACGGGGCCGCTCCTGCAAAAGCCTGTGCTGCCTTTTGTgg                                                                                                                                               >  4:245534/1‑147 (MQ=255)
     ccGCGACAAGGTTGCCAGCGCCAAGATCACTGCGCGTGAGCGGACCCAGTATCTGGATGCCTTGCGCCTGGGCTTGACCCGGTCTTCGTACCTGTCGTCGTAATTGACGGGGCCGCTCCTGCAAAAGCCTGTGCTGCCTTTTGTGGGAg                                                                                                                                            >  6:336577/1‑149 (MQ=255)
       gcgACAAGGTTGCCAGCGCCAAGATCACTGCGCGTGAGCGGACCCAGTATCTGGATGCCTTGCGCCTGGGCTTGACCCGGTCTTCGTACCTGTCGTCGTAATTGACGGGGCCGCTCCTGCAAAAGCCTGTGCTGCCTTTTGTGGGAGCg                                                                                                                                          <  3:139042/149‑1 (MQ=255)
            aaGGTTGCCAGCGCCAAGATCACTGCGCGTGAGCGGACCCAGTATCTGGATGCCTTGCGCCTGGGCTTGACCCGGTCTTCGTACCTGTCGTCGTAATTGACGGGGCCGCTCCTGCAAAAGCCTGTGCTGCCTTTTGTGGGAGCGGCCtt                                                                                                                                     <  8:149215/149‑1 (MQ=255)
              ggTTGGCAGCGCCAAGATCACTGCGCGTGAGGGGAGCCAGTATCTGGATGCCTTGCGCCTGGGCTTGACCCGGTCTTCGTACCTGTCGTCGTAATTGACGGGGCCGCTCCTGCAAAAGCCTGTGCTGCCTTTTGTGGGAGCGGCCTTGc                                                                                                                                   <  3:162002/149‑1 (MQ=255)
                   ccAGCGCCAAGATCACTGCGCGTGAGCGGACCCAGTATCTGGATGCCTTGCGCCTGGGCTTGACCCGGTCTTCGTACCTGTCGTCGTAATTGACGGGGCCGCTCCTGCAAAAGCCTGTGCTGCCTTTTGTGGGAGCGGCCTTGCGTcgc                                                                                                                              <  5:281102/149‑1 (MQ=255)
                    cAGCGCCAAGATCACTGCGCGTGAGCGGACCCAGTATCTGGATGCCTTGCGCCTGGGCTTGACCCGGTCTTCGTACCTGTCGTCGTAATTGACGGGGCCGCTCCTGCAAAAGCCTGTGCTGCCTTTTGTGGGAGCGGCCTTGCGTcgc                                                                                                                              >  6:279891/1‑148 (MQ=255)
                        gcCAAGATCACTGCGCGTGAGCGGACCCAGTATCTGGATGCCTTGCGCCTGGGCTTGACCCGGTCTTCGTACCTGTCGTCGTAATTGACGGGGCCGCTCCTGCAAAAGCCTGTGCTGCCTTTTGTGGGAGCGGCCTTGCGTCGCGatat                                                                                                                         >  2:43390/1‑149 (MQ=255)
                          cAAGTTCACTGCGCGTTAGCGGACCCAGTATCTGGATGCCTTGCGCGTGGGCTTGAGCCGGTCTTCGTACCTGTCGTCGTAATTGACGGGGCCGCTCCTGCAAAAGCCTGTGCTGCCTTTTGTGGGAGCGGCCTTGCGTCGCGATATGa                                                                                                                       <  5:231774/149‑1 (MQ=255)
                            aGATCACTGCGCGTGAGCGGACCCAGTATCTGGATGCCTTGCGCCTGGGCTTGACCCGGTCTTCGTACCTGTCGTCGTAATTGACGGGGCCGCTCCTGCAAAAGCCTGTGCTGCCTTTTGTGGGAGCGGCCTTGCGTCGCGATATGAgg                                                                                                                     <  2:207414/149‑1 (MQ=255)
                             gATCACTGCGCGTGAGCGGACCCAGTATCTGGATGCCTTGCGCCTGGGCTTGACCCGGTCTTCGTACCTGTCGTCGTAATTGACGGGGCCGCTCCTGCAAAAGCCTGTGCTGCCTTTTGTGGGAGCGGCCTTGCGTCGCGATATGAggg                                                                                                                    <  7:279860/149‑1 (MQ=255)
                                  cTGCGCGTGAGCGGACCCAGTATCTGGATGCCTTGCGCCTGGGCTTGACCCGGTCTTCGTACCTGTCGTCGTAATTGACGGGGCCGCTCCTGCAAAAGCCTGTGCTGCCTTTTGTGGGAGCGGCCTTGCGTCGCGATATGAGGGCGAAg                                                                                                               <  5:275530/149‑1 (MQ=255)
                                  cTGCGCGTGAGCGGACCCAGTATCTGGATGCCTTGCGCCTGGGCTTGACCCGGTCTTCGTACCTGTCGTCGTAATTGACGGGGCCGCTCCCACAAAAGGCAGCACAGGCTTTTGCAGGAGCgg                                                                                                                                         >  1:92615/1‑123 (MQ=25)
                                  cTGCGCGTGAGCGGACCCAGTATCTGGATGCCTTGCGCCTGGGCTTGACCCGGTCTTCGTACCTGTCGTCGTAATTGACGGGGCCGCTCCCACAAAAGGCAGCACAGGCTTTTGCAGGAGCgg                                                                                                                                         <  2:92615/123‑1 (MQ=25)
                                   tGCGGGTTAGCGGACCCAGTATCTGGATGGCTTGCGGCTGGGGTTGACCCGGTCTTCGTACCTGTCGTCGTAATTGACGGGGCCGGTCCTGCAAAAGCCTGTGCTGCCTTTTGTGGGAGCGGCCTTGCGTCGCGATATGAGGGCGaa                                                                                                                <  7:253095/147‑1 (MQ=255)
                                         tGAGCGGACCCAGTATCTGGATGCCTTGCGCCTGGGCTTGACCCGGTCTTCGTACCTGTCGTCGTAATTGACGGGGCCGCTCCTGCAAAAGCCTGTGCTGCCTTTTGTGGGAGCGGCCTTGCGTCGCGATATGAGGGCGAAGCCCTCAc                                                                                                        <  7:339570/149‑1 (MQ=255)
                                           aGCGGACCCAGTATCTGGATGCCTTGCGCCTGGGCTTGACCCGGTCTTCGTACCTGTCGTCGTAATTGACGGGGCCGCTCCTGCAAAAGCCTGTGCTGCCTTTTGTGGGAGCGGCCTTGCGTCGCGATATGAGGGCTTCGCCCTCata                                                                                                       <  7:171340/148‑3 (MQ=255)
                                              ggACCGAGTATCTGGATGCCTTGCGCCTGGGCTTGACCCGGTCTTCGTACCTGTCGTCGTAATTGACGGGGCCGCTCCTGCAAAAGCCTGTTCTGCCTTTTGTTGGAGCGGCCTTGCGTCGCGATATGAGGGCGAAGCCCTCACCAGc                                                                                                    <  3:349964/148‑1 (MQ=255)
                                              ggACCCAGTATCTGGATGCCTTGCGCCTGGGCTTGACCCGGTCTTCGTACCTGTCGTCGTAATTGACGGGGCCGCTCCTGCAAAAGCCTGTGCTGCCTTTTGTGGGAGCGGCCTTGCGTCGCGATATGAGGGCTTCGCCCTCatatcgc                                                                                                   <  8:159069/149‑7 (MQ=255)
                                                  ccAGTATCTGGATGCCTTGCGCCTGGGCTTGACCCGGTCTTCGTACCTGTCGTCGTAATTGACGGGGCCGCTCCTGCAAAAGCCTGTGCTGCCTTTTGT‑GGAGCGGCCTTGCGTCGCGATATGAGGGCGAAGCCCTCACCAGCGTcac                                                                                               >  1:273237/1‑148 (MQ=255)
                                                  ccAGTATCTGGATGCCTTGCGCCTGGGCTAGACCCGGTCTTCGTACCTGTCGTCGTAATTGACGGGGCCGCTCCTGCAAAAGCCTGTGCTGCCTTTTGTGGGAGCGGCCTTGCGTCGCGATATGAGGGCGAAGCCCTCACCAGCGTca                                                                                                >  1:214159/1‑148 (MQ=255)
                                                   cAGTATCTGGATGCCTTGCGCCTGGGCTTGACCCGGTCTTCGTACCTGTCGTCGTAATTGACGGGGCCGCTCCTGCAAAAGCCTGTGCTGCCTTTTGTGGGAGCGGCCTTGCGTCGCGATATGAGGGCGAAGCCCTCACCAGCGTctca                                                                                              <  5:336577/149‑1 (MQ=255)
                                                               gCCTTGCGCCTGGGCTTGACCCGGTCTTCGTACCTGTCGTCGTAATTGACGGGGCCGCTCCTGCAAAAGCCTGTGCTGCCTTTTGTGGGAGCGGCCTTGCGTCGCGATATGAGGGCGAAGCCCTCACCAGCGTCACAACCGGTAAatat                                                                                  >  8:211675/1‑149 (MQ=255)
                                                                 cTTGCGCCTGGGCTTGACCCGGTCTTCGTACCTGTCGTCGTAATTGACGGGGCCGCTCCTGCAAAAGCCTGTGCTGCCTTTTGTGGGAGCGGCCTTGCGTCGCGATATGAGGGCGAAGCCCTCACCAGCGTCACAACCGGTAAATATcc                                                                                >  3:66653/1‑149 (MQ=255)
                                                                    gcgcCTGGGCTTGACCCGGTCTTCGTACCTGTCGTCGTAATTGACGGGGCCGCTCCTGCAAAAGCCTGTGCTGCCTTTTGTGGGAGCGGCCTTGCGTCGCGATATGAGGGCGAAGCCCTCACCAGCGTCACAACCGGTAAATATCCg                                                                               <  1:43390/147‑1 (MQ=255)
                                                                           ggCTTGACCCGGTCTTCGTACCTGTCGTCGTAATTGACGGGGCCGCTCCTGCAAAAGCCTGTGCTGCCTTTTGTGGGAGCGGCCTTGCGTCGCGATATGAGGGCGAAGCCCTCACCAGCGTCACAACCGGTAAatat                                                                                  >  1:150070/1‑137 (MQ=255)
                                                                           ggCTTGACCCGGTCTTCGTACCTGTCGTCGTAATTGACGGGGCCGCTCCTGCAAAAGCCTGTGCTGCCTTTTGTGGGAGCGGCCTTGCGTCGCGATATGAGGGCGAAGCCCTCACCAGCGTCACAACCGGTAAatat                                                                                  <  2:150070/137‑1 (MQ=255)
                                                                           ggCTTGACCCGGTCTTCGTACCTGTCGTCGTAATTGACGGGGCCGCTCCTGCAAAAGCCTGTGCTGCCTTTTGTGGGAGCGGCCTTGCGTCGCGATATGAGGGCGAAGCCCTCACCAGCGTCACAACCGGTAAATATCCGGCCCTgg                                                                        <  5:279891/147‑1 (MQ=255)
                                                                           ggCTTGACCCGGTCTTCGTACCTGTCGTCGTAATTGACGGGGCCGCTCCTGCAAAAGCCTGTGCTGCCTTTTGTGGGAGCGGCCGTGCGTCGCGATATGAGGGCGAAGCCCTCACCAGCGTCACAACCGGTAAATATCCGGCCCTGGTc                                                                      <  4:276663/149‑1 (MQ=255)
                                                                            gCTTGACCCGGTCTTCGTACCTGTCGTCGTAATTGACGGGGCCGCTCCTGCAAAAGCCTGTGCTGCCTTTTGTGGGAGCGGCCTTGCGTCGCGATATGAGGGCGAAGCCCTCACCAGCGTCACAACCGGTAAATATCCGGCCCTGGTcc                                                                     <  3:245534/149‑1 (MQ=255)
                                                                                  cccGGTCTTCGTACCTGTCGTCGTAATTGACGGGGCCGCTCCTGCAAAAGCCTGTGCTGCCTTTTGTGGGAGCGGCCTTGCGTCGCGATATGAGGGCGAAGCCCTCACCAGCGTCACAACCGGTAAATATCCGGCCCTGGTCCATCGCa                                                               >  7:14718/1‑149 (MQ=255)
                                                                                  cccGGTCTTCGTACCTGTCGTCGTAATTGACGGGGCCGCTCCTGCAAAAGCCTGTGCTGCCTTTTGTGGGAGCGGCCTTGCGTCGCGATATGAGGGCGAAGCCCTCACCAGCGTCACAACCGGTAAATATCCGGCCCTGGTCCATCGCa                                                               >  2:331916/1‑149 (MQ=255)
                                                                                  cccGGTCTTCGTACCTGTCGTCGTAATTGACGGGGCCGCTCCTGCAAAAGCCTGTGCTGCCTTTTGT‑GGAGCGGCCTTGCGTCGCGATATGAGGGCGAAGCCCTCACCAGCGTCACAACCGGTAAATATCCGGCCCTGGTCCATCGCAc                                                              <  2:273237/149‑1 (MQ=255)
                                                                                   ccGGTCTTCGTACCTGTCGTCGTAATTGACGGGGCCGCTCCTGCAAAAGCCTGTGCTGCCTTTTGTGGGAGc                                                                                                                                           >  7:223823/1‑72 (MQ=255)
                                                                                   ccGGTCTTCGTACCTGTCGTCGTAATTGACGGGGCCGCTCCTGCAAAAGCCTGTGCTGCCTTTTGTGGGAGc                                                                                                                                           <  8:223823/72‑1 (MQ=255)
                                                                                   ccGGTCTTCGTACCTGTCGTCGTAATTGACGGGGCCGCTCCTGCAAAAGCCTGTGCTGCCTTTTGTGGGAGCGGCCTTGCGTCGCGATATGAGGGCGAAGCCCTCACCAGCGTCACAACCGGTAAATATGCGGCCCTGGTCCGTCGcca                                                              >  6:290512/1‑147 (MQ=255)
                                                                                   ccGGTCTTCGTACCTGTCGTCGTAATTGACGGGGCCGCTCCTGCAAAAGCCTGTGCTGCCTTTTGTGGGAGCGGCCTTGCGTCGCGATATGAGGGCGAAGCCCTCACCAGCGTCACAACCGGTAAATATCCGGCCCTGGTCCATCGCAc                                                              <  5:290512/149‑1 (MQ=255)
                                                                                               ccTGTCGTCGTAATTGACGGGGCCGCTCCTGCAAAAGCCCGTGCTGCCTTTTGTGGGAGCGGCCTTGCGTCGCGATATGAGGGCGAAGCCCTCACCAGCGTCACAACCGGTAAATATCCGGCCCTGGTCCATCGCACCAATGCCGCTAc                                                  <  1:147861/149‑1 (MQ=255)
                                                                                                 tGTCGTCGTAATTGACGGGGCCGCTCCTGCAAAAGCCTGTGCTGCCTTTTGTGGGAGCGGCCTTGCGTCGCGATATGAGGGCGAAGCCCTCACCAGCGTCACAACCGGTAAATATCCGGCCCTGGTCCATCGCACCAATGCCGCTACCt                                                >  3:91641/1‑149 (MQ=255)
                                                                                                  gtcgtcGTAATTGACGGGGCCGCTCCTGCAAAAGCCTGTGCTGCCTTTTGTGGGAGCGGCCTTGCGTCGCGATATGAGGGCGAAGCCCTCACCAGCGTCACAACCGGTAAATATCCGGCCCTGGTCCATCGCACCAATGCCGCTACCtt                                               >  3:355428/1‑149 (MQ=255)
                                                                                                   tcgtcgTAATTGACGGGGCCGCTCCTGCAAAAGCCTGTGCTGCCTTTTGTGGGAGCGGCCTTGCGTCGCGATATGAGGGCGAAGCCCTCACCAGCGTCACAACCGGTAAATATCCGGCCCTGGTCCATCGCACCAATGCCGCTACCttt                                              <  1:331916/149‑1 (MQ=255)
                                                                                                     gtcgtAATTGACGGGGCCGCTCCTGCAAAAGCCTGTGCTGCCTTTTGTGGGAGCGGCCTTGCGTCGCGATATGAGGGCGAAGCCCTCACCAGCGTCACAACCGGTAAATATCCGGCCCTGGTCCATCGCACCAATGCCGCTACCTTTg                                             >  2:211350/1‑148 (MQ=255)
                                                                                                        gtAATTGACGGGGCCGCTCCTGCAAAAGCCTGTGCTGCCTTTTGTGGGAGCGGCCTTGCGTCGCGATATGAGGGCGAAGCCCTCACCAGCGTCACAACCGGTAAATATCCGGCCCTGGTCCATCGCACCAATGCCGCTACCTTTGAATg                                         >  5:312207/1‑149 (MQ=255)
                                                                                                           aTTGACGGGGCCGCTCCTGCAAAAGCCTGTGCTGCCTTTTGTGGGAGCGGCCTTGCGTCGCGATATGAGGGCGAAGCCCTCACCAGCGTCACAACCGGTAAATATCCGGCCCTGGTCCATCGCACCAATGCCGCTACCTTTGAATGCaa                                      >  1:310739/1‑149 (MQ=255)
                                                                                                             tGACGGGGCCGCTCCTGCAAAAGCCTGTGCTGCCTTTTGTGGGAGCGGCCTTGCGTCGCGATATGAGGGCGAAGCCCTCACCAGCGTCACAACCGGTAAATATCCGGCCCTGGTCCATCGCACCAATGCCGCTACCTTTGAATGCAATg                                    >  8:7238/1‑149 (MQ=255)
                                                                                                                 ggggCCGCTCCTGCAAAAGCCTGTGCTGCCTTTTGTGGGAGCGGCCTTGCGTCGCGATATGAGGGCGAAGCCCTCACCAGCGTCACAACCGGTAAATATCCGGCCCTGGTCCATCGCACCAATGCCGCTACCTTTGAATGCAATGTTCa                                <  4:355428/149‑1 (MQ=255)
                                                                                                                 ggggCCGCTCCTGCAAAAGCCTGTGCTGCCTTTTGTGGGAGCGGCCTTGCGTCGCGATATGAGGGCGAAGCCCTCACCAGCGTCACAACCGGTAAATATCCGGCCCTGGTCCATCGCACCAATGCCGCTACCTTTGAATGCAATGTTCa                                <  1:249409/149‑1 (MQ=255)
                                                                                                                     ccGCTCCCACAAAAGGCAGCACAGGCTTTTGCAGGAGCGGCCTTGCGTCGCGATATGAGGGCGAAGCCCTCACCAGCGTCACAACCGGTAAATATCCGGCCCTGGTCCATCGCACCAATGCCGCTACCTTTGAATGCAATGtt                                  <  1:152956/143‑1 (MQ=37)
                                                                                                                     ccGCTCCCACAAAAGGCAGCACAGGCTTTTGCAGGAGCGGCCTTGCGTCGCGATATGAGGGCGAAGCCCTCACCAGCGTCACAACCGGTAAATATCCGGCCCTGGTCCATCGCACCAATGCCGCTACCTTTGAATGCAATGtt                                  >  2:152956/1‑143 (MQ=37)
                                                                                                                      cGCTCCTGCAAAAGCCTGTGCTGCCTTTTGTGGGAGCGGCCTTGCGTCGCGATATGAGGGCGAAGCCCTCACCAGCGTCACAACCGGTAAATATCCGGCCCTGGTCCATCGCACCAATGCCGCTACCTTTGAATGCAATGTTCAAag                             >  1:29765/1‑147 (MQ=255)
                                                                                                                         tCCTGCAAAAGCCTGTGCTGCCTTTTGTGGGAGCGGCCTTGCGTCGCGATATGAGGGCGAAGCCCTCACCAGCGTCACAACCGGTAAATATCCGGCCCTGGTCCATCGCACCAATGCCGCTACCTTTGAATGCAATGTTCAAaggag                          >  6:96023/1‑147 (MQ=255)
                                                                                                                         tCCTGCAAAAGCCTGTGCTGCCTTTTGTGGGAGCGGCCTTGCGTCGCGATATGAGGGCGAAGCCCTCACCAGCGTCACAACCGGTAAATATCCGGCCCTGGTCCATCGCACCAATGCCGCTACCTTTGAATGCAATGTTCAAAggagga                        >  2:12191/1‑149 (MQ=255)
                                                                                                                           cTGCAAAAGCCTGTGCTGCCTTTTGTGGGAGCGGCCTTGCGTCGCGATATGAGGGCGAAGCCCTCACCAGCGTCACAACCGGTAAATATCCGGCCCTGGTCCATCGCACCAATGCCGCTACCTTTGAATGCAATGTTCAAAggagga                        <  7:7238/147‑1 (MQ=255)
                                                                                                                           cTGCAAAAGCCTGTGCTGCCTTTTGTGGGAGCGGCCTTGCGTCGCGATATGAGGGCGAAGCCCTCACCAGCGTCACAACCGGTAAATATCCGGCCCTGGTCCATCGCACCAATGCCGCTACCTTTGAATGCAATGTTCAAAGGAGGAc                       <  4:66653/148‑1 (MQ=255)
                                                                                                                           cTGCAAAAGCCTGTGCTGCCTTTTGTGGGAGCGGCCTTGCGTCGCGATATGAGGGCGAAGCCCTCACCAGCGTCACAACCGGTAAATATCCGGCCCTGGTCCATCGCACCAATGCCGCTACCTTTGAATGCAATGTTCAAAGGAGGAc                       >  7:137359/1‑148 (MQ=255)
                                                                                                                           cTGCAAAAGCCTGTGCTGCCTTTTGTGGGAGCGGCCTTGCGTCGCGATATGAGGGCGAAGCCCTCACCAGCGTCACAACCGGTAAATATCCGGCCCTGGTCCATCGCACCAATGCCGCTACCTTTGAATGCAATGTTCAAAGGAGGAc                       >  2:222441/1‑148 (MQ=255)
                                                                                                                             gCAAAAGCCTGTGCTGCCTTTTGTGGGAGCGGCCTTGCGTCGCGATATGAGGGCGAAGCCCTCACCAGCGTCACAACCGGTAAATATCCGGCCCTGGTCCATCGCACCAATGCCGCTACCTTTGAATGCAATGTTCaaa                              <  3:287484/139‑1 (MQ=255)
                                                                                                                             gCAAAAGCCTGTGCTGCCTTTTGTGGGAGCGGCCTTGCGTCGCGATATGAGGGCGAAGCCCTCACCAGCGTCACAACCGGTAAATATCCGGCCCTGGTCCATCGCACCAATGCCGCTACCTTTGAATGCAATGTTCaaa                              >  4:287484/1‑139 (MQ=255)
                                                                                                                                aaaGCCTGTGCTGCCTTTTGTGGGAGCGGCCTTGCGTCGCGATATGAGGGCGAAGCCCTCACCAGCGTCACAACCGGTAAATATCCGGCCCTGGTCCATCGCACCAATGCCGCTACCTTTGAATGCAAt                                     <  3:346061/129‑1 (MQ=255)
                                                                                                                                aaaGCCTGTGCTGCCTTTTGTGGGAGCGGCCTTGCGTCGCGATATGAGGGCGAAGCCCTCACCAGCGTCACAACCGGTAAATATCCGGCCCTGGTCCATCGCACCAATGCCGCTACCTTTGAATGCAAt                                     >  4:346061/1‑129 (MQ=255)
                                                                                                                                aaaGCCTGTGCTGCCTTTTGTGGGAGCGGCCTTGCGTCGCGATATGAGGGCGAAGCCCTCACCAGCGTCACAACCGGTAAATATCCGGCCCTGGTCCATCGCACCAATGCCGCTACCTTTGAATGCAATGTTCAAAGGAGGACGATATg                 >  8:186582/1‑149 (MQ=255)
                                                                                                                                aaaGCCTGTGCTGCCTTTTGTGGGAGCGGCCTTGCGTCGCGATATGAGGGCGAAGCCCTCACCAGCGTCACAACCGGTAAATATCCGGCCCTGGTCCATCGCACCAATGCCGCTACCTTTGAATGCAATGTTCAAAGGAGGACGATATg                 <  8:14718/149‑1 (MQ=255)
                                                                                                                                    ccTGTGCTGCCTTTTGTGGGAGCGGCCTTGCGTCGCGATATGAGGGCGAAGCCCTCACCAGCGTCACAACCGGTAAATATCCGGCCCTGGTCCATCGCACCAATGCCGCTACCTTTGAATGCAATGTTCAAAGGAGGACGATATGCAAc             <  2:214159/149‑1 (MQ=255)
                                                                                                                                    ccTGTGCTGCCTTTTGTGGGAGCGGCCTTGCGTCGCGATATGAGGGCGAAGCCCTCACCAGCGTCACAACCGGTAAATATCCGGCCCTGGTCCATCGCACCAATGCCGCTACCTTTGAATGCAATGTTCAAAGGAGGACGATATGCAAc             >  7:205010/1‑149 (MQ=255)
                                                                                                                                          ctgcCTTTTGTGGGAGCGGCCTTGCGTCGCGATATGAGGGCGAAGCCCTCACCAGCGTCACAACCGGTAAATATCCGGCCCTGGTCCATCGCACCAATGCCGCTACCTTTGAATGCAATGTTCAAAGGAGGACGATATGCAACGTCCAg       >  4:209214/1‑149 (MQ=255)
                                                                                                                                           tgcCTTTTGTGGGAGCGGCCTTGCGTCGCGATATGAGGGCGAAGCCCTCACCAGCGTCACAACCGGTAAATATCCGGCCCTGGTCCATCGCACCAATGCCGCTACCTTTGAATGCAATGTTCAAAGGAGGACGATATGCAACGTCCAgg      >  4:324933/1‑149 (MQ=255)
                                                                                                                                            gcCTTTTGTGGGAGCGGCCTTGCGTCGCGAAATGAGGGCGAAGCCCTCACCAGCGTCACAACCGGTAAATATCCGGCCCTGGTCCATCGCACCAATGCCGCTACCTTTGAATGCAATGTTCAAAGGAGGACGATATGCAACGTCCAGGc     >  1:220332/1‑149 (MQ=255)
                                                                                                                                             ccTTTTGTGGGAGCGGCCTTGCGTCGCGATATGAGGGCGAAGCCCTCACCAGCGTCACAACCGGTAAATATCCGGCCCTGGTCCATCGCACCAATGCCGCTACCTTTGAATGCAATGTTCAAagg                            >  4:142176/1‑125 (MQ=255)
                                                                                                                                              cTTTTGTGGGAGCGGCCTTGCGTCGCGATATGAGGGCGAAGCCCTCACCAGCGTCACAACCGGTAAATATCCGGCCCTGGTCCATCGCACCAATGCCGCTACCTTTGAATGCAATGTTCAAAGGAGGACGATATGCAACGTCCAGGCTc   <  1:222441/149‑1 (MQ=255)
                                                                                                                                                tttGTGGGAGCGGCCTTGCGTCGCGATATGAGGGCGAAGCCCTCACCAGCGTCACAACCGGTAAATATCCGGCCCTGGTCCATCGCACCAATGCCGCTACCTTTGAATGCAATGTTCAAAGGAGGACGATATGCAACGTCCAGGCTcc  >  3:168422/1‑148 (MQ=255)
                                                                                                                                                    |                                                                                                                                               
CATTACCGCGACAAGGTTGCCAGCGCCAAGATCACTGCGCGTGAGCGGACCCAGTATCTGGATGCCTTGCGCCTGGGCTTGACCCGGTCTTCGTACCTGTCGTCGTAATTGACGGGGCCGCTCCTGCAAAAGCCTGTGCTGCCTTTTGTGGGAGCGGCCTTGCGTCGCGATATGAGGGCGAAGCCCTCACCAGCGTCACAACCGGTAAATATCCGGCCCTGGTCCATCGCACCAATGCCGCTACCTTTGAATGCAATGTTCAAAGGAGGACGATATGCAACGTCCAGGCTCC  >  NC_002947/660761‑661052

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