Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,560,602 A→T 5.3% intergenic (‑81/+93) 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,6020AT 5.3% 630.3 / 19.3 222intergenic (‑81/+93)PP_4887/PP_4888hypothetical protein/methyl‑accepting chemotaxis transducer
Reads supporting (aligned to +/- strand):  ref base A (85/116);  new base T (5/7);  total (90/132)
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.89e-02

Only 100 of 223 total aligned reads displayed.

CCGTACGCTGTCACTGTGCAGGAAATTCAGGTCGGCCTGGGCAGCGATTAGTCGGGCGGACATGAGCAACTCCTTTTGTCATGTGTCGCCGACAGCATAGCAACACCGCTCAGGTGTGGGAGCGGCCTTGTGTCGCGAAAGGGCCGCATAGCGGCCCCGGCAATTTCCGCGTGTATGCAACAAACCTGGGGCTGCTATGCAGCCCTTTCGCGACACAAGGCCCACAAAGGACCGCGCAGGCTTCAGGCCCGGAACTGGCCCAGGCTGGCGCGCAGCTGAGCCGCCAGATCATCCAGC  >  NC_002947/5560453‑5560749
                                                                                                                                                      |                                                                                                                                                     
ccGTACGCTGTCACTGTGCAGGAAATTCAGGTCGGCCTGGGCAGCGATTAGTCGGGCGGACATGAGCAACTCCTTTTGTCATGTGTCGCCGACAGCATAGCAACACCGCTCAGGTGTGGGAGCGGCCTTGTGTCGCGAAAGGGCCGCATa                                                                                                                                                       <  1:3057795/150‑1 (MQ=255)
 cGTACGCTGTCACTGTGCAGGAAATTCAGGTCGGCCTGGGCAGCGATTAGTCGGGCGGACATGAGCAACTCCTTTTGTCATGTGTCGCCGACAGCATAGCAACACCGCTCAGGTGTGGGAGCGGCCTTGTGTCGCGAAAGGGCCGCATAg                                                                                                                                                     <  1:6459111/150‑1 (MQ=255)
 cGTACGCTGTCACTGTGCAGGAAATTCAGGTCGGCCTGGGCAGCGATTAGTCGGGCGGACATGAGCAACTCCTTTTGTCATGTGTCGCCGACAGCATAGCAACACCGCTCAGGTGTGGGAGCGGCCTTGTGTCGCGAAAGGGCCGCATAg                                                                                                                                                     >  2:5129665/1‑150 (MQ=255)
   tACGCTGTCACTGTGCAGGAAATTCAGGTCGGCCTGGGCAGCGATTAGTCGGGCGGACATGAGCAACTCCTTTTGTCATGTGTCGCCGACAGCATAGCAACACCGCTCAGGTGTGGGAGCGGCCTTGTGTCGCGAAAGGGCCGCATAGCg                                                                                                                                                   <  1:5895954/150‑1 (MQ=255)
     cGCTGTCACTGTTCAGGAAATTCAGGTCGGCCTGGGCAGCGATTAGTCGGGCGGACATGAGCAACTCCTTTTGTCATGTGTCGCCGACAGCATAGCAACACCGCTCAGGTGTGGGAGCGGCCTTGTGTCGCGAAAGGGCCGCATAGCGGc                                                                                                                                                <  2:4774784/150‑1 (MQ=255)
      gCTGTCACTGTGCAGGAAATTCAGGTCGGCCTGGGCAGCGATTAGTCGGGCGGACATGAGCAACTCCTTTTGTCATGTGTCGCCGACAGCATAGCAACACCGCTCAGGTGTGGGAGCGGCCTTGTGTCGCGAAAGGGCCGCATATCGGcc                                                                                                                                               >  2:356424/1‑150 (MQ=255)
        tGTCACTGTGCAGGAAATTCAGGTCGGCCTGGGCAGCGATTAGTCGGGCGGACATGAGCAACTCCTTTTGTCATGTGTCGCCGACAGCATAGCAACACCGCTCAGGTGTGGGAGCGGCCTTGTGTCGCGAAAGGGCCGCATAGCGGcccc                                                                                                                                             <  1:1976135/150‑1 (MQ=255)
          tCACTGTGCAGGAAATTCAGGTCGGCCTGGGCAGCGATTAGTCGGGCGGACATGAGCAACTCCTTTTGTCATGTGTCGCCGACAGCATAGCAACACCGCTCAGGTGTGGGAGCGGCCTTGTGTCGCGAAAGGGCCGCATAGCGGCCCCgg                                                                                                                                           <  1:6235240/150‑1 (MQ=255)
           cACTGTGCAGGAAATTCAGGTCGGCCTGGGCAGCGATTAGTCGGGCGGACATGAGCAACTCCTTTTGTCATGTGTCGCCGACAGCATAGCAACACCGCTCAGGTGTGGGAGCGGCCTTGTGTCGCGAAAGGGCCGCTATGCGGccctttc                                                                                                                                          <  2:5170932/150‑5 (MQ=255)
            aCTGTGCAGGAAATTCAGGTCGGCCTGGGCAGCGATTAGTCGGGCGGACATGAGCAACTCCTTTTGTCATGTGTCGCCGACAGCATAGCAACACCGCTCAGGTGTGGGAGCGGCCTTTTGTGGCGAAATGGGCGCATAGCGGCCCCGGcc                                                                                                                                         >  2:1384808/1‑149 (MQ=255)
                   aGGAAATTCAGGTCGGCCTGGGCAGCGATTAGTCGGGCGGACATGAGCAACTCCTTTTGTCATGTGTCGCCGACAGCATAGCAACACCGCTCAGGTGTGGGAGCGGCCTTGTTTCGCGGAAGGGCCGCCTAGCCGCGCCGGCAAGTTAcg                                                                                                                                  >  1:3316151/1‑150 (MQ=255)
                   aGGAAATTCAGGTCGGCCTGGGCAGCGATTAGTCGGGCGGACATGAGCAACTCCTTTTGTCATGTGTCGCCGACAGCATAGCAACACCGCTCAGGTGTGGGAGCGGCCTTGTGTCGCGAAAGGGCCGCATAGCGGCCCCGGCatcaacat                                                                                                                                  >  1:5508244/1‑143 (MQ=255)
                   aGGAAATTCAGGTCGGCCTGGGCAGCGATTAGTCGGGCGGACATGAGCAACTCCTTTTGTCATGTGTCGCCGACAGCATAGCAACACCGCTCAGGTGTGGGAGCGGCCTTGTGTCGCGAAAGGGCCGCATAGCGGCCCCGGCAATTTCcg                                                                                                                                  <  2:4383154/150‑1 (MQ=255)
                    ggAAATTCAGGTCGGCCTGGGCAGCGATTAGTCGGGCGGACATGAGCAACTCCTTTTGTCATGTGTCGCCGACAGCATAGCAACACCGCTCAGGTGTGGGAGCGGCCTTGTGTGGCGGAAAGGCCGCAGAGCCGCCCCCGCACATTcccc                                                                                                                                 >  1:450360/1‑148 (MQ=255)
                    ggAAATTCAGGTCGGCCTGGGCAGCGATTAGTCGGGCGGACATGAGCAACTCCTTTTGTCATGTGTCGCCGACAGCATAGCAACACCGCTCAGGTGTGGGAGCGGCCTTGTGTCGCGAAAGGGCCGCATAGCGGCCCCGGCACATTCcgc                                                                                                                                 >  2:1072862/1‑150 (MQ=255)
                    ggAAATTCAGGTCGGCCTGGGCAGCGATTAGTCGGGCGGACATGAGCAACTCCTTTTGTCATGTGTCGCCGACAGCATAGCAACACCGCTCAGGTGTGGGAGCGGCCTTGTGTCGCGAAAGGGCCGCATAGCGGCCCCGGCAATTTCcgc                                                                                                                                 <  1:5113695/150‑1 (MQ=255)
                      aaaTTCAGGTCGGCCTGGGCAGCGATTAGTCGGGCGGACATGAGCAACTCCTTTTGTCATGTGTCGCCGACAGCATAGCAACACCGCTCAGGTGTGGGAGCGGCCTTGTGTCGCGAAAGGGCCGCATAGCGGCCCCGGCAATTTCCGCgt                                                                                                                               <  2:6728863/150‑1 (MQ=255)
                      aaaTTCAGGTCGGCCTGGGCAGCGATTAGTCGGGCGGACATGAGCAACTCCTTTTGTCATGTGTCGCCGACAGCATAGCAACACCGCTCAGGTGTGGGAGCGGCCTTGTGTCGCGAAAGGGCCGCATAGCGGCCCCGGCAATTTCCGCgt                                                                                                                               <  1:3398836/150‑1 (MQ=255)
                         ttCAGGTCGGCCTGGGCAGCGATTAGTCGGGCGGACATGAGCAACTCCTTTTGTCATGTGTCGCCGACAGCATAGCAACACCGCTCAGGTGTGGGAGCGGCCTTGTGTCGCGAAAGGGCCGCATAGCGGCCCCGGCAATTTCCGCGTGTa                                                                                                                            <  2:1629661/150‑1 (MQ=255)
                         ttCAGGTCGGCCTGGGCAGCGATTAGTCGGGCGGACATGAGCAACTCCTTTTGTCATGTGTCGCCGACAGCATAGCAACACCGCTCAGGTGTGGGAGCGGCCTTGTGTCGCGAAAGGGCCGCATAGCGGCCCCGGCAATTTCCGCGTGTa                                                                                                                            <  2:1476142/150‑1 (MQ=255)
                            aGGTCGGCCTGGGCAGCGATTAGTCGGGCGGACATGAGCAACTCCTTTTGTCATGTGTCGCCGACAGCATAGCAACACCGCTCAGGTGTGGGAGCGGCCTTGTGTCGCGAAAGGGCCGCATAGCGGCCCCGGCAATTTCCGCGTGTATGc                                                                                                                         <  1:4619962/150‑1 (MQ=255)
                                  gCCTGGGCAGCGATTAGTCGGGCGGACATGAGCAACTCCTTTTGTCATGTGTCGCCGACAGCATAGCAACACCGCTCAGGTGTGGGAGCGGCCTTGTTGCGACCAAGGGCACCAAAGCGCCCCCGGCCAATTCCGCGAGAAAGCAACcaa                                                                                                                   >  1:3970401/1‑150 (MQ=255)
                                  gCCTGGGCAGCGATTAGTCGGGCGGACATGAGCAACTCCTTTTGTCATGTGTCGCCGACAGCATAGCAACACCGCTCAGGTGTGGGAGCGGCCTTGTGTCGCGAAAGGGCCGCATAGCGGCCCCGGCAATTTCCGCGTGTATGCAACaaa                                                                                                                   <  1:1802601/150‑1 (MQ=255)
                                   ccTGGGCAGCGATTAGTCGGGCGGACATGAGCAACTCCTTTTGTCATGTGTCGCCGACAGCATAGCAACACCGCTCAGGTGTGGGAGCGGCCTTGTGTCGCGGAAGGGCCGCATAGCGGCCCCGGCCATTTCCGCCTGTATGCAACAAAc                                                                                                                  >  2:2249610/1‑150 (MQ=255)
                                   ccTGGGCAGCGATTAGTCGGGCGGACATGAGCAACTCCTTTTGTAATGTGTCGCCGACAGCATAGCAACACCGCTCAGGTGTGGGAGCGGCCTTGTGTCGCGAAAGGGCCGCATAGCGGCCCCGGCAATTTCCGCGTGTATGCAACAAAc                                                                                                                  <  1:320120/150‑1 (MQ=255)
                                    cTGGGCAGCGATTAGTCGGGCGGACATGAGCAACTCCTTTTGTCATGTGTCGCCGACAGCATAGCAACACCGCTCAGGTGTGGGAGCGGCCTTGTGTCGCGAAAGGGCCGCATAGCGGCCCCGGCAATTTCCGCGTGTATGCAACAAAcc                                                                                                                 >  1:4575113/1‑150 (MQ=255)
                                     tGGGCAGCGATTAGTCGGGCGGACATGAGCAACTCCTTTTGTCATGTGTCGCCGACAGCATAGCAACACCGCTCAGGTGTGGGAGCGGCCTTGTGTCGCGAAAGGGCCGCATAGCGGCCCCGGCAATTTCCGCGTGTATGCAACAAACCt                                                                                                                <  2:2243324/150‑1 (MQ=255)
                                      gggCAGCGATTAGTCGGGCGGACATGAGCAACTCCTTTTGTCATGTGTCGCCGACAGCATAGCAACACCGCTCAGGTGTGGGAGCGGCCTTGTGTCGCGAAAGGGCCGCATAGCGGCCCCGGCAATTTCCGCGTGTATGCAACAAACCTg                                                                                                               <  2:6020439/150‑1 (MQ=255)
                                        gCAGCGATTAGTCGGGCGGACATGAGCAACTCCTTTTGTCATGTGTCGCCGACAGCATAGCAACACCGCTCAGGTGTGGGAGCGGCCTTGTGTCGCGAAAGGGCCGCATAGCGGCCCCCGCAATTTACCCGTGTATGCCACAAACCTggg                                                                                                             >  2:2006832/1‑150 (MQ=255)
                                          aGCGATTAGTCGGGCGGACATGAGCAACTCCTTTTGTCATGTGTCGCCGACAGCATAGCAACACCGCTCAGGTGTGGGAGCGGCCTTGTGTCGCGAAAGGGCCGCATAGCGGCCCCGACAATTTCCGCGTGTATGCAACAAACCTGGGgc                                                                                                           >  2:5395590/1‑150 (MQ=255)
                                            cGATTAGTCGGGCGGACATGAGCAACTCCTTTTGTCATGTGTCGCCGACAGCATAGCAACACCGCTCAGGTGTGGGAGCGGCCTTGTGTCGCGAAAGGGCCGCATAGCGGCCCCGGCAATTTCCGCGTGTATGCAACAAACCTGGGgctg                                                                                                         <  2:1225256/150‑1 (MQ=255)
                                                 tgTCGGGCGGACATGAGCAACTCCTTTTGTAATGTGTCGCCGACAGCATAGCAACACCGCTCAGGTGTGGGAGCGGCCTTGTGTCGCGAAAGGGCCGCATAGCGGCCCCGGCAATTTCCGCGTGTATGCAACAAACCTGGGGCTGCTATg                                                                                                    <  1:3684107/149‑1 (MQ=255)
                                                 aGTCGGGCGGACATGAGCAACTCCTTTTGTCATGTGTCGCCGACAGCATAGCAACACCGCTCAGGTGTGGGAGCGGCCTTGTGTCGCGAAAGGGCCGCATAGCGGCCCCGGCAATTTCCGCGTGTATGCAACAAACCTGGGGCTGCTATg                                                                                                    <  1:3536743/150‑1 (MQ=255)
                                                  gTCGGGCGGACATGAGCAACTCCTTTTGTCATGTGTCGCCGACAGCATAGCAACACCGCTCAGGTGTGGGAGCGGCCTTGTGTCGCGAAAGGGCCGCATAGCGGCCCCGGCAATTTCCGCGTGTATGCAACAAACCTGGGGCTGCTATGc                                                                                                   <  1:3356261/150‑1 (MQ=255)
                                                   tCGGGCGGACATGAGCAACTCCTTTTGTCATGTGTCGCCGACAGCATAGCAACACCGCTCAGGTGTGGGAGCGGCCTTGTGTGGCGCAAGGGCCGCCGAGCGGCCCCCGCCCATTGCCCGGGTTTTGAACAAACCAGGGGCGGCGAAGCa                                                                                                  >  2:1136436/1‑150 (MQ=255)
                                                   tCGGGCGGACATGAGCAACTCCTTTTGTCATGTGTCGCCGACAGCATAGCAACACCGCTCAGGTGTGGGAGCGGCCTTGTGTCGCGCAAGGGCCGCATAGCGGCCCCCGCAATTTCCGCGTGTATGCAACAAACCAGGGGCTGCTATGCa                                                                                                  >  2:1567200/1‑150 (MQ=255)
                                                   tCGGGCGGACATGAGCAACTCCTTTTGTCATGTGTCGCCGACAGCATAGCAACACCGCTCAGGTGTGGGAGCGGCCTTGTGTCGCGAAAGGGCCGCATAGCGGCCCCGGCAATTTCCGCGTGTATGCAACAAACCTGGGGCTGCTATGCa                                                                                                  <  2:4913982/150‑1 (MQ=255)
                                                   tCGGGCGGACATGAGCAACTCCTTTTGTCATGTGTCGCCGACAGCATAGCAACACCGCTCAGGTGTGGGAGCGGCCTTGTGTCGCGAAAGGGCCGCAAAGCGGCCCCCGCAAATTCCGCGGGTATGCAACAAACCACGGGCGGCTATGCa                                                                                                  >  2:6708539/1‑150 (MQ=38)
                                                        cGGACATGAGCAACTCCTTTTGTCATGTGTCGCCGACAGCATAGCAACACCGCTCAGGTGTGGGAGCGGCCTTGTGTCGCGCAAGGGCCCCAAAGCGGCCCCCGCAAATTACCCGTGTATGCAACAAAACAGGGGCTGCTAAGCAGcccc                                                                                             >  2:3625962/1‑149 (MQ=37)
                                                        cGGAAATTAGAAACTCCTTTTGTCTTGTGTCGCCTAAAGCATAGAAACACCGCTCTGGTGTGGGTGGGGCCTTTTGTCGCGAATGGGCCGCATAGCGGCCCCGGCAATTTCCGCGTGTATGCAACAAACCTGGGGCTGCATAGCAGcccc                                                                                             <  1:2514105/150‑2 (MQ=255)
                                                           aCATGAGCAACTCCTTTTGTCATGTGTCGCCGACAGCATAGCAACACCGCTCAGGTGTGGGAGCGGCCTTGTGTCGCGAAAGGGCCGCATAGCGGCCCCGGCAATTTCCGCGTGTATGCAACAAACCTGGGGCTGCTATGCAGCCCTTTc                                                                                          <  2:5539987/150‑1 (MQ=255)
                                                             aTGAGCAACTCCTTTTGTCATGTGTCGCCGACAGCATAGCAACACCGCTCAGGTGTGGGAGCGGCCTTGTTTAGCGGAATGGCCCCCAAGCGGCCCCCGCCAATTAAGCCGGTGTTCAACAAAACACGGGCTGCGATGGAACCCattaac                                                                                        >  2:262493/1‑147 (MQ=255)
                                                             aTGAGCAACTCCTTTTGTCATGTGTCGCCGACAGCATAGCAACACCGCTCAGGTGTGGGAGCGGCCTTGTGTCGCGAAAGGGCCGCATAGCGGCCCCGGCAATTTCCGCGTGTATGCAACAAACCTGGGGCTGCTATGCAGCCCTTTcgc                                                                                        <  2:612137/150‑1 (MQ=255)
                                                              tGAGCAACTCCTTTTGTCATGTGTCGCCGACAGCATAGCAACACCGCTCAGGTGTGGGAGCGGCCTTGTGTCGCGAAAGGGCCGCATAGCGGCCCCGGCAATTTCCGCGTGTATGCAACAAACCTGGGGCTGCTATGCAGCCCTTTcgcg                                                                                       <  1:6389049/150‑1 (MQ=255)
                                                                aGCAACTCCTTTTGTCATGTGTCGCCGACAGCATAGCAACACCGCTCAGGTGTGGGAGCGGCCTTGTGTCGCGAAAGGGCCGCATAGCGGCCCCGGCAATTTCCGCGTGTATGCAACAAACCTGGGGCTGCTATGCAGCCCTTTCGCGac                                                                                     <  2:3309346/150‑1 (MQ=255)
                                                                 gCAACTCCTTTTGTCATGTGTCGCCGACAGCATAGCAACACCGCTCAGGTGTGGGAGCGGCCTTGTGTCGCGAAAGGGCCGCATAGCGGCCCCGGCAATTTCCCCCTGTATGCAACAAACCAGGCGCTGCTATGCAGCCCTTTCGCGaca                                                                                    >  2:2821260/1‑150 (MQ=255)
                                                                    aCTCCTTTTGTCATGTGTCGCCGACAGCATAGCAACACCGCTCAGGTGTGGGAGCGGCCTTGTGTCGCGAAAGGGCCGCCTATCGGCCCCGGCAAATTACGCGTGTATTCAACAAACCTGGGGCGGCAATGCACCCCGTTcgccccccaa                                                                                 >  2:5515059/1‑143 (MQ=37)
                                                                               cATGTGTCGCCGACAGCATAGCAACACCGCTCTGGTGTGGGAGCGGCGTTGTGTCGCGATCGGGCCGCATAGCGGCCCCGGCAATTTCCGCGTGTATGCAACAAACCGGGGGCTGCTATGCAGCCCTTTCGCGACACAAGGCCCACAAAg                                                                      <  1:2252439/150‑1 (MQ=38)
                                                                               cATGTGTCGCCGACAGCATAGCAACACCGCTCAGGTGTGGGAGCGGCCTTGTGTCGCGAAAGGGCCGCATAGCGGCCCCGGCAATTTCCGCGTGTATGCAACAAACCTGGGGCTGCTATGCAGCCCTTTCGCGACACAAGGCCCACAAAg                                                                      <  2:2926563/150‑1 (MQ=255)
                                                                                   tgtCGCCGACAGCATAGCAACACCGCTCAGGTGTGGGAGCGGCCTTGTGTCGCGAAAGGGCCGCGTAGCGGCCCCGGCAATTTCCGCGTGTATGCAACAAACCTGGGGCTGCTATGCAGCCCTTTCGCGACACAAGGCCCACAAAGGAcc                                                                  >  1:5700571/1‑150 (MQ=255)
                                                                                   tgtCGCCGACAGCATAGCAACACCGCTCAGGTGTGGGAGCGGCCTTGTGTCGCGAAAGGGCCGCATAGCGGCCCCGGCAATTTCCGCGTGTAt                                                                                                                           >  2:1201239/1‑93 (MQ=255)
                                                                                   tgtCGCCGACAGCATAGCAACACCGCTCAGGTGTGGGAGCGGCCTTGTGTCGCGAAAGGGCCGCATAGCGGCCCCGGCAATTTCCGCGTGTATGCAACAAACCTGGGGCTGCTATGCAGCCCTTTCGCGACACAAGGCCCACAAAGGAcc                                                                  <  2:2507833/150‑1 (MQ=255)
                                                                                       gCCGACAGCATAGCAACACCGCTCAGGTGTGGGAGCGGCCTTGTGTGGCGCAAAGGGCCCAAAGCGGCCCCCGCCAATTTCCCGGTTTTTGAAAAAAACAGGGGGGGGTATTGCGCCCGTTGGCGACACAAGGCCCAAAAAGGACCgcgc                                                              >  2:324876/1‑150 (MQ=255)
                                                                                        ccGACTGCATAGCAACACCGCTCAGGTGTGGGAGCGGCCTTGTGTCGCGAACGGGCCGCATAGCGGCCCCGGCAATTTCCGCGTGTATGCAACAAACCTGGGGCTGCTATGCAGCCCTTTCGCGACACAAGGCCCACAAAGGACCGCGCa                                                             <  1:6064805/150‑1 (MQ=255)
                                                                                             aGCATAGCAACACCGCTCAGGTGTGGGAGCGGCCTTGGGTCGCCCAAGGGCCCCAGACCCGCCCCCGGCAATTGCGCCGGGTTTGCACAAAACGGGGGCGGGGATGCCGCCCCGTGGCCGCCAAAGGCCCCCCAAGGGCCGCCCACGcct                                                        >  2:2426853/1‑148 (MQ=255)
                                                                                                 tAGCAACACCGCTCTGGTGTGGGTGCGGCCTTGTGTCGCTAAAGGGCCGCATAGCGGCCCCGGCAATTTCCGCGTGTATGCAACAAACCTGGGGCTGCTATGCAGCCCTTTCGCGACACAAGGCCCACAAAGGACCGCGCAGGCTTCAgg                                                    <  1:2667420/150‑1 (MQ=255)
                                                                                                    cAACACCGCTCAGGTGTGGGAGCGGCCTTGTGTCGCGAAAGGGCCGCATAGCGGCCCCGGCAATTTCCGCGTGTATGCAACAAACCTGGGGCTGCTATGCAGCCCTTTCGCGACACAAGGCCCACAAAGGACCGCGCAGGCTTCAGGccc                                                 <  2:2690276/150‑1 (MQ=255)
                                                                                                     aaCACCGCTCCGGTGTGGGAGCGGCCTTGTGTCGCGAAAGGGCCGCATAGCGGCCCCGGCAATTTCCGCGTGTATGCAACAAACCTGGGGCTGCTATGCAGCCCTTTCGCGACACAAGGCCCACAAAGGACCGCGCAGGCTTCAGGCCCg                                                <  2:6252267/150‑1 (MQ=255)
                                                                                                        cccGCTCTGGTGTGGGCGCGGCCTTGTGTCGCGAACGGCCCGCATAGCGGCCCCGGCAATTTCCGCGTGTATGCAACAAACCGGGGGCTGCTATGCAGCCCTTTCGCGACACAAGGCCCACAAAGGACCGCGCAGGCTTCAGGCCCGGaa                                             <  1:2858901/149‑1 (MQ=38)
                                                                                                        acCGCTCAGGTGTGGGAGCGGCCTTGTGTGGCGAAATGGCCGCAAAGCGGCCCCCGCAATTTACGCGTGTGTTCAACAAAACAGGGGCGGCTATGCAGCCCGTTAGCGACCCAACGCCCACAAAGGACCCCGCCGGCGTCAGGCCCGGac                                             >  2:4674796/1‑149 (MQ=255)
                                                                                                        acCGCTCAGGTGTGGGAGCGGCCTTGTGTCGCGAAAGGGCCGCATAGCGGCCCCGGCAATTTCCGCGTGTATGCAACAAACCTGGGGCTGCTATGCAGCCCTTTCGTGACACAAGGCCCACAAAGGACCGCGCAGGCTTCAGGCCCGGaa                                             >  2:4002257/1‑150 (MQ=255)
                                                                                                           gcgctggTGTGGGTGGGGCGTTGTGTCGCGAACGGGCCGAATCGCGGCCCCGGCAATTTCCGCGTGTATGCAACAACCCTGGGGCTGCTATGCAGCCCTTTCGCGACACAAGGCCCACAAAGGACCGCGCAGGCTTCAGGCCCGGAACTg                                          <  1:6549717/145‑1 (MQ=38)
                                                                                                           gCTCAGGTGTGGGAGCGGCCTTGTGTCGCGAAAGGGCCGCATAGCGGCCCCCGCAATTTCCGCGGGTATGCAACAAACCTGGGGCTGCTATGCAGCCCGTTCGCGACACAAGGCCCACAAAGGACCGCGCAGGCGTCAGGCCCGGAAcag                                          >  2:4572821/1‑148 (MQ=255)
                                                                                                             gctggTGTGGGAGGGGCCTTGTGGCGCTAATGGGCCGAATCGCGGCCCCGGCAATTTCCGCGTGTATGCAACAACCCTGGGGCTGCTATGCAGCCCTTTCGCGACACAAGGCCCACAAAGGACCGCGCAGGCTTCAGGCCCGGAACTGGc                                        <  1:78903/147‑1 (MQ=38)
                                                                                                              cAGGTGTGGGAGCGGCCTTGTGTCGCGAAAGGGCCGCATAGCGGCCCCGGCAATTTCCGCGTGTATGCAACAAACCTGGGGCTGCTATGCAGCCCTTTCGCGACACAAGGCCCACAAAGGACCGCGCAGGCTTCAGGCCCGGAACTGGcc                                       <  2:6681441/150‑1 (MQ=255)
                                                                                                              cAGGTGTGGGAGCGGCCTTGTGTCGCGAAAGGGCCGCATAGCGGCCCCGGCAATTCCCGCGTGTATGCACCAAACCTGGGGCTGCTATGCAGCCCTTTCGCGACACAAGGCCCACAAAGGACCGCGCAGGCTTCAGGCCCGGAACTGGcc                                       <  2:6859714/150‑1 (MQ=255)
                                                                                                                ggTGTGGGCGCGGCCTTGTGTCGCGAACGGGCCGCATAGCGGCCCCGGCAATTTCCGCGTGTATGCAACAAACCTGGGGCTGCTATGCAGCCCTTTCGCGACACAAGGCCCACAAAGGACCGCGCAGGCTTCAGGCCCGGAACTGGCCCa                                     <  1:4361889/150‑1 (MQ=255)
                                                                                                                  tgtgGGAGCGGCCTTGTGTCGCGAAAGGGCCGCATAGCGGCCCCGGCAATTTCCGCGTGTATGCAACAAACCTGGGGCTGCTATGCAGCCCTTTCGCGACACAAGGCCCACAAAGGACCGCGCAGGCTTCAGGCCCGGAACTGGCCCAgg                                   <  2:450328/150‑1 (MQ=255)
                                                                                                                       gAGCGGCCTTGTGTCGCGAAAGGGCCGCATAGCGGCCCCGGCAATTTCCGCGTGTATGCAACAAACCTGGGGCTGCTATGCAGCCCTTTCGCGACACAAGGCCCACAAAGGACCGCGCAGGCTTCAGGCCCGGAACTGGCCCAGGCTGgc                              <  2:1103185/150‑1 (MQ=255)
                                                                                                                          cGGCCTTGTGTCGCGAAAGGGCCCCAAATCGGCCCCCGCAAATTACGCGTGTATGCAAAAAAACCGGGGCGGCTATTCAGCCCATTAGCGAAAAAAAGCCCACAAAAGAACCCCCAAGCATAAGGCCCCGAACAGGCCCAAGCAGgcgcg                           >  2:4735852/1‑150 (MQ=255)
                                                                                                                            gCCTTGTGTCGCGAAAGGGCCGCATAGCGGCCCCGGCAATTTCCGCGTGTATGCAACAAACCTGGGGCTGCTATGCAGCCCTTTCGCGACACAAGGCCCACAAAGGACCGCGCAGGCTTCAGGCCCGGAACTGGCCCAGGCTGGCGCGCa                         >  2:1368517/1‑150 (MQ=255)
                                                                                                                            gCCTTGTGGCGCGAAAGGGCCGCATAGCGGCCCCGGCAATTTCCGCGTGTATGCAACAAACCTGGGGCTGCTATGCAGCCCTTTCGCGACACAAGGCCCACAAAGGACCGCGCAGGCTTCAGGCCCGGAACTGGCCCAGGCTGGCGCGCa                         <  1:6779870/150‑1 (MQ=255)
                                                                                                                               ttGTGTCGCGAAAGGGCCGCATAGCGGCCCCGGCAATTTCCGCGTGTATGCAACAAACCTGGGGCTGCTATGCAGGCCTTTCGCGACACAAGGCCCACAAAGGACCGCGCAGGCTTCAGGCCCGGAACTGGCCCAGGCTGGCGCGCAGCt                      >  1:2138875/1‑150 (MQ=255)
                                                                                                                               ttGTGTCGCGAAAGGGCCGCATAGCGGCCCCGGCAATTTCCGCGTGTATGCAACAAACCTGGGGCTGCTATGCAGCCCTTTCGCGACACAAGGCCCACAAAGGACCGCGCAGGCTTCAGGCCCGGAACTGGCCCAGGCTGGCGCGCAGCt                      >  2:2764543/1‑150 (MQ=255)
                                                                                                                                 gtgtGGCGAATGGGCCGCATAGCGGCCCCGGCAATTTCCGCGTGTATGCAACAAACCTGGGGCTGCTATGCAGCCCTTTCGCGACACAAGGCCCACAAAGGACCGCGCAGGCTTCAGGCCCGGAACTGGCCCAGGCTGGCGCGCAGCTGa                    <  1:4138737/150‑1 (MQ=255)
                                                                                                                                 gtgtCGCGAACTGGCCGCATAGCGGCCCCGGCAAGTTCCGCGTGTATGCAACAAACCTGGGTCTGCTATGCAGCCCTTTCGCGACACAAGGCCCACAAAGGACCGCGCAGGCTTCAGGCCCGGAACTGGCCCAGGCTGGCGCGCAGCTGa                    <  1:1113230/150‑1 (MQ=255)
                                                                                                                                 gtgtCGCGAAAGGGCCGCATAGCGGCCCCGGCAATTTCCGCGTGTATGCAACAACCCTGGGGCTGCTATGCAGCCCTTTCGCGACACAAGGCCCACAAAGGACCGCGCAGGCTTCAGGCCCGGAACTGGCCCAGGCTGGCGCGCAGCTGa                    <  2:2177909/150‑1 (MQ=255)
                                                                                                                                 gtgtCGCGAAAGGGCCGCATAGCGGCCCCGGCAATTTCCGCGTGTATGCAACAAACCTGGGGCTGCTATGCAGCCCTTTCGCGACACAAGGCCCACAAAGGACCGCGCAGGCTTCAGGCCCGGAACTGGCCCAGGCTGGCGCGCAGCTGa                    >  2:6362866/1‑150 (MQ=255)
                                                                                                                                 gtgtCGCGAAAGGGCCGCATAGCGGCCCCGGCAATTTCCGCGTGTATGCAACAAACCTGGGGCTGCTATGCAGCCCTTTCGCGACACAAGGCCCACAAAGGACCGCGCAGGCTTCAGGCCCGGAACTGGCCCAGGCTGGCGCGCAGCTGa                    <  2:1017027/150‑1 (MQ=255)
                                                                                                                                    tCGCGAAAGGGCCGCATAGCGGCCCCGGCAATTTCCGCGTGTATGCAACAAACCTGGGGCTGCTATGCAGCCCTTTCGCGACACAAGGCCCACAAAGGACCGCGCAGGCTTCAGGCCCGGAACTGGCCCAGGCTGGCGCGCAGCTGAgcc                 >  1:4954874/1‑150 (MQ=255)
                                                                                                                                      gcgAAAGGGCCGCATAGCGGCCCCGGCAATTTCCGCGTGTATGCAACAAACCTGGGGCTGCTATGCAGCCCTTTCGCGACACAAGGCCCACAAAGGACCGCGCAGGCTTCAGGCCCGGAACTGGCCCAGGCTGGCGCGCAGCTGAgccgc               >  2:3901914/1‑150 (MQ=255)
                                                                                                                                      gcgAAAGGGCCGCATAGCGGCCCCGGCAATTTCCGCGTGTATGCAACAAACCTGGGGCTGCTATGCAGCCCTTTCGCGACACAAGGCCCACAAAGGACCGCGCAGGCTTCAGGCCCGGAACTGGCCCAGGCTGGCGCGCAGCTGAgccgc               >  2:1702389/1‑150 (MQ=255)
                                                                                                                                      gcgAAAGGGCCGCATAGCGGCCCCGGCAATTTCCGCGTGTATGCAACAAACCTGGGGCTGCTATGCAGCCCTTTCGCGACACAAGGCCCACAAAGGACCGCGCAGGCTTCAGGCCCGGAACTGGCCCAGGCTGGCGCGCAGCTGAgccgc               >  1:476148/1‑150 (MQ=255)
                                                                                                                                      gcgAAAGGGCCGCATAGCGGCCCCGGCAATTTCCGCGTGTATGCAACAAACCTGGGGCTGCTAGGCAGCCCTTTCGCGACACAAGGCCCACCAAGGACCGCGCAGGCTTCAGGCCCGGAACTGGCCCAGGCTGGCGCGCAGCTGAgccgc               >  1:664072/1‑150 (MQ=255)
                                                                                                                                       cgAAAGGGCCGCATAGCGCCCCCGGCAATTTCCGCGTGTATGCAACAAACCTGGGGCTGCTATGCAGCCCTTTCGCGACACAAGGCCCACAAAGGACCGCGCAGGCTTCAGGCCCGGAACTGGCCCAGGCTGGCGCGCAGCTGAgccgcc              <  1:779228/150‑1 (MQ=255)
                                                                                                                                        gAAAGGGCCGCATAGCGGCCCCGGCAATTTCCGCGTGTATGCAACACACCTGGGGCTGCTATGCAGCCCTTTCGCGACACAAGGCCCACAAAGGACCGCGCAGGCTTCAGGCCCGGAACTGGCCCAGGCTGGCGCGCAGCTGAGCCGCCa             <  2:3802338/150‑1 (MQ=255)
                                                                                                                                          aaaGGGCCGCATAGCGGCCCCGGCAATTTCCGCGTGTATGCAACAAACCTGGGGCTGCTATGCAGCCCTTTCGCGACACAAGGCCCACAAAGGACCGCGCAGGCTTCAGGCCCGGAACTGGCCCAGGCTGGCGCGCAGCTGAGCCGCCAg            >  2:5084538/1‑150 (MQ=255)
                                                                                                                                           cggggCCGCTATGCGGCCCCGGCAATTTCCGCGTGTATGCAACAAACCTGGGGCTGCTATGCAGCCCTTTCGCGACACAAGGCCCACAAAGGACCGCGCAGGCTTCAGGCCCGGAACTGGCCCAGGCTGGCGCGCAGCTGAGCCGCCAGa           >  1:1329387/3‑150 (MQ=255)
                                                                                                                                            ggggCCGCTATGCGGCCCCGGCAATTTCCGCGTGTATGCAACAAACCTGGGGCTGCTATGCAGCCCTTTCGCGACACAAGGCCCACAAAGGACCGCGCAGGCTTCAGGCCCGGAACTGGCCCAGGCTGGCGCGCAGCTGAGCCGCCAGat          >  2:4148349/2‑150 (MQ=255)
                                                                                                                                            ggggCCGCTATGCGGCCCCGGCAATTTCCGCGTGTATGCAACAAACCTGGGGCTGCTATGCAGCCCTTTCGCGACACAAGGCCCACAAAGGACCGCGCAGGCTTCAGGCCCGGAACTGGCCCAGGCTGGCGCGCAGCTGAGCCGCCAGat          >  2:3732735/2‑150 (MQ=255)
                                                                                                                                             gggCCGCTATGCGGCCCCGGCAATTTCCGCGTGTATGCAACAAACCTGGGGCTGCTATGCAGCCCTTTCGCGACACAAGGCCCACAAAGGACCGCGCAGGCTTCAGGCCCGGAACTGGCCCAGGCTGGCGCGCAGCTGAGCCGCCAGatc         >  1:2865585/1‑150 (MQ=255)
                                                                                                                                               gCCGCATAGCGGCCCCGGCAATTTCCGCGTGTATGCAACAAACCTGGGGCTGCTATGCAGCCCTTTCGCGACACAAGGCCCACAAAGGACCGCGCAGGCTTCAGGCCCGGAACTGGCCCAGGCTGGCGCGCAGCTGAGCCGCCAGatcat       <  1:5129666/150‑1 (MQ=255)
                                                                                                                                               gCCGCATAGCGGCCCCGGCAATTTCCGCGTGTATGCAACAAACCTGGGGCTGCTATGCAGCCCTTTCGCGACACAAGGCCCACAAAGGACCGCGCAGGCTTCAGGCCCGGAACTGGCCCAGGCTGGCGCGCAGCTGAGCCGCCAGatcat       >  1:5798254/1‑150 (MQ=255)
                                                                                                                                               gCCGCATAGCGGCCCCGGCAATTTCCGCGTGTATGCAACAAACCTGGGGCTGCTATGCAGCCCTTTCGCGACACAAGGCCCACAAAGGACCGCGCAGGCTTCAGGCCCGGAACTGGCCCAGGCTGGCGCGCAGCTGAGCCGCCAGatcat       >  2:580993/1‑150 (MQ=255)
                                                                                                                                               gCCGCATAGCGGCCCCGGCAATTTCCGCGTGTATGCAACAAACCTGGGGCTGCTATGCAGCCCTTTCGCGACACAAGGCCCACAAAGGACCGCGCAGGCTTCAGGCCCGGAACTGGCCCAGGCTGGCGCGCAGCTGAGCCGCCAGatcat       <  2:4101654/150‑1 (MQ=255)
                                                                                                                                                ccGCATAGCGGCCCCGGCAATTTCCGCGTGTATGCAACAAACCTGGGGCTGCTATGCAGCCCTTTCGCGACACAAGGCCCACAAAGGACCGCGCAGGCTTCAGGCCCGGAACTGGCCCAGGCTGGCGCGCAGCTGAGCCGCCAGatcatc      <  1:1050692/150‑1 (MQ=255)
                                                                                                                                                ccGCATAGCGGCCCCGACAATTTCCGCGTGTATGCAACAAACCTGGGGCTGCTATGCAGCCCTTTCGCGACACAAGGCCCACAAAGGACCGCGCAGGCTTCAGGCCCGGAACTGGCCCAGGCTGGCGCGCAGCTGAGCCGCCAGatcatc      <  1:5395590/150‑1 (MQ=255)
                                                                                                                                                 cGCATAGCGGCCCCGGCAATTTCCGCGTGTATGCAACAAACCTGGGGCTGCTATGCAGCCCTTTCGCGACACAAGGCCCACAAAGGACCGCGCAGGCTTCAGGCCCGGAACTGGCCCAGGCTGGCGCGCAGCTGAGCCGCCAGATCATcc     >  1:3830173/1‑150 (MQ=255)
                                                                                                                                                   cATAGCGGCCCCGGCAATTTCCGCGTGTATGCAACAAACCTGGGGCTGCTATGCAGCCCTTTCGCGACACAAGGCCCACAAAGGACCGCGCAGGCTTCAGGCCCGGAACTGGCCCAGGCTGGCGCGCAGCTGAGCCGCCAGATCATCCAg   <  1:5301844/150‑1 (MQ=255)
                                                                                                                                                    aTAGCGGCCCCGGCAATTTCCGCGTGTATGCAACAAACCTGGGGCTGCTATGCAGCCCTTTCGCGACACAAGGCCCACAAAGGACCGCGCAGGCTTCAGGCCCGGAACTGGCCCAGGCTGGCGCGCAGCTGAGCCGCCAGATCATCCAGc  <  2:2140203/150‑1 (MQ=255)
                                                                                                                                                      |                                                                                                                                                     
CCGTACGCTGTCACTGTGCAGGAAATTCAGGTCGGCCTGGGCAGCGATTAGTCGGGCGGACATGAGCAACTCCTTTTGTCATGTGTCGCCGACAGCATAGCAACACCGCTCAGGTGTGGGAGCGGCCTTGTGTCGCGAAAGGGCCGCATAGCGGCCCCGGCAATTTCCGCGTGTATGCAACAAACCTGGGGCTGCTATGCAGCCCTTTCGCGACACAAGGCCCACAAAGGACCGCGCAGGCTTCAGGCCCGGAACTGGCCCAGGCTGGCGCGCAGCTGAGCCGCCAGATCATCCAGC  >  NC_002947/5560453‑5560749

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: