Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,963,789 A→T 7.7% M1M (ATG→TTG) † yegE → phosphodiesterase domain‑containing diguanylate cyclase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029471,963,7890AT 7.7% 230.9 / 7.1 78M1M (ATG→TTG) †yegEphosphodiesterase domain‑containing diguanylate cyclase
Reads supporting (aligned to +/- strand):  ref base A (37/35);  new base T (4/2);  total (41/37)
Fisher's exact test for biased strand distribution p-value = 6.78e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.79e-01

TGGGAGCGGGCTTGCCCGCGAACACGGGCGAAGCCCGTGCCATCCACCGCGGTGTCTTCTTCGCGGGTAAACCCGCTCCTACAATAGCCAGCGCATACGCTGTACAAATTCCCCGATTCCTTGTCATAAAAATCCCCTCCACAGCGATGGCAAATGGCCTACGCTGTGCTTGTCTGGGTCGATTTACGCCTTATCTGCTCTGTGACAGGGACTCCATGAAGGGACATCGCACACTCGAAGCGCCAAAGTTGCTCGGCATCACCTGGCCTTTCATCGCCGTCGTGGTCC  >  NC_002947/1963643‑1963930
                                                                                                                                                   |                                                                                                                                             
tgGGAGCGGGCTTGCCCGCGAACACGGGCGAAGCCCGTGCCATCCACCGCGGTGTCTTCTTCGCGGGTAAACCCGCTCCTACAATAGCCAGCGCATACGCTGTACAAATTCCCCGATTCCTTGTCATAAAAATCCCCTCCACAGCGAt                                                                                                                                              <  6:110573/148‑1 (MQ=255)
  ggAGCGGGCTTGCCCGCGAACACGGGCGAAGCCCGTGCCATCCACCGCGGTGTCTTCTTCGCGGGTAAACCCGCTCCTACAATAGCCAGCGCATACGCTGTACAAATTCCCCGATTCCTTGTCATAAAAATCCCCTCCACAGCGATgg                                                                                                                                            <  5:73589/148‑1 (MQ=255)
     gCGGGCTTGCCCGCGAACACGGGCGAAGCCCGTGCCATCCACCGCGGTGTCTTCTTCGCGGGTAAACCCGCTCCTACAATAGCCAGCGCATACGCTGTACAAATTCCCCGATTCCTTGTCATAAAAATCCCCTCCACAGCGATGGCaaa                                                                                                                                        >  3:252672/1‑149 (MQ=255)
           ttGCCCGCGAACACGGGCGAAGCCCGTGCCATCCACCGCGGTGTCTTCTTCGCGGGTACACCCGCTCCTACAATAGCCAGCGGATACGCTGTACAGATTCCCCGATTCGTTGTCATAAGGATGCCCTCCGCTGCGATGGGATATGGcgt                                                                                                                                  >  2:320246/1‑147 (MQ=18)
                ggatggcacGGGCTTCGCCCGTGTCATCCACCGCGGTGTCTTCTTCGCGGGTAAACCCGCTCCTACAATAGCCAGCGCATACGCTGTACAAATTCCCCGATTCCTTGTCATAAAAATCCCCTCCACAGCGATGGCAAATGGCCTACGCt                                                                                                                             >  4:303331/7‑149 (MQ=38)
                 gcgAACACGGGCGAAGCCCGTGCCATCCACCGCGGTGTCTTCTTCGCGGGTAAACCCGCTCCTACAATAGCCAGCGCATACGCTGTACAAATTCCCCGATTCCTTGTCATAAAAATCCCCTCCACAGCGATGGCAAATGGCCTACGCtg                                                                                                                            >  6:329542/1‑149 (MQ=255)
                   tggcacGGGCTTCGCCCGTGTCATCCACCGCGGTGTCTTCTTCGCGGGTAAACCCGCTCCTACAATAGCCAGCGCATACGCTGTACAAATTCCCCGATTCCTTGTCATAAAAATCCCCTCCACAGCGATGGCAAATGGCCTACGCtgtg                                                                                                                          <  3:303331/146‑1 (MQ=38)
                   gAACACGGGCGAAGCCCGTGCCATCCACCGCGGTGGCTTCTTCGCGGGTAAACCCGCTCCTACAATAGCCAGCGCATACGCTGTACAAATTCCCCGATTCCTTGTCATAAAAATCCCCTCCACAGCGATGGCAAATGGCCTACGCtgtg                                                                                                                          <  8:201199/149‑1 (MQ=255)
                    aaCACGGGCGAAGCCCGTGCCATCCACCGCGGTGTCTTCTTCGCGGGTAAACCCGCTCCTACAATAGCCAGCGCATACGCTGTACAAATTCCCCGATTCCTTGTCATAAAAATCCCCTCCACAGCGATGGCAAATGGCCTACGCTGTGc                                                                                                                         <  5:106282/149‑1 (MQ=255)
                    aaCACGGGCGAAGCCCGTGCCATCCACCGCGGTGTCTTCTTCGCGGGTAAACCCGCTCCTACAATAGCCAGCGCATACGCTGTACAAATTCCCCGATTCCTTGTCATAAAAATCCCCTCCACAGCGATGGCAAATGGCCTACGCTGTGc                                                                                                                         <  1:214332/149‑1 (MQ=255)
                      cacGGGCGAAGCCCGTGCCATCCACCGCGGTGTCTTCTTCGCGGGTAAACCCGCTCCTACAATAGCCAGCGCATACGCTGTACAAATTCCCCGATTCCTTGTCATAAAAATCCCCTCCACAGCGATGGCAAATGGCCTACGCTGTGCtt                                                                                                                       <  6:254733/149‑1 (MQ=255)
                      cacGGGCGAAGCCCGTGCCATCCACCGAGGTGTCTTATTCGCGGGTAAACCCGCTCCTACAATAGCCAGCGCATACGCTGTACAAATTCCCCGATTCCTTGTCATAAAAATCCCCTCCACAGCGATGGCAAATGGCCTACGCTGTGCtt                                                                                                                       >  2:5239/1‑149 (MQ=255)
                      cacGGGCGAAGCCCGTGCAATCCACCGCGGTGTCTTCTTCGCGGGTAAACCCGCTCCTACAATAGCCAGCGCATACGCTGTACAAATTCCCCGATTCCTTGTCATAAAAATCCCCTCCACAGCGATGGCAAATGGCCTACGCTGTGCtt                                                                                                                       <  1:58805/149‑1 (MQ=255)
                       acGGGCGAAGACCGTGCCATCCACCGCGGTGTCTTCTTAGCGGGTAAACCCGCTCCTACAATAGCCACCGCATACCCTGTACAAATTACCCGAATCATTGTAATAAAAATCCCCACCACATCGATGGCAAAATGCCAACGAAGTGCTTg                                                                                                                      >  6:47778/1‑149 (MQ=21)
                         gggCGAAGCCCGTGCCATCCACCGCGGTGTCTTCTTCGCGGGTAAACCCGCTCCTACAATAGCCAGCGCATACGCTGTACAAATTCCCCGATTCCTTGTCATAAAAATCCCCTCCACAGCGATGGCAAATGGCCTACGCTGTGCTTGTc                                                                                                                    >  6:354720/1‑149 (MQ=255)
                         gggCGAAGCCCGTGCCATCCACCGCGGTGTCTTCTTCGCGGGTAAACCCGCTCCTACAATAGCCAGCGCATACGCTGTACAAATTCCCCGATTCCTTGTCATAAAAATCCCCTCCACAGCGATGGCAAATGGCCTACGCTGTGCTTGTc                                                                                                                    >  3:325194/1‑149 (MQ=255)
                         gggCGAAGCCCGTGCCATCCACCGCGGTGTCTTCTTCGCGGGTAAACCCGCTCCTACAAAAGCCAGCGCATACGCTGTACAAATTCCCCGATTCCTTGTCATAAAAATCCCCTCCACAGCGATGGCAAATGGCCTACGCTGTGCTTGTc                                                                                                                    >  1:298466/1‑149 (MQ=255)
                             gAAGCCCGTGCCATCCACCGCGGTGTCTTCTTCGCGGGTAAACCCGCTCCTACAATAGCCAGCGCATACGCTGTACAAATTCCCCGATTCCTTGTCATAAAAATCCCCTCCACAGCGATGGCAAATGGCCTACGCTGTGCTTGTCTggg                                                                                                                <  7:114742/149‑1 (MQ=255)
                              aaGCCCGTGCCATCCACCGCGGTGTCTTCTTCGCGGGTAAACCCGCTCCTACAATAGCCAGCGCATACGCTGTACAAATTCCCCGATTCCTTGTCATAAAAATCCCCTCCACAGCGATGGCAAATGGCCTACGCTGTGCTTGTCTGGGt                                                                                                               >  8:101624/1‑149 (MQ=255)
                              aaGCCCGTGCCATCCACCGCGGTGTCTTCTTCGCGGGTAAACCCGCTCCTACAATAGCCAGCGCATACGCTGTACAAATTCCCCGATTCCTTGTCATAAAAATCCCCTCCACAGCGATGGCAAATGGCCTACGCTGTGCTTGTCTGGGt                                                                                                               <  3:369110/149‑1 (MQ=255)
                              aaGCCCGTGCCATCCACCGCGGTGTCTTCTTCGCGGGTAAACCCGCTCCTACAATAGCCAGCGCATACGCTGTACAAATTCCCCGATTCCTTGTCATAAAAATCCCCTCCACAGCGATGGCAAATGGCCTACGCTGTGCTTGTCTGGGt                                                                                                               <  2:373639/149‑1 (MQ=255)
                                       ccATCCACCGCGGTGTCTTCTTCGCGGGTAAACCCGCTCCTACAATAGCCAGCGCATACGCTGTACAAATTCCCCGATTCCTTGTCATAAAAATCCCCTCCACAGCGATGGCAAATGGCCTACGCTGTGCTTGTCTGGGTCGATTTAc                                                                                                       <  4:252672/148‑1 (MQ=255)
                                         aTCCACCGCGGTGTCTTCTTCGCGGGTAAACCCGCTTCCTACAATAGCCAGCGCATACGCTGTACAAATTCCCCGATTCCTTGTCATAAAAATCCCCTCCACAGCGATGGCAAATGGCCTACGCTGTGCTTGTCTGGGt                                                                                                               >  1:246993/1‑139 (MQ=255)
                                         aTCCACCGCGGTGTCTTCTTCGCGGGTAAACCCGCTCCTACAATAGCCAGCGCATACGCTGTACAAATTCCCCGATTCCTTGTCATAAAAATCCCCTCCACAGCGATGGCAAATGGCCTACGCTGTGCTTGTCTGGGt                                                                                                               <  2:246993/138‑1 (MQ=255)
                                            cACCGCGGTGTCTTCTTCGCGGGTAAACCCGCTCCTACAATAGCCAGCGCATACGCTGTACAAATTCCCCGATTCCTTGTCATAAAAATCCCCTCCACAGCGATGGCAAATGGCCTACGCTGTGCTTGTCTGGGTCGATTTACGCCtt                                                                                                  <  7:101624/148‑1 (MQ=255)
                                                   gtgtCTTCTTCGCGGGTAAACCCGCTCCTACAATAGCCAGCGCATACGCTGTACAAATTCCCCGATTCCTTGTCATAAAAATCCCCTCCACAGCGATGGCAAATGGCCTACGCTGTGCTTGTCTGGGTCGATTTACGCCTTATCTGctc                                                                                          <  6:91194/149‑1 (MQ=255)
                                                    tgtCTTCTTCGCGGGTAAACCCGCTCCTACAATAGCCAGCGCATACGCTGTACAAATTCCCCGATTCCTTGTCATAAAAATCCCCTCCACAGCGATGGCAAATGGCCTACGCTGTGCTTGTCTGGGTCGATTTACGCCTTATCTGctct                                                                                         >  5:214344/1‑149 (MQ=255)
                                                      tcttctTCGCGGGTAAACCCGCTCCTACAATAGCCAGCGCATACGCTGTACAAATTCCCCGATTCCTTGTCATAAAAATCCCCTCCACAGCGATGGCAAATGGCCTACGCTGTGCTTGTCTGGGTCGATTTACGCCTTATCTGCTCtgt                                                                                       <  4:50969/149‑1 (MQ=255)
                                                        ttcttcGCGGGTAAACCCGCTCCTACAATAGCCAGCGCATACGCTGTACAAATTCCCCGTTTCCTTGTCATAAAAATCCCCTCCACAGCGATGGCAAATGGACTACGCTGTGCTTGTCTGGGTCGATTTACGCCTTATCTGCTCTGTGa                                                                                     >  1:74815/1‑149 (MQ=255)
                                                        ttcttcGCGGGTAAACCCGCTCCTACAATAGCCAGCGCATACGCTGTACAAATTCCCCGATTCCTTGTCATAAAAATCCCCTCCACAGCGATGGCAAATGGCCTACGCTGTGCTTGTCTGGGTCGATTTACGCCTTATCTGctc                                                                                          >  2:142848/1‑144 (MQ=255)
                                                        ttcttcGCGGGTAAACCCGCTCCTACAATAGCCAGCGCATACGCTGTACAAATTCCCCGATTCCTTGTCATAAAAATCCCCTCCACAGCGATGGCAAATGGCCTACGCTGTGCTTGTCTGGGTCGATTTACGCCTTATCTGctc                                                                                          >  5:78296/1‑144 (MQ=255)
                                                        ttcttcGCGGGTAAACCCGCTCCTACAATAGCCAGCGCATACGCTGTACAAATTCCCCGATTCCTTGTCATAAAAATCCCCTCCACAGCGATGGCAAATGGCCTACGCTGTGCTTGTCTGGGTCGATTTACGCCTTATCTGctc                                                                                          <  6:78296/144‑1 (MQ=255)
                                                        ttcttcGCGGGTAAACCCGCTCCTACAATAGCCAGCGCATACGCTGTACAAATTCCCCGATTCCTTGTCATAAAAATCCCCTCCACAGCGATGGCAAATGGCCTACGCTGTGCTTGTCTGGGTCGATTTACGCCTTATCTGctc                                                                                          <  1:142848/144‑1 (MQ=255)
                                                                gggTAAACCCGCTCCTACAATAGCCAGCGCATACGCTGTACAAATTCCCCGATTCCTTGTCATAAAAATCCCCTCCACAGCGATGGCAAATGGCCTACGCTGTGCTTGTCTGGGTCGATTTACGCCTTATCTGCTCTGTGACAGGGACt                                                                             <  5:285138/149‑1 (MQ=255)
                                                                gggTAAACCCGCTCCTACAATAGCCAGCGCATACGCTGTACAAATTCCCCGATTCCTTGTCATAAAAATCCCCTCCACAGCGATGGCAAATGGCCTACGCTGTGCTTGTCTGGGTCGATTTACGCCTTATCTGCTCTGTGACAGGGACt                                                                             >  3:84683/1‑149 (MQ=255)
                                                                  gTAAACCCGCTCCTACAATAGCCAGCGCATACGCTGTACAAATTCCCCGATTCCTTGTCATAAAAATCCCCTCCACAGCGATGGCAAATGGCCTACGCTGTGCTTGTCTGGGTCGATTTACGCCTTATCTGCTCTGTGACAGGGACTcc                                                                           <  6:214344/149‑1 (MQ=255)
                                                                  gTAAACCCGCTCCTACAATAGCCAGCGCATACGCTGTACAAATTCCCCGATTCCTTGTCATAAAAATCCCCTCCACAGCGATGGCAAATGGCCTACGCTGTGCTTGTCTGGGTCGATTTACGCCTTATCTGCTCTGTGACAGGGACTc                                                                            >  1:322725/1‑148 (MQ=255)
                                                                   tAAACCCGCTCCTACAATAGCCAGCGCATACGCTGTACAAATTCCCCGATTCCTTGTCATAAAAATCCCCTCCACAGCGATGGCAAATGGCCTACGCTGTGCTTGTCTGGGTCGATTTACGCCTTATCTGCTCTGTGACAGGGACTCCa                                                                          <  4:84683/149‑1 (MQ=255)
                                                                     aaCCCGCTCCTACAATAGCCAGCGCATACGCTGTACAAATTCCCCGATTCCTTGTCATAAAAATCCCCTCCACAGCGATGGCAAATGGCCTACGCTGTGCTTGTCTGGGTCGATTTACGCCTTATCTGCTCTGTGACAGGGACTCCATg                                                                        <  2:322725/149‑1 (MQ=255)
                                                                       cccGCTCCTACAATAGCCAGCGCATACGCTGTACAAATTCCCCGATTCCTTGTCATAAAAATCCCCTCCACAGCGATGGCAAATGGCCTACGCTGTGCTTGTCTGGGTCGATTTACGc                                                                                                     >  8:94710/1‑118 (MQ=255)
                                                                       cccGCTCCTACAATAGCCAGCGCATACGCTGTACAAATTCCCCGATTCCTTGTCATAAAAATCCCCTCCACAGCGATGGCAAATGGCCTACGCTGTGCTTGTCTGGGTCGATTTACGCCTTATCTGCTCTGTGACAGGGACTCCACGaa                                                                      <  5:354720/149‑1 (MQ=255)
                                                                       cccGCTCCTACAATAGCCAGCGCATACGCTGAACAAATTCCCCGATTCCTTGTCATAAAAATCCCCTCCACAGCGATGGCAAATGGCCTACGCAGTGCTTGTCTGGGTCGATTTACGCCTTATCTGCTCGGTGACAGGGACTCCATg                                                                        >  8:72521/1‑147 (MQ=255)
                                                                          gCTCCTACAATAGCCAGCGCATACGCTGTACAAATTCCCCGATTCCTTGTCATAAAAATCCCCTCCACAGCGTTGGCAAATGGCCTACGCTGTGCTTGTCTGGGTCGATTTACGCCTTATCTGCTCTGTGACAGGGACTCCATg                                                                        <  6:19494/144‑1 (MQ=255)
                                                                          gCTCCTACAATAGCCAGCGCATACGCTGTACAAATTCCCCGATTCCTTGTCATAAAAATCCCCTCCACAGCGTTGGCAAATGGCCTACGCTGTGCTTGTCTGGGTCGATTTACGCCTTATCTGCTCTGTGACAGGGACTCCATg                                                                        >  5:19494/1‑144 (MQ=255)
                                                                                tACAATAGCCAGCGCATACGCTGTACAAATTCCCCGATTCCTTGTCATAAAAATCCCCTCCACAGCGATGGCAAATGGCCTACGCTGTGCTTGTCTGGGTCGATTTACGCCTTATCTGCTCTGTGACAGGGACTCCATGAAGGGACATc                                                              >  3:170961/1‑149 (MQ=255)
                                                                                 aCAATAGCCAGCGCATACGCTGTACAAATTCCCCGATTCCTTGTCATAAAAATCCCCTCCACAGCGATGGCAAATGGCCTACGCTGTGCTTGTCTGGGTCGATTTACGCCTTATCTGCTCTGTGACAGGGACTCCATGAAGGGACATc                                                              >  8:219254/1‑148 (MQ=255)
                                                                                  cAATAGCCAGCGCATACGCTGTACAAATTCCCCGATTCCTTGTCATAAAAATCCCCTCCACAGCGATGGCAAATGGCCTACGCTGTGCTTGTCTGGGTCGATTTACGCCTTATCTGCTCTGTGACAGGGACTCCATGAAGGGACATc                                                              <  4:179421/147‑1 (MQ=255)
                                                                                   aaTAGCCAGCGCATACGCTGTACAAATTCCCCGATTCCTTGTCATAAAAATCCCCTCCACAGCGATGGCAAATGGCCTACGCTGTGCTTGTCTGGGTCGATTTACGCCTTATCTGCTCTGTGACAGGGACTCCATGAAGGGACATCGca                                                           >  3:371295/1‑149 (MQ=255)
                                                                                         cAGCGCATACGCTGTACAAATTCCCCGATTCCTTGTCATAAAAATCCCCTCCACAGCGATGGCAAATGGCCTACGCTGTGCTTGTCTGGGTCGATTTACGCCTTATCTGCTCTGTGACAGGGACTCCATGAAGGGACATCGCACACTc                                                      >  3:83897/1‑148 (MQ=255)
                                                                                         cAGCGCATACGCTGTACAAATTCACCGATTCCTTGTCATAAAAATCACCTCCACAGCGATAGCAAATGGCCTAAGCTGTGATTGTCTGGGTCGATTTACGCCTTATCTGCTATGTGACAGGGACTCAATGAAGGGACAACGCACACTc                                                      >  8:235198/1‑148 (MQ=255)
                                                                                           gcgcACACGCTGTACAAATTCCCCGATTCCTTGTCATAAAAATCCCCTCCACAGCGATGGCAAATGGCCTACGCTGTGCTTGTCTGGGTCGATTTACGCCTTATCTGCTCTGTGACAGGGACTCCATGAAGGGACATCGCACACTCGaa                                                   <  8:332624/149‑1 (MQ=255)
                                                                                             gcATACGCTGTACAAATTCCCCGATTCCTTGTCATAAAAATCCCCTCCACAGCGATGGCAAATGGCCTACGCTGTGCTTGTCTGGGTCGATTTACGCCTTATCTGCTCTGTGACAGGGACTCCATGAAGGGACATCGCACACTCGAAgc                                                 <  2:298466/149‑1 (MQ=255)
                                                                                               aTACGCTGTACAAATTCCCCGTTTCCTTGTCATAAAAATCCCCTCCACAGCGATGGCAAATGGCCTACGCTGTGCTTGTCTGGGTCGATTTACGCCTTATCTGCTCTGTGACAGGGACTCCATGAAGGGACATCGCACACTCGAAgcgc                                               <  2:74815/149‑1 (MQ=255)
                                                                                                tACGCTGTACAAATTCCCCGATTCCTTGTCATAAAAATCCCCTCCACAGCGATGGCAAATAGCCTACGCTGTGCTTGTCTGGGTCGATTTACGCCTTATCTGCTCTGTGACAGGGACTCCATGAAGGGACATCGCACACTCGAAgcgc                                               >  1:109533/1‑148 (MQ=255)
                                                                                                   gCTGTACAAATTCCCCGATTCCTTGTCATAAAAATCCCCTCCACAGCGATGGCAAATGGCCTACGCTGTGCTTGTCTGGGTCGATTTACGCCTTATCTGCTCTGTGACAGGGACTCCATGAAGGGACATCGCACACTCGAAGCGCCaaa                                           >  4:303565/1‑149 (MQ=255)
                                                                                                     tGTACAAATTCCCCGATTCCTTGTCATAAAAATCCCCTCCACAGCGATGGCAAATGGCCTACGCTGTGCTTGTCTGGGTCGATTTACGCCTTATCTGCTCTGTGACAGGGACTCCATGAAGGGACATCGCACACTCGAAGCGCCAAAGt                                         >  7:336856/1‑149 (MQ=255)
                                                                                                        aCAAATTCCCCGATTCCTTGTCATAAAAATCCCCTCCACAGCGATGGCAAATGGCCTACGCTGTGCTTGTCTGGGTCGATTTACGCCTTATCTGCTCTGTGACAGGGACTCCATGAAGGGACATCGCACACTCGAAGCGCCAAAGTTg                                       <  5:47778/148‑1 (MQ=255)
                                                                                                        aCAAATTCCCCGATTCCTTGTCATAAAAATCCCCTCCACAGCGATGGCAAATGGCCTACGCTGTGCTTGTCTGGGTCGATTTACGCCTTATCTGCTCTGTGACAGGGACTCCATGAAGGGACATCGCACACTCGAAGCGCCAAAGTTGc                                      >  2:108147/1‑149 (MQ=255)
                                                                                                          aaaTTCCCCGATTCCTTGTCATAAAAATCCCCTCCACAGCGATGGCAAATGGCCTACGCTGTGCTTGTCTGGGTCGATTTACGCCTTATCTGCTCTGTGACAGGGACTCCATGAAGGGACATCGCACACTCGAAGCGCCAAAGTTGCTc                                    >  3:85526/1‑149 (MQ=255)
                                                                                                           aaTTCCCCGATTCCTTGTCATAAGAATCCCCTCCACAGCGATGGCAAATGGCCTACGCTGTGCTTGTCTGGGTCGATTTACGCCTTATCTGCTCTGTGACAGGGACTCCATGAAGGGACATCGCACACTCGAAGCGCCAAAGTTGCTCg                                   <  3:182944/149‑1 (MQ=255)
                                                                                                           aaTTCCCCGATTCCTTGTCATAAAAATCCCCTCCACAGCGATGGCAAATGGCCTACGCTGTGCTTGTCTGGGTCGATTTACGCCTTATCTGCTCTGTGACAGGGACTCCATGAAGGGACATCGCACACTCGAAGCGCCAAAGTTGCTCg                                   <  3:303565/149‑1 (MQ=255)
                                                                                                           aaTTCCCCGATTCCTTGTCATAAAAATCCCCTCCACAGCGATGGCAAATGGCCTACGCTGTGCTTGTCTGGGTCGATTTACGCCTTATCTGCTCTGTGACAGGGACTCCATGAAGGGACATCGCACACTCGAAGCGCCAAAGTTGCTCg                                   <  7:219254/149‑1 (MQ=255)
                                                                                                                 ccGATTCCTTGTCATAAAAATCCCCTCCACAGCGATGGCAAATGGCCTACGCTGTGCTTGTCTGGGTCGATTTACGCCTTATCTGCTCTGTGACAGGGACTCCATGAAGGGACATCGCACACTCGAAGCGCCAAAGTTGCTCGGCATCa                             >  7:261103/1‑149 (MQ=255)
                                                                                                                      tCCTTGTCATAAAAATCCCCTCCACAGCGATGGCAAATGGCCTACGCTGTGCTTGTCTGGGTCGATTTACGCCTTATCTGCTCTGTGACAGGGACTCCATGAAGGGACATCGCACACTCGAAGCGCCAAAGTTGCTCGGCATCACCTgg                        >  3:222933/1‑149 (MQ=255)
                                                                                                                         ttGTCATAAAAATCCCCTCCACAGCGATGGCAAATGGCCTACGCTGTGCTTGTCTGGGTCGATTTACGCCTTATCTGCTCTGTGACAGGGACTCCATGAAGGGACATCGCACACTCGAAGCGCCAAAGTTGCTCGGAATCCCCTGGCCt                     >  6:88694/1‑149 (MQ=255)
                                                                                                                           gTCATAAAAATCCCCTCCACAGCGATGGCAAATGGCCTACGCTGTGCTTGTCTGGGTCGATTTACGCCTTATCTGCTCTGTGACAGGGACTCCATGAAGGGACATCGCACACTCGAAGCGCCAAAGTTGCTCGGCATCACCTGGCCt                     <  2:109533/147‑1 (MQ=255)
                                                                                                                            tCATAAAAATCCCCTCCACAGCGATGGCAAATGGCCTACGCTGTGCTTGTCTGGGTCGATTTACGCCTTATCTGCTCTGTGACAGGGACTCCATGAAGGGACATCGCACACTCGAAGCGCCAAAGTTGCTCGGCATCACCTGGCCTTTc                  <  4:222933/149‑1 (MQ=255)
                                                                                                                             cATAAAAATCCCCTCCACAGCGATGGCAAATGGCCTACGCTGTGCTTGTCTGGGTCGATTTACGCCTTATCTGCTCTGTGACAGGGACTCCATGAAGGGACATCGCACACTCGAAGCGCCAAAGTTGCTCGGCATCACCTGGCCTTTCa                 >  2:244687/1‑149 (MQ=255)
                                                                                                                              aTAAAAATCCCCTCCACAGCGATGGCAAATGGCCTACGCTGTGCTTGTCTGGGTCGATTTACGCCTTATCTGCTCTGTGACAGGGACTCCATGAAGGGACATCGCACACTCGAAGCGCCAAAGTTGCTCGGCATCACCTGGCCTTTCAt                >  3:166935/1‑149 (MQ=255)
                                                                                                                                         cTCCACAGCGATGGCAAATGGCCTACGCTGTGCTTGTCTGGGTCGATTTACGCCTTATCTGCTCTGTGACAGGGACTCCATGAAGGGACATCGCACACTCGCAGCGCCAAAGTTGCTCGGCATCACCTGGCCTTTCATCGCCGTCGTgg     >  5:112623/1‑149 (MQ=255)
                                                                                                                                         cTCCACAGCGATGGCAAATGGCCTACGCTGTGCTTGTCTGGGTCGATTTACGCCTTATCTGCTCTGTGACAGGGACTCCATGAAGGGACATCGCACACTCGAAGCGCCAAAGTTGCTCGGCATCACCTGGCCTTTCATCGCCGTCGTgg     <  2:231021/149‑1 (MQ=255)
                                                                                                                                         cTCCACAGCGATGGCAAATGGCCTACGCTGTGCTTGTCTGGGTCGATTTACGCCTTATCTGCGCTGTGCCAGGGACTCCATGAAGGGACATCGCAAACTCGAAGCGCCAAAGTTGCTCGGCATCACCTGGCCTTTCATCGCCGTCGTgg     >  5:117774/1‑149 (MQ=255)
                                                                                                                                         cTCCACAGCGATGGCAAATGGCCTACGCTGTGCTTGGCTGGGTCGATTTACGCCTTATCTGCTCTGTGACAGGGACTCCATGAAGGGACATCGCACACTCGAAGCGCCAAAGTTGCTCGGCATCACCTGGCCTTTCATCGCCGTCGTgg     <  8:336856/149‑1 (MQ=255)
                                                                                                                                         aagcacaGCGTAGGCCATTTGCCATCGCTGTGCTTGTCTGGGTCGATTTACGCCTTATCTGCTCTGTGACAGGGACTCCATGAAGGGACATCGCACACTCGAAGCGCCAAAGTTGCTCGGCATCACCTGGCCTTTCATCGCCGTCGTgg     >  4:122000/4‑149 (MQ=255)
                                                                                                                                         aagcacaGCGTAGGCCATTTGCCATCGCTGTGCTTGTCTGGGTCGATTTACGCCTTATCTGCTCTGTGACAGGGACTCCATGAAGGGACATCGCACACTCGAAGCGCCAAAGTTGCTCGGCATCACCTGGCCTTTCATCGCCGTCGTgg     >  1:332231/4‑149 (MQ=255)
                                                                                                                                            cacaGCGTAGGCCATTTGCCATCGCTGTGCTTGTCTGGGTCGATTTACGCCTTATCTGCTCTGTGACAGGGACTCCATGAAGGGACATCGCACACTCGAAGCGCCAAAGTTGCTCGGCATCACCTGGc                       >  2:109910/1‑128 (MQ=255)
                                                                                                                                            cacaGCGTAGGCCATTTGCCATCGCTGTGCTTGTCTGGGTCGATTTACGCCTTATCTGCTCTGTGACAGGGACTCCATGAAGGGACATCGCACACTCGAAGCGCCAAAGTTGCTCGGCATCACCTGGc                       <  1:109910/128‑1 (MQ=255)
                                                                                                                                            cacaGCGATGGCAAATGGCCTACGCTGTGCTTGTCTGGGTCGATTTACGCCTTATATGCTCTGTGACAGGGACTCCATGAAGGGACATCGAACACTCGAAGCGCCAAAGTTGCTCGGCATCACCTGGCCTTTCATCGCCGTCGTGGTcc  >  1:153216/1‑149 (MQ=255)
                                                                                                                                                   |                                                                                                                                             
TGGGAGCGGGCTTGCCCGCGAACACGGGCGAAGCCCGTGCCATCCACCGCGGTGTCTTCTTCGCGGGTAAACCCGCTCCTACAATAGCCAGCGCATACGCTGTACAAATTCCCCGATTCCTTGTCATAAAAATCCCCTCCACAGCGATGGCAAATGGCCTACGCTGTGCTTGTCTGGGTCGATTTACGCCTTATCTGCTCTGTGACAGGGACTCCATGAAGGGACATCGCACACTCGAAGCGCCAAAGTTGCTCGGCATCACCTGGCCTTTCATCGCCGTCGTGGTCC  >  NC_002947/1963643‑1963930

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: