Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,930,187 G→C 9.1% L10L (CTC→CTG cheA ← chemotaxis histidine kinase CheA

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,930,1870GC 9.1% 271.5 / 15.0 88L10L (CTC→CTGcheAchemotaxis 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 = 8.32e-01

AGCTGAAGGAAGCCGGCGCCCCCTTTTACAGTGTGGAAACCGCGAAAGATCGCATTGAGCAGGTCGGCATCGTCAGGTCGGCTTTCCAGCTCGACCAGTTGCTCGGACAATTGCTCAAGTATTTCGCCGGCCTCTACCAGGAAATCCTGGAGGATTTCTTCATCGGCGCCGAAGCTCATCAAACGTGCTCCTTAAAAACCAAGGCTGGACAGCAGATCATCGACATCGTCCTGACCGGACACGACGTCTTCACGCTTATCGGCATGAATCTGCGGACCTTCACCCCGAGTCGGATGTTT  >  NC_002947/4930038‑4930336
                                                                                                                                                     |                                                                                                                                                     
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                                                                                                                                               >  2:257570/1‑150 (MQ=255)
        gaagCCGGCGCCCCCTTTTACAGTGTGGAAACCGCGAAAGATCGCATTGAGCAGGTCGGCATCGTCAGGTCGGCTTTCCAGCTCGACCAGTTGCTCGGACAATTGCTCAAGTATTTCGCCGGCCTCTACCAGGAAATCCTGGAGGATttc                                                                                                                                               <  1:618094/150‑1 (MQ=255)
          agCCGGCGCCCCCTTTTACAGTGTGGAAACCGCGAAAGATCGCATTGAGCAGGTCGGCATCGTCAGGTCGGCTTTCCAGCTCGACCAGTTGCTCGGACAATTGCTCAAGTATTTCGCCGGCCTCTACCAGGAAATCCTCCAGGATttcct                                                                                                                                             <  2:151794/150‑3 (MQ=255)
           gCCGGCGCCCCCTTTTACAGTGTGGAAACCGCGAAAGATCGCATTGAGCAGGTCGGCATCGTCAGGTCGGCTTTCCAGCTCGACCAGTTGCTCGGACAATTGCTCAAGTATTTCGCCGGCCTCTACCAGGAAATCCTCCAGGATttcctg                                                                                                                                            <  2:1414040/150‑4 (MQ=255)
           gCCGGCGCCCCCTTTTACAGTGTGGAAACCGCGAAAGATCGCATTGAGCAGGTCGGCATCGTCAGGTCGGCTTTCCAGCTCGACCAGTTGCTCGGACAATTGCTCAAGTATTTCGCCGGCCTCTACCAGGAAATCCTCCAGGATttcctg                                                                                                                                            <  2:930952/150‑4 (MQ=255)
            ccGGCGCCCCCTTTTACAGTGTGGAAACCGCGAAAGATCGCATTGAGCAGGTCGGCATCGTCAGGTCGGCTTTCCAGCTCGACCAGTTGCTCGGACAATTGCTCAAGTATTTCGCCGGCCTCTACCAGGAAATCCTggagga                                                                                                                                                   >  2:444737/1‑142 (MQ=255)
            ccGGCGCCCCCTTTTACAGTGTGGAAACCGCGAAAGATCGCATTGAGCAGGTCGGCATCGTCAGGTCGGCTTTCCAGCTCGACCAGTTGCTCGGACAATTGCTCAAGTATTTCGCCGGCCTCTACCAGGAAATCCTggagga                                                                                                                                                   <  1:444737/142‑1 (MQ=255)
              ggCGCCCCCTTTTACAGTGTGGAAACCGCGAAAGATCGCATTGAGCAGGTCGGCATCGTCAGGTCGGCTTTCCAGCTCGACCAGTTGCTCGGACAATTGCTCAAGTATTTCGCCGGCCTCTACCAGGAAATCCTGGAGGATTTCTTCATc                                                                                                                                         <  1:2359181/150‑1 (MQ=255)
              ggCGCCCCCTTTTACAGTGTGGAAACCGCGAAAGATCGCATTGAGCAGGTCGGCATCGTCAGGTCGGCTTTCCAGCTCGACCAGTTGCTCGGACAATTGCTCAAGTATTTCGCCGGCCTCTACCAGGAAATCCTGGAGGATTTCTTCATc                                                                                                                                         <  2:1032480/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                                                                                                                          >  1:2017181/1‑126 (MQ=255)
                                                     aTTGAGCAGGTCGGCATCGTCAGGTCGGCTTTCCAGCTCGACCAGTTGCTCGGACAATTGCTCAAGTATTTCGCCGGCCTCTACCAGGAAATCCTGGAGGATTTCTTCATCGGCGCCGAAGCtcat                                                                                                                          <  2:2017181/126‑1 (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:415885/150‑1 (MQ=255)
                                                                                gCTTTCCAGCTCGACCAGTTGCTCGGACAATTGCTCAAGTATTTCGCCGGCCTCTACCAGGAAATCCTGGAGGATTTCTTCATCGGCGCCGAAGCTCATCAAACGTGCTCCTTAAAAACCAAGGCTGGACAGCAGATCATCGACATCGTc                                                                       <  2:798302/150‑1 (MQ=255)
                                                                                gCTTTCCAGCTCGACCAGTTGCTCGGACAATTGCTCAAGTATTTCGCCGGCCTCTACCAGGAAATCCTGGAGGATTTCTTCATCGGCGCCGAAGCTCATCAAACGTGCTCCTTAAAAACCAAGGCTGGACAGCAGATCATCGACATCGTc                                                                       <  1:257570/150‑1 (MQ=255)
                                                                                gCTTTCCAGCTCGACCAGTTGCTCGGACAATTGCTCAAGTATTTCGCCGGCCTCTACCAGGAAATCCTGGAGGATTTCTTCATCGGCGCCGAAGCTCATCAAACGTGCTCCTTAAAAACCAAGGCTGGACAGCAGATCATCGACATCGTc                                                                       <  1:1048135/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                                                           >  1:503515/1‑145 (MQ=255)
                                                                                                 gTTGCTCGGACAATTGCTCAAGTATTTCGCCGGCCTCTACCAGGAAATCCTGGAGGATTTCTTCATCGGCGCCGAAGCTCATCAAACGTGCTCCTTAAAAACCAAGGCTGGACAGCAGATCATCGACATCGTCCTGACCGGacac                                                           <  2:503515/145‑1 (MQ=255)
                                                                                                 gTTGCTCGGACAATTGCTCAAGTATTTCGCCGGCCTCTACCAGGAAATCCTGGAGGATTTCTTCATCGGCGCCGAAGCTCATCAAACGTGCTCCTTAAAAACCAAGGCTGGACAGCAGATCATCGACATCGTCCTGACCGGACACGACGt                                                      >  2:1048454/1‑150 (MQ=255)
                                                                                                      tCGGACAATTGCTCAAGTATTTCGCCGGCCTCTACCAGGAAATCCTGGAGGATTTCTTCATCGGCGCCGAAGCTCATCGAACGTGCTCCTTAAAAACCAAGGCTGGa                                                                                            >  2:1205487/1‑107 (MQ=255)
                                                                                                      tCGGACAATTGCTCAAGTATTTCGCCGGCCTCTACCAGGAAATCCTGGAGGATTTCTTCATCGGCGCCGAAGCTCATCGAACGTGCTCCTTAAAAACCAAGGCTGGa                                                                                            <  1:1205487/107‑1 (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                                                                                                                                                  >  1:2173234/1‑39 (MQ=255)
                                                                                                                    aaGTATTTCGCCGGCCTCTACCAGGAAATCCTGGAGGAt                                                                                                                                                  <  2:2173234/39‑1 (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                                                 <  2:1366381/118‑1 (MQ=255)
                                                                                                                                      tACCAGGAAATCCTGGAGGATTTCTTCATCGGCGCCGAAGCTCATCAAACGTGCTCCTTAAAAACCAAGGCTGGACAGCAGATCATCGACATCGTCCTGACCGGACACGACGTCTTCa                                                 >  1:1366381/1‑118 (MQ=255)
                                                                                                                                        ccAGGAAATCCTGGAGGATTTCTTCATCGGCGCCGAAGCTCATCAAACGTGCTCCTTAAAAACCAAGGCTGGACAGCAGATCATCGACATCGTCCTGACCGGACACGACGTCTTCACGCTTATCGGCATGAATCTGCGGACCTTCAcccc               >  1:2268143/1‑150 (MQ=255)
                                                                                                                                         cAGGAAATCCTGGAGGATTTCTTCATCGGCGCCGAAGCTCATCAAACGTGCTCCTTAAAAACCAAGGCTGGACAGCAGATCATCGACATCGTCCTGACCGGACACGACGTCTTCACGCTTATCGGc                                      >  2:1157747/1‑126 (MQ=255)
                                                                                                                                         cAGGAAATCCTGGAGGATTTCTTCATCGGCGCCGAAGCTCATCAAACGTGCTCCTTAAAAACCAAGGCTGGACAGCAGATCATCGACATCGTCCTGACCGGACACGACGTCTTCACGCTTATCGGc                                      <  1:1157747/126‑1 (MQ=255)
                                                                                                                                         cAGGAAATCCTGGAGGATTTCTTCATCGGCGCCGAAGCTCATCAAACGTGCTCCTTAAAAACCAAGGCTGGACAGCAGATCATCGACATCGTCCTGACCGGACACGACGTCTTCACGCTTATCGGCATGAATCTGCGGACCTTCACCCCg              >  1:2394840/1‑150 (MQ=255)
                                                                                                                                         cAGGAAATCCTGGAGGATTTCTTCATCGGCGCCGAAGCTCATCAAACGTGCTCCTTAAAAACCAAGGCTGGACAGCAGATCATCGACATCGTCCTGACCGGACACGACGTCTTCACGCTTATCGGCATGAATCTGCGGACCTTCACCCCg              >  1:2267429/1‑150 (MQ=255)
                                                                                                                                           ggAAATCCTGGAGGATTTCTTCATCGGCGCCGAAGCTCATCAAACGTGCTCCTTAAAAACCAAGGCTGGACAGCAGATCATCGACATCGTCCTGACCGGACACGACGTCTTCACGCTTATCGGCATGAAt                                <  1:1711715/130‑1 (MQ=255)
                                                                                                                                           ggAAATCCTGGAGGATTTCTTCATCGGCGCCGAAGCTCATCAAACGTGCTCCTTAAAAACCAAGGCTGGACAGCAGATCATCGACATCGTCCTGACCGGACACGACGTCTTCACGCTTATCGGCATGAAt                                >  2:1711715/1‑130 (MQ=255)
                                                                                                                                           agAAATCCTCCAGGATTTCTTCATCGGCGCCGAAGCTCATCAAACGTGCTCCTTAAAAACCAAGGCTGGACAGCAGATCATCGACATCGTc                                                                       <  2:1970040/90‑1 (MQ=255)
                                                                                                                                           agAAATCCTCCAGGATTTCTTCATCGGCGCCGAAGCTCATCAAACGTGCTCCTTAAAAACCAAGGCTGGACAGCAGATCATCGACATCGTc                                                                       >  1:1970040/2‑91 (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)
                                                                                                                                                     gaggaTTTCTTCATCGGCGCCGAAGCTCATCAAACGTTTTTCTTAAAAACCAAGGCTAGACCGCAGATCATCGACATCGTCCTGACCGGACACGCGGTCTTCACCCTTTTCGGCTTGAATCTGCGCACCTTCACCGCTAGTCGGATGttt  <  2:531898/150‑1 (MQ=255)
                                                                                                                                                     gaggaTTTCTTCATCGGCGCCGAAGCTCATCAAACGTGCTCCTTAAAAACCAAgg                                                                                                 <  1:2377342/55‑1 (MQ=255)
                                                                                                                                                     gaggaTTTCTTCATCGGCGCCGAAGCTCATCAAACGTGCTCCTTAAAAACCAAgg                                                                                                 >  2:2377342/1‑55 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
AGCTGAAGGAAGCCGGCGCCCCCTTTTACAGTGTGGAAACCGCGAAAGATCGCATTGAGCAGGTCGGCATCGTCAGGTCGGCTTTCCAGCTCGACCAGTTGCTCGGACAATTGCTCAAGTATTTCGCCGGCCTCTACCAGGAAATCCTGGAGGATTTCTTCATCGGCGCCGAAGCTCATCAAACGTGCTCCTTAAAAACCAAGGCTGGACAGCAGATCATCGACATCGTCCTGACCGGACACGACGTCTTCACGCTTATCGGCATGAATCTGCGGACCTTCACCCCGAGTCGGATGTTT  >  NC_002947/4930038‑4930336

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: