Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,983,221 T→A 5.8% M1289K (ATG→AAG)  PP_5242 → GAF domain‑containing sensor histidine kinase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,983,2210TA 5.8% 208.5 / 2.2 69M1289K (ATG→AAG) PP_5242GAF domain‑containing sensor histidine kinase
Reads supporting (aligned to +/- strand):  ref base T (32/33);  new base A (2/2);  total (34/35)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.55e-01

AGACGTTGCATCCGTTTTTGCGTACCCAGCCGATTCAGGAGACCTATCGCGCCACCCACCAGGCCAGGCTGGACCTCGAAGCCGGCATCGAAGCCGCGCGGGCGCTGTCCGAAGAAGTGCTGCTCGAAGGCCTGATCGAAACCCTGATGGGCCACCTGACCCAGCACTCCGGTGCCGACCATGGCGCGCTGCTGATTGTAAGCGGCGCCGAATTCCAGATGGCCGCCCTGGCCTACATCGACGATGCCGGGCTGCACGTGAGCATGGACAACTGCCAGAAGTTCATGACCC  >  NC_002947/5983074‑5983364
                                                                                                                                                   |                                                                                                                                               
agaCGTTGCATCCGTTTTTGCGTACCCAGCCGATTCAGGAGACCTATCGCGCCACCCACCAGGCCAGGCTGGACCTCGAAGCCGGCATCGAAGCCGCGCGGGCGCTGTCCGAAGAAGTGCTGCTCGAAGGCCTGATCGAAACCCTGATg                                                                                                                                                >  1:261105/1‑149 (MQ=255)
         aTCCGTTTTTGCGTACCCAGCCGATTCAGGAGACCTATCGCGCCACCCACCAGGCCAGGCTGGACCTCGAAGCCGGCATCGAAGCCGCGCGGGCGCTGTCCGAAGAAGTGCTGCTCGAAGGCCTGATCGAAACCCTGATGGGCCACCt                                                                                                                                        >  6:370067/1‑148 (MQ=255)
          tCCGTTTTTGCGTACCCAGCCGATTCAGGAGACCTATCGCGCCACCCACCAGGCCAGGCTGGACCTCGAAGCCGGCATCGAAGCCGCGCGGGCGCTGTCCGAAGAAGTGCTGCTCGAAGGCCTGATCGAACCCCTGATGGGCCACCTGa                                                                                                                                      >  4:252118/1‑149 (MQ=255)
              tttttGCGTACCCAGCCGATTCAGGAGACCTATCGCGCCACCCACCAGGCCAGGCTGGACCTCGAAGCCGGCATCGAAGCCGCGCGGGCGCTGTCCGAAGAAGTGCTGCTCGAAGGCCTGATCGAAACCCTGATGGGCCACCTGACCCa                                                                                                                                  <  1:281303/149‑1 (MQ=255)
              tttttCCGTCCCCCGCGTATTCCGGTGACGTTTCGGGCCACCCACAAGGCATGGCTGTACCTCCAAGCCGGCATGGACGCCGCGGGGGGGCCGTCCGAAGAAGTGCTGCTCGAAGGCCTGATCGAAACCCTGATGGGCCACCTGACCCa                                                                                                                                  <  1:293958/149‑1 (MQ=255)
                        cccAGCCGATTCAGGAGACCTATCGCGCCACCCACCAGGCCAGGCTGGACCTCGAAGCCGGCATCGAAGCCGCGCGGGCGCTGTCCGAAGAAGTGCTGCTCGAAGGCCTGATCGAAACCCTGATGGGCCACCTGACCCAGCACTCCGGt                                                                                                                        <  8:149353/149‑1 (MQ=255)
                         ccAGCCGATTCAGGAGACCTATCGCGCCACCCACCAGGCCAGGCTGGACCTCGAAGCCGGCATCGAAGCCGCGCGGGCGCTGTCCGAAGAAGTGCTGCTCGAAGGCCTGATCGAAACCCTGATGGGCCACCTGACCCAGCACTCCGGTg                                                                                                                       >  5:308153/1‑149 (MQ=255)
                         ccAGCCGATTCAGGAGACCTATCGCGCCACCCACCAGGCCAGGCTGGACCTCGAAGCCGGCATCGAAGCCGCGCGGGCGCTGTCCGAAGAAGTGCTGCTCGAAGGCCTGATCGAAACCCTGATGGGCCACCTGACCCAGCACTCCGGTg                                                                                                                       >  3:87793/1‑149 (MQ=255)
                           aGCCGATTCAGGAGACCTATCGCGCCACCCACCAGGCCAGGCTGGACCTCGAAGCCGGCATCGAAGCCGCGCGGGCGCTGTCCGAAGAAGTGCTGCTCGAAGGCCTGATCGAAACCCTGATGGGCCACCTGACCCAGCACTCCGGTGcc                                                                                                                     >  7:301984/1‑149 (MQ=255)
                              cGATTCAGGAGACCTATCGCGCCACCCACCAGGCCAGGCTGGACCTCGAAGCCGGCATCGAAGCCGCGCGGGCGCTGTCCGAAGAAGTGCTGCTCGAAGGCCTGATCGAAACCCTGATGGGCCACCTGACCCAGCACTCCGGTGCCGAc                                                                                                                  <  7:320787/149‑1 (MQ=255)
                               gATTCAGGAGACCTATCGCGCCACCCACCAGGCCAGGCTGGACCTCGAAGCCGGCATCGAAGCCGCGCGGGCGCTGTCCGAAGAAGTGCTGCTCGAAGGCCTGATCGAAACCCTGATGGGCCACCTGACCCAGCACTCCGGTGCCGAcc                                                                                                                 <  3:274482/149‑1 (MQ=255)
                                     ggAGACCTATCGCGCCACCCACCAGGCCAGGCTGGACCTCGAAGCCGGCATCGAAGCCGCGCGGGCGCTGTCCGAAGAAGTGCTGCTCGAAGGCCTGATCGAAACCCTGATGGGCCACCTGAccc                                                                                                                                   <  6:184250/125‑1 (MQ=255)
                                     ggAGACCTATCGCGCCACCCACCAGGCCAGGCTGGACCTCGAAGCCGGCATCGAAGCCGCGCGGGCGCTGTCCGAAGAAGTGCTGCTCGAAGGCCTGATCGAAACCCTGATGGGCCACCTGAccc                                                                                                                                   >  5:184250/1‑125 (MQ=255)
                                     ggAGACCTATCGCGCCACCCACCAGGCCAGGCTGGACCTCGAAGCCGGCATCGAAGCCGCGCGGGCGCTGTCCGAAGAAGTGCTGCTCGAAGGCCTGATCGAAACCCTGATGGGCCACCTGACCCAGCACTCCGGTGCCGACCATGgcg                                                                                                           >  3:134849/1‑149 (MQ=255)
                                            tATCGCGCCACCCACCAGGCCAGGCTGGACCTCGAAGCCGGCATCGAAGCCGCGCGGGCGCTGTCCGAAGAAGTGCTGCTCGAAGGCCTGATCGAAACCCTGATGGGCCACCTGACCCAGCACTCCGGTGCCGACCATGGCGCgctgct                                                                                                    >  5:354044/1‑149 (MQ=255)
                                            tATCGCGCCACCCACCAGGCCAGGCTGGACCTCGAAGCCGGCATCGAAGCCGCGCGGGCGCTGTCCGAAGAAGTGCTGCTCGAAGGCCTGATCGAAACCCTGATGGGCCACCTGACCCAGCACTCCGGTGCCGACCATGGCGCgctgct                                                                                                    >  2:240537/1‑149 (MQ=255)
                                              tCGCGCCACCCACCAGGCCAGGCTGGACCTCGAAGCCGGCATCGAAGCCGCGCGGGCGCTGTCCGAAGAAGTGCTGCTCGAAGGCCTGATCGAAACCCTGATGGGCCACCTGACCCAGCACTCCGGTGCCGACCATGGCGCGCTGCTGa                                                                                                  <  6:85943/149‑1 (MQ=255)
                                                  gcCACCCACCAGGCCAGGCTGGACCTCGAAGCCGGCATCGAAGCCGCGCGGGCGCTGTCCGAAGAAGTGCTGCTCGAAGGCCTGATCGAAACCCTGATGGGCCACCTGACCCAGCACTCCGGTGCCGACCATGGCGCGCTGCTGATTGt                                                                                              <  6:308153/149‑1 (MQ=255)
                                                      cccaccAGGCCAGGCTGGACCTCGAAGCCGGCATCGAAGCCGCGCGGGCGCTGTCCGAAGAAGTGCTGCTCGAAGGCCTGATCGAAACCCTGATGGGCCACCTGACCCAGCACTCCGGTGCCGACCATGGCGCGCTGCTGATTGTAAgc                                                                                          >  7:84925/1‑149 (MQ=255)
                                                      cccaccAGGCCAGGCTGGACCTCGAAGCCGGCATCCAAGCCGCGCGGGCGCTGTCCGAAGAAGTGCTGCTCCAAGGCCTGATCGAAACCCTGCTGGGCCACCTGACCCAGCACTCCGGTGCCGACCATGGCGCGCGGCTGATTCTAAgc                                                                                          >  7:79661/1‑149 (MQ=255)
                                                       ccaccaGGCCAGGCTGGACCTCGAAGCCGGCATCGAAGCCGCGCGGGCGCTGTCCGAAGAAGTGCTGCTCGAAGGCCTGATCGAAACCCTGATGGGCCACCTGACCCAGCACTCCGGTGCCGACCATGGCGCGCTGCTGATTGTAAgcg                                                                                         <  8:301984/149‑1 (MQ=255)
                                                         accaGGCCAGGCTGGACCTCGAAGCCGGCATCGAAGCCGCGCGGGCGCTGTCCGAAGAAGTGCTGCTCGAAGGCCTGATCGAAACCCTGATGGGCCACCTGACCCAGCACTCCGGTGCCGACCGTGGCGCGCTGCTGATTGTAAgcggc                                                                                       >  8:36885/1‑149 (MQ=255)
                                                            aggccaggcTGGACCTCGAAGCCGGCATCGAAGCCGCGCGGGCGCTGTCCGAAGAAGTGCTGCTCGAAGGCCTGATCGAAACCCTGATGGGCCACCTGACCCAGCACTCCGGTGCCGACCATGGCGCGCTGCTGATTGTAAGCGGCGcc                                                                                    <  2:261105/149‑1 (MQ=255)
                                                                caggcTGGACCTCGAAGCCGGCATCGAAGCCGCGCGGGCGCTGTCCGAAGAAGTGCTGCTCGAAGGCCTGATCGAAACCCTGATGGGCCACCTGACCCAGCACTCCGGTGCCGACCATGGCGCGCTGCTGATTGTAAGCGGCGCCGAAt                                                                                >  8:197598/1‑149 (MQ=255)
                                                                 aggcTGGACCTCGAAGCCGGCATCGAAGCCGCGCGGGCGCTGTCCGAAGAAGTGCTGCTCGAAGGCCTGATCGAAACCCTGATGGGCCACCTGACCCAGCACTCCGGTGCCGACCATGGCGCGCTGCTGATTGTAAGCGGCGCCGAAtt                                                                               <  7:197598/149‑1 (MQ=255)
                                                                    cTGGACCTCGAAGCCGGCATCGAAGCCGCGCGGGCGCTGTCCGAAGAAGTGCTGCTCGAAGGCCTGATCGAAACCCTGATGGGCCACCTGACCCAGCACTCCGGTGCCGACCATGGCGCGCTGCTGATTGTAAGCGGCGCCGAATTCCa                                                                            <  4:87793/149‑1 (MQ=255)
                                                                     tGGACCTCGAAGCCGGCATCGAAGCCGCGCGGGCGCTGTCCGAAGAAGTGCTGCTCGAAGGCCTGATCGAAACCCTGATGGGc                                                                                                                                             >  8:121200/1‑83 (MQ=255)
                                                                     tGGACCTCGAAGCCGGCATCGAAGCCGCGCGGGCGCTGTCCGAAGAAGTGCTGCTCGAAGGCCTGATCGAAACCCTGATGGGc                                                                                                                                             <  7:121200/83‑1 (MQ=255)
                                                                      ggACCTCGAAGCCGGCATCGAAGCCGCGCGGGCGCTGTCCGAAGAAGTGCTGCTCGAAGGCCTGATCGAAACCCTGATGGGCCACCTGACCCAGCACTCCGGTGCCGACCATGGCGCGCTGCTGATTGTAAGCGGCGcc                                                                                    >  4:283427/1‑139 (MQ=255)
                                                                      ggACCTCGAAGCCGGCATCGAAGCCGCGCGGGCGCTGTCCGAAGAAGTGCTGCTCGAAGGCCTGATCGAAACCCTGATGGGCCACCTGACCCAGCACTCCGGTGCCGACCATGGCGCGCTGCTGATTGTAAGCGGCGcc                                                                                    <  3:283427/139‑1 (MQ=255)
                                                                       gACCTCGAAGCCGGCATCGAAGCCGCGCGGGCGCTGTCCGAAGAAGTGCTGCTCGAAGGCCTGATCGAAACCCTGATGGGCCACCTGACCCAGCACTCCGGTGCCGACCATGGCGCGCTGCTGATTGTAAGCGGCGCCGAATTCCAGa                                                                          <  2:5893/148‑1 (MQ=255)
                                                                           tCGAAGCCGGCATCGAAGCCGCGCGGGCGCTGTCCGAAGAAGTGCTGCTCGAAGGCCTGATCGAAACCCTGATGGGCCACCTGACCCAGCACTCCGGTGCCGACCATGGCGCGCTGCTGATTGTAAGCGGCGCCGAATTCCAGATGgcc                                                                     >  4:108499/1‑149 (MQ=255)
                                                                            cGAAGCCGGCATCGAAGCCGCGCGGGCGCTGTCCGAAGAAGTGCTGCTCGAAGGCCTGATCGAAACCCTGATGGGCCACCTGACCCAGCACTCAGGTGCCGACCATGGCGCGCTGCTGATTGTAAGCGGCGCCGAATTCCAGATGgccg                                                                    >  4:196676/1‑149 (MQ=255)
                                                                              aaGCCGGCATCGAAGCCGCGCGGGCGCTGTCCGAAGAAGTGCTGCTCGAAGGCCTGATCGAAACCCTGATGGGCCACCTGACCCAGCACTCCGGTGCCGACCGTGGCGCGCTGCTGATTGTAAGCGGCGCCGAATTCCAGATGgccgcc                                                                  <  7:36885/149‑1 (MQ=255)
                                                                                  ctgCATCGAAGCCGCGCGGGCGCTGTCCGAAGAAGTGCTGCTCGAAGGCCTGATCGAAACCCTGATGGGCCACCTGACCCAGCACTCCGGTGCCGACCATGGCGCGCTGCTGATTGTAAGCGGCGCCGAATTCCAGATGGCCGCCCTgg                                                              <  3:91371/147‑1 (MQ=255)
                                                                                  cGGCATCGAAGCCGCGCGGGCGCTGTCCGAAGAAGTGCTGCTCGAAGGCCTGATCGAAACCCTGATGGGCCACCTGACCCAGCACTCCGGTGCCGACCATGGCGCGCTGCTGATTGTAAGCGGCGCCGAATTCCAGATGGCCGCCCTgg                                                              >  7:145055/1‑149 (MQ=255)
                                                                                    gCATCGAAGCCGCGCGGGCGCTGTCCGAAGAAGTGCTGCTCGAAGGCCTGATCGAAACCCTGATGGGCCACCTGACCCAGCACTCCGGTGCCGACCATGGCGCGCTGCTGATTGTAAGCGGCGCCGAATTCCAGATGGCCGCCCTGGcc                                                            >  3:3875/1‑149 (MQ=255)
                                                                                      aTCGAAGCCGCGCGGGCGCTGTCCGAAGAAGTGCTGCTCGAAGGCCTGATCGAAACCCTGATGGGCCACCTGACCCAGCACTCCGGTGCCGACCATGGCGCGCTGCTGATTGTAAGCGGCGCCGAATTCCAGATGGCCGCCCTGGCCt                                                           >  6:3675/1‑148 (MQ=255)
                                                                                          aaGCCGCGCGGGCGCTGTCCGAAGAAGTGCTGCTCGAAGGCCTGATCGAAACCCTGATGGGCCACCTGACCCAGCACTCCGGTGCCGACCATGGCGCGCTGCTGATTGTAAGCGGCGCCGAATTCCAGATGGCCGCCCTGGCCTACATc                                                      >  2:194324/1‑149 (MQ=255)
                                                                                           aGCCTCGCGGGCGCTGTCCGAAGAAGTGCTGCTCGAAGGCCTGATCGAAACCCTGATGGGCCACCTGACCCAGCACTCCGGTGCCGACCATGGCGCgctgct                                                                                                    <  4:247330/102‑1 (MQ=255)
                                                                                           aGCCGCGCGGGCGCTGTCCGAAGAAGTGCTGCTCGAAGGCCTGATCGAAACCCTGATGGGCCACCTGACCCAGCACTCCGGTGCCGACCATGGCGCgctgct                                                                                                    >  3:247330/1‑102 (MQ=255)
                                                                                           aGCCGCGCGGGCGCTGTCCGAAGAAGTGCTGCTCGAAGGCCTGATCGAAACCCTGATGGGCCACCTGACCCAGCACTCCGGTGCCGACCATGGCGCGCTGCTGATTGTAAGCGGCGCCGAATTCCAGATGGCCGCCCTGGCCTACATcg                                                     <  2:289633/149‑1 (MQ=255)
                                                                                                cgcgGGCGCTGTCCGAAGAAGTGCTGCTCGAAGGCCTGATCGAAACCCTGATGGGCCACCTGACCCAGCACTCCGGTGCCGACCATGGCGCGCTGCTGATTGTAAGCGGCGCCGAATTCCAGATGGCCGCCCTGGCCTACATCGACGAt                                                <  1:240537/149‑1 (MQ=255)
                                                                                                       gcggTCCGAAGATGTGCTGCTCGAAGGCCTGATCGAAACCCTGATGGGCCACCTGACCCAGCACTCCGGTGCCGACCATGGCGCGCTGCTGATTGTAAGCGGCGCCGAATTCCAGATGGCCGCCCTGGCCTACATCGACGATGCCGGGc                                         <  8:79661/146‑1 (MQ=255)
                                                                                                       gcTGTCCGAAGAAGTGCTGCTCGAAGGCCTGAGCGAAACCCTGATGGGCCACCTGACCCAGCACTCCGGTGCCGACCATGGCGCGCTGCTGATTGTAAGCGGCGCCGAATTCCAGATGGCCGCCCTGGCCTACATCGACGATGCCGGGc                                         <  8:84925/149‑1 (MQ=255)
                                                                                                           tCCGAAGAAGTGCTGCTCGAAGGCCTGATCGAAACCCTGATGGGCCACCTGACCCAGCACTCCGGTGCCGACCATGGCGCGCTGCTGATTGTAAGCGGCGCCGAATTCCAGATGGCCGCCCTGGCCTACATCGACGATGCCGGGCTGc                                      <  7:212937/148‑1 (MQ=255)
                                                                                                              gaagaaGTGCTGCTCGAAGGCCTGATCGAAACCCTGATGGGCCACCTGACCCAGCACTCCGGTGCCGACCATGGCGCGCTGCTGATTGTAAGCGGCGCCGAATTCCAGATGGCCGCCCTGGCCTACATCGACGATGCCGGGCTGCACGt                                  >  6:346145/1‑149 (MQ=255)
                                                                                                              gaagaaGTGCTGCTCGAAGGCCTGATCGAAACCCTGATGGGCCACCTGACCCAGCACTCCGGTGCCGACCATGGCGCGCTGCTGATTGTAAGCGGCGCCGAATTCCAGATGGCCGCCCTGGCCTACATCGACGATGCCGGGCTGCACGt                                  <  8:145055/149‑1 (MQ=255)
                                                                                                               aagaagTGCTGCTCGAAGGCCTGATCGAAACCCTGATGGGCCACCTGACCCAGCACTCCGGTGCCGACCATGGCGCGCTGCTGATTGTAAGCGGCGCCGAATTCCAGATGGCCGCCCTGGCCTACATCGACGATGCCGGGCTGCAc                                    <  1:26410/146‑1 (MQ=255)
                                                                                                               aagaagTGCTGCTCGAAGGCCTGATCGAAACCCTGATGGGCCACCTGACCCAGCACTCCGGTGCCGACCATGGCGCGCTGCTGATTGTAAGCGGCGCCGAATTCCAGATGGCCGCCCTGGCCTACATCGACGATGCCGGGCTGCAc                                    >  2:26410/1‑146 (MQ=255)
                                                                                                               aagaagTGCTGCTCGAAGGCCTGATCGAAACCCTGATGGGCCACCTGACCCAGCACTCCGGTGCCGACCATGGCGCGCTGCTGATTGTAAGCGGCGCCGAATTCCAGATGGCCGCCCTGGCCTACATCGACGATGCCGGGCTGCACGt                                  >  3:82796/1‑148 (MQ=255)
                                                                                                               aagaagTGCTGCTCGAAGGCCTAATCGAAACCCTGATGGGCCACCTGACCCAGTACTCCGGTGCCGACCATGGCGCGCTGCTGATTGTAAGCGGCGCCGAATTCCAGATGGCCGCCCTGGCCTACATCGACGATGCCGGGCTGCACGt                                  >  4:44090/1‑148 (MQ=255)
                                                                                                                   agTGCTGCTCGAAGGCCTGATCGAAACCCTGATGGGCCACCTGACCCAGCACTCCGGTGCCGACCATGGCGCGCTGCTGATTGTAAGCGGCGCCGAATTCCAGATGGCCGCCCTGGCCTACATCGACGATGCCGGGCTGCACGTGAGCa                             <  3:108499/149‑1 (MQ=255)
                                                                                                                     tgctgcTCGAAGGCCTGATCGAAACCCTGATGGGCCACCTGACCCAGCACTCCGGTGCCGACCATGGCGCGCTGCTGATTGTAAGCGGCGCCGAATTCCAGATGGCCGCCCTGGCCTACATCGACGATGCCGGGCTGCACGTGAGCATg                           <  4:3875/149‑1 (MQ=255)
                                                                                                                     tgctgcTCGAAGGCCTGATCGAAACCCTGATGGGCCACCTGACCCAGCACTCCGGTGCCGACCATGGCGCGCTGCTGATTGTAAGCGGCGCCGAATTCCAGATGGCCGCCCTGGCCTACATCGACGATGCCGGGCTGCACGTGAGCATg                           <  4:134849/149‑1 (MQ=255)
                                                                                                                             gAAGGCCTGATCGAAACCCTGATGGGCCACCTGACCCAGCACTCCGGTGCCGACCATGGCGCGCTGCTGATTGTAAGCGGCGCCGAATTCCAGATGGCCGCCCTGGCCTACATCGACGATGCCGGGCTGCACGTGAGCATGGACAACTg                   <  3:196676/149‑1 (MQ=255)
                                                                                                                                ggCCTGATCGAAACCCTGATGGGCCACCTGACCCAGCACTCCGGTGCCGACCATGGCGCGCTGCTGATTGTAAGCGGCGCCGAATTCCAGATGGCCGCCCTGGCCTACATCGACGATGCCGGGCTGCACGTGAGCATGGACAACTGCCa                <  5:3675/149‑1 (MQ=255)
                                                                                                                                    tGATCGAAACCCTGAAGGGCCACCTGACCCAGCACTCCGGTGCCGACCATGGCGCGCTGCTGATTGTAAGCGGCGCCGAATTCCAGATGGCCGCCCTGGCCTACATCGACGATGcc                                             >  1:1118/1‑116 (MQ=255)
                                                                                                                                    tGATCGAAACCCTGAAGGGCCACCTGACCCAGCACTCCGGTGCCGACCATGGCGCGCTGCTGATTGTAAGCGGCGCCGAATTCCAGATGGCCGCCCTGGCCTACATCGACGATGcc                                             >  1:7338/1‑116 (MQ=255)
                                                                                                                                    tGATCGAAACCCTGAAGGGCCACCTGACCCAGCACTCCGGTGCCGACCATGGCGCGCTGCTGATTGTAAGCGGCGCCGAATTCCAGATGGCCGCCCTGGCCTACATCGACGATGcc                                             <  2:1118/116‑1 (MQ=255)
                                                                                                                                    tGATCGAAACCCTGAAGGGCCACCTGACCCAGCACTCCGGTGCCGACCATGGCGCGCTGCTGATTGTAAGCGGCGCCGAATTCCAGATGGCCGCCCTGGCCTACATCGACGATGcc                                             <  2:7338/116‑1 (MQ=255)
                                                                                                                                       tCGAAACCCTGATGGGCCACCTGACCCAGCACTCCGGTGCCGACCATGGCGCGCTGCTGATTGTAAGCGGCGCCGAATTCCAGATGGCCGCCCTgg                                                              <  3:11833/96‑1 (MQ=255)
                                                                                                                                       tCGAAACCCTGATGGGCCACCTGACCCAGCACTCCGGTGCCGACCATGGCGCGCTGCTGATTGTAAGCGGCGCCGAATTCCAGATGGCCGCCCTgg                                                              >  4:11833/1‑96 (MQ=255)
                                                                                                                                         gAAACCCTGATGGGCCACCTGACCCAGCACTCCGGTGCCGACCATGGCGCGCTGCTGATTGTAAGCGGCGCCGAATTCCAGATGGCCGCCCTGGCCTACATCGACGATGCCGGGCTGCACGTGAGCGTGGACAACTGCCAGAAGTTCAt       <  1:194324/149‑1 (MQ=255)
                                                                                                                                          aaaCCCTGATGGGCCACCTGACCCAGCACTCCGGTGCCGACCATGGCGCGCTGCTGATTGTAAGCGGCGCCGAATTCCAGATGGCCGCCCTGGCCTACATCGACGATGCCGGGCTGCACGTGAGCATGGACAACTGCCAGAAGTTCAt       >  3:279636/1‑148 (MQ=255)
                                                                                                                                           aaCCCTGATGGGCCACCTGACCCAGCACTCCGGTGCCGACCATGGCGCGCTGCTGATTGTAAGCGGCGCCGAATTCCAGATGGCCGCCCTGGCCTACATCGACGATGCCGGGCTGCACGTGAGCATGGACAACTGCCAGAAGTTCATGa     >  2:35826/1‑149 (MQ=255)
                                                                                                                                              cgtGATCGGCAACCTGACCCAGCACTCCGGTGCCGACCATGTCGCGCTGCTGATTGTAAGCGGCGCCGAATTCCAGATGGCCGCCCTGGCCTACATCGACGATGCCGGGCTGCACGTGAGCATGGACAACTGCCAGAAGTTCATGAccc  <  4:279636/147‑1 (MQ=255)
                                                                                                                                              ccTGATGGGCCACCTGACCCAGCACTCCGGTGCCGACCATGGCGCGCTGCTGATTGTAAGCGGCGCCGAATTCCAGATGGCCGCCCTGGCCTACATCGACGATGCCGGGCTGCACGTGAGCATGGACAACTGCCAGAAGTTCATGAcc   >  4:16342/1‑148 (MQ=255)
                                                                                                                                              ccTGATGGGCCACCTGACCCAGCACTCCGGTGCCGACCATGGCGCGCTGCTGATTGTAAGCGGCGCCGAATTCCAGATGGCCGCCCTGGCCTACATCGACGATGCCGGGCTGCACGTGAGCATGGACAACTGCCAGAAGTTCATGAcc   >  1:249190/1‑148 (MQ=255)
                                                                                                                                                   |                                                                                                                                               
AGACGTTGCATCCGTTTTTGCGTACCCAGCCGATTCAGGAGACCTATCGCGCCACCCACCAGGCCAGGCTGGACCTCGAAGCCGGCATCGAAGCCGCGCGGGCGCTGTCCGAAGAAGTGCTGCTCGAAGGCCTGATCGAAACCCTGATGGGCCACCTGACCCAGCACTCCGGTGCCGACCATGGCGCGCTGCTGATTGTAAGCGGCGCCGAATTCCAGATGGCCGCCCTGGCCTACATCGACGATGCCGGGCTGCACGTGAGCATGGACAACTGCCAGAAGTTCATGACCC  >  NC_002947/5983074‑5983364

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