Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,930,186 G→C 9.1% Q11E (CAG→GAG)  cheA ← chemotaxis histidine kinase CheA

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,930,1860GC 9.1% 270.8 / 14.9 88Q11E (CAG→GAG) cheAchemotaxis histidine kinase CheA
Reads supporting (aligned to +/- strand):  ref base G (39/41);  new base C (3/5);  total (42/46)
Fisher's exact test for biased strand distribution p-value = 7.16e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.21e-01

CAGCTGAAGGAAGCCGGCGCCCCCTTTTACAGTGTGGAAACCGCGAAAGATCGCATTGAGCAGGTCGGCATCGTCAGGTCGGCTTTCCAGCTCGACCAGTTGCTCGGACAATTGCTCAAGTATTTCGCCGGCCTCTACCAGGAAATCCTGGAGGATTTCTTCATCGGCGCCGAAGCTCATCAAACGTGCTCCTTAAAAACCAAGGCTGGACAGCAGATCATCGACATCGTCCTGACCGGACACGACGTCTTCACGCTTATCGGCATGAATCTGCGGACCTTCACCCCGAGTCGGATGT  >  NC_002947/4930037‑4930334
                                                                                                                                                     |                                                                                                                                                    
cAGCTGAAGGAAGCCGGCGCCCCCTTTTACAGTGTGGAAACCGCGAAAGATCGCATTGAGCAGGTCGGCATCGTCAGGTCGGCTTTCCAGCTCGACCAGTTGCTCGGACAATTGCTCAAGTATTTCGCCGGCCTCTACCAGGAAATCCTg                                                                                                                                                      <  1:2346072/150‑1 (MQ=255)
 aGCTGAAGGAAGCCGGCGCCCCCTTTTACAGTGTGGAAACCGCGAAAGATCGCATTGAGCAGGTCGGCATCGTCAGGTCGGCTTTCCAGCTCGACCAGTTGCTCGGACAATTGCTCAAGTATTTCGCCTGCCTCTACCAGGAAATCCTgg                                                                                                                                                     >  1:797181/1‑150 (MQ=255)
  gCTGAAGGAAGCCGGCGCCCCCTTTTACAGTGTGGAAACCGCGAAAGATCGCATTGAGCAGGTCGGCATCGTCAGGTCGGCTTTCCAGCTCGACCAGTTGCTCGGACAATTGCTCAAGTATTTCGCCGGCCTCTACCAGGAAATCCTgga                                                                                                                                                    <  1:1582686/150‑1 (MQ=255)
   cTGAAGGAAGCCGGCGCCCCCTTTTACAGTGTGGAAACCGCGAAAGATCGCATTGAGCAGGTCGGCATCGTCAGGTCGGCTTTCCAGCTCGACCAGTTGCTCGGACAATTGCTCAAGTATTTCGCCGGCCTCTACCAGGAAATCCTggag                                                                                                                                                   <  1:821166/150‑1 (MQ=255)
     gaaggaagCCGGCGCCCCCTTTTACAGTGTGGAAACCGCGAAAGATCGCATTGAGCAGGTCGGCATCGTCAGGTCGGCTTTCCAGCTCGACCAGTTGCTCGGACAATTGCTCAAGTATTTCGCCGGCCTCTACCAGGAAATCCTggagga                                                                                                                                                 <  2:1177279/150‑1 (MQ=255)
        ggaagCCGGCGCCCCCTTTTACAGTGTGGAAACCGCGAAAGATCGCATTGAGCAGGTCGGCATCGTCAGGTCGGCTTTCCAGCTCGACCAGTTGCTCGGACAATTGCTCAAGTATTTCGCCGGCCTCTACCAGGAAATCCTGGAGGAttt                                                                                                                                              >  2:295201/1‑150 (MQ=255)
        ggaagCCGGCGCCCCCTTTTACAGTGTGGAAACCGCGAAAGATCGCATTGAGCAGGTCGGCATCGTCAGGTCGGCTTTCCAGCTCGACCAGTTGCTCGGACAATTGCTCAAGTATTTCGCCGGCCTCTACCAGGAAATCCTCCAGGAttt                                                                                                                                              >  1:1414040/1‑150 (MQ=255)
         gaagCCGGCGCCCCCTTTTACAGTGTGGAAACCGCGAAAGATCGCATTGAGCAGGTCGGCATCGTCAGGTCGGCTTTCCAGCTCGACCAGTTGCTCGGACAATTGCTCAAGTATTTCGCCGGCCTCTACCAGGAAATCCTGGAGGATttc                                                                                                                                             <  1:618094/150‑1 (MQ=255)
         gaagCCGGCGCCCCCTTTTACAGTGTGGAAACCGCGAAAGATCGCATTGAGCAGGTCGGCATCGTCAGGTCGGCTTTCCAGCTCGACCAGTTGCTCGGACAATTGCTCAAGTATTTCGCCGGCCTCTACCAGGAAATCCTGGAGGATttc                                                                                                                                             >  2:257570/1‑150 (MQ=255)
           agCCGGCGCCCCCTTTTACAGTGTGGAAACCGCGAAAGATCGCATTGAGCAGGTCGGCATCGTCAGGTCGGCTTTCCAGCTCGACCAGTTGCTCGGACAATTGCTCAAGTATTTCGCCGGCCTCTACCAGGAAATCCTCCAGGATttcct                                                                                                                                           <  2:151794/150‑3 (MQ=255)
            gCCGGCGCCCCCTTTTACAGTGTGGAAACCGCGAAAGATCGCATTGAGCAGGTCGGCATCGTCAGGTCGGCTTTCCAGCTCGACCAGTTGCTCGGACAATTGCTCAAGTATTTCGCCGGCCTCTACCAGGAAATCCTCCAGGATttcctg                                                                                                                                          <  2:930952/150‑4 (MQ=255)
            gCCGGCGCCCCCTTTTACAGTGTGGAAACCGCGAAAGATCGCATTGAGCAGGTCGGCATCGTCAGGTCGGCTTTCCAGCTCGACCAGTTGCTCGGACAATTGCTCAAGTATTTCGCCGGCCTCTACCAGGAAATCCTCCAGGATttcctg                                                                                                                                          <  2:1414040/150‑4 (MQ=255)
             ccGGCGCCCCCTTTTACAGTGTGGAAACCGCGAAAGATCGCATTGAGCAGGTCGGCATCGTCAGGTCGGCTTTCCAGCTCGACCAGTTGCTCGGACAATTGCTCAAGTATTTCGCCGGCCTCTACCAGGAAATCCTggagga                                                                                                                                                 >  2:444737/1‑142 (MQ=255)
             ccGGCGCCCCCTTTTACAGTGTGGAAACCGCGAAAGATCGCATTGAGCAGGTCGGCATCGTCAGGTCGGCTTTCCAGCTCGACCAGTTGCTCGGACAATTGCTCAAGTATTTCGCCGGCCTCTACCAGGAAATCCTggagga                                                                                                                                                 <  1:444737/142‑1 (MQ=255)
               ggCGCCCCCTTTTACAGTGTGGAAACCGCGAAAGATCGCATTGAGCAGGTCGGCATCGTCAGGTCGGCTTTCCAGCTCGACCAGTTGCTCGGACAATTGCTCAAGTATTTCGCCGGCCTCTACCAGGAAATCCTGGAGGATTTCTTCATc                                                                                                                                       <  2:1032480/150‑1 (MQ=255)
               ggCGCCCCCTTTTACAGTGTGGAAACCGCGAAAGATCGCATTGAGCAGGTCGGCATCGTCAGGTCGGCTTTCCAGCTCGACCAGTTGCTCGGACAATTGCTCAAGTATTTCGCCGGCCTCTACCAGGAAATCCTGGAGGATTTCTTCATc                                                                                                                                       <  1:2359181/150‑1 (MQ=255)
                  gcCCCCTTTTACAGTGTGGAAACCGCGAAAGATCGCATTGAGCAGGTCGGCATCGTCAGGTCGGCTTTCCAGCTCGACCAGTTGCTCGGACAATTGCTCAAGTATTTCGCCGGCCTCTACCAGGAAATCCTGGAGGATTTCTTCATCGgc                                                                                                                                    <  2:1307800/150‑1 (MQ=255)
                    ccccTTTTACAGTGTGGAAACCGCGAAAGATCGCATTGAGCAGGTCGGCATCGTCAGGTCGGCTTTCCAGCTCGACCAGTTGCTCGGACAATTGCTCAAGTATTTCGCCGGCCTCTACCAGGAAATCCTCCAGGATttcct                                                                                                                                           <  1:1607073/141‑3 (MQ=255)
                    ccccTTTTACAGTGTGGAAACCGCGAAAGATCGCATTGAGCAGGTCGGCATCGTCAGGTCGGCTTTCCAGCTCGACCAGTTGCTCGGACAATTGCTCAAGTATTTCGCCGGCCTCTACCAGGAAATCCTCCAGGATttcct                                                                                                                                           >  2:1607073/1‑139 (MQ=255)
                      ccTTTTACAGTGTGGCAACCGCGAAAGATCGCATTGAGCAGGGCGGCGTCGTCAGGGCGGCTTTCCAGCTCGACCAGTTGCTCGGACAATTGCTCAAGTATTTCGCCGGCCTCTACCAGGAAATCCTGGAGGATTTCTTCATCGGCGCCg                                                                                                                                >  2:399747/1‑150 (MQ=255)
                      ccTTTTACAGTGTGGAAACCGCGAAAGATCGCATTGAGCAGGTCGGCATCGTCAGGTCGGCTTTCCAGCTCGACCAGTTGCTCGGACAATTGCTCAAGTATTTCGCCGGCCTCTACCAGGAAATCCTGGAGGATTTCTTCATCGGCGCCg                                                                                                                                >  2:978763/1‑150 (MQ=255)
                                  tgGAAACCGCGAAAGATCGCATTGAGCAGGTCGGCATCGTCAGGTCGGCTTTCCAGCTCGACCAGTTGCTCGGACAATTGCTCAAGTATTTCGCCGGCCTCTACCAGGAAATCCTGGAGGATTTCTTCATCGGCGCCGAAGCTCATCaaa                                                                                                                    <  1:1763982/150‑1 (MQ=255)
                                   ggAAACCGCGAAAGATCGCATTGAGCAGGTCGGCATCGTCAGGTCGGCTTTCCAGCTCGACCAGTTGCTCGGACAATTGCTCAAGTATTTCGCCGGCCTCTACCAGGAAATCCTGGAGGATTTCTTCATCGGCGCCGAAGCTCATCAAAc                                                                                                                   <  1:1695219/150‑1 (MQ=255)
                                     aaaCCGCGAAAGATCGCATTGAGCAGGTCGGCATCGTCAGGTCGGCTTTCCAGCTCGACCAGTTGCTCGGACAATTGCTCAAGTATTTCGCCGGCCTCTACCAGGAAATCCTGGAGGAttt                                                                                                                                              <  2:121589/121‑1 (MQ=255)
                                     aaaCCGCGAAAGATCGCATTGAGCAGGTCGGCATCGTCAGGTCGGCTTTCCAGCTCGACCAGTTGCTCGGACAATTGCTCAAGTATTTCGCCGGCCTCTACCAGGAAATCCTGGAGGAttt                                                                                                                                              >  1:121589/1‑121 (MQ=255)
                                       aCCGCGAAAGATCGCATTGAGCAGGTCGGCATCGTCAGGTCGGCTTTCCAGCTCGACCAGTTGCTCGGACAATTGATCTATCCGCACCTCTGCCTCTACCAGGAAATCCTGGCGGTTCTTTTCCCCCGCCTTCAACCTCATCAAACGTGc                                                                                                               >  1:655976/1‑150 (MQ=255)
                                         cgctAAATATCGCATTGAGCAGGTCGGCATCGTCAGGTCGGCTTTCCAGCTCGACCAGTTGCTCGGACAATTGCTCAAGTATTTCGCCGGCCTCTACCAGGAAATCCTGGAGGATTTCTTCATCGGCGCCGAAGCTCATCAAACGTGCTc                                                                                                             <  1:500011/150‑1 (MQ=255)
                                          gcgAAAGATCGCATTGAGCAGGTCGGCATCGTCAGGTCGGCTTTCCAGCTCGACCAGTTGCTCGGACAATTGCTCAAGTATTTCGCCGGCCTCTACCAGGAAATCCTGGAGGATTTCTTCATCGGCGCCGAAGCTCATCAAACGTGCTcc                                                                                                            >  2:1048135/1‑150 (MQ=255)
                                              aaGATCGCATTGAGCAGGTCGGCATCGTCAGGTCGGCTTTCCAGCTCGACCAGTTGCTCGGACAATTGCTCAAGTATTTCGCCTGCCTCTACCAGGAAATCCTGGAGGATTTCTTCATCGGCGCCGAAGCTCATCAAACGTGCTCCTTaa                                                                                                        <  2:797181/150‑1 (MQ=255)
                                                  tCGCATTGAGCAGGTCGGCATCGTCAGGTCGGCTTTCCAGCTCGACCAGTTGCTCGGACAATTGCTCAAGTATTTCGCCGGCCTCTACCAGGAAATCCTGGAGGATTTCTTCATCGGCGCCGAAGCTCATCAAAc                                                                                                                   >  2:770908/1‑135 (MQ=255)
                                                  tCGCATTGAGCAGGTCGGCATCGTCAGGTCGGCTTTCCAGCTCGACCAGTTGCTCGGACAATTGCTCAAGTATTTCGCCGGCCTCTACCAGGAAATCCTGGAGGATTTCTTCATCGGCGCCGAAGCTCATCAAAc                                                                                                                   <  1:770908/135‑1 (MQ=255)
                                                  tCGCATTGAGCAGGTCGGCATCGTCAGGTCGGCTTTCCAGCTCGACCAGTTGCTCGGACAATTGCTCAAGTATTTCGCCGGCCTCTACCAGGAAATCCTGGAGGATTTCTTCATCGGCGCCGAAGCTCATCAAACGTGCTCCTTAAAAAc                                                                                                    >  1:886562/1‑150 (MQ=255)
                                                     cATTGAGCAGGTCGGCATCGTCATGTCGGCTTTCCAGCTCGACCAGTTTCTCGGACAATTGCTCAAGTATTTCGCCGGCCTCTACCAGGAAATCCTGGAGGATTTCTTCATCGGCGCCGAAGCTCATCAAACGTGCTCCTTAAAAACCaa                                                                                                 <  1:181319/150‑1 (MQ=255)
                                                      aTTGAGCAGGTCGGCATCGTCAGGTCGGCTTTCCAGCTCGACCAGTTGCTCGGACAATTGCTCAAGTATTTCGCCGGCCTCTACCAGGAAATCCTGGAGGATTTCTTCATCGGCGCCGAAGCtcat                                                                                                                        <  2:2017181/126‑1 (MQ=255)
                                                      aTTGAGCAGGTCGGCATCGTCAGGTCGGCTTTCCAGCTCGACCAGTTGCTCGGACAATTGCTCAAGTATTTCGCCGGCCTCTACCAGGAAATCCTGGAGGATTTCTTCATCGGCGCCGAAGCtcat                                                                                                                        >  1:2017181/1‑126 (MQ=255)
                                                        tGAGCAGGTCGGCATCGTCAGGTCGGCTTTCCAGCTCGACCAGTTGCTCGGACAATTGCTCAAGTATTTCGCCGGCCTCTACCAGGAAATCCTGGAGGATTTCTTCATCGGCGCCGAAGCTCATCAAACGTGCTCCTTAAAAACCacttc                                                                                              <  1:1667538/150‑5 (MQ=255)
                                                         gAGCAGGTCGGCATCGTCAGGTCGGCTTTCCAGCTCGACCAGTTGCTCGGACAATTGCTCAAGTATTTCGCCGGCCTCTACCAGGAAATCCTGGAGGATTTCTTCATCGGCGCCGAAGCTCATCAAACGTGCTCCTTAAAAACCAAGGCt                                                                                             <  2:158959/150‑1 (MQ=255)
                                                         gAGCAGGTCGGCATCGTCAGGTCGGCTTTCCAGCTCGACCAGTTGCTCGGACAATTGCTAAAGTATTTGGGCGGCCTCTACCAGGAAATCCTGGAGTATTTCTGCATCGGCTCCGAATCTCATCAAACGTGCTCCTTAAAAACCGAGGCt                                                                                             <  2:882213/150‑1 (MQ=255)
                                                         gAGCAGGTCGGCACCGTCAGGTCGGCTTTCCAGCTCCACCAGTTGCTCGGACAATTGCTCAAGTATTTCGCCGGCCTTTACCATGAAATCCTGGAGGATTTCTTCATCGGCGCCGAAGCTCATCAAACGTGCTCCTTAAAAACCAAGGCt                                                                                             <  2:1332346/150‑1 (MQ=255)
                                                               gTCGGCATCGTCAGGTCGGCTTTCCAGCTCGACCAGTTGCTCGGACAATTGCTCAAGTATTTCGCCGGCCTCTACCAGGAAATCCTGGAGGATTTCTTCATCGGCGCCGAAGCTCATCAAACGTGCTCCTTAAAAACCAAGGCTGGAcag                                                                                       <  1:334762/150‑1 (MQ=255)
                                                                 cGGCATCGTCAGGTCGGCTTTCCAGCTCGACCAGTTGCTCGGACAATTGCTCAAGTATTTCGCCGGCCTCTACCAGGAAATCCTGGAGGATTTCTTCATCGGCGCCGAAGCTCATCAAACGTGCTCCTTAAAAACCAAGGCTGGAcagca                                                                                     <  2:911861/150‑1 (MQ=255)
                                                                   gCATCGTCAGGTCGGCTTTCCAGCTCGACCGGTTGCTCGGACAATTGCTCAAGTATTTCGCCGGCCGCTACCAGGAAATCCTGGAGGATTTCTTCATCGGCGCCGAAGCTCATCAAACGTGCTCCTTAAAAACCAAGGCTGGACAGCAGa                                                                                   >  1:1818100/1‑150 (MQ=255)
                                                                        gTCAGGTCGGCTTTCCAGCTCGACCAGTTGCTCGGACAATTGCTCAAGTATTTCGCCGGCCTCTACCAGGAAATCCTGGAGGATTTCTTCATCGGCGCCGAAGCTCATCAAACGTGCTCCTTAAAAACCAAGGCTGGACAGCAGatcatc                                                                              <  1:207469/150‑1 (MQ=255)
                                                                          cAGGTCGGCTTTCCAGCTCGACCAGTTGCTCGGACAATTGCTCAAGTATTTCGCCGGCCTCTACCAGGAAATCCTGGAGGATTTCTTCATCGGCGCCGAAGCTCATCAAACGTGCTCCTTAAAAACCAAGGCTGGACAGCAGATCATCGa                                                                            <  1:978763/150‑1 (MQ=255)
                                                                              tCGGCTTTCCAGCTCGACCAGTTGCTCGGACAATTGCTCAAGTATTTCGCCGGCCTCTACCAGGAAATCCTGGAGGATTTCTTCATCGGCGCCGAAGCTCATCAAACGTGCTCCTTAAAAACCAAGGCTGGACAGCAGATCATCGACATc                                                                        >  2:1427516/1‑150 (MQ=255)
                                                                                 gCTTTCCAGCTCGACCAGTTGCTCGGACAATTGCTCAAGTATTTCGCCGGCCTCTACCAGGAAATCCTGGAGGATTTCTTCATCGGCGCCGAAGCTCATCAAACGTGCTCCTTAAAAACCAAGGCTGGACAGCAGATCATCGACATCGTc                                                                     <  1:1048135/150‑1 (MQ=255)
                                                                                 gCTTTCCAGCTCGACCAGTTGCTCGGACAATTGCTCAAGTATTTCGCCGGCCTCTACCAGGAAATCCTGGAGGATTTCTTCATCGGCGCCGAAGCTCATCAAACGTGCTCCTTAAAAACCAAGGCTGGACAGCAGATCATCGACATCGTc                                                                     <  1:415885/150‑1 (MQ=255)
                                                                                 gCTTTCCAGCTCGACCAGTTGCTCGGACAATTGCTCAAGTATTTCGCCGGCCTCTACCAGGAAATCCTGGAGGATTTCTTCATCGGCGCCGAAGCTCATCAAACGTGCTCCTTAAAAACCAAGGCTGGACAGCAGATCATCGACATCGTc                                                                     <  1:257570/150‑1 (MQ=255)
                                                                                 gCTTTCCAGCTCGACCAGTTGCTCGGACAATTGCTCAAGTATTTCGCCGGCCTCTACCAGGAAATCCTGGAGGATTTCTTCATCGGCGCCGAAGCTCATCAAACGTGCTCCTTAAAAACCAAGGCTGGACAGCAGATCATCGACATCGTc                                                                     <  2:798302/150‑1 (MQ=255)
                                                                                   tttCCAGCTCGACCAGTTGCTCGGACAATTGCTCAAGTATTTCGCCGGCCTCTACCAGGAAATCCTGGAGGATTTCTTCATCGGCGCCGAAGCTCATCAAACGTGCTCCTTAAAAACCAAGGCTGGACAGCAGATCATCGACATCGTCCt                                                                   >  1:1989538/1‑150 (MQ=255)
                                                                                       cAGCTCGACCAGTTGCTCGGACAATTGCTCAAGTATTTCGCCGGCCTCTACCAGGAAATCCTGGAGGATTTCTTCATCGGCGCCGAAGCTCATCAAACGTGCTCCTTAAAAACCAAGGCTGGACAGCAGATCATCGACATCGTCCTGAcc                                                               >  2:2400487/1‑150 (MQ=255)
                                                                                         gCTCGACCAGTTGCTCGGACAATTGCTCAAGTATTTCGCCGGCCTCTACCAGGAAATCCTGGAGGATTTCTTCATCGGCGCCGAAGCTCATCAAACGTGCTCCTTAAAAACCAAGGCTGGACAGCAGATCATCGACATCGTCCTGACCgg                                                             >  2:2250794/1‑150 (MQ=255)
                                                                                            cGACCAGTTGCTCGGACAATTGCTCAAGTATTTCGCCGGCCTCTACCAGGAAATCCTGGAGGATTTCTTCATCGGCGCCGAAGCTCATCAAACGTGCTCCTTAAAAACCAAGGCTGGACAGCAGATCATCGACATCGTCCTGACCGGaca                                                          <  2:1989538/150‑1 (MQ=255)
                                                                                               ccAGTTGCTCGGACAATTGCTCAAGTAGGTCGCCGGCCTCTACCAGGAAATCCTGGAGGATTTCTTCATCGGCGCCGAAGCTAATCAAACGTGCTTCCTCGCAACCACGCAGTGCCACCAGAACATCGACATCGTCCTGACCGGACacga                                                       >  2:530275/1‑150 (MQ=255)
                                                                                                 aGTTGCTCGGACAATTGCTCAAGTATTTCGCCGGCCTCTACCAGGAAATCCTGGAGGATTTCTTCATCAGCGCCGAAGCTCATCAAACGTGCTCCTTAAAAACCAAGGCTGGACAGCAGa                                                                                   <  1:512155/120‑1 (MQ=255)
                                                                                                 aGTTGCTCGGACAATTGCTCAAGTATTTCGCCGGCCTCTACCAGGAAATCCTGGAGGATTTCTTCATCAGCGCCGAAGCTCATCAAACGTGCTCCTTAAAAACCAAGGCTGGACAGCAGa                                                                                   >  2:512155/1‑120 (MQ=255)
                                                                                                  gTTGCTCGGACAATTGCTCAAGTATTTCGCCGGCCTCTACCAGGAAATCCTGGAGGATTTCTTCATCGGCGCCGAAGCTCATCAAACGTGCTCCTTAAAAACCAAGGCTGGACAGCAGATCATCGACATCGTCCTGACCGGacac                                                         <  2:503515/145‑1 (MQ=255)
                                                                                                  gTTGCTCGGACAATTGCTCAAGTATTTCGCCGGCCTCTACCAGGAAATCCTGGAGGATTTCTTCATCGGCGCCGAAGCTCATCAAACGTGCTCCTTAAAAACCAAGGCTGGACAGCAGATCATCGACATCGTCCTGACCGGacac                                                         >  1:503515/1‑145 (MQ=255)
                                                                                                  gTTGCTCGGACAATTGCTCAAGTATTTCGCCGGCCTCTACCAGGAAATCCTGGAGGATTTCTTCATCGGCGCCGAAGCTCATCAAACGTGCTCCTTAAAAACCAAGGCTGGACAGCAGATCATCGACATCGTCCTGACCGGACACGACGt                                                    >  2:1048454/1‑150 (MQ=255)
                                                                                                       tCGGACAATTGCTCAAGTATTTCGCCGGCCTCTACCAGGAAATCCTGGAGGATTTCTTCATCGGCGCCGAAGCTCATCGAACGTGCTCCTTAAAAACCAAGGCTGGa                                                                                          <  1:1205487/107‑1 (MQ=255)
                                                                                                       tCGGACAATTGCTCAAGTATTTCGCCGGCCTCTACCAGGAAATCCTGGAGGATTTCTTCATCGGCGCCGAAGCTCATCGAACGTGCTCCTTAAAAACCAAGGCTGGa                                                                                          >  2:1205487/1‑107 (MQ=255)
                                                                                                       tCGGACAATTGCTCAAGTATTTCGCCGGCCTCTACCAGGAAATCCTGGAGGATTTCTTCATCGGCGCCGAAGCTCATCAAACGTGCTCCTTAAAAACCAAGCCTGGAGAGAAGATGATCGACATCGTCCTGACCGGACACGACGTCTTCa                                               >  2:1004082/1‑150 (MQ=255)
                                                                                                             aaTTGCTCAAGTATTTCGCCGGCCTCTACCAGGAAATCCTGGAGGATTTCTTCATCGGCGCCGAAGCTCATCAAACGTGCTCCTTAAAAACCAAGGCTGGACAGCAGATCATCGACATCGTCCTGACCGGACACGACGTCTTCACGCTTa                                         >  1:509780/1‑150 (MQ=255)
                                                                                                                tGCTCAAGTATTTCGCCGGCCTCTACCAGGAAATCCTGGAGGATTTCTTGATCGGCGCCGAAGCTCATCAAACGTGCTCCTTAAAAACCAAGGCTGGACAGCAGATCATCGACATCGTCCTGACCGGACACGACGTCTTCACGCTTATCg                                      <  1:2357865/150‑1 (MQ=255)
                                                                                                                tGCTCAAGTATTTCGCCGGCCTCTACCAGGAAATCCTGGAGGATTTCTTCATCGGCGCCGAAGCTCATCAAACGTGCTCCTTAAAAACCAAGGCTGGACAGCAGATCATCGACATCGTCCTGACCGGACACGACGTCTTCACGCTTATCg                                      <  1:2250794/150‑1 (MQ=255)
                                                                                                                   tCAAGTATTTCGCCGGCCTCTACCAGGAAATCCTGGAGGATTTCTTCATCTGCGCCGAAGCTCATCAAACGTGCTCCTTAAAAACCAAGGCTGGACAGCAGATCATCGACATCGTCCTGACCGGACACGACGTCTTCACGCTTATCGGCa                                   >  2:946745/1‑150 (MQ=255)
                                                                                                                   tCAAGTATTTCGCCGGCCTCTACCAGGAAATCCTGGAGGATTTCTTCATCGGCGCCGAAGCTCATCAAACGTGCTCCTTAAAAACCAAGGCTGGACAGCAGATCATCGACATCGTCCTGACCGGACACGACGTCTTCACGCTTATCGGCa                                   >  1:672565/1‑150 (MQ=255)
                                                                                                                    cAAGTATTTCGCCGGCCTCTACCAGGAAATCCTGGAGGATTTCTTCATCGGCGCCGAAGCTCATCAAACGTGCTCCTTAAAAACCAAGGCTGGACAGCAGATCATCGACATCGTCCTGACCGGACACGACGTCTTCACGCTTATCGGCAt                                  <  1:1427516/150‑1 (MQ=255)
                                                                                                                     aaGTATTTCGCCGGCCTCTACCAGGAAATCCTGGAGGAt                                                                                                                                                <  2:2173234/39‑1 (MQ=255)
                                                                                                                     aaGTATTTCGCCGGCCTCTACCAGGAAATCCTGGAGGAt                                                                                                                                                >  1:2173234/1‑39 (MQ=255)
                                                                                                                       gTATTTCGCCGGCCTCTACCAGGAAATCCTGGAGGATTTCTTCATCGGCGCCGAAGCTCATCAAACGTGCTCCTTAAAAAc                                                                                                    <  2:1094896/81‑1 (MQ=255)
                                                                                                                       gTATTTCGCCGGCCTCTACCAGGAAATCCTGGAGGATTTCTTCATCGGCGCCGAAGCTCATCAAACGTGCTCCTTAAAAAc                                                                                                    >  1:1094896/1‑81 (MQ=255)
                                                                                                                       gTATTTCGCCGGCCTCTACCAGGAAATCCTGGAGGATTTCTTCATCGGCGCCGAAGCTCATCAAACGTGCTCCTTAAAAACCAAGGCTGGACAGCAGATCATCGACATCGTCCTGACCGGACACGACGTCTTCACGCTTATCGGCATGaa                               >  1:836663/1‑150 (MQ=255)
                                                                                                                          tttCGCCGGCCTCTACCAGGAAATCCTGGAGGATTTCTTCATCGGCGCCGAAGCTCATCAAACGTGCTCCTTAAAAACCAAGGCTGGACAGCAGATCATCGACATCGTCCTGACCGGACACGACGTCTTCACGCTTATCGGCATGAATCt                            >  2:2298461/1‑150 (MQ=255)
                                                                                                                             cGCCGGCCTCTACCAGGAAATCCTGGAGGATTTCTTCATCGGCGCCGAAGCTCATCAAACGTGCTCCTTAAAAACCAAGGCTGGACAGCAGATCATCGACATCGTCCTGACCGGACACGACGTCTTCACGCTTATCGGCATGAATCTGCg                         <  1:1048454/150‑1 (MQ=255)
                                                                                                                                  gCCTCTACCAGGAAATCCTGGAGGATTTCTTCATCGGCGCCGAAGCTCATCAAACGTGCTCCTTAAAAACCAAGGCTg                                                                                            >  2:2295937/1‑78 (MQ=255)
                                                                                                                                  gCCTCTACCAGGAAATCCTGGAGGATTTCTTCATCGGCGCCGAAGCTCATCAAACGTGCTCCTTAAAAACCAAGGCTg                                                                                            <  1:2295937/78‑1 (MQ=255)
                                                                                                                                   ccTCTACCAGGAAATCCTGGAGGATTTCTTCATCGGCGCCGAAGCTCATCAAACGTGCTCCTTAAAAACCAAGGCTGGACAGCAGATCATCGACATCGTCCTGACCGGACACGACGTCTTCACGCTTATCGGCATGAATCTGCGGACCtt                   >  1:1425565/1‑150 (MQ=255)
                                                                                                                                     tctACCAGGAAATCCTGGAGGATTTCTTCATCGGCGCCGAAGCTCATCAAACGTGCTCCTTAAAAACCAAGGCTg                                                                                            <  1:1324962/75‑1 (MQ=255)
                                                                                                                                     tctACCAGGAAATCCTGGAGGATTTCTTCATCGGCGCCGAAGCTCATCAAACGTGCTCCTTAAAAACCAAGGCTg                                                                                            >  2:1324962/1‑75 (MQ=255)
                                                                                                                                       tACCAGGAAATCCTGGAGGATTTCTTCATCGGCGCCGAAGCTCATCAAACGTGCTCCTTAAAAACCAAGGCTGGACAGCAGATCATCGACATCGTCCTGACCGGACACGACGTCTTCa                                               >  1:1366381/1‑118 (MQ=255)
                                                                                                                                       tACCAGGAAATCCTGGAGGATTTCTTCATCGGCGCCGAAGCTCATCAAACGTGCTCCTTAAAAACCAAGGCTGGACAGCAGATCATCGACATCGTCCTGACCGGACACGACGTCTTCa                                               <  2:1366381/118‑1 (MQ=255)
                                                                                                                                         ccAGGAAATCCTGGAGGATTTCTTCATCGGCGCCGAAGCTCATCAAACGTGCTCCTTAAAAACCAAGGCTGGACAGCAGATCATCGACATCGTCCTGACCGGACACGACGTCTTCACGCTTATCGGCATGAATCTGCGGACCTTCAcccc             >  1:2268143/1‑150 (MQ=255)
                                                                                                                                          cAGGAAATCCTGGAGGATTTCTTCATCGGCGCCGAAGCTCATCAAACGTGCTCCTTAAAAACCAAGGCTGGACAGCAGATCATCGACATCGTCCTGACCGGACACGACGTCTTCACGCTTATCGGc                                    <  1:1157747/126‑1 (MQ=255)
                                                                                                                                          cAGGAAATCCTGGAGGATTTCTTCATCGGCGCCGAAGCTCATCAAACGTGCTCCTTAAAAACCAAGGCTGGACAGCAGATCATCGACATCGTCCTGACCGGACACGACGTCTTCACGCTTATCGGc                                    >  2:1157747/1‑126 (MQ=255)
                                                                                                                                          cAGGAAATCCTGGAGGATTTCTTCATCGGCGCCGAAGCTCATCAAACGTGCTCCTTAAAAACCAAGGCTGGACAGCAGATCATCGACATCGTCCTGACCGGACACGACGTCTTCACGCTTATCGGCATGAATCTGCGGACCTTCACCCCg            >  1:2394840/1‑150 (MQ=255)
                                                                                                                                          cAGGAAATCCTGGAGGATTTCTTCATCGGCGCCGAAGCTCATCAAACGTGCTCCTTAAAAACCAAGGCTGGACAGCAGATCATCGACATCGTCCTGACCGGACACGACGTCTTCACGCTTATCGGCATGAATCTGCGGACCTTCACCCCg            >  1:2267429/1‑150 (MQ=255)
                                                                                                                                            ggAAATCCTGGAGGATTTCTTCATCGGCGCCGAAGCTCATCAAACGTGCTCCTTAAAAACCAAGGCTGGACAGCAGATCATCGACATCGTCCTGACCGGACACGACGTCTTCACGCTTATCGGCATGAAt                              >  2:1711715/1‑130 (MQ=255)
                                                                                                                                            ggAAATCCTGGAGGATTTCTTCATCGGCGCCGAAGCTCATCAAACGTGCTCCTTAAAAACCAAGGCTGGACAGCAGATCATCGACATCGTCCTGACCGGACACGACGTCTTCACGCTTATCGGCATGAAt                              <  1:1711715/130‑1 (MQ=255)
                                                                                                                                            agAAATCCTCCAGGATTTCTTCATCGGCGCCGAAGCTCATCAAACGTGCTCCTTAAAAACCAAGGCTGGACAGCAGATCATCGACATCGTc                                                                     >  1:1970040/2‑91 (MQ=255)
                                                                                                                                            agAAATCCTCCAGGATTTCTTCATCGGCGCCGAAGCTCATCAAACGTGCTCCTTAAAAACCAAGGCTGGACAGCAGATCATCGACATCGTc                                                                     <  2:1970040/90‑1 (MQ=255)
                                                                                                                                                  ccTGGAGGATTTCTTCATCGGCGCCGAAGCTCATCAAACGTGCTCCTTAAAAACCAAGGCTGGACAGCAGATCATCGACATCGTCCTGACCGGACACGACGTCTTCACGCTTATCGGCATGAATCTGCGGACCTTCACCCCGAGTCGGAt    <  2:886562/150‑1 (MQ=255)
                                                                                                                                                    tGGAGGATTTCTTCATCGGCGCCGAAGCTCATCAAACGTGCTCCTTAAAAACCAAGGCTGGACAGCAGATCATCGACATCGTCCTGACCGGACACGACGTCt                                                  <  2:737594/102‑1 (MQ=255)
                                                                                                                                                    tGGAGGATTTCTTCATCGGCGCCGAAGCTCATCAAACGTGCTCCTTAAAAACCAAGGCTGGACAGCAGATCATCGACATCGTCCTGACCGGACACGACGTCt                                                  >  1:737594/1‑102 (MQ=255)
                                                                                                                                                    tGGAGGATTTCTTCATCGGCGCCGAAGCTCATCAAACGTGCTCCTTAAAAACCAAGGCTGGACAGCAGATCATCGACATCGTCCTGACCGGACACGACGTCTTCACGCTTATCGGCATGAATCTGCGGACCTTCACCCCGAGTCGGATGt  >  1:2042903/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                    
CAGCTGAAGGAAGCCGGCGCCCCCTTTTACAGTGTGGAAACCGCGAAAGATCGCATTGAGCAGGTCGGCATCGTCAGGTCGGCTTTCCAGCTCGACCAGTTGCTCGGACAATTGCTCAAGTATTTCGCCGGCCTCTACCAGGAAATCCTGGAGGATTTCTTCATCGGCGCCGAAGCTCATCAAACGTGCTCCTTAAAAACCAAGGCTGGACAGCAGATCATCGACATCGTCCTGACCGGACACGACGTCTTCACGCTTATCGGCATGAATCTGCGGACCTTCACCCCGAGTCGGATGT  >  NC_002947/4930037‑4930334

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