Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,560,601 T→A 7.6% intergenic (‑80/+94) PP_4887 ← / ← PP_4888 hypothetical protein/methyl‑accepting chemotaxis transducer

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,560,6010TA 7.6% 573.8 / 8.6 210intergenic (‑80/+94)PP_4887/PP_4888hypothetical protein/methyl‑accepting chemotaxis transducer
Reads supporting (aligned to +/- strand):  ref base T (69/114);  new base A (17/4);  total (91/119)
Fisher's exact test for biased strand distribution p-value = 2.50e-04
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 214 total aligned reads displayed.

CCGTACGCTGTCACTGTGCAGGAAATTCAGGTCGGCCTGGGCAGCGATTAGTCGGGCGGACATGAGCAACTCCTTTTGTCATGTGTCGCCGACAGCATAGCAACACCGCTCAGGTGTGGGAGCGGCCTTGTGTCGCGAAAGGGCCGCATAGCGGCCCCGGCAATTTCCGCGTGTATGCAACAAACCTGGGGCTGCTATGCAGCCCTTTCGCGACACAAGGCCCACAAAGGACCGCGCAGGCTTCAGGCCCGGAACTGGCCCAGGCTGGCGCGCAGCTGAGCCGCCAGATCATCCAGCA  >  NC_002947/5560453‑5560750
                                                                                                                                                    |                                                                                                                                                      
   tACGCTGTCACTGTGCAGGAAATTCAGGTCGGCCTGGGCAGCGATTAGTCGGGCGGACATGAGCAACTCCTTTTGTCATGTGTCGCCGACAGCATAGCAACACCGCTCAGGTGTGGGAGCGGCCTTGTGTCGCGACAGGGCCGCATAGCg                                                                                                                                                    >  2:3680914/1‑150 (MQ=255)
   tACGCTGTCACTGTGCAGGAAATTCAGGTCGGCCTGGGCAGCGATTAGTCGGGCGGACATGAGCAACTCCTTTTGTCATGTGTCGCCGACAGCATAGCAACACCGCTCAGGTGTGGGAGCGGCCTTGTGTCGCGAAAGGGCCGCATAGCg                                                                                                                                                    >  2:1842824/1‑150 (MQ=255)
       cTGTCACTGTGCAGGAAATTCAGGTCGGCCTGGGCAGCGATTAGTCGGGCGGACATGAGCCACTCCTTTTGTCATGTGTCGCCGACAGCATAGCAACACCGCTCAGGTGTGGGAGCGGCCTTGTGTCGCCAAAGGGCCCCCAAGCGGccc                                                                                                                                               >  2:970912/1‑150 (MQ=255)
       cTGTCACTGTGCAGGAAATTCAGGTCGGCCTGGGCAGCGATTAGTCGGGCGGACATGAGCAACTCCTTTTGTCATGTGTCGCCGACAGCATAGCAACACCGCTCAGGTGTGGGAGCGGCCTTGTGTCGCGAAAGGGCCGCATAGCGGccc                                                                                                                                               <  1:2607697/150‑1 (MQ=255)
           cACTGTGCAGGAAATTCAGGTCGGCCTGGGCAGCGATTAGTCGGGCGGACATGAGCAACTCCTTTTGTCATGTGTCGCCGACAGCATAGCAACACCGCTCAGGTGTGGGAGCGGCCTTGTGGCGCGAAAGGGCCCCCTAGCGGCCCCGGc                                                                                                                                           >  2:1363396/1‑150 (MQ=255)
            aCTGTGCAGGAAATTCAGGTCGGCCTGGGCAGCGATTAGTCGGGCGGACATGAGCAACTCCTTTTGTCATGTGTCGCCGACAGCATAGCAACACCGCTCAGGTGTGGGAGCGGCCTTGTGTCGCGAAAGGGCCGCTATGCGGccctttcg                                                                                                                                          <  2:1303857/150‑6 (MQ=255)
            aCTGTGCAGGAAATTCAGGTCGGCCTGGGCAGCGATTAGTCGGGCGGACATGAGCAACTCCTTTTGTCATGTGTCGCCGACAGCATAGCAACACCGCTCAGGTGTGGGAGCGGCCTTGTGTCGCGAAAGGGCCGCAGAGCGGCCCCGGcc                                                                                                                                          >  2:837513/1‑149 (MQ=255)
              tgtgCAGGAAATTCAGGTCGGCCTGGGCAGCGATTAGTCGGGCGGACATGAGCAACTCCTTTTGTCATGTGTCGCCGACAGCATAGCAACACCGCTCAGGTGTGGGAGCGGCCTTGTGTCGCGCAAGGGCCGCATAGCGGCCCCCGccaa                                                                                                                                        >  2:1908978/1‑147 (MQ=255)
              tgtgCAGGAAATTCAGGTCGGCCTGGGCAGCGATTAGTCGGGCGGACATGAGCAACTCCTTTTGTCATGTGTCGCCGACAGCATAGCAACACCGCTCAGGTGTGGGAGCGGCCTTGTGTCGCGAAAGGGCCGCCAAGCGGCCCCGGccag                                                                                                                                        >  2:3455507/1‑147 (MQ=255)
                 gCAGGAAATTCAGGTCGGCCTGGGCAGCGATTAGTCGGGCGGACATGAGCAACTCCTTTTGTCATGTGTCGCCGACAGCATAGCAACACCGCTCAGGTGTGGGAGCGGCCTTGTGTCGCGAAAGGGCCGCATAGCGGCCCCGGCAATTTc                                                                                                                                     <  1:421105/150‑1 (MQ=255)
                  gaGGTAATTAAGGTCGGCCTGGGCATGGTTTATTCGGGCGGACATGAGCAACTCCTTTTGTCATGTGTCGCCGACAGCATAGCAACACCGCTCAGGTGTGGGAGCGGCCTTGTGTCGCGAAAGGGCCGCATAGCGGCCCCGGCAATTTcc                                                                                                                                    <  2:3915347/149‑1 (MQ=38)
                  cAGGAAATTCAGGTCGGCCTGGGCAGCGATTAGTCGGGCGGACATGAGCAACTCCTTTTGTCATGTGTCGCCGACAGCATAGCAACACCGCTCAGGTGTGGGAGCGGCCTTGTGTCGCGGAAGGGCCGCAAAGCGGCCCCGGCacggtcc                                                                                                                                    >  2:3147406/1‑144 (MQ=255)
                  cAGGAAATTCAGGTCGGCCTGGGCAGCGATTAGTCGGGCGGACATGAGCAACTCCTTTTGTCATGTGTCGCCGACAGCATAGCAACACCGCTCAGGTGTGGGAGCGGCCTTGTGTCGCGACAGGGCCGCCTAGCGGCCCCGGcccattgc                                                                                                                                    >  1:5470046/1‑143 (MQ=255)
                  cAGGAAATTCAGGTCGGCCTGGGCAGCGATTAGTCGGGCGGACATGAGCAACTCCTTTTGTCATGTGTCGCCGACAGCATAGCAACACCGCTCAGGTGTGGGAGCGGCCTTGTGTCGCGAAAGGGCCGCATAGCGGCCCCGGCAATTTcc                                                                                                                                    >  2:808304/1‑150 (MQ=255)
                   aGGAAATTCAGGTCGGCCTGGGCAGCGATTAGTCGGGCGGACATGAGCAACTCCTTTTGTCATGTGTCGCCGACAGCATAGCAACACCGCTCAGGTGTGGGAGCGGCCTTGTGTAGCGGAAGGGCCGCATAGCGGCCCCCGCAAattacc                                                                                                                                   >  1:2795043/1‑147 (MQ=255)
                     gAAATTCAGGTCGGCCTGGGCAGCGATTAGTCGGGCGGACATGAGCAACTCCTTTTGTCATGTGTCGCCGACAGCATAGCAACACCGCTCAGGTGTGGGAGCGGCCTTGTGTCGCGAAAGGGCCGCATAGCGGCCCCGGCACTTTCcgcg                                                                                                                                 >  2:4475749/1‑150 (MQ=255)
                     gAAATTCAGGTCGGCCTGGGCAGCGATTAGTCGGGCGGACATGAGCAACTCCTTTTGTCATGTGTCGCCGACAGCATAGCAACACCGCTCAGGTGTGGGAGCGGCCTTGTGTCGCGAAAGGGCCCCATAGCGGCCCCCGCAAattacccc                                                                                                                                 >  1:2711370/1‑145 (MQ=255)
                      aaaTTCAGGTCGGCCTGGGCAGCGATTAGTCGGGCGGACATGAGCAACTCCTTTTGTCATGTGTCGCCGACAGCATAGCAACACCGCTCAGGTGTGGGAGCGGCCTTGTGTCGCGAAAGGGCCGCATAGCGGCCCCGGCAATTTCCGCgt                                                                                                                                <  2:1575377/150‑1 (MQ=255)
                        aTTCAGGTCGGCCTGGGCAGCGATTAGTCGGGCGGACATGAGCGACTCCTTTTGTCATGTGTCGCCGACAGCATAGCAACACCGCTCAGGTGTGGGAGCGGCCTTGTGTCGCGAAAGGGCCGCATAGCGGCCCCGGCAATTTCCGCgtgt                                                                                                                              <  1:3705129/150‑1 (MQ=255)
                          tCAGGTCGGCCTGGGCAGCGATTAGTCGGGCGGACGTGAGCAACTCCTTTTGTCATGTGTCGCCGACAGCATAGCAACACCGCTCAGGTGTGGGAGCGGCCTTGTGTCGCGAAAGGGCCGCATAGCGGCCCCGGCAATTTCCGCGTGTAt                                                                                                                            <  2:1636356/150‑1 (MQ=255)
                           cAGGTCGGCCTGGGCAGCGATTAGTCGGGCGGACATGAGCAACTCCTTTTGTCATGTGTCGCCGACAGCATAGCAACACCGCTCAGGTGTGGGAGCGGCCTTGTGTCGCGAAAGGGCCGCATAGCGGCCCCGGCAATTTCCGCGTGTATg                                                                                                                           <  2:2627094/150‑1 (MQ=255)
                                   ccTGGGCAGCGATTAGTCGGGCGGACATGAGCAACTCCTTTTGTCATGTGTCGCCGACAGCATAGCAACACCGCTCAGGTGTGGGAGCGGCCTTGTGTCGCGAAAGGGCCGCATAGCGGCCCCGGCAATTTCCGCGTGTATGCAACAAAc                                                                                                                   <  1:3244978/150‑1 (MQ=255)
                                     tGGGCAGCGATTAGTCGGGCGGACATGAGCAACTCCTTTTGTCATGTGTCGCCGACAGCATAGCAACACCGCTCAGGTGTGGGAGCGGCCTTGTGTCGCGAAAGGGCCGCATAGCGGCCCCGGCAATTTCCGCGTGTATGCAACAAACCt                                                                                                                 <  1:5446048/150‑1 (MQ=255)
                                        gCAGCGATTAGTCGGGCGGACATGAGCAACTCCTTTTGTCTTGTGTCGCCGACAGCATAGCAACACCGCTCAGGTGTGGGAGCGGCCTTGTGTCGCGAATGGGCCGCATAGCGGCCCCGGCAATTTCCGCGTGTATGCAACAAACCTggg                                                                                                              <  1:3509804/150‑1 (MQ=255)
                                          aGCGATTAGTCGGGCGGACATGAGCAACTCCTTTTGTCATGTGTCGCCGACAGCATAGCAACACCGCTCAGGTGTGGGAGCGGCCTTGTGTCGCGAAAGGGCCGCCTAGCGGCCCCCGCCAATTCCGCGTGTATGCAACAAAACGGGGgc                                                                                                            >  2:2475123/1‑150 (MQ=39)
                                                   tCGGGCGGACATGAGCAACTCCTTTTGTCATGTGTCGCCGACAGCATAGCAACACCGCTCAGGTGTGGGAGCGGCCTTGTGTCGCGAAAGGGCCGCATAGCGGCCCCGGCAATTTCCGCGTGTATGCAACAAACCTGGGGCTGCTATGCa                                                                                                   <  1:125500/150‑1 (MQ=255)
                                                   tCGGGCGGACATGAGCAACTCCTTTTGTCATGTGTCGCCGACAGCATAGCAACACCGCTCAGGTGTGGGAGCGGCCTTGTGTCGCGAAAGGGCCGCATAGCGGCCCCGGCAATTTCCGCGTGTATGCAACAAACCTGGGGCTGCTATGCa                                                                                                   <  2:348747/150‑1 (MQ=255)
                                                   tCGGGCGGACATGAGCAACTCCTTTTGTCATGTGTCGCCGACAGCATAGCAACACCGCTCAGGTGTGGGAGCGGCCTTGTGTCGCGAAAGGGCCGCATAGCGGCCCCGGCAATTTCCGCGTGTATGCAACAAACCTGGGGCTGCTATGCa                                                                                                   <  1:2152578/150‑1 (MQ=255)
                                                     gggCGGACATGAGCAACTCCTTTTGTCATGTGTCGCCGACAGCATAGCAACACCGCTCAGGTGTGGGAGCGGCCTTGTGTGGCGAAAAGGCCGCATAGCGGCCCCCGCCAATTACGCCGGGATGGAACAAACCTGGGGCGGCGATGCAGc                                                                                                 >  2:2559116/1‑150 (MQ=255)
                                                     atgatctggcggctcagctgcgcgccagcctgggccagttccgggcctgaagcctgcgcggtcctttgtgGGCCTTGTGTCGCGAAAGGGCCGCATAGCGGCCCCGGCAATTTCCGCGTGTATGCAACAAACCTGGGGCTGCTATGCAGc                                                                                                 >  1:1878004‑M1/71‑150 (MQ=255)
                                                     atgatctggcggctcagctgcgcgccagcctgggccagttccgggcctgaagcctgcgcggtcctttgtgGGCCTTGTGTCGCGAAAGGGCCGCATAGCGGCCCCGGCAATTTCCGCGTGTATGCAACAAACCTGGGGCTGCTATGCAGc                                                                                                 >  2:4300409‑M1/71‑150 (MQ=255)
                                                         ggACATGAGCAACTCCTTTTGTCATGTGTCGCCGACAGCATAGCAACACCGCTCAGGTGTGGGAGCGGCCTTGTGTCGCGAAAGGGCCGCATAGCGGCCCCGGCAATTTCCGCGTGTATGCAACAAACCTGGGGCTGCATAGCAGcccca                                                                                             <  2:5507657/150‑3 (MQ=255)
                                                          gACATGAGCAACTCCTTTTGTCATGTGTCGCCGACAGCATAGCAACACCGCTCAGGTGTGGGAGCGGCCTTGTGTCGCGAAAGGGCCGCATAGCGGCCCCGGCAATTTCCGCGTGTATGCAACAAACCTGGGGCTGCATAGCAGccccag                                                                                            <  2:4291113/150‑4 (MQ=255)
                                                           tggcggctcagctgcgcgccagcctgggccagttccgggcctgaagcctgcgcggtcctttgtgGGCCTTGTGTCGCGAAAGGGCCGCATAGCGGCCCCGGCAATTTCCGCGTGTATGCAACAAACCTGGGGCTGCTATGCAGCCCTTTc                                                                                           >  1:214458‑M1/65‑150 (MQ=255)
                                                               gAGCAACTCCTTTTGTCATGTGTCGCCGACAGCATAGCAACACCGCTCAGGTGTGGGTGCGGCCTTGTGTCGCGAAAGGGCCGCATAGCGGCCCCGGCAATTTCCGCGTGTATGCAACAAACCTGGGGCTGCTATGCAGCCCTTTCGCGa                                                                                       <  2:2798507/150‑1 (MQ=255)
                                                               gAGCAACTCCTTTTGTCATGTGTCGCCGACAGCATAGCAACACCGCTCAGGTGTGGGAGCGGCCTTGTGTCGCGAAAGTGCCGCAAAGCGGCCCCCGCAAATTCCGCCTGTATGCAACAAAACTGGGGCTGCTAAGCAGCCCATTCGCGa                                                                                       >  2:3044119/1‑150 (MQ=38)
                                                               gAGCAACTCCTTTTGTCATGTGTCGCCGACAGCATAGCAACACCGCTCAGGTGTGGGAGCGGCCTTGTGTCGCGAAAGGGCCGCATAGCGGCCCCGGCAATTTCCGCGTGTATGCAACAAACCTGGGGCTGCTATGCAGCCCTTTCGCGa                                                                                       <  1:4072943/150‑1 (MQ=255)
                                                               gAGCAACTCCTTTTGTCATGTGTCGCCGACAGCATAGCAACACCGCTCAGGTGTGGGAGCGGCCTTGTGTCGCGAAAGGGCCGCATAGCGGCCCCGGCAATTTCCGCGTGTATGCAACAAACCTGGGGCTGCTATGCAGCCCTTTCGCGa                                                                                       <  2:4823136/150‑1 (MQ=255)
                                                                aGCAACTCCTTTTGTCATGTGTCGCCGACAGCATAGCAACACCGCTCAGGTGTGGGAGCGGCCTTGTGTCGCGGGAGGGCCGCGTTTCGGCCCCGGCACTTTCCGGGTTTATGGCACAAACCTGGGGGTGGTATGGAGGGCTTTggcggg                                                                                      >  1:849042/1‑148 (MQ=255)
                                                                  cAACTCCTTTTGTCATGTGTCGCCGACAGCATAGCAACACCGCTCAGGTGTGGGAGCGGCCTTGTGTGGGGTAAGGGGCGCATAGGGGGCCCGGGAATTTTCGGGTGTTTGGGACAAAGCTGGGGGTGGTTTTGAGGCCTTTAGGGacac                                                                                    >  2:1511851/1‑150 (MQ=255)
                                                                  cAACTCCTTTTGTCATGTGTCGCCGACAGCATAGCAACACCGCTCAGGTGTGGGAGCGGCCTTGTGTCGCGAAAGGGCCGCATAGCGGCCCCGGCAATTTCCGCGTGTATGCAACAAACCTGGGGCTGCTATGCAGCCCTTTCGCGacac                                                                                    <  2:261550/150‑1 (MQ=255)
                                                                    aCTCCTTTTGTCATGTGTCGCCGACAGCATAGCAACACCGCTCAGGTGTGGGAGCGGCCTTGTTTCGCGAAAGGGCCGCATAGCCGCCCCCGCAAATTCCGCGTGTATTCAACAAAACAGGGGCTGCTATGCAGCCCAGTCGCGACACaa                                                                                  >  2:2404687/1‑150 (MQ=255)
                                                                    aCTCCTTTTGTCATGTGTCGCCGACAGCATAGCAACACCGCTCAGGTGTGGGAGCGGCCTTGTGTCGCGAAAGGGCCGCATAGCGGCCCCGGCAATTTCCGCGTGTATGCAACAAACCTGGGGCTGCTATGCAGCCCTTTCGCGACACaa                                                                                  <  1:3071341/150‑1 (MQ=255)
                                                                      tCCTTTTGTCATGTGTCGCCGACAGCATAGCAACACCGCTCAGGTGTGGGAGCGGCCTTGTGTCGCGAAAGGGCCGCAAAGCGGCCCCCGCAATTTACGCGGGTATGCAACAAACCTGGGGCTGCAATGCAGCCCTTTCGCGACACAAgg                                                                                >  2:2301900/1‑150 (MQ=37)
                                                                      tCCTTTTGTCATGTGTCGCCGACAGCATAGCAACACCGCTCAGGTGTGGGAGCGGCCGTGTGTGGCGAAAGGGCCGCCTAGCGGCCCCCGCCAATTGCCCGGGGATGCAACAAACCAGGGGGGGCGATTCAGCCCGGTAGCGCCACaacg                                                                                >  2:1116397/1‑148 (MQ=255)
                                                                       ccTTTTGTCATGTGTCGCCGACAGCATAGCAACACCGCTCTGGTGTGGGAGCGGCCTTGTGTCGCGAAGGGGCCGCATCGCGGCCCCGGCAATTTCCGCGTGTATGCAACAAACCTGGGGCTGCTATGCAGCCCTTTCGCGACACAAGGc                                                                               <  1:725246/150‑1 (MQ=255)
                                                                       ccTTTTGTCATGTGTCGCCGACAGAATAGCAACACCGCTCTGGTGTGGGAGCGGCTTTGTGTCGCGAAAGGGCCGCATAGCGGCCCCGGCAATTTCCGCTTTTATGCAACAAACCTGGGGCTGCTTTGCAGCCCTTTCGCGACACAAGGc                                                                               <  1:1836924/150‑1 (MQ=37)
                                                                         ttttGTCATGTGTCGCCGACAGCATAGCAACACCGCTCAGGTGTGGGAGCGGCCTTGTGTCGCGAAAGGGCCGCATAGCGGcccc                                                                                                                                              <  1:5246920/85‑1 (MQ=255)
                                                                          tttGTCATGTGTCGCCGACAGCATAGCAACACCGCTCAGGTGTGGGAGCGGCCTTGTGTGGCGAAAGGGCCCCCAAGCGGCCCCCGCCAATTCCGCGGGTGTGCAACAAACCCGGGGCGGCGATGCCGCCCGTTCGCGACACACAGCCca                                                                            >  1:3749648/1‑150 (MQ=255)
                                                                              tCATGTGTCGCCGACAGCATAGCAACACCGCTCAGGTGTGGGAGCGGCCTTGTGTCGCGAAAGGGCCGCATAGCGGCCCCGGCCAGTTCCGCGGGTATGCAACAAACCTGGGGCTGCGATGCAGCCCGTTCGCGACACAAGGCCcacaca                                                                        >  2:5151470/1‑148 (MQ=38)
                                                                               cttgtgTCGCCGACAGCATAGCAACACGGCTCTGGTGTGGGAGGGGCGTTGTGTCGCGAATGGGCCGCATAGCGGCCCCGGCAATTTCCGCGTGTATGCAACAAACGTGGGGCTGCTATGCAGCCCTTTCGCGACACAAGGCCCACAAAg                                                                       <  1:788911/148‑1 (MQ=38)
                                                                               cATGTGTCGCCGACAGCATAGCAACACCGCTCAGGTGTGGGAGCGGCCTTGTGTCGCGAATGGGCCGCATCGCGGCCCCGGCAATTTCCGCGTGTATGCAACAAACCTGGGGCTGCTATGCAGCCCTTTCGCGACACAAGGCCCACAAAg                                                                       <  2:2082159/150‑1 (MQ=255)
                                                                               cATGTGTCGCCGACAGAATAGAAACCCCGCTCTGGTGTGGGAGGGGCCTTGGGTCGCCAACGGGCCGAATTGCGGCCCCGGCAATTTCCGCGTGTATGCAACAAACCTGGGGCTGCTATGCAGCCCTTTCGCGACACAAGGCCCACAAAg                                                                       <  1:5304877/150‑1 (MQ=38)
                                                                                aTGTGTCGCCGACAGCATAGCAACACCGCTCAGGTGTGGGAGCGGCCTTGTGTGGCGAAAGTGGCGCATATCGGCCCCGGCAATTTGAGGGTGTTTTGAACAAAACTGGGGGTGCTATGAATGCCTTTGGCGAAACAAGGCCCACAAAgg                                                                      >  2:5612162/1‑150 (MQ=255)
                                                                                aTGTGTCGCCGACAGCATAGCAACACCGCTCAGGTGTGGGAGCGGCCTTGTGTCGCGCAAGGGCCGCATAGCGGCCCCGGCCAATTCCGCGGGTATGCAACAAACCAGGGGCTGCGATGCAGCCCGTTGGCGACACAAGGCCCACAAAgg                                                                      >  2:4939371/1‑150 (MQ=38)
                                                                                aTGTGTCGCCGACAGCATAGCAACACCGCTCAGGTGTGGGAGCGGCCTTGGGGCGCGCAAGGGCCCCAAAGCGGCCCCCGCCAGTTCCCCGTGTATGCAACCAACCGGGGGCGGCGCTTCCGCCCGGTCGCGACCCAACGCCcacccagg                                                                      >  2:334966/1‑145 (MQ=255)
                                                                                   tgtCGCCGACAGCATAGCAACACCGCTCAGGTGTGGGAGCGGCCTTGTGTCGCGAAAGGGCCGCATAGCGGCCCCGGCAATTTCCGCGTGTATGCAACAAACCTGGGGCTGCTATGCAGCCCTTTCGCGACACAAGGCCCACAAAGGAcc                                                                   <  1:1818153/150‑1 (MQ=255)
                                                                                        ccGACAGCATAGCAACACCGCTCAGGTGTGGGAGCGGCCTTGTGTCGCGAAAGGGCCGCATAGCGGCCCCGGCAATTTCCGCGTGTATGCAACAAACCTGGGGCTGCTATGCAGCCCTTTCGCGACACAAGGCCCACAAAGGACCGCGCa                                                              <  2:4427033/150‑1 (MQ=255)
                                                                                         cGACAGCATAGCAACACCGCTCAGGTGTGGGAGCGGCCTTGTGTCGCGAAAGGGCCGCATAGCGGCCCCGGCAATTTCCGCGTGTATGCAACAAACCTGGGGCTGCTATGCAGCCCTTTCGCGACACAAGGCCCACAAAGGACCGCGCAg                                                             <  2:4804749/150‑1 (MQ=255)
                                                                                          gACAGCATAGCAACACCGCTCAGGTGTGGGAGCGGCCTTGTGTGGCGGAAGGGGCGCCAAGCGGCCCCCGCCAATTTCCCGGGTGTGGAACAAAACGGGGGCGGCGATGGACCCCCTTCGCGGCCCCAGGCCCACAAAAGACCccgcccg                                                            >  2:2515081/1‑147 (MQ=255)
                                                                                              gCATAGCACCCCCGCTCGGGTGTGGGCGCGGCGTTGTGTCGCGAACGGGCCGCATAGCGGCCCCGGCAATTTCCGCGTGTATGCAACAAACCTGGGGCTGCTATGCAGCCCTTTCGCGACACAAGGCCCACAAAGGACCGCGCAGGCTTc                                                        <  1:3367689/150‑1 (MQ=38)
                                                                                              gCATAGCAACACCGCTCAGGTGTGGGAGCGGCCTTGTGTCGCGTAAGGGCCGCATAGCGGCCCCCGCAAATTGCGCGTGTATGCCACAAAACAGGGGCTGCGATGCAGCCCATTCGCGACACAAGGCCCACCAACGACCGCGCCGGcgga                                                        >  2:4904586/1‑147 (MQ=37)
                                                                                              gCATAGCAACACCGCTCAGGTGTGGGAGCGGCCTTGTGTCGCGAAAGGGCCGCATAGCGGCCCCGGCAAATTCCCCGTGTATGCAACAAACCTGGGGCTGCTATGCACCCCTTTCGCGACACAAGGCCCCCAAAGGACCGCCCAGGCTTc                                                        >  2:3594626/1‑150 (MQ=38)
                                                                                               cATAGCAACACCGCTCAGGTGTGGGAGCGGCCTTGTGTCGCGGAATGGCCGCCTAGCGGCCCCCGCCCATTCCCCGGGTATGCCACAAACCGGGGGCTGCTATGCAGCCCGTTAGCCCCACAACGCCCCCAAACGACCCCGCCGGCGTCa                                                       >  2:3091344/1‑150 (MQ=255)
                                                                                                aTAGCACCACCGCTCCGGTGTGGGAGCGGCCTTGTGTCGCGAACGGGCCGCATAGCGGCCCCGGCAATTTCCGCGTGTATGCAACAAACCTGGGGCTGCTATGCAGCCCTTTCGCGACACAAGGCCCACAAAGGACCGCGCAGGCTTCAg                                                      <  2:5353449/150‑1 (MQ=255)
                                                                                                 tAGCAACACCGCTCAGGTGTGGGAGCGGCCTTGTGTCGCGAAAGGGCCGCATAGCGGCCCCGGCAATTTCCGCGTGTATGCCACAAACCGGGGGCTGCTCTGCAGCCCGTTCGCGACACAAGGCCCACAAAGGACCGCGCCCGCTTCAgg                                                     >  2:1278323/1‑150 (MQ=38)
                                                                                                 aaGCAACACCGCTCAGGTGTGGGAGCGGCCTTGTGTCTCGAAAGTGACGCATAGCGGCCCCCGCAAATTACGCCTGTAATCAACAAAACTGGGGCTGCTAATCAACCCATGCGCGACACAAGGACCACAAATGACCCCGCAGGcgactgg                                                     >  2:4665765/2‑144 (MQ=255)
                                                                                                     aaCACCGCTCAGGTGTGGGAGCGGCCTTGTGTCGCGAAAGGGCCGCATAGCGGCCCCGGCAATTTCCGCGTGTATGCAACAAACCTGGGGCTGCTATGCAGCCCTTTCGCGACACAAGGCCCACAAAGGACCGCGCAGGCTTCAGGCCCg                                                 <  1:5601943/150‑1 (MQ=255)
                                                                                                        acCGCTCAGGTGTGGGAGCGGCCTTGTGTCGCGAAAGGGCCGCATAGCGGCCCCGGCCATTTACGCGTGTATGCAACAAACCGGGGGCTGCTATGCAGCCCATTCGCGACACAAGGCCCACAAAGGACCGCGCAGGCTTCAGGCCCGGaa                                              >  2:893011/1‑150 (MQ=255)
                                                                                                            cTCAGGTGTGGGAGCGGCCTTGTGTCGCGAAAGGGCCGCATAGCGGCCCCGGCAATTTCCGCGTGTATGCAACAAACCTGGGGCTGCTATGCAGCCCTTTCGCGACACAAGGCCCACAAAGGACCGCGCAGGCTTCAGGCCCGGAACTgg                                          <  2:1760604/150‑1 (MQ=255)
                                                                                                             tCAGGTGTGGGAGCGGCCTTGTGTCGCGAAAGGGCCGCATAGCGGCCCCGGCAATTTCCGCGTGTATGCAACAAACCTGGGGCTGCTAt                                                                                                      >  1:1355284/1‑89 (MQ=255)
                                                                                                                 gtgtGGGAGCGGCCTTGTGTCGCGACAGGGCCGCATAGCGGCCCCGGCAATTTCCGCGTGTATGCAACCAACCTGGGGCTGCTATGCAGCCCTTTCGCGACACAAGGCCCACAAAGGACCGCGCAGGCTTCAGGCCCGGAACTGGCCCAg                                     <  2:4685274/150‑1 (MQ=255)
                                                                                                                             CCTTGTGTCGCGAAAGGGCCGCATAGCGGCCCCGGCAATTTCCGCGTGTATGCAACAAACCTGGGGCTGCTATGCAGCCCTTTCGCGACACAAGGCCCACAAAGGACCGCGCAGGCTTCAGGCCCGGAACTGGCCCAGGCTGGCGCGCAg                         <  1:4498613‑M1/150‑1 (MQ=255)
                                                                                                                             CCTTGTGTCGCGAAAGGGCCGCATAGCGGCCCCGGCAATTTCCGCGTGTATGCAACAAACCTGGGGCTGCTATGCAGCCCTTTCGCGACACAAGGCCCACAAAGGACCGCGCAGGCTTCAGGCCCGGAACTGGCCCAGGCTGGCGCGCAg                         <  2:5105975‑M1/150‑1 (MQ=255)
                                                                                                                                 GTGTCGCGAAAGGGCCGCATAGCGGCCCCGGCAATTTCCGCGTGTATGCAACAAACCTGGGGCTGCTATGCAGCCCTTTCGCGACACAAGGCCCACAAAGGACCGCGCAGGCTTCAGGCCCGGAACTGGCCCAGGCTGGCGCGCAGCTGa                     >  2:5404847‑M1/1‑150 (MQ=255)
                                                                                                                                 GTGTCGCGAAAGGGCCGCATAGCGGCCCCGGCAATTTCCGCGTGTATGCAACAAACCTGGGGCTGCTATGCAGCCCTTTCGCGACACAAGGCCCACAAAGGACCGCGCAGGCTTCAGGCCCGGAACTGGCCCAGGCTGGCGCGCAGCTGa                     <  1:2237516‑M1/150‑1 (MQ=255)
                                                                                                                                 GTGTCGCGAAAGGGCCGCATAGCGGCCCCGGCAATTTCCGCGTGTATGCAACAAACCTGGGGCTGCTATGCAGCCCTTTCGCGACACAAGGCCCACAAAGGACCGCGCAGGCTTCAGGCCCGGAACTGGCCCAGGCTGGCGCGCAGCTGa                     >  2:1686440‑M1/1‑150 (MQ=255)
                                                                                                                                   GTCGCGAAAGGGGCGCATAGCGGCCCCGGCAATTTCCGCGTGTATGCAACAAACCTGGGGCTGCTATGCAGCCCTTTCGCGACACAAGGCCCACAAAGGACCGCGCAGGCTTCTGGCCCGGAACTGGCCCAGGCTGGCGCGCAGCTGAgc                   <  2:3794616‑M1/150‑1 (MQ=255)
                                                                                                                                    TCGCGAAAGGGCCGCATAGCGGCCCCGGCAATTTCCGCGTGTATGCAACAAACCTGGGGCTGCTATGCAGCCCTTTCGCGACACAAGGCCCACAAAGGACCGCGCAGGCTTCAGGCCCGGAACTGGCCCAGGCTGGCGCGCAGCTGAgcc                  <  1:4061887‑M1/150‑1 (MQ=255)
                                                                                                                                    TCGCGAAAGGGCCGCATAGCGGCCCCGGCAATTTCCGCGTGTATGCAACAAACCTGGGGCTGCTATGCAGCCCTTTCGCGACACAAGGCCCACAAAGGACCGCGCAGGCTTCAGGCCCGGAACTGGCCCAGGCTGGCGCGCAGCTGAgcc                  <  1:2485379‑M1/150‑1 (MQ=255)
                                                                                                                                      GCGAAAGGGCCGCATAGCGGCCCCGGCAATTTCCGCGTGTATGCAACAAACCTGGGGCTGCTATGCAGCCCTTTCGCGACACAAGGCCCACAAAGGACCGCGCAGGCTTCAGGCCCGGAACTGGCCCAGGCTGGCGCGCAGCTGAgccgc                >  2:782371‑M1/1‑150 (MQ=255)
                                                                                                                                      GCGAAAGGGCCGCATAGCGGCCCCGGCAATTTCCGCGTGTATGCAACAAACCTGGGGCTGCTATGCAGCCCTTTCGCGACACAAGGCCCACAAAGGACCGCGCAGGCTTCAGGCCCGGAACTGGCCCAGGCTGGCGCGCAGCTGAgccgc                >  2:5293604‑M1/1‑150 (MQ=255)
                                                                                                                                      GCGAAAGGGCCGCATAGCGGCCCCGGCAATTTCCGCGTGTATGCAACAAACCTGGGGCTGCTATGCAGCCCTTTCGCGACACAAGGCCCACAAAGGACCGCGCAGGCTTCAGGCCCGGAACTGGCCCAGGCTGGCGCGCAGCTGAgccgc                <  2:4327584‑M1/150‑1 (MQ=255)
                                                                                                                                        gaccGGGCCGCATAGCGGCCCCGGCAATTTCCGCGTGTATGCAACACACCTGGGGCTGCTATGCAGCCCTTTCGCGACACAAGGCCCACAAAGGACCGCGCAGGCTTCAGGCCCGGAACTGGCCCAGGCTGGCGCGCAGCTGAGCCGCCa              >  1:5253850‑M1/5‑150 (MQ=255)
                                                                                                                                        gacaGGGCCGCATCGCGGCCCCGGCAATTTCCGCGTGTATGCAACAAACCTGGGGCTGCTATGCAGCCCTTTCGCGACACAAGGCCCACAAAGGACCGCGCAGGCTTCAGGCCCGGAACTGGCCCAGGCTGGCGCGCAGCTGAGCCGCCa              <  1:4178918/147‑1 (MQ=255)
                                                                                                                                        GAACGGGCCTCATCGCGGCCCCGGCAATTTCCGCGTGTATGCAACAAACCTGGGGCTGCTATGCAGCCCTTTCGCGACACAAGGCCCACAAAGGACCGCGCAGGCTTCAGGCCCGGAACTGGCCCAGGCTGGCGCGCAGCTGAGCCGCCa              >  1:4872932‑M1/1‑150 (MQ=255)
                                                                                                                                        GAAAGGGCCGCCTCGCGGCCCCGGCAATTACCGCGTGTATGCAACAACCCGGGGGCTGCTATGCAGCCCTTTCGCGAAACAAGGCCCACAAAGGACCGCGCAGGCTTCAGGCCCGGAACTGGCCCAGGCTGGCGCGCAGCTGAGCCGCCa              >  1:970912‑M1/1‑150 (MQ=255)
                                                                                                                                        GAAAGGGCCGCATAGCGGCCCCGGCAATTTCCGCGTGTATGCAACAAACCTGGGGCTGCTATGCAGCCCTTTCGCGACACAAGGCCCACAAAGGACCGCGCAGGCTTCAGGCCCGGAACTGGCCCAGGCTGGCGCGCAGCTGAGCCGCCa              <  1:1703534‑M1/150‑1 (MQ=255)
                                                                                                                                        GAAAGGGCCGCATAGCGGCCCCGGCAATTTCCGCGTGTATGCAACAAACCTGGGGCTGCTATGCAGCCCTTTCGCGACACAAGGCCCACAAAGGACCGCGCAGGCTTCAGGCCCGGAACTGGCCCAGGCTGGCGCGCAGCTGAGCCGCCa              >  1:2559172‑M1/1‑150 (MQ=255)
                                                                                                                                         AAAGGGCCGCATAGCGGCCCCGGCAATTTCCGCGTGTATGCAACAAACCTGGGGCTGCTATGCAGCCCTTTCGCGACACAAGGCCCACAAAGGACCGCGCAGGCTTCAGGCCCGGAACTGGCCCAGGCTGGCGCGCAGCTGAGCCGCCAg             <  1:1483454‑M1/150‑1 (MQ=255)
                                                                                                                                           AGGGCCGCATAGCGGCCCCGGCAATTTCCGCGTGTATGCAACAAACCTGGGGCTGCTATGCAGCCCTTTCGCGACACAAGGCCCACAAAGGACCGCGCAGGCTTCAGGCCCGGAACTGGCCCAGGCTGGCGCGCAGCTGAGCCGCCAGat           <  2:4589101‑M1/150‑1 (MQ=255)
                                                                                                                                           AGGGCCGCATAGCGGCCCCGGCAATTTCCGCGTGTATGCAACAAACCTGGGGCTGCTATGCAGCCCTTTCGCGACACAAGGCCCACAAAGGACCGCGCAGGCTTCAGGCCCGGAACTGGCCCAGGCTGGCGCGCAGCTGAGCCGCCAGat           <  2:4215967‑M1/150‑1 (MQ=255)
                                                                                                                                            GGGCCGCTATGCGGCCCCGGCAATTTCCGCGTGTATGCAACAAACCTGGGGCTGCTATGCAGCCCTTTCGCGACACAAGGCCCACAAAGGACCGCGCAGGCTTCAGGCCCGGAACTGGCCCAGGCTGGCGCGCAGCTGAGCCGCCAGatc          <  1:404098‑M1/150‑1 (MQ=255)
                                                                                                                                             GGCCGCATAGCGGCCCCGGCAATTTCCGCGTGTATGCAACAAACCTGGGGCTGCTATGCAGCCCTTTCGCGACACAAGGCCCACAAAGGACCGCGCAGGCTTCAGGCCCGGAACTGGCCCAGGCTGGCGCGCAGCTGAGTCGCCAGatca         <  1:1130272‑M1/150‑1 (MQ=255)
                                                                                                                                                cGCATAGCGGCCCCGGCAATTTCCGCGTGTATGCAACAAACCTGGGGCTGCTATGCAGCCCTTTCGCGACACAAGGCCCACAAAGGACCGCGCAGGCTTCAGGCCCGGAACTGGc                                         <  1:3524330/115‑1 (MQ=255)
                                                                                                                                                cGCATAGCGGCCCCGGCAATTTCCGCGTGTATGCAACAAACCTGGGGCTGCTATGCAGCCCTTTCGCGACACAAGGCCCACAAAGGACCGCGCAGGCTTCAGGCCCGGAACTGGCCCAGGCTGGCGCGCAGCTGAGCCGCCAGATCATcc      >  2:2840192/1‑150 (MQ=255)
                                                                                                                                                 gCATAGCGGCCCCGGCAATTTCCGCGTGTATGCAACAAACCTGGGGCTGCTATGCAGCCCTTTCGCGACACAAGGCCCACAAAGGACCGCGCAGGCTTCAGGCCCGGAACTGGCCCAGGCTGGCGCGCAGCTGAGCCGCCAGATCATCCa     >  2:5458497/1‑150 (MQ=255)
                                                                                                                                                    tAGCGGCCCCGGCAATTTCCGCGTGTATGCAACAAACCTGGGGCTGCTATGCAGCCCTTTCGCGACACAAGGCCCACAAAGGACCGCGCAGGCTTCAGGCCCGGAACTGGCCCAGGCTGGCGCGCAGCTGAGCCGCCAGATCATCCAGCa  <  2:4435789/150‑1 (MQ=255)
                                                                                                                                                    tAGCGGCCCCGGCAATTTCCGCGTGTATGCAACAAACCTGGGGCTGCTATGCAGCCCTTTCGCGACACAAGGCCCACAAAGGACCGCGCAGGCTTCAGGCCCGGAACTGGCCCAGGCTGGCGCGCAGCTGAGCCGCCAGATCATCCAGCa  <  2:1599745/150‑1 (MQ=255)
                                                                                                                                                    tAGCGGCCCCGGCAATTTCCGCGTGTATGCAACAAACCTGGGGCTGCTATGCAGCCCTTTCGCGACACAAGGCCCACAAAGGACCGCGCAGGCTTCAGGCCCGGAACTGGCCCAGGCTGGCGCGCAGCTGAGCCGCCAGATCATCCAGCa  <  2:645836/150‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                      
CCGTACGCTGTCACTGTGCAGGAAATTCAGGTCGGCCTGGGCAGCGATTAGTCGGGCGGACATGAGCAACTCCTTTTGTCATGTGTCGCCGACAGCATAGCAACACCGCTCAGGTGTGGGAGCGGCCTTGTGTCGCGAAAGGGCCGCATAGCGGCCCCGGCAATTTCCGCGTGTATGCAACAAACCTGGGGCTGCTATGCAGCCCTTTCGCGACACAAGGCCCACAAAGGACCGCGCAGGCTTCAGGCCCGGAACTGGCCCAGGCTGGCGCGCAGCTGAGCCGCCAGATCATCCAGCA  >  NC_002947/5560453‑5560750

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 11 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: