Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,930,142 G→C 5.5% 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,1420GC 5.5% 226.4 / 2.6 72L10L (CTC→CTG) cheAchemotaxis histidine kinase CheA
Reads supporting (aligned to +/- strand):  ref base G (37/31);  new base C (2/2);  total (39/33)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.49e-01

GCTGAAGGAAGCCGGCGCCCCCTTTTACAGTGTGGAAACCGCGAAAGATCGCATTGAGCAGGTCGGCATCGTCAGGTCGGCTTTCCAGCTCGACCAGTTGCTCGGACAATTGCTCAAGTATTTCGCCGGCCTCTACCAGGAAATCCTGGAGGATTTCTTCATCGGCGCCGAAGCTCATCAAACGTGCTCCTTAAAAACCAAGGCTGGACAGCAGATCATCGACATCGTCCTGACCGGACACGACGTCTTCACGCTTATCGGCATGAATCTGCGGACCTTCACCCCGAGTCGGATGTT  >  NC_002947/4929994‑4930290
                                                                                                                                                    |                                                                                                                                                    
gCTGAAGGAAGCCGGCGCCCCCTTTTACAGTGTGGAAACCGCGAAAGATCGCATTGAGCAGGTCGGCATCGTCAGGTCGGCTTTCCAGCTCGACCAGTTGCTCGGACAATTGCTCAAGTATTTCGCCGGCCTCTACCAGGAAATCCTgg                                                                                                                                                      <  8:54927/149‑1 (MQ=255)
 cTGAAGGAAGCCGGCGCCCCCTTTTACAGTGTGGAAACCGCGAAAGATCGCATTGAGCAGGTCGGCATCGTCAGGTCGGCTTTCCAGCTCGACCAGTTGCTCGGACAATTGCTCAAGTATTTCGCCGGCCTCTACCAGGAAATCCTgga                                                                                                                                                     <  5:174310/149‑1 (MQ=255)
 cTGAAGGAAGCCGCCGCCCCCTTTTACACTGTGGAAACCGCGAAAGATCGCATTGAGCAGGTCGGCAGCGTCAGGTCGGCTTTCCAGCTCGACCAGTTGCTCGGACAATTGCTCAAGTATTTCGCCGGCCTCTACCAGGAAATCCTgga                                                                                                                                                     <  2:108195/149‑1 (MQ=255)
  tGAAGGAAGCCGGCGCCCCCTTTTACAGTGTGGAAACCGCGAAAGATCGCATTGAGCAGGTCGGCATCGTCAGGTCGGCTTTCCAGCTCGACCAGTTGCTCGGACAATTGCTCAAGTATTTCGCCGGCCTCTACCAGGAAATCCTggag                                                                                                                                                    >  3:306620/1‑149 (MQ=255)
   gaaggaagCCGGCGCCCCCTTTTACAGTGTGGAAACCGCGAAAGATCGCATTGAGCAGGTCGGCATCGTCAGGTCGGCTTTCCAGCTCGACCAGTTGCTCGGACAATTGCTCAAGTATTTCGCCGGCCTCTACCAGGAAATCCTggagg                                                                                                                                                   >  6:32483/1‑149 (MQ=255)
       gaagCCGGCGCCTCCTTTTACAGTGTGGAAACCGCGAAAGATCGCATTGAGCAGGTCGGCATCGTCAGGTCGGCTTTCCAGCTCGACCAGTTGCTCGGACAATTGCTCAAGTATTTCGCCGGCCTCTACCAGGAAATCCTGGAGGAttt                                                                                                                                               >  3:71635/1‑149 (MQ=255)
        aagCCGGCGCCCCCTTTTACAGTGTGGAAACCGCGAAAGATCGCATTGAGCAGGTCGGCATCGTCAGGTCGGCTTTCCAGCTCGACCAGTTGCTCGGACAATTGCTCAAGTATTTCGCCGGCCTCTACCAGGAAATCCTGGAGGATttc                                                                                                                                              <  5:32483/149‑1 (MQ=255)
         agCCGGCGCCCCCTTTTACAGTGTGGAAACCGCGAAAGGTCGCATTGAGCAGGTCGGCATCGTCAGGTCGGCTTTCCAGCTCGACCAGTTGCTCGGACAATTGCTCAAGTATTTCGCCGGCCTCTACCAGGAAATCCTGGAGGATttct                                                                                                                                             <  3:17071/149‑1 (MQ=255)
          gCCGGCGCCCCCTTTTACAGTGTGGAAACCGCGAAAGATCGCATTGAGCAGGTCGGCATCGTCAGGTCGGCTTTCCAGCTCGACCAGTTGCTCGGACAATTGCTCAAGTATTTCGCCGGCCTCTACCAGGAAATCCTCCAGGATttcct                                                                                                                                            <  5:379275/149‑3 (MQ=255)
                   cccTTTTACAGTGTGGAAACCGCGAAAGATCGCATTGAGCAGGTCGGCATCGTCAGGTCGGCTTTCCAGCTCGACCAGTTGCTCGGACAATTGCTCAAGTATTTCGCCGGCCTCTACCAGGAAATCCTGGAGGATTTCTTCATCGgcgc                                                                                                                                   >  3:114436/1‑149 (MQ=255)
                   cccTTTTACAGTGTGGAAACCGCGAAAGATCGCATTGAGCAGGTCGGCATCGTCAGGTCGGCTTTCCAGCTAGACCAGTTGCTCGGACAATTGCTCAAGTATTTCGCCGGCCTCTACCAGGAAATGCTGGAGGATTTCTTCATCGgcgc                                                                                                                                   >  7:277059/1‑149 (MQ=255)
                      ttttACAGTGTGGAAACCGCGAAAGATCGCATTGAGCAGGTCGGCATCGTCAGGTCGGCTTTCCAGCTCGACCAGTTGCTCGGACAATTGCTCAAGTATTTCGCCGGCCTCTACCAGGAAATCCTCCAGGATttcct                                                                                                                                            >  1:59053/1‑135 (MQ=255)
                      ttttACAGTGTGGAAACCGCGAAAGATCGCATTGAGCAGGTCGGCATCGTCAGGTCGGCTTTCCAGCTCGACCAGTTGCTCGGACAATTGCTCAAGTATTTCGCCGGCCTCTACCAGGAAATCCTCCAGGATttcct                                                                                                                                            <  2:59053/137‑3 (MQ=255)
                          aCAGTGTGGAAACCGCGAAAGATCGCATTGAGCAGGTCGGCATCGTCAGGTCGGCTTTCCAGCTCGACCAGTTGCTCGGACAATTGCTCAAGTATTTCGCCGGCCTCTACCAGGAAATCCTGGAGGATTTCTTCATCGGCGCCGAAGc                                                                                                                             <  2:135883/148‑1 (MQ=255)
                              tgtgGAAACCGCGAAAGATCGCATTGAGCAGGTCGGCATCGTCAGGTCGGCTTTCCAGCTCGACCAGTTGCTCGGACAATTGCTCAAGTATTTCGCCGGCCTCTACCAGGAAATCCTGGAGGATTTCTTCATCGGCGCCGAAGCtcat                                                                                                                         >  1:260047/1‑148 (MQ=255)
                                  gAAACCGCGAAAGATCGCATTGAGCAGGTCGGCATCGTCAGGTCGGCTTTCCAGCTCGACCAGTTGCTCGGACAATTGCTCAAGTATTTCGCCGGCCTCTACCAGGAAATCCTGGAGGATTTCTTCATCGGCGCCGAAGCTCATCAAAc                                                                                                                    >  4:88663/1‑149 (MQ=255)
                                       cgcgAAAGATCGCATTGAGCAGGTCGGCATCGTCAGGTCGGCTTTCCAGCTCGACCAGTTGCTCGGACAATTGCTCAAGTATTTCGCCGGCCTCTACCAGGAAATCCTGGAGGATTTCTTCATCGGCGCCGAAGCTCATCAAACGTGCt                                                                                                               >  5:86664/1‑149 (MQ=255)
                                       cgcgAAAGATCGCATTGAGCAGGTCGGCATCGTCAGGTCGGCTTTCCAGCTCGACCAGTTGCTCGGACAATTGCTCAAGTATTTCGCCGGCCTCTACCAGGAAATCCTGGAGGATTTCTTCATCGGCGCAGAAGCTCTTCAAACGTTCt                                                                                                               >  8:145632/1‑149 (MQ=255)
                                            aaGATCGCATTGAGCAGGTCGGCATCGTCAGGTCGGCTTTCCAGCTCGACCAGTTGCTCGGACAATTGCTCAAGTATTTCGCCGGCCTCTACCAGGAAATCCTGGAGGATTTCTTCATCGGCGCCGAAGCTCATCAAACGTGCTCCTTa                                                                                                          <  6:280922/149‑1 (MQ=255)
                                                tCGCATTGAGCAGGTCGGCATCGTCAGGTCGGCTTTCCAGCTCGACCAGTTGCTCGGACAATTGCTCAAGTATTTCGCCGGCCTCTACCAGGAAATCCTGGAGGATTTCTTCATCGGCGCCGAAGCTCATCAAACGTGCTCCTTaaaaa                                                                                                      <  8:206380/149‑1 (MQ=255)
                                                           atgTCGGCATCGTCAGGTCGGCTTTCCAGCTCGACCAGTTGCTCGGACAATTGCTCAAGTATTTCGCCGGCCTCTACCAGGAAATCCTGGAGGATTTCTTCATCGGCGCCGAAGCTCATCAAACGTGCTCCTTAAAAACCAAGGCTGGa                                                                                           <  3:113057/147‑1 (MQ=255)
                                                            ggTCGGCATCGTCAGGTCGGCTTTCCAGCTCGACCAGTTGCTCGGACAATTGCTCAAGTATTTCGCCGGCCTCTACCAGGAAATCCTGGAGGATTTCTTCATCGGCGCCGAAGCTCATCAAACGTGCTCCTTAAAAACCAAGGCTGGAc                                                                                          <  2:359188/149‑1 (MQ=255)
                                                            ggTCGGCATCGTCAGGTCGGCTTTCCAGCTCGACCAGTTGCTCGGACAATTGCTCAAGTATTTCGCCGGCCTCTACCAGGAAATCCTGGAGGATTTCTTCATCGGCGCCGAAGCTCATCAAACGTGCTCCTTAAAAACCAAGGCTGGAc                                                                                          <  2:260047/149‑1 (MQ=255)
                                                                ggCATCGTCAGGTCGGCTTTCCAGCTCGACCAGTTGCTCGGACAATTGCTCAAGTATTTCGCCGGCCTCTACCAGGAAATCCTGGAGGATTTCTTCATCGGCGCCGAAGCTCATCAAACGTGCTCCTTAAAAACCAAGGCTGGAcagca                                                                                      <  3:88663/149‑1 (MQ=255)
                                                                ggCATCGTCAGGTCGGCTTTCCAGCTCGACCAGTTGCTCGGACAATTGCTCAAGTATTTCGCCGGCCTCTACCAGGAAATCCTGGAGGATTTCTTCATCGGCGCCGAAGCTCATCAAACGTGCTCCTTAAAAACCAAGGCTGGAcagca                                                                                      <  4:114436/149‑1 (MQ=255)
                                                                      gTCAGGTCGGCTTTCCAGCTCGACCAGTTGCTCGGACAATTGCTCAAGTATTTCGCCGGCCTCTACCAGGAAATCCTGGAGGATTTCTTCATCGGCGCCGAAGCTCATCAAACGTGCTCCTTAAAAACCAAGGCTGGACAGCAGatca                                                                                 <  3:159162/148‑1 (MQ=255)
                                                                        cAGGTCGGCTTTCCAGCTCGACCAGTTGCTCGGACAATTGCTCAAGTATTTCGCCGGCCTCTACCAGGAAATCCTGGAGGATTTCTTCATCGGCGCCGAAGCTCATCAAACGTGCTCCTTAAAAACCAAGGCTGGACAGCAGatcatc                                                                               >  2:306730/1‑148 (MQ=255)
                                                                           gTCGGCTTTCCAGCTCGACCAGTTGCTCGGACAATTGCTCAAGTATTTCGCCGGCCTCTACCAGGAAATCCTGGAGGATTTCTTCATCGGCGCCGAAGCTCATCAAACGTGCTCCTTAAAAACCAAGGCTGGACAGCAGat                                                                                   >  2:306321/1‑141 (MQ=255)
                                                                           gTCGGCTTTCCAGCTCGACCAGTTGCTCGGACAATTGCTCAAGTAATTCGCCGGCCTCTACCAGGAAATCCTGGAGGATTTCTTCATCGGCGCCGAAGCTCATCAAACGTGCTCCTTAAAAACCAAGGCTGGACAGCAGat                                                                                   <  1:306321/141‑1 (MQ=255)
                                                                             cGGCTTTCCAGCTCGACCAGTTGCTCGGACAATTGCTCAAGTATTTCGCCGGCCTCTACCAGGAAATCCTGGAGGATTTCTTCATCGGCGCCGAAGCTCATCAAACGTGCTCCTTAAAAACCAAGGCTGGACAGCAGatcatc                                                                               >  2:190904/1‑143 (MQ=255)
                                                                             cGGCTTTCCAGCTCGACCAGTTGCTCGGACAATTGCTCAAGTATTTCGCCGGCCTCTACCAGGAAATCCTGGAGGATTTCTTCATCGGCGCCGAAGCTCATCAAACGTGCTCCTTAAAAACCAAGGCTGGACAGCAGatcatc                                                                               <  1:190904/143‑1 (MQ=255)
                                                                              ggCTTTCCAGCTCGACCAGTTGCTCGGACAATTGCTCAAGTATTTCGCCGGCCTCTACCAGGAAATCCTGGAGGATTTCTTCATCGGCGCCGAAGCTCATCAAACGTGCTCCTTAAAAACCAAGGCTGGACAGCAGATCATCGACATCg                                                                        >  8:23107/1‑149 (MQ=255)
                                                                               gCTTTCCAGCTCGACCAGTTGCTCGGACAATTGCTCAAGTATTTCGCCGGCCTCTACCAGGAAATCCTGGAGGATTTCTTCATCGGCGCCGAAGCTCATCAAACGTGCTCCTTAAAAACCAAGGCTGGACAGCAGATCATCGACATCg                                                                        >  4:119976/1‑148 (MQ=255)
                                                                                cTTTCCAGCTCGACCAGTTGCTCGGACAATTGCTCAAGTATTTCGCCGGCCTCTACCAGGAAATCCTGGAGGATTTCTTCATCGGCGCCGAAGCTCATCAAACGTGCTCCTTAAAAACCAAGGCTGGACAGCAGATCATCGACATCGTc                                                                      <  6:163278/149‑1 (MQ=255)
                                                                                     cAGCTCGACCAGTTGCTCGGACAATTGCTCAAGTATTTCGCCGGCCTCTACCAGGAAATCCTGGAGGATTTCTTCATCGGCGCCGAAGCTCATCAAACGTGCTCCTTAAAAACCAAGGCTGGACAGCAGATCATCGACATCGTCCTGa                                                                  <  4:71635/148‑1 (MQ=255)
                                                                                     cAGCTCGACCAGTTGCTCGGACAATTGCTCAAGTATTTCGCCGGCCTCTACCAGGAAATCCTGGAGGATTTCTTCATCGGCGCCGAAGCTCATCAAACGTGCTCCTTAAAAACCAAGGCTGGACAGCAGATCATCGACATCGTCCTGAc                                                                 <  3:119976/149‑1 (MQ=255)
                                                                                           gaacAGTTGCTCGGACAATTGCTCAAGTATTTCGCCTGCCTCTACCAGGAAATCCTGGAGGATTTATTCATCGGCGCCGAAGCTCATCAAACGTGCTCCTTAAAAACCAAGGCTGGAGAGAAGATAATCTCATTAGTCATGACCGGac                                                            >  8:160831/4‑148 (MQ=255)
                                                                                              cAGTTGCTCGGACAATTGCTCAAGTATTTCGCCGGCCTCTACCAGGAAATCCTGGAGGATTTCTTCATCGGCGCCGAAGCTCATCAAACGTGCTCCTTAAAAACCAAGGCTGGACAGCAGATCATCGACATCGTCCTGACCGGACacga                                                        <  4:306620/149‑1 (MQ=255)
                                                                                               aGTTGCTCGGACAATTGCTCAAGTATTTCGCCGGCCTCTACCAGGAAATCCTGGAGGATTTCTTCATCGGCGCCGAAGCTCATCAAACGTGCTCCTTAAAAACCAAGGCTGGACAGCAGATCATCGACATCGTCCTGACCGGACacgac                                                       <  5:246584/149‑1 (MQ=255)
                                                                                                    cTCGGACAATTGCTCAAGTATTTCGCCGGCCTCTACCAGGAAATCCTGGAGGATTTCTTCATCGGCGCCGAAGCTCATCAAACGTGCTCCTTAAAAACCAAGGCTGGACAGCAGATCATCg                                                                              >  8:79856/1‑121 (MQ=255)
                                                                                                    cTCGGACAATTGCTCAAGTATTTCGCCGGCCTCTACCAGGAAATCCTGGAGGATTTCTTCATCGGCGCCGAAGCTCATCAAACGTGCTCCTTAAAAACCAAGGCTGGACAGCAGATCATCg                                                                              <  7:79856/121‑1 (MQ=255)
                                                                                                     tCGGACAATTGCTCAAGTATTTCGCCGGCCTCTACCAGGAAATCCTGGAGGATTTCTTCATCGGCGCCGAAGCTCATCAAACGTGCTCCTTAAAAACCAAGGCTGGACAGCAGATCATCGACATCGTCCTGACCGGac                                                            <  6:288259/138‑1 (MQ=255)
                                                                                                     tCGGACAATTGCTCAAGTATTTCGCCGGCCTCTACCAGGAAATCCTGGAGGATTTCTTCATCGGCGCCGAAGCTCATCAAACGTGCTCCTTAAAAACCAAGGCTGGACAGCAGATCATCGACATCGTCCTGACCGGac                                                            >  5:288259/1‑138 (MQ=255)
                                                                                                       tgACAATTGCTCAAGTATTTCGCCGGCCTCTACCAGGAAATCCTGGAGGATTTCTTCATCGGCGCCGAAGCTCATCAAACGTGCTCCTTAAAAACCAAGGCTGGACAGCAGATCATCGACATCGTCCTGACCGGACACGACGTCTTc                                                 <  6:86664/146‑1 (MQ=255)
                                                                                                        gACAATTGCTCAAGTATTTCGACGGCCTCTACCAGGATATCCAGGAGGAGGTCTTCATCGGCGCCGAAGCTCTTCAAGGGTGCTCCGAAAACAGCAAGGGGGGAGGGGCGAGCCTCGCGCTCGTCCTGACGGGACAGGACGTCGTcgcg                                              >  8:133548/1‑149 (MQ=255)
                                                                                                         aCAATTGCTCAAGTATTTCGCCGGCCTCTACCAGGAAATCCTGGAGGATTTCTTCATCGGCGCCGAAGCTCATCAAACGTGCTCCTTAAAAAGCAAGGCTGGACAGCAGATCATCGACATCGTCCTGACCGGACACTACGTCTTCACGc                                             >  3:45880/1‑149 (MQ=255)
                                                                                                          cAATTGCTCAAGTATTTCGCCGGCCTCTACCAGGAAATCCTGGAGGATTTCTTCATCGGCGCCGAAGCTCATCAAACGTGCTCCTTAAAAACCAAGGCTGGACAGCAGATCATCGACATCGTCCTGACCGGACACGACGTCTTCACGc                                             >  1:133780/1‑148 (MQ=255)
                                                                                                          cAATTGCTCAAGTATTTCGCCGGCCTCTACCAGGAAATCCTGGAGGATTTCTTCATCGGCGCCGAAGCTCATCAAACGTGCTCCTTAAAAACCAAGGCTGGACAGCAGATCATCGACATCGTCCTGACCGGACACGACGTCTTCACGc                                             <  2:133780/148‑1 (MQ=255)
                                                                                                           aaTTGCTCAAGTATTTCGCCGGCCTCTACCAGGAAATCCTGGAGGATTTCTTCATCGGCGCCGAAGCTCATCAAACGTGCTCCTTAAAAACCAAGGCTGGACAGCAGATCATCGACATCGTCCTGACCGGACACGACGTCTTCACGCt                                            >  6:340851/1‑148 (MQ=255)
                                                                                                             ttGCTCAAGTATTTCGCCGGCCTCTACCAGGAAATCCTGGAGGATTTCTTCATCGGCGCCGAAGCTCATCAAACGTGCTCCTTAAAAACCAAGGCTGGACAGCAGATCATCGACATCGTCCTGACCGGACACGACGTCTTCACGCTTAt                                         <  5:174920/149‑1 (MQ=255)
                                                                                                                 tCAAGTATTTCGCCGGCCTCTACCAGGAAATCCTGGAGGATTTCTTCATCGGCGCCGAAGCTCATCAAACGTGCTCCTTAAAAACCAAGGCTGGACAGCAGATCATCGACATCGTCCTGACCGGACACGACGTc                                                    <  4:100354/134‑1 (MQ=255)
                                                                                                                 tCAAGTATTTCGCCGGCCTCTACCAGGAAATCCTGGAGGATTTCTTCATCGGCGCCGAAGCTCATCAAACGTGCTCCTTAAAAACCAAGGCTGGACAGCAGATCATCGACATCGTCCTGACCGGACACGACGTc                                                    >  3:100354/1‑134 (MQ=255)
                                                                                                                 tCAAGTATTTCGCCGGCCTCTACCAGGAAATCCTGGAGGATTTCTTCATCGGCGCCGAAGCTCATCAAACGTGCTCCTTAAAAACCAAGGCTGGACAGCAGATCATCGACATCGTCCTGACCGGACACGACGTCTTCACGCTTATCgg                                      >  7:298903/1‑148 (MQ=255)
                                                                                                                 tCAAGTATTTCGCCGGCCTCTACCAGGAAATCCTGGAGGATTTCTTCATCGGCGCCGAAGCTCATCAAACGTGCTCCTTAAAAACCAAGGCTGGACAGCAGATCATCGACATCGTCCTGACCGGACACGACGTCTTCACGCTTATCgg                                      >  4:71621/1‑148 (MQ=255)
                                                                                                                 tCAAGTATTTCGCCGGCCTCTACCAGGAAATCCTGGAGGATTTCTTCATCGGCGCCGAAGCTCATCAAACGTGATCCTTAAAAACCAAGGCTGGACAGCAGATCATCGACATCGTCCTGACCGGACACGACTTCTTCACGCTTATCgg                                      >  1:157535/1‑148 (MQ=255)
                                                                                                                   aaGTATTTCGCCGGCCTCTACCAGGAAATCCTGGAGGATTTCTTCATCGGCGCCGAAGCTCATCAAACGTGCTCCTTAAAAACCAAgg                                                                                                >  4:215925/1‑88 (MQ=255)
                                                                                                                   aaGTATTTCGCCGGCCTCTACCAGGAAATCCTGGAGGATTTCTTCATCGGCGCCGAAGCTCATCAAACGTGCTCCTTAAAAACCAAgg                                                                                                <  3:215925/88‑1 (MQ=255)
                                                                                                                   aaGTATTTCGCCGGCCTCTACCAGGAAATCCTGGAGGATTTCTTCATCGGCGCCGAAGCTCATCAAACGTGCTCCTTAAAAACCAAGGCTGGACAGCAGATCATCGACATCGTCCTGACCGGACACGACGTCTTCACGCTTATCGGCAt                                   >  5:224375/1‑149 (MQ=255)
                                                                                                                     gTATTTCGCCGGCCTCTACCAGGAAATCCTGGAGGATTTCTTCATCGGCGCCGAAGCTCATCAAACGTGCTCCTTAAAAACCAAGGCTGGACAGCAGATCATCGACATCGTCCTGACCGGACACGACGTCTTCACGCTTATCGGCATg                                  >  4:24410/1‑148 (MQ=255)
                                                                                                                      tATTTCGCCGGCCTCTACCAGGAAATCCTGGAGGATTTCTTCATCGGCGCCGAAGCTCATCAAACGTGCTCCTTAAAAACCAAGGCTGGAcagc                                                                                       >  8:248932/1‑94 (MQ=255)
                                                                                                                      tATTTCGCCGGCCTCTACCAGGAAATCCTGGAGGATTTCTTCATCGGCGCCGAAGCTCATCAAACGTGCTCCTTAAAAACCAAGGCTGGAcagc                                                                                       <  7:244008/94‑1 (MQ=255)
                                                                                                                      tATTTCGCCGGCCTCTACCAGGAAATCCTGGAGGATTTCTTCATCGGCGCCGAAGCTCATCAAACGTGCTCCTTAAAAACCAAGGCTGGAcagc                                                                                       <  7:248932/94‑1 (MQ=255)
                                                                                                                      tATTTCGCCGGCCTCTACCAGGAAATCCTGGAGGATTTCTTCATCGGCGCCGAAGCTCATCAAACGTGCTCCTTAAAAACCAAGGCTGGAcagc                                                                                       >  8:244008/1‑94 (MQ=255)
                                                                                                                        tttCGCCGGCCTCTACCAGGAAATCCTGGAGGATTTCTTCATCGGCGCCGAAGCTCATCAAACGTGCTCCTTAAAAACCAAGGCTGGACAGCAGATCATCGACATCGTCCTGACCGGACACGACGTCTTCACGCTTATCGGCATGAATc                              >  3:17476/1‑149 (MQ=255)
                                                                                                                        tttCGCCGGCCTCTACCAGGAAATCCTGGAGGATTTCTTCATCGGCGCCGAAGCTCATCAAACGTGCTCCTTAAAAACCAAGGCTGGACAGCAGATCATCGACATCGTCCTGACCGGACACGACGTCTTCACGCTTATCGGCATGAATc                              >  3:342493/1‑149 (MQ=255)
                                                                                                                              cGGCATCTACCAGGAAATCCTGGAGGATTTCTTCATCGGCGCCGAAGCTCATCAAACGTGCTCCTTAAAAACCAAGGCTGGACAGCAGATCATCGACATCTTCCTGACCTGACactac                                                       >  2:265740/1‑118 (MQ=255)
                                                                                                                                  ctctACCAGGAAATCCTGGAGGATTTCTTCATCGGCGCCGAAGCTCATCAAACGTGCTCCTTAAAAACCAAGGCTGGACAGCAGATCATCGACATCGTCCTGACCGGACACGACGTCTTCACGCTTATCGGCATGAATCTGCGGACCtt                    >  2:232355/1‑149 (MQ=255)
                                                                                                                                   tctACCAGGAAATCCTGGAGGATTTCTTCATCGGCGCCGAAGCTCATCAAACGTGCTCCTTAAAAACCAAGGCTGGACAGCTGATCATCGACATCGTCCTGACCGGACACGACGTCtt                                                  >  4:371087/1‑118 (MQ=255)
                                                                                                                                   tctACCAGGAAATCCTGGAGGATTTCTTCATCGGCGCCGAAGCTCATCAAACGTGCTCCTTAAAAACCAAGGCTGGACAGCTGATCATCGACATCGTCCTGACCGGACACGACGTCtt                                                  <  3:371087/118‑1 (MQ=255)
                                                                                                                                     tACCAGGAAATCCTGGAGGATTTCTTCATCGGCGCCGAAGCTCATCAAACGTGCTCCTTAAAAACCAAGGCTGGACAGCAGATCATCGACATCGTCCTGACCGGACACGACGTCTTCACGCTTATCGGCATGAATCTGCGGACCTTc                   <  4:17476/147‑1 (MQ=255)
                                                                                                                                     tACCAGGAAATCCTGGAGGATTTCTTCATCGGCGCCGAAGCTCATCAAACGTGCTCCTTAAAAACCAAGGCTGGACAGCAGATCATCGACATCGTCCTGACCGGACACGACGTCTTCACGCTTATCGGCATGAATCTGCGGACCTTCAc                 >  5:281657/1‑149 (MQ=255)
                                                                                                                                     tACCAGGAAATCCTGGAGGATTTCTTCATCGGCGCCGAAGCTCATCAAACGTGCTCCTTAAAAACCAAGGCTGGACAGCAGATCATCAACATCGTCCTGACCGGACACGACGTCTTCACGCTTATCGGCATGTATCTGCGGACCTTc                   <  6:335073/147‑1 (MQ=255)
                                                                                                                                       tgaagAAATCCTCCAGGATTTCTTCATCGGCGCCGAAGCTCATCAAACGTGCTCCTTAAAAACCAAGGCTGGACAGCAGATCATCGACATCGTCCTGACCGGACACGACGTCTTCACGCTTATCGGCATGAATCTGCGGACCTTCAccc               >  2:312818/5‑149 (MQ=255)
                                                                                                                                                 cTGGAGGATTTCTTCATCGGCGCCGAAGCTCATCAAACGTGCTCCTTAAAAACCAAGGCTGGACAGCAGATCATCGACATCGTCCTGACCGGACACGACGTCTTCACGCTTATCGGCATGAATCTGCGGACCTTCACCCCGAGTCGGAt     <  1:163903/149‑1 (MQ=255)
                                                                                                                                                  tGGAGGATTTCTTCATCGGCGCCGAAGCTCATCAAACGTGCTCCTTAAAAACCAAGGCTGGACAGCAGATCATCGACATCGTCCTGACCGGACACGACGTCTTCACGCTTATCGGCATGAATCTGCGGACCTTCACCCCGAGTCGGATg    >  5:110107/1‑149 (MQ=255)
                                                                                                                                                  tGGAGGATTTCTTCATCGGCGCCGAAGCTCATCAAACGTGCTCCTTAAAAACCAAGGCTGGACAGCAGATCATCGACATCGTCCTGACCGGACACGACGTCTTCACGCTTATCGGCATGAATCTGCGGACCTTCACCCCGAGTCGGATg    >  1:118707/1‑149 (MQ=255)
                                                                                                                                                  tGGAGGATTTCTTCATCGGCGCCGAAGCTCATCAAACGTGCTCCTTAAAAACCAAGGCTGGACAGCAGATCATCGACATCGTCCTGACCGGACACGACGTCTTCACGCTTAGCGGCATGAATATGCGGACCTTCACCCCGAGTCgtaag    >  4:357030/1‑145 (MQ=255)
                                                                                                                                                   ggaggaTTTCTTCATCGGCGCCGAAGCTCATCAAACGTGCTCCTTAAAAACCAAGGCTGGACAGCAGATCATCGACATCGTCCTGACCGGACACGACGTCTTCACGCTTATCGGCATGAATCTGCGGACCTTCACCCCGAGTCGGATGt   >  7:304576/1‑149 (MQ=255)
                                                                                                                                                    gaggaTTTCTTCATCGGCGCCGAAGCTCATCAAACGTGCTCCTTAAAAACCAAGGCTGGACAGCAGATCATCGACATCGTCCTGACCGGACACGACGTCTTCACGCTTATCGGCATGAATCTGCGGACCTTCACCCCGAGTCGGATGtt  <  7:23107/149‑1 (MQ=255)
                                                                                                                                                    gaggaTTTCTTCATCGGCGCCGAAGCTCATCAAACGTGCTCCTTAAAAACCAAGGCTGGACAGCAGATCATCGACATCGTCCTGACCGGACACGACGTCTTCACGCTTATCGGCATGAATCTGCGGACCTTCACCCCGAGTCGGATGtt  <  6:224375/149‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                    
GCTGAAGGAAGCCGGCGCCCCCTTTTACAGTGTGGAAACCGCGAAAGATCGCATTGAGCAGGTCGGCATCGTCAGGTCGGCTTTCCAGCTCGACCAGTTGCTCGGACAATTGCTCAAGTATTTCGCCGGCCTCTACCAGGAAATCCTGGAGGATTTCTTCATCGGCGCCGAAGCTCATCAAACGTGCTCCTTAAAAACCAAGGCTGGACAGCAGATCATCGACATCGTCCTGACCGGACACGACGTCTTCACGCTTATCGGCATGAATCTGCGGACCTTCACCCCGAGTCGGATGTT  >  NC_002947/4929994‑4930290

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