Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,110,625 T→A 5.0% intergenic (+38/+47) PP_3612 → / ← PP_3613 TonB‑dependent receptor/L‑sorbosone dehydrogenase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,110,6250TA 5.0% 1060.2 / 26.2 320intergenic (+38/+47)PP_3612/PP_3613TonB‑dependent receptor/L‑sorbosone dehydrogenase
Reads supporting (aligned to +/- strand):  ref base T (166/137);  new base A (12/5);  total (178/142)
Fisher's exact test for biased strand distribution p-value = 2.22e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 330 total aligned reads displayed.

CTGAAGAACGTGTTCGACCGTGACTACTATGCCTCGGCACATGGCTCGGTCGACCTGATCACCCCAGGGGCACCGCGCACGCTGGAAGCCAGCTTCAACTACCGCTTCTGACCCACGCGCAATGCCTGTCAGTCAGCAATTGGGCCGCTCTGCGGCCCATCGCCAGCAAGCTTTGACGATTTGATCAGATAGCATTCAGCGATTACGCGTCACTCGCCACACCGTATTCGCCAGGTCATCGGCAATGATCAACGCCCCACGCGGGTCCACCGTCACGCCAACCGGTCGCCCCCGGGTC  >  NC_002947/4110477‑4110774
                                                                                                                                                    |                                                                                                                                                      
 tGAAGAACGTGTTCGACCGTGACTACTATGCCTCGGCACATGGCTCGGTCGACCTGATCACCCCAGGGGCACCGCGCACGCTGGAAGCCAGCTTCAACTACCGCTTCTGACCCACGCGCAATGCCTGTCAGTCAGCAATTGGGCCGctct                                                                                                                                                      >  2:4585892/1‑150 (MQ=255)
  gaagaaCGTGTTCGACCGTGACTACTATGCCTCGGCACATGGCTCGGTCGACCTGATCACCCCAGGGGCACCGCGCACGCTGGAAGCCAGCTTCAACTACCGCTTCTGACCCACGCGCAATGCCTGTCAGTCAGCAATTGGGCCGCTCTg                                                                                                                                                     >  1:3439625/1‑150 (MQ=255)
   aagaaCGTGTTCGACCGTGACTACTATGCCTCGGCACATGGCTCGGTCGACCTGATCACCCCAGGGGCACCGCGCACGCTGGAAGCCAGCTTCAACTACCGCTTCTGACCCACGCGCAATGCCTGTCAGTCAGCAATTGGGCCGCTCTGc                                                                                                                                                    <  2:7384259/150‑1 (MQ=255)
    agaaCGTGTTCGACCGTGACTACTATGCCTCGGCACATGGCTCGGTCGACCTGATCACCCCAGGGGCACCGCGCACGCTGGAAGCCAGCTTCAACTACCGCTTCTGACCCACGCGCAATGCCTGTCAGTCAGCAATTGGGCCGCTCTGCg                                                                                                                                                   <  2:878632/150‑1 (MQ=255)
      aaCGTGTTCGACCGTGACTACTATGCCTCGGCACATGGCTCGGTCGACCTGATCACCCCAGGGGCACCGCGCACGCTGGAAGCCAGCTTCAACTACCGCTTCTGACCCACGCGCAATGCCTGTCAGTCAGCAATTGGGCCGCTCTGCGGc                                                                                                                                                 <  1:1850545/150‑1 (MQ=255)
        cGTGTTCGACCGTGACTACTATGCCTCGGCACATGGCTCGGTCGACCTGATCACCCCAGGGGCACCGCGCACGCTGGAAGCCAGCTTCAACTACCGCTTCTGACCCACGCGCAATGCCTGTCAGTCAGCAATTGGGCCGCTCTGCGGccc                                                                                                                                               >  1:4628021/1‑150 (MQ=255)
           gtTCGACCGTGACTACTATGCCTCGGCACATGGCTCGGTCGACCTGATCACCCCAGGGGCACCGCGCACGCTGGAAGCCAGCTTCAACTACCGCTTCTGACCCACGCGCAATGCCTGTCAGTCAGCAATTGGGCCGCTCTGCGGCCCATc                                                                                                                                           <  1:5081156/150‑1 (MQ=255)
            ttCGACCGTGACTACTATGCCTCGGCACATGGCTCGGTCGACCTGATCACCCCAGGGGCACCGCGCACGCTGGAAGCCAGCTTCAACTACCGCTTCTGACCCACGCGCAATGCCTGTCAGTCAGCAATTGGGCCGCTCTGCGGCCAATCg                                                                                                                                          >  2:4665628/1‑150 (MQ=255)
             tCGACCGTGACTACTATGCCTCGGCACATGGCTCGGTCGACCTGATCACCCCAGGGGCACCGCGCACGCTGGAAGCCAGCTTCAACTACCGCTTCTGACCCACGCGCAATGCCTGTCAGTCAGCAATTGGGCCGCAGAGCGGCCCaattg                                                                                                                                         <  2:6073349/150‑5 (MQ=255)
                   gTGACTACTATGCCTCGGCACATGGCTCGGTCGACCTGATCACCCCAGGGGCACCGCGCACGCTGGAAGCCAGCTTCAACTACCGCTTCTGACCCACGCGCAATGCCTGTCAGTCAGCAATTGGGCCGCTCTGCGGCCCATCGCCAGCaa                                                                                                                                   >  2:1808177/1‑150 (MQ=255)
                        tactaTGCCTCGGCACATGGCTCGGTCGACCTGATCACCCCAGGGGCACCGCGCACGCTGGAAGCCAGCTTCAACTACCGCTTCTGACCCACGCGCAATGCCTGTCAGTCAGCAATTGGGCCGCTCTGCGGCCCATCGCCAGCAAGCttt                                                                                                                              >  1:6752837/1‑150 (MQ=255)
                        tactaTGCCTCGGCACATGGCTCGGTCGACCTGATCACCCCAGGGGCACCGCGCACGCTGGAAGCCAGCTTCAACTACCGCTTCTGACCCACGCGCAATGCCTGTCAGTCAGCAATTGGGCCGCTCTGCGGCCCATCGCCAGCAAGCttt                                                                                                                              >  1:3143077/1‑150 (MQ=255)
                          ctaTGCCTCGGCACATGGCTCGGTCGACCTGATCACCCCAGGGGCACCGCGCACGCTGGAAGCCAGCTTCAACTACCGCTTCTGACCCACGCGCAATGCCTGTCAGTCAGCAATTGGGCCGCTCTGCGGCCCATCGCCAGCAAGCTTTGa                                                                                                                            <  2:1261591/150‑1 (MQ=255)
                          ctaTGCCTCGGCACATGGCTCGGTCGACCTGATCACCCCAGGGGCACCGCGCACGCTGGAAGCCAGCTTCAACTACCGCTTCTGACCCACGCGCAATGCCTGTCAGTCAGCAATTGGGCCGCTCTGCGGCCCATCGCCAGCAAGCTTTGa                                                                                                                            <  1:5412391/150‑1 (MQ=255)
                          ctaTGCCTCGGCACATGGCTCGGTCGACCTGATCACCCCAGGGGCACCGCGCACGCTGGAAGCCAGCTTCAACTACCGCTTCTGACCCACGCGCAATGCCTGTCAGTCAGCAATTGGGCCGCTCTGCGGCCCATCGCCAGCAAGCTTTGa                                                                                                                            >  2:3552606/1‑150 (MQ=255)
                          ctaTGCCTCGGCACATGGCTCGGTCGACCTGATCACCCCAGGGGCACCGCGCACGCTGGAAGCCAGCTTCAACTACCGCTTCTGACCCACGCGCAATGCCTGTCAGTCAGCAATTGGGCCGCTCTGCGGCCCATCGCCAGCAAGCTTTGa                                                                                                                            >  2:1923913/1‑150 (MQ=255)
                          ctaTGCCTCGGCACATGGCTCGGTCGACCTGATCACCCCAGGGGCACCGCGCACGCTGGAAGCCAGCTTCAACTACCGCTTCTGACCCACGCGCAATGCCTGTCAGTCAGCAATTGGGCCGCTCTGCGGCCCATCGCCAGCAAGCTTTGa                                                                                                                            >  2:5411624/1‑150 (MQ=255)
                           taTGCCTCGGCACATGGCTCGGTCGACCTGATCACCCCAGGGGCACCGCGCACGCTGGAAGCCAGCTTCAACTACCGCTTCTGACCCACGCGCAATGCCTGTCAGTCAGCAATTGGGCCGCTCTGCGGCCCATCGCCAGCAAGCTTTGAc                                                                                                                           >  1:6968935/1‑150 (MQ=255)
                              gCCTCGGCACATGGCTCGGTCGACCTGATCACCCCAGGGGCACCGCGCACGCTGGAAGCCAGCTTCAACTACCGCTTCTGACCCACGCGCAATGCCTGTCAGTCAGCAATTGGGCCGCTCTGCGGCCCATCGCCAGCAAGCTTTGACGAt                                                                                                                        >  2:4886732/1‑150 (MQ=255)
                                cTCGGCACATGGCTCGGTCGACCTGATCACCCCAGGGGCACCGCGCACGCTGGAAGCCAGCTTCAACTACCGCTTCTGACCCACGCGCAATGCCTGTCAGTCAGCAATTGGGCCGCTCTGCGGCCCATCGCCAGCAAGCTTTGACGAttt                                                                                                                      >  1:5497266/1‑150 (MQ=255)
                                  cGGCACATGGCTCGGTCGACCTGATCACCCCAGGGGCACCGCGCACGCTGGAAGCCAGCTTCAACTACCGCTTCTGACCCACGCGCAATGCCTGTCAGTCAGCAATTGGGCCGCTCTGCGGCCCATCGCCAGCAAGCTTTGACGATTTGa                                                                                                                    >  2:3481853/1‑150 (MQ=255)
                                    gCACATGGCTCGGTCGACCTGATCACCCCAGGGGCACCGCGCACGCTGGAAGCCAGCTTCAACTACCGCTTCTGACCCACGCGCAATGCCTGTCAGTCAGCAATTGGGCCGCTCTGCGGCCCATCGCCAGCAAGCTTTGACGATTTGATc                                                                                                                  <  2:7008908/150‑1 (MQ=255)
                                        aTGGCTCGGTCGACCTGATCACCCCAGGGGCACCGCGCACGCTGGAAGCCAGCTTCAACTACCGCTTCTGACCCACGCGCAATGCCTGTCAGTCAGCAATTGGGCCGCTCTGCGGCCCATCGCCAGCAAGCTTTGACGATTTGATCAGAt                                                                                                              <  2:1666969/150‑1 (MQ=255)
                                         tGGCTCGGTCGACCTGATCACCCCAGGGGCACCGCGCACGCTGGAAGCCAGCTTCAACTACCGCTTCTGACCCACGCGCAATGCCTGTCAGTCAGCAATTGGGCCGCTCTGCGGCCCATCGCCAGCAAGCTTTGACGATTTGATCAGATa                                                                                                             <  2:2220563/150‑1 (MQ=255)
                                          ggCTCGGTCGACCTGATCACCCCAGGGGCACCGCGCACGCTGGAAGCCAGCTTCAACTACCGCTTCTGACCCACGCGCAATGCCTGTCAGTCAGCAATTGGGCCGCTCTGCGGCCCATCGCCAGCAAGCTTTGACGATTTGATCAGATAg                                                                                                            >  2:5068131/1‑150 (MQ=255)
                                             tCGGTCGACCTGATCACCCCAGGGGCACCGCGCACGCTGGAAGCCAGCTTCAACTACCGCTTCTGACCCACGCGCAATGCCTGTCAGTCAGCAATTGGGCCGCTCTGCGGCCCATCGCCAGCAAGCTTTGACGATTTGATCAGATAGCAt                                                                                                         >  2:2000973/1‑150 (MQ=255)
                                             tCGGTCGACCTGATCACCCCAGGGGCACCGCGCACGCTGGAAGCCAGCTTCAACTACCGCTTCTGACCCACGCGCAATGCCTGTCAGTCAGCAATTGGGCCGCTCTGCGGCCCATCGCCAGCAAGCTTTGACGATTTGATCAGATAGCAt                                                                                                         >  1:1017301/1‑150 (MQ=255)
                                             tCGGTCGACCTGATCACCCCAGGGGCACCGCGCACGCTGGAAGCCAGCTTCAACTACCGCTTCTGACCCACGCGCAATGCCTGTCAGTCAGCAATTGGGCCGCTCTGCGGCCCATCGCCAGCAAGCTTTGACGATTTGATCAGATAGCAt                                                                                                         >  1:5743010/1‑150 (MQ=255)
                                                    aCCTGATCACCCCCGGGGCACCGCGCACGCTGGAAGCCAGCTTCAACTACCGCTTCTGACCCACGCGCAATGCCTGTCAGTCAGCAATTGGGCCGCTCTGCGGCCCATCGCCAGCAAGCTTTGACGATTTGATCAGATAGCATTCAGCGa                                                                                                  <  1:610333/150‑1 (MQ=255)
                                                      cTGATCACCCCAGGGGCACCGCGCACGCTGGAAGCCAGCTTCAACTACCGCTTCTGACCCACGCGCAATGCCTGTCAGTCAGCAATTGGGCCGCTCTGCGGCCCATCGCCAGCAAGCTTTGACGATTTGATCAGATAGCATTCAGCGAtt                                                                                                >  1:2721891/1‑150 (MQ=255)
                                                          tCACCCCAGGGGCACCGCGCACGCTGGAAGCCAGCTTCAACTACCGCTTCTGACCCACGCGCAATGCCTGTCAGTCAGCAATTGGGCCGCTCTGCGGCCCATCGCCAGCAAGCTTTGACGATTTGATCAGATAGCATTCAGCGATTAcgc                                                                                            <  1:2281766/150‑1 (MQ=255)
                                                              cccAGGGGCACCGCGCACGCTGGAAGCCAGCTTCAACTACCGCTTCTGACCCACGCGCAATGCCTGTCAGTCAGCAATTGGGCCGCTCTGCGGCCCATCGCCAGCAAGCTTTGACGATTTGATCAGATAGCATTCAGCGATTACGCGTCa                                                                                        >  1:4488218/1‑150 (MQ=255)
                                                              cccAGGGGCACCGCGCACGCTGGAAGCCAGCTTCAACTACCGCTTCTGACACACGCGCAATGCCTGTCAGTCAGCAATTGGGCCGCTCTGCGGCCCATCGCCAGCAAGCTTTGACGATTTGATCAGATAGCATTCAGCGATTACGCGTCa                                                                                        >  2:4746730/1‑150 (MQ=255)
                                                               ccAGGGGCACCGCGCACGCTGGAAGCCAGCTTCAACTACCGCTTCTGACCCACGCGCAATGCCTGTCAGTCAGCAATTGGGCCGCTCTGCGGCCCATCGCCAGCAAGCTTTGACGATTTGATCAGATAGCATTCAGCGATTACGCGTCAc                                                                                       >  1:2657692/1‑150 (MQ=255)
                                                                cAGGGGCACCGCGCACGCTGGAAGCCAGCTTCAACTACCGCTTCTGACCCACGCGCAATGCCTGTCAGTCAGCAATTGGGCCGCTCTGCGGCCCATCGCCAGCAAGCTTTGACGATTTGATCAGATAGCATTCAGCGATTACGCGTCACt                                                                                      >  1:3637381/1‑150 (MQ=255)
                                                                  ggggCACCGCGCACGCTGGAAGCCAGCTTCAACTACCGCTTCTGACCCACGCGCAATGCCTGTCAGTCAGCAATTGGGCCGCAGAGCGGCCCaatt                                                                                                                                          <  2:144399/96‑4 (MQ=255)
                                                                     gCACCGCGCACGCTGGAAGCCAGCTTCAACTACCGCTTCTGACCCACGCGCAATGCCTGTCAGTCAGCAATTGGGCCGCTCTGCGGCCCATCGCCAGCAAGCTTTGACGATTTGATCAGATAGCATTCAGCGATTACGCGTCACTCGCca                                                                                 <  2:5237778/150‑1 (MQ=255)
                                                                      cACCGCGCACGCTGGAAGCCAGCTTCAACTACCGCTTCTGACCCACGCGCAATGCCTGTCAGTCAGCAATTGGGCCGCTCTGCGGCCCATCGCCAGCAAGCTTTGACGATTTGATCAGATAGCATTCAGCGATTACGCGTCACTcaccac                                                                                >  2:4459259/1‑150 (MQ=255)
                                                                            gcACGCTGGAAGCCAGCTTCAACTACCGCTTCTGACCCACGCGCAATGCCTGTCAGTCAGCAATTGGGCCGCTCTGCGGCCCATCGCCAGCAAGCTTTGACGATTTGATCAGATAGCATTCAGCGATTACGCGTCACTCGCcacattgtc                                                                          <  2:5600856/150‑6 (MQ=255)
                                                                             cACGCTGGAAGCCAGCTTCAACTACCGCTTCTGACCCACGCGCAATGCCTGTCAGTCAGCAATTGGGCCGCTCTGCGGCCCATCGCCAGCAAGCTTTGACGATTTg                                                                                                                     <  1:4902485/106‑1 (MQ=255)
                                                                               cGCTGGAAGCCAGCTTCAACTACCGCTTCTGACCCACGCGCAATGCCTGTCAGTCAGCAATTGGGCCGCTCTGCGGCCCATCGCCAGCAAGCTTTGACGATTTGATCAGATAGCATTCAGCGATTACGCGTCACTCGCCACACCGTATTc                                                                       <  2:4246518/150‑1 (MQ=255)
                                                                               cGCTGGAAGCCAGCTTCAACTACCGCTTCTGACCCACGCGCAATGCCTGTCAGTCAGCAATTGGGCCGCTCTGCGGCCCATCGCCAGCAAGCTTTGACGATTTGATCAGATAGCATTCAGCGATTACGCGTCACTCGCCACACCGTATTc                                                                       <  1:5411624/150‑1 (MQ=255)
                                                                                 cTGGAAGCCAGCTTCAACTACCGCTTCTGACCCACGCGCAATGCCTGTCAGTCAGCAATTGGGCCGCTCTGCGGCCCATCGCCAGCAAGCTTTGACGATTTGATCAGATAGCATTCAGCGATTACGCGTCACTCGCCACACCGTATTCGc                                                                     <  2:3760608/150‑1 (MQ=255)
                                                                                 cTGGAAGCCAGCTTCAACTACCGCTTCTGACCCACGCGCAATGCCTGTCAGTCAGCAATTGGGCCGCTCTGCGGCCCATCGCCAGCAAGCTTTGACGATTTGATCAGATAGCATTCAGCGATTACGCGTCACTCGCCACACCGTATTCGc                                                                     <  2:1037654/150‑1 (MQ=255)
                                                                                   ggAAGCCAGCTTCAACTACCGCTTCTGACCCACGCGCAATGCCTGTCAGTCAGCAATTGGGCCGGTCTGCGGCCCATCGCCAGCAAGCTTTGACGATTTGATCAGATAGCATTCAGCGATTACGCGTCACTCGCCACACCGTATTCGCCa                                                                   <  1:4094328/150‑1 (MQ=255)
                                                                                   ggAAGCCAGCTTCAACTACCGCTTCTGACCCACGCGCAATGCCTGTCAGTCAGCAATTGGGCCGCTCTGCGGCCCATCGCCAGCAAGCTTTGACGATTTGATCAGATAGCATTCAGCGATTACGCGTCACTCGCCACACCGTATTCGCCa                                                                   <  1:771245/150‑1 (MQ=255)
                                                                                        ccAGCTTCAACTACCGCTTCTGACCCACGCGCAATGCCTGTCAGTCAGCAATTGGGCCGCTCTGCGGCCCATCGCCAGCAAGCTTTGACGATTTGATCAGATAGCATTCAGCGATTACGCGTCACTCGCCACACCGTATTCGCCAGGTCa                                                              >  2:6565154/1‑150 (MQ=255)
                                                                                          atcTTAAACTACCGCTTCTCACCCACCCGAAATGCCTTTCATTCAGAAATTGGGCCGCTCTGCGGCCCATCGCCAGCAAGCTTTGACGATTTGATCAGATAGCATTCAGCGATTACGCGTCACTCGCCACACCGTATTCGCCAGGTCATc                                                            <  1:878004/148‑1 (MQ=255)
                                                                                          aGCTTCAACTACCGCTTCTGACCCACGCGCAATGCCTGTCAGTCAGCAATTGGGCCGCTCTGCGGCCCATCGCCAGCAAGCTTTGACGATTTGATCAGATAGCATTCAGCGATTACGCGTCACTCGCCACACCGTATTCGCCAGGTCATc                                                            <  2:4253262/150‑1 (MQ=255)
                                                                                             ttCAACTACCGCTTCTGACCCACGCGCAATGCCTGTCAGTCAGCAATTGGGCCGCTCTGCGGCCCATCGCCAGCAAGCTTTGACGATTTGATCAGATAGCATTCAGCGATTACGCGTCACTCGCCACACCGTATTCGCCAGGTCATCGGc                                                         >  2:3077696/1‑150 (MQ=255)
                                                                                                aaCTACCGCTTCTGACCCACGCGCAATGCCTGTCAGTCAGCAATTGGGCCGCTCTGCGGCCCATCGCCAGCAAGCTTTGACGATTTGATCAGATAGCATTCAGCGATTACGCGTCACTCGCCACACCGTATTCGCCAGGTCATCGGCAAt                                                      >  1:1595655/1‑150 (MQ=255)
                                                                                                aaCTACCGCTTCTGACCCACGCGCAATGCCTGTCAGTCAGCAATTGGGCCGCTCTGCGGCCCATCGCCAGCAAGCTTTGACGATTTGATCAGATAGCATTCAGCGATTACGCGTCACTCGCCACACCGTATTCGCCAGGTCATCGGCAAt                                                      >  1:3959428/1‑150 (MQ=255)
                                                                                                  cTACCGCTTCTGACCCACGCGCAATGCCTGTCAGTCAGCAATTGGGCCGCTCTGCGGCCCATCGCCAGCAAGCTTTGACGATTTGATCAGATAGCATTCAGCGATTACGCGTCACTCGCCACACCGTATTCGCCAGGTCATCGGCAATGa                                                    >  1:404336/1‑150 (MQ=255)
                                                                                                   tACCGCTTCTGACCCACGCGCAATGCCTGTCAGTCAGCAATTGGGCCGCTCTGCGGCCCATCGCCAGCAAGCTTTGACGATTTGATCAGATAGCATTCAGCGATTACGCGTCACTCGCCACACCGTATTCGCCAGGTCATCGGCAATGAt                                                   >  1:1429311/1‑150 (MQ=255)
                                                                                                   tACCGCTTCTGACCCACGCGCAATGCCTGTCAGTCAGCAATTGGGCCGCTCTGCGGCCCATCGCCAGCAAGCTTTGACGATTTGATCAGATAGCATTCAGCGATTACGCGTCACTCGCCACACCGTATTCGCCAGGTCATCGGCAATGAt                                                   >  2:6450512/1‑150 (MQ=255)
                                                                                                    aCCGCTTCTGACCCACGCGCAATGCCTGTCAGTCAGCAATTGGGCCGCTCTGCGGCCCATCGCCAGCAAGCTTTGACGATTTGATCAGATAGCATTCAGCGATTACGCGTCACTCGCCACACCGTATTCGCCAGGTCATCGGCAATGATc                                                  >  1:4191613/1‑150 (MQ=255)
                                                                                                     ccGCTTCTGACCCACGCGCAATGCCTGTCAGTCAGCAATTGGGCCGCTCTGCGGCCCATCGCCAGCAAGCTTTGACGATTTGATCAGATAGCATTCAGCGATTACGCGTCACTCGCCACACCGTATTCGCCAGGTCATCGGCAATGATCa                                                 >  1:191519/1‑150 (MQ=255)
                                                                                                             gACCCACGCGCAATGCCTGTCAGTCAGCAATTGGGCCGCTCTGCGGCCCATCGCCAGCAAGCTTTGACGATTTGATCa                                                                                                                 <  1:6517171/78‑1 (MQ=255)
                                                                                                              aCCCACGCGCAATGCCTGTCAGTCAGCAATTGGGCCGCTCTGCGGCCCATCGCCAGCAAGCTTTGACGATTTGATCAGATAGCATTCAGCGATTACGCGTCACTCGCCACACCGTATTCGCCAGGTCATCGGCAATGATCAACGCCCCAc                                        >  1:6626297/1‑150 (MQ=255)
                                                                                                              aCCCACGCGCAATGCCTGTCAGTCAGCAATTGGGCCGCTCTGCGGCCCATCGCCAGCAAGCTTTGACGATTTGATCAGATAGCATTCAGCGATTACGCGTCACTCGCCACACCGTATTCGCCAGGTCATCGGCAATGATCAACGCCCCAc                                        >  2:5726156/1‑150 (MQ=255)
                                                                                                                ccACGCGCAATGCCTGTCAGTCAGCAATTGGGCCGCTCTGCGGCCCATCGCCAGCAAGCTTTGACGATTTGATCAGATAGCATTCAGCGATTACGCGTCACTCGCCACACCGTATTCGCCAGGTCATCGGCAATGATCAACGCCCCAcgc                                      >  1:2579795/1‑150 (MQ=255)
                                                                                                                ccACGCGCAATGCCTGTCAGTCAGCAATTGGGCCGCTCTGCGGCCCATCGCCAGCAAGCTTTGACGATTTGATCAGATAGCATTCAGCGATTACGCGTCACTCGCCACACCGTATTCGCCAGGTCATCGGCAATGATCAACGCCCCAcgc                                      >  1:4280780/1‑150 (MQ=255)
                                                                                                                  aCGCGCAATGCCTGTCAGTCAGCAATTGGGCCGCTCTGCGGCCCATCGCCAGCAAGCTTTGACGATTTGATCAGATAGCATTCAGCGATTACGCGTCACTCGCCACACCGTATTCGCCAGGTCATCGGCAATGATCAACGCCCCACGCgg                                    >  2:2147500/1‑150 (MQ=255)
                                                                                                                      gcAATGCCTGTCAGTCAGCAATTGGGCCGCTCTGCGGCCCATCGCCAGCAAGCTTTGACGATTTGATCAGATAGCATTCAGCGATTACGCGTCACTCGCCACACCGTATTCGCCAGGTCATCGGCAATGATCAACGCCCCACGCGGGTcc                                >  2:7183173/1‑150 (MQ=255)
                                                                                                                        aaTGCCTGTCAGTCAGCAATTGGGCCGCTCTGCGGCCCATCGCCAGCAAGCTTTGACGATTTGATCAGATAGCATTCAGCGATTACGCGTCACTCGCCACACCGTATTCGCCAGGTCATCGGCAATGATCAACGCCCCACGCGGGTCCAc                              >  1:1253351/1‑150 (MQ=255)
                                                                                                                        aaTGCCTGTCAGTCAGCAATTGGGCCGCTCTGCGGCCCATCGCCAGCAAGCTTTGACGATTTGATCAGATAGCATTCAGCGATTACGCGTCACTCGCCACACCGTATTCGCCAGGTCATCGGCAATGATCAACGCCCCACGCGGGTCCAc                              >  2:4813746/1‑150 (MQ=255)
                                                                                                                           gCCTGTCAGTCAGCAATTGGGCCGCTCTGCGGCCCATCGCCAGCAAGCTTTGACGATTTGATCAGATAGCATTCAGCGATTACGCGTCACTCGCCACACCGTATTCGCCAGGTCATCGGCAATGa                                                    <  2:5888566/125‑1 (MQ=255)
                                                                                                                             cTGTCAGTCAGCAATTGGGCCGCTCTGCGGCCCATCGCCAGCAAGCTTTGAGGATTTGATCAGATAGCATTCAGCGATTACGCGTCACTCGCCACACCGTATTCGCCAGGTCATCGGCAATGATCAACGCCCCACGCGGGTCCACCGTCa                         >  2:774897/1‑150 (MQ=255)
                                                                                                                             cTGTCAGTCAGCAATTGGGCCGCTCTGCGGCCCATCGCCAGCAAGCTTTGACGATTTGATCAGATAGCATTCAGCGATTACGCGTCACTCGCCACACCGTATTCGCCAGGTCATCGGCAATGATCAACGCCCCACGCGGGTCCACCGTCa                         <  1:1978715/150‑1 (MQ=255)
                                                                                                                             cTGTCAGTCAGCAAATGGGCCGCTCTGCGGCACATCGCCAGCAAGCTTTGACGATTTGATCAGATAGCATACAGCGATTACGCGTCACTCGCCACACCGTATTCGCCAGGTCATCGGCAATGATCAACGCCCCACGCGGGTCCACCGTCa                         >  1:4787948/1‑150 (MQ=255)
                                                                                                                              tGTCAGTCAGCAATTGGGCCGCTCTGCGGCCCATCGCCAGCAAGCTTTGACGATTTGATCAGATAGCATTCAGCGATTACGCGTCACTCGCCACACCGTATTCGCCAGGTCATCGGCAATGATCAACGCCCCACGCGGGTCCACCGTCAc                        >  1:1673587/1‑150 (MQ=255)
                                                                                                                               gtcagtcaGCAATTGGGCCGCTCTGCGGCCCATCGCCAGCAAGCTTTGACGATTTGATCAGATAGCATTCAGCGATTACGCGTCACTCGCCACACCGTATTCGCCAGGTCATCGGCAATGATCAACGCCCCACGCGGGTCCACCGTCACg                       >  1:7201635/1‑150 (MQ=255)
                                                                                                                                tcagtcagCAATTGGGCCGCTCTGCGGCCCATCGCCAGCAAGCTTTGACGATTTGATCAGATAGCATTCAGCGATTACGCGTCACTCGCCACACCGTATTCGCCAGGTCATCGGCAATGATCAACGCCCCACGCGGGTCCACCGTCACGc                      <  2:7266494/150‑1 (MQ=255)
                                                                                                                                tcagtcagCAATTGGGCCGCTCTGCGGCCCATCGCCAGCAAGCTTTGACGATTTGATCAGATAGCATTCAGCGATTACGCGTCACTCGCCACACCGTATTCGCCAGGTCATCGGCAATGATCAACGCCCCACGCGGGTCCACCGTCACGc                      >  1:7008419/1‑150 (MQ=255)
                                                                                                                                 cagtcagCAATTGGGCCGCTCTGCGGCCCATCGCCAGCAAGCTTTGACGATTTGATCAGATAGCATTCAGCGATTACGCGTCACTCGCCACACCGTATTCGCCAGGTCATCGGCAATGATCAACGCCCCACGCGGGTCCACCGTCACGcc                     <  1:5146732/150‑1 (MQ=255)
                                                                                                                                  agtcagCAATTGGGCCGCTCTGCGGCCCATCGCCAGCAAGCTTTGACGATTTGATCAGATAGCATTCAGCGATTACGCGTCACTCGCCACACCGTATTCGCCAGGTCATCGGCAATGATCAACGCCCCACGCGGGTCCACCGTCACGCCa                    >  2:2751418/1‑150 (MQ=255)
                                                                                                                                    tcagCAATTGGGCCGCTCTGCGGCCCATCGCCAGCAAGCTTTGACGATTTGATCAGATAGCATTCAGCGATTACGCGTCACTCGCCACACCGTATTCGCCAGGTCATCGGCAATGATCAACGCCCCACGCGGGTCCACCGTCACGCCAAc                  >  2:3171162/1‑150 (MQ=255)
                                                                                                                                    tcagCAATTGGGCCGCTCTGCGGCCCATCGCCAGCAAGCTTTGACGATTTGATCAGATAGCATTCAGCGATTACGCGTCACTCGCCACACCGTATTCGCCAGGTCATCGGCAATGATCAACGCCCCACGCGGGTCCACCGTCACGCCAAc                  <  2:426338/150‑1 (MQ=255)
                                                                                                                                      tggcgatGGGCCGCAGAGCGGCCCGATTGCCCGCACGCTTGGACGATTTGATCAGTTAGCATTCAGCGATTACGCGTCACTCGCCACACCGGATTCGCCAGGTCATCGGCACTGATCAACGCCCCACGCGGGGCCACCGTCACGCCAAcc                 >  2:5346368/7‑150 (MQ=255)
                                                                                                                                      agCAATTGGGCCGCTCTGCGGCCCATCGCCAGCAAGCTTTGACGATTTGATCAGATAGCATTCAGCGATTACGCGTCACTCGCCACACCGTATTCGCCAGGTCATCGGCAATGATCAACGCCCCACGCGGGTCCACCGTCACGCCAACCg                >  2:6851963/1‑150 (MQ=255)
                                                                                                                                      agCAATTGGGCCGCTCTGCGGCCCATCGCCAGCAAGCTTTGACGATTTGATCAGATAGCATTCAGCGATTACGCGTCACTCGCCACACCGTATTCGCCAGGTCATCGGCAATGATCAACGCCCCACGCGGGTCCACCGTCACGCCAACCg                <  2:6201665/150‑1 (MQ=255)
                                                                                                                                      agCAATTGGGCCGCTCTGCGGCCCATCGCCAGCAAGCTTTGACGATTTGATCAGATAGCATTCAGCGATTACGCGTCACTCGCCACACCGTATTCGCCAGGTCATCGGCAATGATCAACGCCCCACGCGGGTCCACCGTCACGCCAACCg                <  2:5497266/150‑1 (MQ=255)
                                                                                                                                       ggcgatGGGCCGCAGAGCGGCTCATCGCCAGCAAGCTTTGACGATTTGATCAGATAGCATTCAGCGATTACGCGTCACTCGCCACACCGTATTCGCCAGGTCATCGGCAATGATCAACGCCCCACGCGGGTCCACCGTCACGCCAACCgg               >  2:2446809/6‑150 (MQ=255)
                                                                                                                                        cAATTGGGCCGCTCTGCGGCCCATCGCCAGCAAGCTTTGACGATTTGATCAGATAGCATTCAGCGATTACGCGTCACTCGCCACACCGTATTCGCCAGGTCATCGGCAATGATCAACGCCCCACGCGGGTCCACCGTCACGCCAACCGGt              <  2:2548725/150‑1 (MQ=255)
                                                                                                                                        cAATTGGGCCGCTCTGCGGCCCATCGCCAGCAAGCTTTGACGATTTGATCAGATAGCATTCAGCGATTACGCGTCACTCGCCACACCGTATTCGCCAGGTCATCGGCAATGATCAACGCCCCACGCGGGTCCACCGTCACGCCAACCGGt              >  1:4423057/1‑150 (MQ=255)
                                                                                                                                         aaTTGGGCCGCTCTGCGGCCCATCGCCAGCAAGCTTTGACGATTTGATCAGATAGCATTCAGCGATTACGCGTCACTCGCCACACCGTATTCGCCAGGTCATCGGCAATGATCAACGCCCCACGCGGGTCCACCGTCACGCCAACCGGTc             >  2:2282740/1‑150 (MQ=255)
                                                                                                                                         aaTTGGGCCGCTCTGCGGCCCATCGCCAGCAAGCTTTGACGATTTGATCAGATAGCATTCAGCGATTACGCGTCACTCGCCACACCGTATTCGCCAGGTCATCGGCAATGATCAACGCCCCACGCGGGTCCACCGTCACGCCAACCGGTc             <  2:1313851/150‑1 (MQ=255)
                                                                                                                                          gatGGGCCGCAGAGCGGCCCATCGCCAGCAAGCTTTGACGATTTGATCAGATAGCATTCAGCGATTACGCGTCACTCGCCACACCGTATTCGCCAGGTCATCGGCAATGATCAACGCCCCACGCGGGTCCACCGTCACGCCAACCGGTCg            >  2:1377601/3‑150 (MQ=255)
                                                                                                                                          gatGGGCCGCAGAGCGGCCCATCGCCAGCAAGCTTTGACGATTTGATCAGATAGCATTCAGCGATTACGCGTCACTCGCCACACCGTATTCGCCAGGTCATCGGCAATGATCAACGCCCCACGCGGGTCCACCGTCACGCCAACCGGTCg            >  1:6354791/3‑150 (MQ=255)
                                                                                                                                          aTTGGGCCGCTCTGCGGCCCATCGCCAGCAAGCTTTGACGATTTGATCAGATAGCATTCAGCGATTACGCGTCACTCGCCACACCGTATTCGCCAGGTCATCGGCAATGATCAACGCCCCACGCGGGTCCACCGTCACGCCAACCGGTCg            >  1:2972486/1‑150 (MQ=255)
                                                                                                                                            tGGGCCGCAGAGCGGCCCATCGCCAGCAAGCTTTGACGATTTGATCAGATAGCATTCAGCGATTACGCGTCACTCGCCACACCGTATTCGCCAGGTCATCGGCAATGATCAACGCCCCACGCGGGTCCACCGTCACGCCAACCGGTCGcc          >  2:4049192/1‑150 (MQ=255)
                                                                                                                                              ggCCGCTCTGCGGCCCATCGCCAGCAAGCTTTGACGATTTGATCAGATAGCATTCAGCGATTACGCGTCACTCGCCACACCGTATTCGCCAGGTCATCGGCAATGATCAACGCCCCACGCGGGTCCACCGTCACGCCAACCGGTCGcccc        >  1:2173136/1‑150 (MQ=255)
                                                                                                                                               gCCGCTCTGCGGCCCATCGCCAGCAAGCTTTGACGATTTGATCAGATAGCATTCAGCGATTACGCGTCACTCGCCAGACCGTATTCGCCAGGTCATCGGCAATGATCAACGCCCCACGCGGGTCCACCGTCACGCCAACCGGTCGccccc       >  1:6174732/1‑150 (MQ=255)
                                                                                                                                               gCCGCTCTGCGGCCCATCGCCAGCAAGCTTTGACGATTTGATCAGATAGCATTCAGCGATTACGCGTCACTCGCCACACCGTATTCGCCAGGTCATCGGCAATGATCAACGCCCCACGCGGGTCCACCGTCACGCCAACCGGTCGccccc       >  1:3880444/1‑150 (MQ=255)
                                                                                                                                               gCCGCTCTGCGGCCCATCGCCAGCAAGCTTTGACGATTTGATCAGATAGCATTCAGCGATTACGCGTCACTCGCCACACCGTATTCGCCAGGTCATCGGCAATGATCAACGCCCCACGCGGGTCCACCGTCACGCCAACCGGTCGccccc       >  1:4752079/1‑150 (MQ=255)
                                                                                                                                                 cGCTCTGCGGCCCATCGCCAGCAAGCTTTGACGATTTGATCAGATAGCATTCAGCGATTACGCGTCACTCGCCACACCGTATTCGCCAGGTCATCGGCAATGATCAACGCCCCACGCGGGTCCACCGTCACGCCAACCGGTCGCCCCCgg     <  2:6722506/150‑1 (MQ=255)
                                                                                                                                                 cGCTCTGCGGCCCATCGCCAGCAAGCTTTGACGATTTGATCAGATAGCATTCAGCGATTACGCGTCACTCGCCACACCGTATTCGCCAGGTCATCGGCAATGATCAACGCCCCACGCGGGTCCACCGTCACGCCAACCGGTCGCCCCCgg     <  1:5982963/150‑1 (MQ=255)
                                                                                                                                                  gCTCTGCGGCCCATCGCCAGCAAGCTTTGACGATTTGATCAGATAGCATTCAGCGATTACGCGTCACTCGCCACACCGTATTCGCCAGGTCATCGGCAATGATCAACGCCCCACGCGGGTCCACCGTCACGCCAACCGGTCGCCCCCggg    >  1:3174501/1‑150 (MQ=255)
                                                                                                                                                    tctGCGGCCCATCGCCAGCAAGCTTTGACGATTTGATCAGATAGCATTCAGCGATTACGCGTCACTCGCCACACCGTATTCGCCAGGTCATCGGCAATGATCAACGCCCCACGCGGGTCCACCGTCACGCCAACCGGTCGCCCCCGGGTc  >  1:4984354/1‑150 (MQ=255)
                                                                                                                                                    tctGCGGCCCATCGCCAGCAAGCTTTGACGATTTGATCAGATAGCATTCAGCGATTACGCGTCACTCGCCACACCGTATTCGCCAGGTCATCGGCAATGATCAACGCCCCACGCGGGTCCACCGTCACGCCAACCGGTCGCCCCCGGGTc  <  1:5315128/150‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                      
CTGAAGAACGTGTTCGACCGTGACTACTATGCCTCGGCACATGGCTCGGTCGACCTGATCACCCCAGGGGCACCGCGCACGCTGGAAGCCAGCTTCAACTACCGCTTCTGACCCACGCGCAATGCCTGTCAGTCAGCAATTGGGCCGCTCTGCGGCCCATCGCCAGCAAGCTTTGACGATTTGATCAGATAGCATTCAGCGATTACGCGTCACTCGCCACACCGTATTCGCCAGGTCATCGGCAATGATCAACGCCCCACGCGGGTCCACCGTCACGCCAACCGGTCGCCCCCGGGTC  >  NC_002947/4110477‑4110774

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