Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,378,260 A→T 6.7% intergenic (‑37/‑156) grpE ← / → recN heat shock protein GrpE/DNA repair/recombination protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,378,2600AT 6.7% 208.8 / 2.6 69intergenic (‑37/‑156)grpE/recNheat shock protein GrpE/DNA repair/recombination protein
Reads supporting (aligned to +/- strand):  ref base A (27/37);  new base T (2/3);  total (29/40)
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.89e-01

TGCGGCCAGCTGTTCCTCGAGCTCCAGAACACGGGTATCAGCTGCATTGCCTGCACCGGTTTCTTCCACGTTAAGGTCTTTCTCGTTCAGCTGCTCATCAGCCATGGGGTCTCTCCTGCGCAATTTTGTGGACTGCGAGCCAGGCTCGCCTTGGATGTCGGCTATATGGGGCCGGAAAAACCAGGTTCAAGGGGCAAGGTAGTTTTCAAGACTGCAACATAATCTGCGAAGGGGCATTGGCAGGCCTGAAAAAAGTACTGTATAAATAACCAGACCTGATTTCCGGGA  >  NC_002947/5378119‑5378406
                                                                                                                                              |                                                                                                                                                  
tGCGGCCAGCTGTTCCTCGAGCTCCAGAACACGGGTATCAGCTGCATTGCCTGCACCGGTTTCTTCCACGTTAAGGTCTTTCTCGTTCAGCTGCTCATCAGCCATGGGGTCTCTCCTGCGCAATTTTGTGGACTGCGAGCCAGGCTCGc                                                                                                                                             >  6:41786/1‑149 (MQ=255)
    ccccgCTGTCCCTCCAGCACCATAACAGGGGTATCCGCTGAATTGCCTTCACCGTCCTTTTCCCCCCCCCCAGTCTTTCTCGTTCAGCTGCTCATCAGCCATGGGGTCTCTCCTGCGCAATTTTGTGGACTGCGAGCCTGGCTCGcagt                                                                                                                                          <  1:130848/145‑4 (MQ=255)
     ccAGCTGTTCCTCGAGCTCCAGAACACGGGTATCAGCTGCATTGCCTGCACCGGTTTCTTCCACGTTAAGGTCTTTCTCGTTCAGCTGCTCATCAGCCATGGGGTCTCTCCTGCGCAATTTTGTGGACTGCGAGCCTGGCTCGcagtcc                                                                                                                                        <  6:62277/149‑6 (MQ=255)
         cTGTTCCTCGAGCTCCAGAACACGGGTATCAGCTGCATTGCCTGCACCGGTTTCTTCCACGTTAAGGTCTTTCTCGTTCAGCTGCTCATCAGCCATGGGGTCTCTCCTGCGCAATTTTGTGGACTGCGAGCCAGGCTCGCCTTGGATGt                                                                                                                                    <  3:79682/149‑1 (MQ=255)
          tGTTCCTCGAGCTCCAGAACACGGGTATCAGCTGCATTGCCTGCACCGGTTTCTTCCACGTTAAGGTCTTTCTCGTTCAGCTGCTCATCAGCCATGGGGTCTCTCCTGCGCAATTTTGTGGACTGCGAGCCAGGCTCGCCTTGGATGTc                                                                                                                                   >  4:129072/1‑149 (MQ=255)
            ttCCTCGAGCTCCAGAACACGGGTATCAGCTGCATTGCCTGCACCGGTTTCTTCCACGTTAAGGTCTTTCTCGTTCAGCTGCTCATCAGCCATGGGGTCTCTCCTGCGCAATTTTGTGGACTGCGAGCCAGGCTCGCCTTGGATGTCgg                                                                                                                                 <  8:233624/149‑1 (MQ=255)
            ttCCTCGAGCTCCAGAACACGGGTATCAGCTGCATTGCCTGCACCGGTTTCTTCCACGTTAAGGTCTTTCTCGTTCAGCTGCTCATCAGCCATGGGGTCTCTCCTGCGCAATTTTGTGGACTGCGAGCCAGGCTCGCCTTGGATGTCg                                                                                                                                  <  2:4149/148‑1 (MQ=255)
            ttCCTCGAGCTCCAGAACACGGGTATCAGCTGCATTGCCTGCACCGGTTTCTTCCACGGTAAGGTCTTTCTCGTTCAGCTGCTCATCAGCCATGGGGTCTCTCCTGCGCAATTTTGTGGACTGCGAGCCAGGCTCGCCTTGGATGTc                                                                                                                                   <  7:281997/147‑1 (MQ=255)
              ccTCGAGCTCCAGAACACGGGTATCAGCTGCATTGCCTGCACCGGTTTCTTCCACGTTAAGGTCTTTCTCGTTCAGCTGCTCATCAGCCATGGGGTCTCTCCTGCGCAATTTTGTGGACTGCGAGCCAGGCTCGCCTTGGATGTCGGCt                                                                                                                               >  4:259556/1‑149 (MQ=255)
                 cGAGCTCCAGAACACGGGTATCAGCTGCATTGCCTGCACCGGTTTCTTCCACGTTAAGGTCTTTCTCGTTCAGCTGCTCATCAGCCATGGGGTCTCTCCTGCGCAATTTTGTGGACTGCGAGCCAGGCTCGCCTTGGATGTCGTCtata                                                                                                                            >  8:134338/1‑149 (MQ=255)
                   aGCTCCAGAACACGGGTATCAGCTGCATTGCCTGCACCGGTTTCTTCCACGTTAAGGTCTTTCTCGTTCAGCTGCTCATCAGCCATGGGGTCTCTCCTGCGCAATTTTGTGGACTGCGAGCCAGGCTCGCCTTGGATGTCGGCTATATg                                                                                                                          >  1:69804/1‑149 (MQ=255)
                    gCTCCAGAACACGGGTATCAGCTGCATTGCCTGCACCGGTTTCTTCCACGTTAAGGTCTTTCTCGTTCAGCTGCTCATCAGCCATGGGGTCTCTCCTGCGCAATTTTGTGGACTGCGAGCCAGGCTCGCCTTGGATGTCGGCTATATgg                                                                                                                         <  3:129072/149‑1 (MQ=255)
                    gCTCCAGAACACGGGTATCAGCTGCATTGCCTGCACCGGTTTCTTCCACGTTAAGGTCTTTCTCGTTCAGCTGCTCATCAGCCATGGGGTCTCTCCTGCGCAATTTTGTGGACTGCGAGCCAGGCTCGCCTTGGATGTCGGCTATATgg                                                                                                                         <  6:189768/149‑1 (MQ=255)
                      tCCAGAACACGGGTATCAGCTGCCTTGGCTGCACCGGTTTCTTCCACGTTAACGTCTTTCTCGTTCGGCTGCTCATCAGCCATGGGGTCTCTCCTGCGGAATTTTGTGGACTGCGAGCCAGGCTCGCCTTGGATGTGGGGTATATgggg                                                                                                                       >  6:217482/1‑149 (MQ=255)
                      tCCAGAACACGGGTATCAGCTGCATTGCCTGCACCGGTTTCTTCCACGTTAAGGTCTTTCTCGTTCAGCTGCTCATCAGCCATGGGGTCTCTCCTGCGCAATTTTGTGGACTGCGAGCCAGGCTCGCCTTGGATGTCGGCTATATgggg                                                                                                                       >  1:12731/1‑149 (MQ=255)
                        cAGAACACGGGTATCAGCTGCATTGCCTGCACCGGTTTCTTCCACGTTAAGGTCTTTCTCGTTCAGCTGCTCATCAGCCATGGGGTCTCTCCTGCGCAATTTTGTGGACTGCGAGCCAGGCTCGCCTTGGATGTCGGCTATATGGGGcc                                                                                                                     >  5:316373/1‑149 (MQ=255)
                         aGAACACGGGTATCAGCTGCATTGCCTGCACCGGTTTCTTCCACGTTAAGGTCTTTCTCGTTCAGCTGCTCATCAGCCATGGGGTCTCTCCTGCGCAATTTTGTGGACTGCGAGCCAGGCTCGCCTTGGATGTCGGCTATATGGGGCCg                                                                                                                    <  4:365841/149‑1 (MQ=255)
                              acGGGTATCAGCTGCATTGCCTGCTCCGGTTTCTTCCACGTTAAGGTCTTTCTCGTTCAGCTGCTCATCAGCCATGGGGTCTCTCCTGCGCAATTTTGTGGACTGCGAGCCAGGCTCGCCTTGGATGTCGGCTATATGGGGCCGGaaaa                                                                                                               <  7:134338/149‑1 (MQ=255)
                                gggTATCAGCTGCATTGCCTGCACCGGTTTCTTCCACGTTAAGGTCTTTCTCGTTCAGCTGCTCATCAGCCATGGGGTCTCTCCTGCGCAATTTTGTGGACTGCGAGCCAGGCTCGCCTTGGATGTCGGCTATATGGGGCCGGAAAAAc                                                                                                             <  1:7592/149‑1 (MQ=255)
                                 ggTATCAGCTGCATTGCCTGCACCGGTTTCTTCCACGTTAAGGTCTTTCTCGTTCAGCTGCTCATCAGCCATGGGGTCTCTCCTGCGCAATTTTGTGGACTGCGAGCCAGGCTCGCCTTGGATGTCGGCTATATGGGGCCGGAAAAAc                                                                                                             >  6:399789/1‑148 (MQ=255)
                                   tATCAGCTGCATTGCCTGCACCGGTTTCTTCCACGTTAAGGTCTTTCTCGTTCAGCTGCTCATCAGCCATGGGGTCTCTCCTGCGCAATTTTGTGGACTGCGAGCCAGGCTCGCCTTGGATGTCGGCTATATGGGGCCGGAAAAACCAg                                                                                                          <  6:316373/149‑1 (MQ=255)
                                   tATCAGCTGCATTGCCTGCACCGGTTTCTTCCACGTTAAGGTCTTTCTCGTTCAGCTGCTCATCAGCCATGGGGTCTCTCCTGCGCAATTTTGTGGACTGCGAGCCAGGCTCGCCTTGGATGTCGGCTATATGGGGCCGGAAAAACCAg                                                                                                          <  4:1065/149‑1 (MQ=255)
                                     tCAGCTGCATTGCCTGCACCGGTTTCTTCCACGTTAAGGTCTTTCTCGTTCAGCTGCTCATCAGCCATGGGGTCTCTCCTGCGCAATTTTGTGGACTGCGAGCCAGGCTCGCCTTGGATGTCGGCTATATGGGGCCGGAAAAACCAg                                                                                                          <  4:375273/147‑1 (MQ=255)
                                     tCAGCTGCATTGCCTGCACCGGTTTCTTCCACGTTAAGGTCTTTCTCGTTCAGCTGCTCATCAGCCATGGGGTCTCTCCTGCGCAATTTTGTGGACTGCGAGCCAGGCTCGCCTTGGATGTCGGCTATATGGGGCCGGAAAAACCAg                                                                                                          >  3:375273/1‑147 (MQ=255)
                                       aacTGCATTGCCTGCACCGGTTTCTTCCACGTTAAGGTCTTTCTCGTTCAGCTGCTCATCAGCCATGGGGTCTCTCCTGCTCAATTTTGTGGACTGCGAGCCAGGCTCGCCTTGGa                                                                                                                                       <  8:303089/114‑1 (MQ=255)
                                       aGCTGCATTGCCTGCACCGGTTTCTTCCACGTTAAGGTCTTTCTCGTTCAGCTGCTCATCAGCCATGGGGTCTCTCCTGCGCAATTTTGTGGACTGCGAGCCAGGCTCGCCTTGGATGTCGGCTATATGGGGCCGGAAAAACCAGGtt                                                                                                       >  6:408749/1‑148 (MQ=255)
                                        gCTGCATTGCCTGCACCGGTTTCTTCCACGTTAAGGTCTTTCTCGTTCAGCTGCTCATCAGCCATGGGGTCTCTCCTGCGCAATTTTGTGGACTGCGAGCCAGGCTCGCCTTGGATGTCGGCTATATGGGGCCGGAAAAACCAGGTTc                                                                                                      >  6:375136/1‑148 (MQ=255)
                                                 ccTGCACCGGTTTCTTCCACGTTAAGGTCTTTCTCGTTCAGCTGCTCATCAGCCATGGGGTCTCTCCTGCGCAATTTTGTGGACTGCGAGCCAGGCTCGCCTTGGATGTCGGCTATATGGGGCCGGAAAAACCAGGTTCAAGGGGCAAg                                                                                            <  1:92217/149‑1 (MQ=255)
                                                     cACCGGTTTCTTCCACGTTAAGGTCTTTCTCGTTCAGCTGCTCATCAGCCATGGGGTCTCTCCTGCGCAATTTTGTGGACTGCGAGCCAGGCTCGCCTTGGATGTCGGCTATATGGGGCCGGAAAAACCAGGTTCAAGGGGCAAGGt                                                                                          <  5:399789/147‑1 (MQ=255)
                                                     cACCGGTTTCTTCCACGTTAAGGTCTTTCTCGTTCAGCTGCTCATCAGCCATGGGGTCTCTCCTGCGCAATTTTGTGGACTGCGAGCCAGGCTCGCCTTGGATGTCGGCTATATGGGGCCGGAAAAACCAGGTTCAAGGGGCAAGGt                                                                                          <  1:268952/147‑1 (MQ=255)
                                                        cGGTTTCTTCCACGTTAAGGTCTTTCTCGTTCAGCTGCTCATCAGCCATGGGGTCTCTCCTGCGCAATTTTGTGGACTGCGAGCCAGGCTCGCCTTGGATGTCGGCTATATGGGGCCGGAAAAACCAGGTTCAAGGGGCAAGGTAGttt                                                                                     >  4:222705/1‑149 (MQ=255)
                                                         ggTTTCTTCCACGTTAAGGTCTTTCTCGTTCAGCTGCTCATCAGCCATGGGGTCTCTCCTGCGCAATTTTGTGGACTGCGAGCCAGGCTCGCCTTGGATGTCGGCTATATGGGGCCGGAAAAACCAGGTTCAAGGGGCAAGGTAGtttt                                                                                    <  4:180281/149‑1 (MQ=255)
                                                              tcttcCACGTTAAGGTCTTTCTCGTTCAGCTGCTCATCAGCCATGGGGTCTCTCCTGCGCAATTTTGTGGACTGCGAGCCAGGCTCGCCTTGGATGTCGGCTATATGGGGCCGGAAAAACCAGGTTCAAGGGGCAAGGTAGTTTTCaa                                                                                 <  5:342052/148‑1 (MQ=255)
                                                              tcttcCACGTTAAGGTCTTTCTCGTTCAGCTGCTCATCAGCCATGGGGTCTCTCCTGCGCAATTTTGTGGACTGCGAGCCAGGCTCGCCTTGGATGTCGGCTATATGGGGCCGGAAAAACCAGGTTCAAGGGGCAAGGTAGTTTTCaa                                                                                 <  5:41786/148‑1 (MQ=255)
                                                                ttcCACGTTAAGGTCTTTCTCGTTCAGCTGCTCATCAGCCATGGGGTCTCTCCTGCGCAATTTTGTGGACTGCGAGCCAGGCTCGCCTTGGATGTCGGCTATATGGGGCCGGAAAAACCAGGTTCAAGGGGCAAg                                                                                            >  6:382157/1‑135 (MQ=255)
                                                                ttcCACGTTAAGGTCTTTCTCGTTCAGCTGCTCATCAGCCATGGGGTCTCTCCTGCGCAATTTTGTGGACTGCGAGCCAGGCTCGCCTTGGATGTCGGCTATATGGGGCCGGAAAAACCAGGTTCAAGGGGCAAg                                                                                            <  5:382157/135‑1 (MQ=255)
                                                                ttcCACGTTAAGGTCTTTCTCGTTCAGCTGCTCATCAGCCATGGGGTCTCTCCTGCGCAATTTTGTGGACTGCGAGCCAGGCTCGCCTTGGATGTCGGCTATATGGGGCCGGAAAAACCAGGTTCAAGGGGCAAGGTa                                                                                         <  5:5235/138‑1 (MQ=255)
                                                                ttcCACGTTAAGGTCTTTCTCGTTCAGCTGCTCATCAGCCATGGGGTCTCTCCTGCGCAATTTTGTGGACTGCGAGCCAGGCTCGCCTTGGATGTCGGCTATATGGGGCCGGAAAAACCAGGTTCAAGGGGCAAGGTa                                                                                         >  6:5235/1‑138 (MQ=255)
                                                                  ccACGTTAAGGTCTTTCTCGTTCAGCTGCTCATCAGCCATGGGGTCTCTCCTGCGCAATTTTGTGGACTGCGAGCCAGGCTCGCCTTGGATGTCGGCTATATGGGGCCGGAAAAACCAGGTTCAAGGGGCAAGGTAGTTTTCAAGa                                                                               <  3:259556/146‑1 (MQ=255)
                                                                   cACGTTAAGGTCTTTCTCGTTCAGCTGCTCATCAGCCATGGGGTCTCTCCTGCGCAATTTTGTGGACTGCGAGCCAGGCTCGCCTTGGATGTCGGCTATATGGGGCCGGAAAAACCAGGTTCAAGGGGCAAGGTAGTTTTCAAGACTGc                                                                           <  3:19504/149‑1 (MQ=255)
                                                                     cGTTAAGGTCTTTCTCGTTCAGCTGCTCATCAGCCATGGGGGCTCTCCTGCGCAATTTTGTGGACTGCGAGCCAGGCTCGCCTTGGATGTCGGCTATATGGGGCCGGAAAAACCAGGTTCAAGGGGCAAGGTAGTTTTCAAGACTGCaa                                                                         <  6:369504/149‑1 (MQ=255)
                                                                      gTTAAGGTCTTTCTCGTTCAGCTGCTCATCAGCCATGGGGTCTCTCCTGCGCAATTTTGTGGACTGCGAGCCAGGCTCGCCTTGGATGTCGGCTATATGGGGCCGGAAAAACCAGGTTCAAGGGGCAAGGTAGTTTTCAAGACTGCAAc                                                                        <  5:408749/149‑1 (MQ=255)
                                                                           ggTCTTTCTCGTTCAGCTGCTCATCAGCCATGGGGTCTCTCCTGCGCAATTTTGTGGACTGCGAGCCAGGCTCGCCTTGGATGTCGGCTATATGGGGCCGGAAAAACCAGGTTCAAGGGGCAAGGTAGTTTTCAAGACTGCAACATAAt                                                                   <  3:261021/149‑1 (MQ=255)
                                                                            gTCTTTCTCGTTCAGCTGCTCATCAGCCATGGGGTCTCTCCTGCGCAATTTTGTGGACTGCGAGCCAGGCTCGCCTTGGATGTCGGCTATATGGGGCCGGAAAAACCAGGTTCAAGGGGCAAGGTAGTTTTCAAGACTGCAACATAATc                                                                  >  5:248891/1‑149 (MQ=255)
                                                                            gTCTTTCTCGTTCAGCTGCTCATCAGCCATGGGGTCTCTCCTGCGCAATTTTGTGGACTGCGAGCCAGGCTCGCCTTGGATGTCGGCTATATGGGGCCGGAAAAACCAGGTTCAAGGGGCAAGGTAGTTTTCAAGACTGCAACATAATc                                                                  <  3:222705/149‑1 (MQ=255)
                                                                                    cGTTCAGCTGCTCATCAGCCATGGGGTCTCTCCTGCGCAATTTTGTGGACTGCGAGCCAGGCTCGCCTTGGATGTCGGCTATATGGGGCCGGAAAAACCAGGTTCAAGGGGCAAGGTAGTTTTCAAGACTGCAACATAATCTGCGAAgg                                                          >  6:353498/1‑149 (MQ=255)
                                                                                     gTTCAGCTGCTCATCAGCCATGGGGTCTCTCCTGCGCAATTTTGTGGACTGCGAGCCAGGCTCGCCTTGGATGTCGGCTATATGGGGCCGGAAAAACCAGGTTCAAGGGGCAAGGTAGTTTTCAAGAc                                                                              <  8:380591/128‑1 (MQ=255)
                                                                                     gTTCAGCTGCTCATCAGCCATGGGGTCTCTCCTGCGCAATTTTGTGGACTGCGAGCCAGGCTCGCCTTGGATGTCGGCTATATGGGGCCGGAAAAACCAGGTTCAAGGGGCAAGGTAGTTTTCAAGAc                                                                              >  7:380591/1‑128 (MQ=255)
                                                                                          gctgctCATCAGCCATGGGGTCTCTCCTGCGCAATTTTGTGGACTGCGAGCCAGGCTCGCCTTGGATGTCGGCTATATGGGGCCGGAAAAACCAGGTTCAAGGGGCAAGGTAGTTTTCAAGACTGCAACATAATCTGCGAAGGGGCAtt                                                    <  4:282127/149‑1 (MQ=255)
                                                                                            tgctCATCAGCCATGGGGTCTCTCCTGCGCAATTTTGTGGACTGCGAGCCAGGCTCGCCTTGGATGTCGGCTATATGGGGCCGGAAAAACCAGGTTCAAGGGGCAAGGTAGTTTTCAAGACTGCAACATAATCTGCGAAGGGGCATTgg                                                  >  2:64806/1‑149 (MQ=255)
                                                                                                      ccATGGGGTCTCTCCTGCGCAATTTTGTGGACTGCGAGCCAGGCTCGCCTTGGATGTCGGCTATATGGGGCCGGAAAAACCAGGTTCAAGGGGCAAGGTAGTTTTCAAGACTGCAACATAATCTGCGAAGGGGCATTGGCAGGCCTGa                                         >  7:231755/1‑148 (MQ=255)
                                                                                                      ccATGGGGTCTCTCCTGCGCAATTTTGTGGACTGCGAGCCAGGCTCGCCTTGGATGTCGGCTATATGGGGCCGGAAAAACCAGGTTCAAGGGGCAAGGTAGTTTTCAAGACTGCAACATAATCTGCGAAGGGGCATTGGCAGGCCTGa                                         >  5:356321/1‑148 (MQ=255)
                                                                                                        aTGGGGTCTCTCCTGCGCAATTTTGTGGACTGCGAGCCAGGCTCGCCTTGGATGTCGGCTATATGGGGCCGGAAAAACCAGGTTCAAGGGGCAAGGTAGTTTTCAAGACTGCAACATAATCTGCGAAGGGGCATTGGCAGGCCTGaaaa                                      <  2:12731/149‑1 (MQ=255)
                                                                                                        aTGGGGTCTCTCCTGCGCAATTTTGTGGACTGCGAGCCAGGCTCGCCTTGGATGTCGGCTATATGGGGCCGGAAAAACCAGGTTCAAGGGGCAAGGTAGTTTTCAAGACTGCAACATAATCTGCGAAGGGGCATTGGCAGGCCTGaaa                                       >  2:224934/1‑148 (MQ=255)
                                                                                                           gggTCTCTCCTGCGCAATTTTGTGGACTGCGAGCCAGGCTCGCCTTGGATGTCGGCTATATGGGGCCGGAAAAACCAGGTTCAAGGGGCAAGGTAGTTTTCAAGACTGCAACATAATCTGCGAAGGGGCATTGGCAGGCCTGAAAAAAg                                   <  5:375136/149‑1 (MQ=255)
                                                                                                             gTCTCTCCTGCGCAATTTTGTGGACTGCGAGCCAGGCTCGCCTTGGATGTCGGCTATATGGGGCCGGAAAAACCAGGTTCAAGGGGCAAGGTAGTTTTCAAGACTGCAACATAATCTGCGAAGGGGCATTGGCAGGCCTGAAAAAAGt                                  >  5:111616/1‑148 (MQ=255)
                                                                                                                      gcgcAATTTTGTGGACTGCGAGCCAGGCTCGCCTTGGATGTCGGCTATATGGGGCCGGAAAAACCAGGTTCAAGGGGCAAGGTAGTTTTCAAGACTGCAACATAATCTGCGAAGGGGCATTGGCAGGCCTGAAAAAAGTACTGTataaa                        >  4:271428/1‑149 (MQ=255)
                                                                                                                            ttttGTGGACTGCGAGCCAGGCTCGCCTTGGATGTCGGCTATATGGGGCCGGAAAAACCAGGTTCAAGGGGCAAGGTAGTTTTCAAGACTGCAACATAATCTGCGAAGGGGCATTGGCAGGCCTGAAAAAAGTACTGTATAAATAACCa                  <  6:248891/149‑1 (MQ=255)
                                                                                                                             tttGTGGACTGCGAGCCAGGCTCGCCTTGGATGTCGGCTATATGGGGCCGGAAAAACCAGGTTCAAGGGGCAAGGTAGTTTTCAAGACTGCAACATAATCTGCGAAGGGGCATTGGCAGGCCTGAAAAAAGTACTGTATAAATAAc                    >  2:241818/1‑146 (MQ=255)
                                                                                                                                   gACTGCGAGCCAGGCTCGCCTTGGATGTCCGCTTTATGGGTCAGGAAAAACCAGTTGCAAGGGGCAAGGTAGTTTTCTAGACTGCAACATAATCTTCGAAGGGGCTTTGGCTGGCCTGAAAAAAGTACTGTataa                         <  6:361777/135‑1 (MQ=255)
                                                                                                                                   gACAGCGAGCCAGGCTCGCCTTGGATGTCTGCTAAATGGGGCCGGAAAAACCAGGTTAAAGGGGCAAGGTAGTTTTCAAGACTGCAACATAATCTGCGAAGGGGCATTGGCAGGCCTGAATAAAGTACTGTataa                         >  5:361777/1‑135 (MQ=255)
                                                                                                                                     cTGCGAGCCAGGCTCGCCTTGGATGTCGGCTATATGGGGCCGGAAAAACCAGGTTCAAGGGGCAAGGTAGTTTTCAAGACTGCAACATAATCTGCGAAGGGGCATTGGCAGGCCTGAAAAAAGTACTGTATAAATAACCAGACCTGAtt         <  1:64806/149‑1 (MQ=255)
                                                                                                                                      tGCGAGCCAGGCTCGCCTTGGATGTCGGCTATATGGGGCCGGAAAAACCAGGTTCAAGGGGCAAGGTAGTTTTCAAGACTGCAACATAATCTGCGAAGGGGCATTGGCAGGCCTGAAAAAAGTACTGTATAAATAACCAGACCTGAttt        >  2:220270/1‑149 (MQ=255)
                                                                                                                                         gAGCCTGGCTCGCCTTGGATGTCGGCTATATGGGGCCGGAAAAACCAGGTTCAAGGGGCAAGGTAGTTTTCAAGACTGCTACATAATCTGCGAAGGGGCATTGGCAGGCCTGAAAAAAGTACTGTATAAATTACCAGACCTGATTTCCg     >  2:151411/1‑149 (MQ=255)
                                                                                                                                         gAGCCTGGCTCGCCTTGGATGTCGGCTATATGGGGCCGGAAAAACCAGGTTCAAGGGGCAAGGTAGTTTTCAAGACTGCAACATAATCTGCGAAGGGGCATTGGCAGGCCt                                           >  1:293299/1‑111 (MQ=255)
                                                                                                                                         gAGCCTGGCTCGCCTTGGATGTCGGCTATATGGGGCCGGAAAAACCAGGTTCAAGGGGCAAGGTAGTTTTCAAGACTGCAACATAATCTGCGAAGGGGCATTGGCAGGCCt                                           <  2:293299/111‑1 (MQ=255)
                                                                                                                                          aGCCAGGCTCGCCTTGGATGTCGTCTATATGGGTCCGGAAAAACCAGGTTCAAGGGGCAAGGTAGTTTTCAAGACTGCAACATAATCTGCGAAGGGGCATTGGCAGGCCTGAAAAAAGTACTGTATAAATAACCAGACCTGATTTCCgg    <  5:217482/149‑1 (MQ=255)
                                                                                                                                            ccAGGCTCGCCTTGGATGTCGGCTATATGGGGCCGGAAAAACCAGGTTCAAGGGGCAAGGTAGTTTTCAAGACTGCAACATAATCTGCGAAGGGGCATTGGCAGGCCTGAAAAAAGTACTGTATAAATAACCAGACCTGATTTCCGGGa  <  6:356321/149‑1 (MQ=255)
                                                                                                                                             cAGGCTCGCCTTGGATGTCGGCTATATGGGGCCGGAAAAACCAGGTTCAAGGGGCAAGGTAGTTTTCAAGACTGCAACATAATCTGCGAAGGGGCATTGGCAGGCCTGAAAAAAGTACTGTATAAATAACCAGACCTGATTTCCGGGa  <  6:111616/148‑1 (MQ=255)
                                                                                                                                              |                                                                                                                                                  
TGCGGCCAGCTGTTCCTCGAGCTCCAGAACACGGGTATCAGCTGCATTGCCTGCACCGGTTTCTTCCACGTTAAGGTCTTTCTCGTTCAGCTGCTCATCAGCCATGGGGTCTCTCCTGCGCAATTTTGTGGACTGCGAGCCAGGCTCGCCTTGGATGTCGGCTATATGGGGCCGGAAAAACCAGGTTCAAGGGGCAAGGTAGTTTTCAAGACTGCAACATAATCTGCGAAGGGGCATTGGCAGGCCTGAAAAAAGTACTGTATAAATAACCAGACCTGATTTCCGGGA  >  NC_002947/5378119‑5378406

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: