Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 2,246,242 A→G 8.3% intergenic (+31/+60) ibpA → / ← PP_1983 small heat shock protein IbpA/PAS/PAC/GAF sensor‑containing diguanylate cyclase/phosphodiesterase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029472,246,2420AG 8.3% 203.0 / 6.9 72intergenic (+31/+60)ibpA/PP_1983small heat shock protein IbpA/PAS/PAC/GAF sensor‑containing diguanylate cyclase/phosphodiesterase
Reads supporting (aligned to +/- strand):  ref base A (35/31);  new base G (3/3);  total (38/34)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.90e-01

AAGGCCGCCGGCCTGGCCAATGGCCTGCTCAGCATCGACCTGCTGCGCATCGTGCCTGAAGAAGCCAAGGCCAAGCGCATTCCGATCAATGGCGAAAAACCAGCGCTGAACTGATCGAGCTGCAACGAAAAAAAGGGCACCTCCAGAGGTGCCCTTTTTCGTGCTTCTCACGCAGTCCCGCGTAGGTGCGCGGCCGTCATGCGATCAGTAATCTGCCGGCGCCTCCAGCAACATCTGCCGAAACTCCGGTAACGGTAATGGCCGGCTGTGCAGATAACCCTGGT  >  NC_002947/2246098‑2246381
                                                                                                                                                |                                                                                                                                             
aaGGCCGCCGGCCTGGCCAATGGCCTGCTCAGCATCGACCTGCTGCGCATCGTGCCTGAAGAAGCCAAGGCCAAGCGCATTCCGATCAATGGCGAAAAACCAGCGCTGAACTGATCGAGCTGCAACGAAAAAAAGGGCACCTCCagag                                                                                                                                           <  6:247780/148‑1 (MQ=255)
   gccgccGGCCTGGCCAATGGCCTGCTCAGCATCGACCTGCTGCGCATCGTGCCTGAAGAAGCCAAGGCCAAGCGCATTCCGATCAATGGCGAAAACCCAGCGCTGAACTGATCGAGCTGCAACGAAAAAAAGGGCACCTCCAGAGGTGc                                                                                                                                       >  3:41163/1‑149 (MQ=255)
   gccgccGGCCTGGCCAATGGCCTGCTCAGCATCGACCTGCTGCGCATCGTGCCTGAAGAAGCCAAGGCCAAGCGCATTCCGATCAATGGCGAAAAACCAGCGCTGAACTGATCGAGCTGCAACGAAAAAAAGGGCACCTCCAGAGGTGc                                                                                                                                       <  3:4454/149‑1 (MQ=255)
     cgcctGCCTGGCCAATGGCCTGCTCAGCATCGACCTGATGCGCATCGTGCCTGAAGAAGCCAAGGCCAAGCGCATTCCGATCAAAGGCGAAAAACCAGCGCTGAACTGATCGAGCTGCAACGAAAAAAAGGGCACCTCCAGAGGTGccc                                                                                                                                     <  5:169090/149‑1 (MQ=255)
     cgccgGCCTGGCCAATGGCCTGCTCAGCATCGACCTGCTGCGCATCGTGCCTGAAGAAGCCAAGGCCAAGCGCATTCCGATCAATGGCGAAAAACCAGCGCTGAACTGATCGAGCTGCAACGAAAAAAAGGGCACCTCCAGAGGTGccc                                                                                                                                     <  7:134011/149‑1 (MQ=255)
     cgccgGCCTGGCCAATGGCCTGCTCAGCATCGACCTGCTGCGCATCGTGCCTGAAGAAGCCAAGGCCAAGCGCATTCCGATCAATGGCGAAAAACCAGCGCTGAACTGATCGAGCTGCAACGAAAAAAAGGGCACCTCCAGAGGTGcc                                                                                                                                      >  1:154947/1‑148 (MQ=255)
     cgccgGCCTGGCCAATGGCCTGCTCAGCATCGACCTGCTGCGCATCGTGCCTGAAGAAGCCAAGGCCAAGCGCATTCCGATCAATGGCGAAAAACCAGCGCTGAACTGATCGAGCTGCAACGAAAAAAAGGGCACCTCCAGAGGTGcc                                                                                                                                      <  2:154947/148‑1 (MQ=255)
       ccgGCCTGGCCAATGGCCTGCTCAGCATCGACCTGCTGCGCATCGTGCCTGAAGAAGCCAAGGCCAAGCGCATTCCGATCAATGGCGAAAAACCAGCGCTGAACTGATCGAGCTGCAACGAAAAAAAGGGCACCTCCAGAGGTGCCCtt                                                                                                                                   >  3:896/1‑149 (MQ=255)
           ccTGGCCAATGGCCTGCTCAGCATCGACCTGCTGCGCATCGTGCCTGAAGAAGCCAAGGCCAAGCGCATTCCGATCAATGGCGAAAAACCAGCGCTGAACTGATCGAGCTGCAACGAAAAAAAGGGCACCTCCAGAGGTGCCCTTTTTc                                                                                                                               >  3:188578/1‑149 (MQ=255)
              ggCCAATGGCCTGCTCAGCATCGACCTGCTGCGCATCGTGCCTGAAGAAGCCAAGGCCAAGCGCATTCCGATCAATGGCGAAAAACCAGCGCTGAACTGATCGAGCTGCAACGAAAAAAAGGGCACCTCCAGAGGTGCCCTTTTTCGt                                                                                                                             >  2:222293/1‑148 (MQ=255)
               gCCAATGGCCTGCTCAGCATCGACCTGCTGCGCATCGTGCCTGAAGAAGCCAAGGCCAAGCGCATTCCGATCAATGGCGAAAAACCAGCGCTGAACTGATCGAGCTGCAACGAAAAAAAGGGCACCTCCAGAGGTGCCCTTTTTCGTg                                                                                                                            <  4:896/148‑1 (MQ=255)
               gCCAATGGCCTGCTCAGCATCGACCTGCTGCGCATCGTGCCTGAAGAAGCCAAGGCCAAGCGCATTCCGATCAATGGCGAAAAACCAGCGCTGAACTGATCGAGCTGCAACGAAAAAAAGGGCACCTCCAGAGGTGCCCTTTTTCGTGc                                                                                                                           >  3:58348/1‑149 (MQ=255)
                 cAATGGCCTGCTCAGCATCGACCTGCTGCGCATCGTGCCTGAAGAAGCCAAGGCCAAGCGCATTCCGATCAATGGCGAAAAACCAGCGCTGAACTGATCGAGCTGCAACGAAAAAAAGGGCACCTCCAGAGGTGCCCTTTTTCGTGCtt                                                                                                                         <  2:297878/149‑1 (MQ=255)
                  aaTGGCCTGCTCAGCATCGACCTGCTGCGCATCGTGCCTGAAGAAGCCAAGGCCAAGCGCATTCCGATCAATGGCGAAAAACCAGCGCTGAACTGATCGAGCTGCAACGAAAAAAAGGGCACCTCCAGAGGTGCCCTTTTTCGTg                                                                                                                            <  8:244693/145‑1 (MQ=255)
                  aaTGGCCTGCTCAGCATCGACCTGCTGCGCATCGTGCCTGAAGAAGCCAAGGCCAAGCGCATTCCGATCAATGGCGAAAAACCAGCGCTGAACTGATCGAGCTGCAACGAAAAAAAGGGCACCTCCAGAGGTGCCCTTTTTCGTg                                                                                                                            >  7:244693/1‑145 (MQ=255)
                        cTGCTCAGCATCGACCTGCTGCGCATCGTGCCTGAAGAAGCCAAGGCCAAGCGCATTCCGATCAATGGCGAAAAACCAGCGCTGAACTGATCGAGCTGCAACGAAAAAAAGGGCACCTCCAGAGGTGCCCTTTTTCGTGCTTCTCACg                                                                                                                   >  4:104868/1‑148 (MQ=255)
                        cTGCTCAGCATCGACCTGCTGCGCATCGTGCCTGAAGAAGCCAAGGCCAAGCGCATTCCGATCAATGGCGAAAAACCAGCGCTGAACTGATCGAGCTGCAACGAAAAAAAGGGCACCTCCAGAGGTGCCCTTTTTCGTGCTTCTCACGc                                                                                                                  <  8:168588/149‑1 (MQ=255)
                         tGCTCAGCATCGACCTGCTGCGCATCGTGCCTGAAGAAGCCAAGGCCAAGCGCATTCCGATCAATGGCGAAAAACCAGCGCTGAACTGATCGAGCTGCAACGAAAAAAAGGGCACCTCTGGAGGTGccc                                                                                                                                     <  3:126864/129‑1 (MQ=255)
                         tGCTCAGCATCGACCTGCTGCGCATCGTGCCTGAAGAAGCCAAGGCCAAGCGCATTCCGATCAATGGCGAAAAACCAGCGCTGAACTGATCGAGCTGCAACGAAAAAAAGGGCACCTCTGGAGGTGccc                                                                                                                                     >  4:126864/1‑129 (MQ=255)
                         tGCTCAGCATCGACCTGCTGCGCATCGTGCCTGAAGAAGCCAAGGCCAAGCGCATTCCGATCAATGGCGAAAAACCAGCGCTGAACTGATCGAGCTGCAACGAAAAAAAGGGCACCTCCAGAGGTGCCCTTTTTCGTGCTTCTCACg                                                                                                                   >  4:131506/1‑147 (MQ=255)
                              aGCATCAACCTGCTGCGCATCGTGCCTGAAGAAGCCAAGGCCAAGCGCATTCCGATCAATGGCGAAAAACCAGCGCTGAACTGATCGAGCTGCAACGAAAAAAAGGGCACCTCCAGAGGTGCCCTTTTTCGTGCTTCTCACGCAGTc                                                                                                              >  7:249030/1‑147 (MQ=255)
                                    gACCTGCTGCGCATCGTGCCTGAAGAAGCCAAGGCCAAGCGCATTCCGATCAATGGCGAAAAACCAGCGCTGAACTGATCGAGCTGCAACGAAAAAAAGGGCACCTCCAGAGGTGCCCTTTTTCGTGCTTCTCACGCAGTCCCGCGTAg                                                                                                     <  1:272983/149‑1 (MQ=255)
                                    gACCTGCTGCGCATCGTGCCTGAAGAAGCCAAGGCCAAGCGCATTCCGATCAATGGCGAAAAACCAGCGCTGAACTGATCGAGCTGCAACGAAAAAAAGGGCACCTCCAGAGGTGCCCTTTTTCGTGCTTCTCACGCAGTCCCGCGTAg                                                                                                     <  6:54947/149‑1 (MQ=255)
                                    gACCTGCTGCGCATCGTGCCTGAAGAAGCCAAGGCCAAGCGCATTCCGATCAATGGCGAAAAACCAGCGCTGAACTGATCGAGCTGCAACGAAAAAAAGGGCACCTCCAGAGGTGCCCTTTTTCGTGCTTCTCACGCAGTCCCGCGTAg                                                                                                     <  4:58348/149‑1 (MQ=255)
                                     aCCTGCTGCGCATCGTGCCTGAAGAAGCCAAGGCCAAGCGCATTCCGATCAATGGCGAAAAACCAGCGCTGAACTGATCGAGCTGCAACGAAAAAAAGGGCACCTCTGGAGGTGCCCtttttttcgt                                                                                                                           >  3:26466/1‑122 (MQ=255)
                                     aCCTGCTGCGCATCGTGCCTGAAGAAGCCAAGGCCAAGCGCATTCCGATCAATGGCGAAAAACCAGCGCTGAACTGATCGAGCTGCAACGAAAAAAAGGGCACCTCTGGAGGTGCCCtttttttcgt                                                                                                                           <  4:26466/127‑6 (MQ=255)
                                     aCCTGCTGCGCATCGTGCCTGAAGAAGCCAAGGCCAAGCGCATTCCGATCAATGGCGAAAAACCAGCGCTGAACTGATCGAGCTGCAACGAAAAAAAGGGCACCTCCagag                                                                                                                                           <  7:107342/111‑1 (MQ=255)
                                     aCCTGCTGCGCATCGTGCCTGAAGAAGCCAAGGCCAAGCGCATTCCGATCAATGGCGAAAAACCAGCGCTGAACTGATCGAGCTGCAACGAAAAAAAGGGCACCTCCagag                                                                                                                                           >  8:107342/1‑111 (MQ=255)
                                      ccTGCTGCGCATCGTGCCTGAAGAAGCCAAGGCCAAGCGCATTCCGATCAATGGCGAAAAACCAGCGCTGAACTGATCGATCTGCAACGAAAAAAAGGGCACCTCCAGAGGTGCCCTTTTTCGTGCTTCTCACGCAGTCCCGCGTAGGt                                                                                                   <  8:14595/149‑1 (MQ=255)
                                      ccTGCTGCGCATCGTGCCTGAAGAAGCCAAGGCCAAGCGCATTCCGATCAATGGCGAAAAACCAGCGCTGAACTGATCGAGCTGCAACGAAAAAAAGGGCACCTCCAGAGGTGCCCTTTTTCGTGCTTCTCACGCAGTCCCGCGTa                                                                                                      <  2:143006/146‑1 (MQ=255)
                                      ccTGCTGCGCATCGTGCCTGAAGAAGCCAAGGCCAAGCGCATTCCGATCAATGGCGAAAAACCAGCGCTGAACTGATCGAGCTGCAACGAAAAAAAGGGCACCTCCAGAGGTGCCCTTTTTCGTGCTTCTCACGCAGTCCCGCGTa                                                                                                      >  1:143006/1‑146 (MQ=255)
                                         gctgcGCATCGTGCCTGAAGAAGCCAAGGCCAAGCGCATTCCGATCAATGGCGAAAAACCAGCGCTGAACTGATCGAGCTGCAACGAAAAAAAGGGCACCTCCAGAGGTGCCCTTTTTCGTGCTTCTCACGCAGTccc                                                                                                            >  7:280660/1‑138 (MQ=255)
                                         gctgcGCATCGTGCCTGAAGAAGCCAAGGCCAAGCGCATTCCGATCAATGGCGAAAAACCAGCGCTGAACTGATCGAGCTGCAACGAAAAAAAGGGCACCTCCAGAGGTGCCCTTTTTCGTGCTTCTCACGCAGTccc                                                                                                            <  8:280660/138‑1 (MQ=255)
                                          ctgcGCATCGTGCCTGAAGAAGCCAAGGCCAAGCGCATTCCGATCAATGGCGAAAAACCAGCGCTGAACTGATCGAGCTGCAACGAAAAAAAGGGCACCTCCAGAGGTGCCCTTTTTCGTGCTTCTCACGCAGTCCCGCGTAGGTgcg                                                                                                <  4:188578/148‑1 (MQ=255)
                                            gcgcATCGTGCCTGAAGAAGCCAAGGCCAAGCGCATTCCGATCAATGGCGAAAAACCAGCGCTGAACTGATCGAGCTGCAACGAAAAAAAGGGCACCTCCAGAGGTGCCCTTTTTCGTGCTTCTCACGCAGTCCCGCGTAGGTGCGCgg                                                                                             >  1:162094/1‑149 (MQ=255)
                                             cgcATCGTGCCTGAAGAAGCCAAGGCCAAGCGCATTCCGATCAATGGCGAAAAACCAGCGCTGAACTGATCGAGCTGCAACGAAAAAAAGGGCACCTCCACGAGGTGCCCTTTTTCGTGCTTCTCACGCAGTCCCGCGTAGGTgcgc                                                                                               <  1:222293/147‑1 (MQ=255)
                                             cgcATCGTGCCTGAAGAAGCCAAGGCCAAGCGCATTCCGATCAATGGCGAAAAACCAGCGCTGAACTGATCGAGCTGCAACGAAAAAAAGGGCACCTCCAGAGGTGCCCTTTTTCGTGCTTCTCACGCAGTCCCGCGTAGGTGCGCgg                                                                                             <  1:236448/148‑1 (MQ=255)
                                                aTCGTGCCTGAAGAAGCCAAGGCCAAGCGCATTCCGATCAATGGCGAAAAACCAGCGCTGAACTGATCGAGCTGCAACGAAAAAAAGGGCACCTCCAGAGGTGCCCTTTTTCGTGCTTCTCACGCAGTCCCGCGTAGGTGCGCGGCCGt                                                                                         >  5:282978/1‑149 (MQ=255)
                                                aTCGTGCCTGAAGAAGCCAAGGCCAAGCGCATTCCGATCAATGGCGAAAAACCAGCGCTGAACTGATCGAGCTGCAACGAAAAAAAGGGCACCTCCAGAGGTGCCCTTTTTCGTGCTTCTCACGCAGTCCCGCGTAGGTGCGCGGCCGt                                                                                         >  5:184573/1‑149 (MQ=255)
                                                 tCGTTCCTGAAGAAGCCAAGTCCAAGCGCATTCCGATCAATCGCGAAAAACCAGCGCTGAACTGATCGAGCTGTAACGAAAAAAAGGGCACCTCCAGAGGTGCCCTTTTTCGTGCTTCTCACGCAGTCCCGCGTAGGTGCGCGGCCGTc                                                                                        <  3:104/149‑1 (MQ=255)
                                                      ccTGAAGAAGCCAAGGCCAAGCGCATTCCGATCAATGGCGAAAAACCAGCGCTGAACTGATCGAGCTGCAACGAAAAAAAGGGCACCTCCaga                                                                                                                                            <  5:88192/93‑1 (MQ=255)
                                                      ccTGAAGAAGCCAAGGCCAAGCGCATTCCGATCAATGGCGAAAAACCAGCGCTGAACTGATCGAGCTGCAACGAAAAAAAGGGCACCTCCaga                                                                                                                                            >  6:88192/1‑93 (MQ=255)
                                                      ccTGAAGAAGCCAAGGCCAAGCGCATTCCGATCAATGGCGAAAAACCAGCGCTGAACTGATCGAGCTGCAACGAAAAAAAGGGCACCTCCAGAGGTGCCCTTTTTCGTg                                                                                                                            <  2:165300/109‑1 (MQ=255)
                                                      ccTGAAGAAGCCAAGGCCAAGCGCATGCCGATCAATGGCGAAAAACCAGCGCTGAACTGATCGAGCTGCAACGAAAAAAAGGGCACCTCCAGAGGTGCCCTTTTTCGTg                                                                                                                            >  1:165300/1‑109 (MQ=255)
                                                         gaagaaGCCAAGGCCAAGCGCATTCCGATCAATGGCGAAAAACCAGCGCTGAACTGATCGAGCTGCAACGAAAAAAAGGGCACCTCCAGAGGTGCCCTTTTTCGTGCTTCTCACGCAGTCCCGCGTAGGTGCGCGGCCGTCATGCGAt                                                                                 >  4:109755/1‑148 (MQ=255)
                                                         gaagaaGCCAAGGCCAAGCGCATTCCGATCAATGGCGAAAAACCAGCGCTGAACTGATCGAGCTGCAACGAAAAAAAGGGCACCTCCAGAGGTGCCCTTTTTCGTGCTTCTCACGCAGTCCCGCGTAGGTGCGCGGCCGTCATGCGAt                                                                                 >  2:35010/1‑148 (MQ=255)
                                                          aagaagCCAAGGCCAAGCGCATTCCGATCAATGGCGAAAAACCAGCGCTGAACTGATCGAGCTGCAACGAAAAAAAGGGCACCTCCAGAGGTGCCCTTTTTCGTGCTTCTCACGCAGTCCCGCGTAGGTGCGCGGCCGTCATGCGATc                                                                                <  7:7753/148‑1 (MQ=255)
                                                          aagaagCCAAGGCCAAGCGCATTCCGATCAATGGCGAAAAACCAGCGCTGAACTGATCGAGCTGCAACGAAAAAAAGGGCACCTCCAGAGGTGCCCTTTTTCGTGCTTCTCACGCAGTCCCGCGTAGGTGCGCGGCCGTCATGCGATc                                                                                <  2:162094/148‑1 (MQ=255)
                                                          aagaagCCAAGGCCAAGCGCATTCCGATCAATGGCGAAAAACCAGCGCTGAACTGATCGAGCTGCAACGAAAAAAAGGGCACCTCCAGAGGTGCCCTTTTTCGTGCTTCTCACGCAGTCCCGCGTAGGTGCGCGGCCGTCATGCGATCa                                                                               <  3:104868/149‑1 (MQ=255)
                                                                ccaaggccaagCGCATTCCGATCAATGGCGAAAAACCAGCGCTGAACTGATCGAGCTGCAACGAAAAAAAGGGCACCTCCAGAGGTGCCCTTTTTCGTGCTTCTCACGCAGTCCCGCGTAGGTGCGCGGCCGTCATGCGATCAGTAATc                                                                         >  4:282582/1‑149 (MQ=255)
                                                                       caagCGCATTCCGATCAATGGCGAAAAACCAGCGCTGAACTGATCGAGCTGCAACGAAAAAAAGGGCACCTCCAGAGGTGCCCTTTTTCGTGCTTCTCACGCAGTCCCGCGTAGGTGCGCGGCCGTCATGCGATCAGTAATCTGCCGgc                                                                  <  6:184573/149‑1 (MQ=255)
                                                                        aagCGCATTCCGATCAATGGCGAAAAACCAGCGCTGAACTGATCGAGCTGCAACGAAAAAAAGGGCACCTCCAGAGGTGCCCTTTTTCGTGCTTCTCACGCAGTCCCGCGTAGGTGCGCGGCCGTCATGCGATCAGTAATCTGCCGgcg                                                                 >  8:160348/1‑149 (MQ=255)
                                                                             cATTCCGATCAATGGCGAAAAACCAGCGCTGAACTGATCGAGCTGCAACGAAAAAAAGGGCACCTCCAGAGGTGCCCTTTTTCGTGCTTCTCACGCAGTCCCGCGTAGGTGCGCGGCCGTCATGCGATCAGTAATCTGCCGGCGCCTcc                                                            <  3:131506/149‑1 (MQ=255)
                                                                                  cGATCAATGGCGAAAAACCAGCGCTGAACTGATCGAGCTGCAACGAAAAAAAGGGCACCTCCAGAGGTGCCCTTTTTCGTGCTTCTCACGCAGTCCCGCGTAGGTGCGCGGCCGTCATGCGATCAGTAATCTGCCGGCGCCTCCAGCaa                                                       >  2:81488/1‑149 (MQ=255)
                                                                                         tGGCGAAAAACCAGCGCTGAACTGATCGAGCTGCAACGAAAAAAAGGGCACCTCCAGAGGTGCCCTTTTTCGTGCTTCTCACGCAGTCCCGCGTAGGTGAGCGGCCGTCATGCGATCAGTAATCTGCCGGCGCCTCCAGCAACATCTGc                                                >  1:166237/1‑149 (MQ=255)
                                                                                         tGGCGAAAAACCAGCGCTGAACTCATCGAGCTGCAACGAAAAAAAGGGCACCTCCATAGGTGCCCTTTTTCGTGCTTCTCACGCAGTCCAGCGTAGGTGCGCGGCCGTCATGCGATCAGTAATCTGCCGGCGCCTCCAGCAACATCTGc                                                >  7:54484/1‑149 (MQ=255)
                                                                                            cGAAAAACCAGCGCTGAACTGATCGAGCTGCAACGAAAAAAAGGGCACCTCCAGAGGTGCCCTTTTTCGTGCTTCTCACGCAGTCCCGCGTAGGTGCGCGGCCGTCATGCGATCAGTAATCTGCCGGCGCCTCCAGCAACATCTGCCg                                              >  1:161530/1‑148 (MQ=255)
                                                                                                    cAGCGCTGAACTGATCGAGCTGCAACGAAAAAAAGGGCACCTCCAGAGGTGCCCTTTTTCGTGCTTCTCACGCAGTCCCGCGTAGGTGCGCGGCCGTCATGCGATCAGTAATCTGCCGGCGCCTCCAGCAACATCTGCCGAAACTCcgg                                     >  3:8746/1‑149 (MQ=255)
                                                                                                      gcgcTGAACTGATCGAGCTGCAACGAAAAAAAGGGCACCTCCAGAGGTGCCCTTTTTCGTGCTTCTCACGCAGTCCCGCGTAGGTGCGCGGCCGTCATGCGATCAGTAATCTGCCGGCGCCTCCAGCAACATCTGCCGAAACTCcgg                                     >  1:272839/1‑147 (MQ=255)
                                                                                                      gcgcTGAACTGATCGAGCTGCAACGAAAAAAAGGGCACCTCCAGAGGTGCCCTTTTTCGTGCTTCTCACGCAGTCCCGCGTAGGTGCGCGGCCGTCATGCGATCAGTAATCTGCCGGCGCCTCCAGCAACATCTGCCGAAACTCcgg                                     <  2:272839/147‑1 (MQ=255)
                                                                                                       cgcTGAACTGATCGAGCTGCAACGAAAAAAAGGGCACCTCCAGAGGTGCCCTTTTTCGTGCTTCTCACGCAGTCCCGCGTAGGTGCGCGGCCGTCATGCGATCAGTAATCTGCCGGCGCCTCCAGCAACATCTGCCGAAACTCcg                                      <  3:248597/145‑1 (MQ=255)
                                                                                                       cgcTGAACTGATCGAGCTGCAACGAAAAAAAGGGCACCTCCAGAGGTGCCCTTTTTCGTGCTTCTCACGCAGTCCCGCGTAGGTGCGCGGCCGTCATGCGATCAGTAATCTGCCGGCGCCTCCAGCAACATCTGCCGAAACTCcg                                      >  4:248597/1‑145 (MQ=255)
                                                                                                               tGATCGAGCTGCAACGAAAAAAAGGGCACCTCCAGAGGTGCCCTTTTTCGTGCTTCTCACGCAGTCCCGCGTAGGTGCGCGGCCGTCATGCGATCAGTAATCTGCCGGCGCCTCCAGCAACATCTGCCGAAACTCcggtaacggta                             >  1:285522/1‑146 (MQ=255)
                                                                                                               tGATCGAGCTGCAACGAAAAAAAGGGCACCTCCAGAGGTGCCCTTTTTCGTGCTTCTCACGCAGTCCCGCGTAGGTGCGCGGCCGTCATGCGATCAGTAATCTGCCGGCGCCTCCAGCAACATCTGCCGAAACTCcggtaacggta                             <  2:285522/146‑1 (MQ=255)
                                                                                                               tGATCGAGCTGCAACGAAAAAAAGGGCACCTCCAGAGGTGCCCTTTTTCGTGCTTCTCACGCAGTCCCGCGTAGGTGCGCGGCCGTCATGCGATCAGTAATCTGCCGGCGCCTCCAGCAACATCTGCCGAAACTCCGGTAACGGTAATg                          >  3:187473/1‑149 (MQ=255)
                                                                                                                  tCGAGCTGCAACGAAAAAAAGGGCACCTCCAGAGGTGCCCTTTTTCGTGCTTCTCACGCAGTCCCGCGTAGGTGCGCGGCCGTCATGCGATCAGTAATCTGCCGGCGCCTCCAGCAACATCTGc                                                >  6:53566/1‑124 (MQ=255)
                                                                                                                    gAGCTGCAACGAAAAAAAGGGCACCTCCAGAGGTGCCCTTTTTCGTGCTTCTCACGCAGTCCCGCGTAGGTGCGCGGCCGTCATGCGATCAGTAATCTGCCGGCGCCTc                                                             >  7:232537/1‑109 (MQ=255)
                                                                                                                    gAGCTGCAACGAAAAAAAGGGCACCTCCAGAGGTGCCCTTTTTCGTGCTTCTCACGCAGTCCCGCGTAGGTGCGCGGCCGTCATGCGATCAGTAATCTGCCGGCGCCTc                                                             <  8:232537/109‑1 (MQ=255)
                                                                                                                            aCGAAAAAAAGGGCACCTCCAGAGGTGCCCTTTTTCGTGCTTCTCACGCAGTCCCGCGTAGGTGCGCGGCCGTCATGCGATCAGTAATCTGCCGGCGCCTCCAGCAACATCTGCCGAAACTCCGGTAACGGTAATGGCCGGCtgtg                >  5:205957/1‑146 (MQ=255)
                                                                                                                             cGAAAAAAAGGGCACCTCCAGAGGTGCCCTTTTTCGTGCTTCTCACGCAGTCCCGCGTAGGTGCGCGGCCGTCATGCGATCAGTAATCTGCCGGCGCCTCCAGCAACATCTGCCGAAACTCCGGTAACGGTAATGGCCGGCTGTGCAGa            <  3:282582/149‑1 (MQ=255)
                                                                                                                                aaaaaaGGGCACCTCCAGAGGTGCCCTTTTTCGTGCTTCTCACGCAGTCCCGCGTAGGTGCGCGGCCGTCATGCGATCAGTAATCTGCCGGCGCCTCCAGCAACATCTGCCGAAACTCCGGTAACGGTAATGGCCGGCTGTGCAGATaa         >  1:199928/1‑149 (MQ=255)
                                                                                                                                 aaaaaGGGCACCTCCAGAGGTGCCCTTTTTCGTGCTTCTCACGCAGTCCCGCGTAGGTGCGCGGCCGTCATGCGATCAGTAATCTGCCGGCGCCTCCAGCAACATCTGCCGAAACTCCGGTAACGGTAATGGCCGGCTGTGCAGATAAc        >  7:167425/1‑149 (MQ=255)
                                                                                                                                  aaaaGGGCACCTCCAGAGGTGCCCTTTTTCGTGCTTCTCACGCAGTCCCGCGTAGGTGCGCGGCCGTCATGCGATCAGTAATCTGCCGGCGCCTCCAGCAACATCTGCCGAAACTCCGGTAACGGTAATGGCCGGCTGTGCAGATAAcc       <  2:166237/149‑1 (MQ=255)
                                                                                                                                      gggCACCTCTGGAGGTGCCCTTTTTCGTGCTTCTCACGCAGTCCCGGCGTAGGTGCGCGGCCGTCATGCGATCAGTAATCTGCCGGCGCATCCAGCAACATCTGCCGAAACTCCGGTAACGGTAATGGCCGGCTGTGCAGATACCCCta    >  1:218401/1‑148 (MQ=255)
                                                                                                                                      gggCACCTCTGGAGGTGCCCTTTTTCGTGCTTCTCACGCAGTCCCGCGTAGGTGCGCGGCCGTCATGCGATCAGTAATCTGCCGGCGCCTCCAGCAACATCTGCCGAAACTCCGGTAACGGTAATGGCCGGCTGTGCAGATAAccc      <  2:218401/146‑1 (MQ=255)
                                                                                                                                      gggCACCTCCAGAGGTGCCCTTTTTCGTGCTTCTCACGCAGTCCCGCGTAGGTGCGCGGCCGTCATGCGATCAGTAATCTGCCGGCGCCTCCAGCAACATCTGCCGAAACTCCGGTAACGGTAATGACCGGCTGTGCAGATAACCCTgg   <  1:35010/149‑1 (MQ=255)
                                                                                                                                        gCACCTCCAGAGGTGCCCTTTTTCGTGCTTCTCACGCAGTCCCGCGTAGGTGCGCGGCCGTCATGCGATCAGTAATCTGCCGGCGCCTCCAGCAACATCTGCCGAAATAACCGTAACGCTAAATGCCGGCATTGCAGAAAACCCAGGt  >  1:257088/1‑148 (MQ=255)
                                                                                                                                                |                                                                                                                                             
AAGGCCGCCGGCCTGGCCAATGGCCTGCTCAGCATCGACCTGCTGCGCATCGTGCCTGAAGAAGCCAAGGCCAAGCGCATTCCGATCAATGGCGAAAAACCAGCGCTGAACTGATCGAGCTGCAACGAAAAAAAGGGCACCTCCAGAGGTGCCCTTTTTCGTGCTTCTCACGCAGTCCCGCGTAGGTGCGCGGCCGTCATGCGATCAGTAATCTGCCGGCGCCTCCAGCAACATCTGCCGAAACTCCGGTAACGGTAATGGCCGGCTGTGCAGATAACCCTGGT  >  NC_002947/2246098‑2246381

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: