Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,606,450 A→T 6.0% intergenic (‑158/+53) PP_1407 ← / ← PP_t64 sulfate transporter/tRNA‑Gly

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029471,606,4500AT 6.0% 320.7 / 6.8 97intergenic (‑158/+53)PP_1407/PP_t64sulfate transporter/tRNA‑Gly
Reads supporting (aligned to +/- strand):  ref base A (43/48);  new base T (4/2);  total (47/50)
Fisher's exact test for biased strand distribution p-value = 4.26e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.99e-01

CTTTCGAATAGTTTCACGAATCCAGTGAGCATATCCGCCACAAGCCGTTTCGCTTGCGCACGCAAAGGCCCGATGCGCATTGAAGCGGATCGGTGGGTAATTGGCCATGGGCCAGGTCAAGGCAGGATGAAAAACGAAAAAAGGAGCCAAGGCTCCTTTTGACTGCGATGACGTTGCTGGAACGCCATGCGGGACTGCGATGTGGAGCGGGTAGAGGGAATCGAACCCTCAACTAAAGCTTGGGAAGCTTTCGTTTTGCCACTAAACTATACCCGCATCTGAGCTCGACTTTTTACCA  >  NC_002947/1606301‑1606598
                                                                                                                                                     |                                                                                                                                                    
cTTTCGAATAGTTTCACGAATCCAGTGAGCATATCCGCCACAAGCCGTTTCGCTTGCGCACGCAAAGGCCCGATGCGCATTGAAGCGGATCGGTGGGTAATTGGCCATGGGCCAGGTCAAGGCAGGATGAAAAACGAAAAAAGGAGCCaa                                                                                                                                                      >  1:1157645/1‑150 (MQ=255)
 tttCGAATAGTTTCACGAATCCAGTGAGCATATCCGCCACAAGCCGTTTCGCTTGCGCACGCAAAGGCCCGATGCGCATTGAAGCGGATCGGTGGGTAATTGGCCATGGGCCAGGTCAAGGCAGGATGAAAAACGAAAAAAGGAGCCAAg                                                                                                                                                     >  1:2033627/1‑150 (MQ=255)
 tttCGAATAGTTTCACGAATCCAGTGAGCATATCCGCCACAAGCCGTTTCGCTTGCGCACGCAAAGGCCCGATGCGCATTGAAGCGGATCGGTGGGTAATTGGCCATGGGCCAGGTCAAGGCAGGATGAAAAACGAAAAAAGGAGCCAAg                                                                                                                                                     <  1:1310374/150‑1 (MQ=255)
  ttCGAATAGTTTCACGAATCCAGTGAGCATATCCGCCACAAGCCGCTTCGCTTGCGCACGCAAAGGCCCGAGGCGCATTGAAGCGGATCGGTGGGTAATTGGCCATGGGCCAGGTCAAGGCAGGATGAAAAACGAAAAAAGGAGCCAAgg                                                                                                                                                    <  1:264740/150‑1 (MQ=255)
     gAATAGTTTCACGAATCCAGTGAGCATATCCGCCACAAGCCGTTTCGCTTGCGCACGCAAAGGCCCGATGCGCATTGAAGCGGATCGGTGGGTAATTGGCCATGGGCCAGGTCAAGGCAGGATGAAAAACGAAAAAAGGAGCCAAGGCTc                                                                                                                                                 >  2:2075730/1‑150 (MQ=255)
       aTAGTTTCACGAATCCAGTGAGCATATCCGCCACAAGCCGTTTCGCTTGCGCACGCAAAGGCCCGATGCGCATTGAAGCGGATCGGTGGGTAATTGGCCATGGGCCAGGTCAAGGCAGGATGAAAAACGAAAAAAGGAGCCAAGGCTCCt                                                                                                                                               >  2:1525888/1‑150 (MQ=255)
         aGTTTCACGAATCCAGTGAGCATATCCGCCACAAGCCGTTTCGCTTGCGCACGCAAAGGCCCGATGCGCATTGAAGCGGATCGGTGGGTAATTGGCCATGGGCCAGGTCAAGGCAGGATGAAAAACGAAAAAAGGAGCCAAGGCTCCttt                                                                                                                                             <  1:1824081/150‑1 (MQ=255)
          gTTTCACGAATCCAGTGAGCATATCCGCCACAAGCCGTTTCGCTTGCGCACGCAAAGGCCCGATGCGCATTGAAGCGGATCGGTGGGTAATTGGCCATGGGCCAGGTCAAGGCAGGATGAAAAACGAAAAAAGGAGCCTTGGCTCCtttt                                                                                                                                            <  1:2036350/150‑1 (MQ=255)
          gTTTCACGAATCCAGTGAGCATATCCGCCACAAGCCGTTTCGCTTGCGCACGCAAAGGCCCGATGCGCATTGAAGCGGATCGGTGGGTAATTGGCCATGGGCCAGGTCAAGGCAGGATGAAAAACGAAAAAAGGAGCCAAGGCTCCtttt                                                                                                                                            <  1:1525888/150‑1 (MQ=255)
           tttCACGAATCCAGTGAGCATATCCGCCACAAGCCGTTTCGCTTGCGCACGCAAAGGCCCGATGCGCATTGAAGCGGATCGGTGGGTAATTGGCCATGGGCCAGGTCAAGGCAGGATGAAAAACGAAAAAAGGAGCCAAGGCTCCTTTTg                                                                                                                                           >  2:150206/1‑150 (MQ=255)
               aCGAATCCAGTGAGCATATCCGCCACAAGCCGTTTCGCTTGCGCACGCAAAGGCCCGATGCGCATTGAAGCGGATCGGTGGGTAATTGGCCATGGGCCAGGTCAAGGCAGGATGAAAAACGAAAAAAGGAGCCAAGGCTCCTTTTGACTg                                                                                                                                       >  2:1878039/1‑150 (MQ=255)
                cGAATCCAGTGAGCATATCCGCCACAAGCCGTTTCGCTTGCGCACGCAAAGGCCCGATGCGCATTGAAGCGGATCGGTGGGTAATTGGCCATGGGCCAGGTCAAGGCAGGATGAAAAACGAAAAAAGGAGCCAAGGCTCCTTTTGACTGc                                                                                                                                      <  2:1986147/150‑1 (MQ=255)
                   aTCCAGTGAGCATATCCGCCACAAGCCGTTTCGCTTGCGCACGCAAAGGCCCGATGCGCATTGAAGCGGATCGGTGGGTAATTGGCCATGGGCCAGGTCAAGGCAGGATGAAAAACGAAAAAAGGAGCCAAGGCTCCTTTTGACTGCGAt                                                                                                                                   >  1:2198653/1‑150 (MQ=255)
                      cAGTGAGCATATCCGCCACAAGCCGTTTCGCTTGCGCACGCAAAGGCCCGATGCGCATTGAAGCGGATCGGTGGGTAATTGGCCATGGGCCAGGTCAAGGCAGGATGAAAAACGAAAAAAGGAGCCAAGGCTCCTTTTGACTGCGATGAc                                                                                                                                >  1:1651362/1‑150 (MQ=255)
                       aGTGAGCATATCCGCCACAAGCCGTTTCGCTTGCGCACGCAAAGGCCCGATGCGCATTGAAGCGGATCGGTGGGTAATTGGCCATAGGCCAGGTCAAGGCAGGATGAAAAACGAAAAAAGGAGCCAAGGCTCCTTTTGACTGCGATGACg                                                                                                                               >  2:104302/1‑150 (MQ=255)
                            gCATATCCGCCACAAGCCGTTTCGCTTGCGCACGCAAAGGCCCGATGCGCATTGAAGCGGATCGGTGGGTAATTGGCCATGGGCCAGGTCAAGGCAGGATGAAAAACGAAAAAAGGAGCCAAGGCTCCTTTTGACTGCGATGACGTTGCt                                                                                                                          <  1:1345027/150‑1 (MQ=255)
                            gCATATCCGCCACAAGCCGTTTCGCTTGCGCACGCAAAGGCCCGATGCGCATTGAAGCGGATCGGTGGGTAATTGGCCATGGGCCAGGTCAAGGCAGGATGAAAAACGAAAAAAGGAGCCAAGGCTCCTTTTGACTGCGATGACGTTGCt                                                                                                                          <  2:2198653/150‑1 (MQ=255)
                                    gCCACAAGCCGTTTCGCTTGCGCACGCAAAGGCCCGATGCGCATTGAAGCGGATCGGTGGGTCATTGACCATGGGCCAGGTCAAGGCAGGATGAAAAACGAAAAAAGGAGCCAAGGCTCCTTTTGACTGCGATGACGTTGCTGGAACGcc                                                                                                                  >  2:2026115/1‑150 (MQ=255)
                                          aGCCGTTTCGCTTGCGCACGCAAAGGCCCGATGCGCATTGAAGCGGATCGGTGGGTAATTGGCCATGGGCCAGGTCAAGGCAGGATGAAAAACGAAAAAAGGAGCCAAGGCTCCTTTTGACTGCGATGACGTTGCTGGAACGCCATGCgg                                                                                                            >  1:777099/1‑150 (MQ=255)
                                           gCCGTTTCGCTTGCGCACGCAAAGGCCCGATGCGCATTGAAGCGGATCGGTGGGTAATTGGCCATGGGCCAGGTCAAGGCAGGATGAAAAACGAAAAAAGGAGCCAAGGCTCCTTTTGACTGCGATGACGTTGCTGGAACGCCATGCggg                                                                                                           >  1:1848554/1‑150 (MQ=255)
                                                     ttGCGCACGCAAAGGCCCGATGCGCATTGAAGCGGATCGGTGGGTAATTGGCCATGGGCCAGGTCAAGGCAGGATGAAAAACGAAAAAAGGAGCCAAGGCTCCTTTTGACTGCGATGACGTTGCTGGAACGCCATGCGGGACTGCGAtgt                                                                                                 <  2:1157645/150‑1 (MQ=255)
                                                     ttGCGCACGCAAAGGCCCGATGCGCATTGAAGCGGATCGGTGGGTAATTGGCCATGGGCCAGGTCAAGGCAGGATGAAAAACGAAAAAAGGAGCCAAGGCTCCTTTTGACTGCGATGACGTTGCTGGAACGCCATGCGGGACTGCGAtgt                                                                                                 >  2:828598/1‑150 (MQ=255)
                                                     ttGCGCACGCAAAGGCCCGATGCGCATTGAAGCGGATCGGTGGGTAATTGGCCATGGGCCAGGTCAAGGCAGGATGAAAAACGAAAAAAGGAGCCAAGGCTCCTTTTGACTGCGATGACGTTGCTGGAACGCCATGCGGGACTGCGAtgt                                                                                                 <  2:2271429/150‑1 (MQ=255)
                                                     ttGCGCACGCAAAGGCCCGATGCGCATTGAAGCGGATCGGTGGGGAATTGGCCATGGGCCAGGTCAAGGCAGGATGAAAAACGAAAAAAGGAGCCAAGGCTCCTTTTGACTGCGATGACGTTGCTGGAACGCCATGGGGGACTGCGAtgt                                                                                                 >  1:1927748/1‑150 (MQ=255)
                                                      tGCGCACGCAAAGGCCCGATGCGCATTGAAGCGGATCGGTGGGTAATTGGCCATGGGCCAGGTCAAGGCAGGATGAAAAACGAAAAAAGGAGCCAAGGCTCCTTTTGACTGCGATGACGTTGCTGGAACGCCATGCGGGACTGCGAtgtg                                                                                                <  1:1500407/150‑1 (MQ=255)
                                                             gcaAAGGCCCGATGCGCATTGAAGCGGATCGGTGGGTAATTGGCCATGGGCCAGGTCAAGGCAGGATGAAAAACGAAAAAAGGAGCCAAGGCTCCTTTTGACTGCGATGACGTTGCTGGAACGCCATGCGGGACTGCGATGTGGAGCggg                                                                                         <  2:2398361/150‑1 (MQ=255)
                                                             gcaAAGGCCCGATGCGCATTGAAGCGGATCGGTGGGTAATTGGCCATGGGCCAGGTCAAGGCAGGATGAAAAACGAAAAAAGGAGCCAAGGCTCCTTTTGACTGCGATGACGTTGCTGGAACGCCATGCGGGACTGCGATGTGGAGCggg                                                                                         >  2:1225002/1‑150 (MQ=255)
                                                                  ggCCCGATGCGCATTGAAGCGGATCGGTGGGTAATTGGCCATGGGCCAGGTCAAGGCAGGATGAAAAACGAAAAAAGGAGCCAAGGCTCCTTTTGACTGCGATGACGTTGCTGGAACGCCATGCGGGACTGCGATGTGGAGCGGGTagag                                                                                    <  1:1878039/150‑1 (MQ=255)
                                                                      cGATGCGCATTGAAGCGGATCGGTGGGTAATTGGCCATGGGCCAGGTCAAGGCAGGATGAAAAACGAAAAAAGGAGCCAAGGCTCCTTTTGACTGCGATGACGTTGCTGGAACGCCATGCGGGACTGCGATGTGGAGCGGGTAGAGGGaa                                                                                <  1:388387/150‑1 (MQ=255)
                                                                      cGATGCGCATTGAAGCGGATCGGTGGGTAATTGGCCATGGGCCAGGTCAAGGCAGGATGAAAAACGAAAAAAGGAGCCAAGGCTCCTTTTGACTGCGATGACGTTGCTGGAACGCCATGCGGGACTGCGATGTGGAGCGGGTAGAGGGaa                                                                                >  2:1017788/1‑150 (MQ=255)
                                                                       gATGCGCATTGAAGCGGATCGGTGGGTAATTGGCCATGGGCCAGGTCAAGGCAGGATGAAAAACGAAAAAAGGAGCCAAGGCTCCTTTTGACTGCGATGACGTTGCTGGAACGCCATGCGGGACTGCGATGTGGAGCGGGTAGAGGGAAt                                                                               >  2:1145534/1‑150 (MQ=255)
                                                                       gATGCGCATTGAAGCGGATCGGTGGGTAATTGGCCATGGGCCAGGTCAAGGCAGGATGAAAAACGAAAAAAGGAGCCAAGGCTCCTTTTGACTGCGATGACGTTGCTGGAACGCCATGCGGGACTGCGATGTGGAGCGGGTAGAGGGAAt                                                                               <  1:1145534/150‑1 (MQ=255)
                                                                        aTGCGCATTGAAGCGGATCGGTGGGTAATTGGCCATGGGCCAGGTCAAGGCAGGATGAAAAACGAAAAAAGGAGCCAAGGCTCCTTTTGACTGCGATGACGTTGCTGGAACGCCATGCGGGACTGCGATGTGGAGCGGGTAGAGGGAATc                                                                              <  2:1744703/150‑1 (MQ=255)
                                                                           cgcATTGAAGCGGATCGGTGGGTAATTGGCCATGGGCCAGGTCAAGGCAGGATGAAAAACGAAAAAAGGAGCCAAGGCTCCTTTTGACTGCGATGACGTTGCTGGAACGCCATGCGGGACTGCGATGTGGAGCGGGTAGAGGGAATCGaa                                                                           <  2:351016/150‑1 (MQ=255)
                                                                            gcATTGAAGCGGATCGTTGGGTAATTGGCCATGGGCCAGGTCAAGGCAGGATGAAAAACGAAAAAAGGAGCCAAGGCTCCTTTTGACTGCGATGACGTTGCTGGAACGCCATGCGGGACTGCGATGTGGAGCGGGTAGAGGGAATCGAAc                                                                          <  2:551392/150‑1 (MQ=255)
                                                                            gcATTGAAGCGGATCGGTGGGTAATTGGCCATGGGCCAGGTCAAGGCAGGATGAAAAACGAAAAAAGGAGCCTTGGCTCCtttttt                                                                                                                                          <  2:1459616/86‑3 (MQ=39)
                                                                            gcATTGAAGCGGATCGGTGGGTAATTGGCCATGGGCCAGGTCAAGGCAGGATGAAAAACGAAAAAAGGAGCCTTGGCTCCtttttt                                                                                                                                          >  1:1459616/1‑84 (MQ=39)
                                                                                  aaGCGGATCGGTGGGTAATTGGCCATGGGCCAGGTCAAGGCAGGATGAAAAACGAAAAAAGGAGCCAAGGCTCCTTTTGACTGCGATGACGTTGCTGGAACGCCATGCGGGACTGCGATGTGGAGCGGGTAGAGGGAATc                                                                              >  2:456774/1‑140 (MQ=255)
                                                                                  aaGCGGATCGGTGGGTAATTGGCCATGGGCCAGGTCAAGGCAGGATGAAAAACGAAAAAAGGAGCCAAGGCTCCTTTTGACTGCGATGACGTTGCTGGAACGCCATGCGGGACTGCGATGTGGAGCGGGTAGAGGGAATc                                                                              <  1:456774/140‑1 (MQ=255)
                                                                                  aaGCGGATCGGTGGGTAATTGGCCATGGGCCAGGTCAAGGCAGGATGAAAAACGAAAAAAGGAGCCAAGGCTCCTTTTGACTGCGATGACGTTGCTGGAACGCCATGCGGGACTGCGATGTGGAGCGGGTAGAGGGAATCGAACCCTCaa                                                                    <  2:1848554/150‑1 (MQ=255)
                                                                                    gCGGATCGGTGGGTAATTGGCCATGGGCCAGGTCAAGGCAGGATGAAAAACGAAAAAAGGAGCCAAGGCTCCTTTTGACTGCGATGACGTTGCTGGAACGCCATGCGGGACTGCTATGTGGAGCGGGTAGAGGGAATCGAACCCTCAACt                                                                  >  1:1001883/1‑150 (MQ=255)
                                                                                      ggATCGGTGGGTAATTGGCCATGGGCCAGGTCAAGGCAGGATGAAAAACGAAAAAAGGAGCCAAGGCTCCTTTTGACTGCGATGACGTTGCTGGAACGCCATGCGGGACTGCGATGTGGAGCGGGTaga                                                                                     >  2:595749/1‑129 (MQ=255)
                                                                                      ggATCGGTGGGTAATTGGCCATGGGCCAGGTCAAGGCAGGATGAAAAACGAAAAAAGGAGCCAAGGCTCCTTTTGACTGCGATGACGTTGCTGGAACGCCATGCGGGACTGCGATGTGGAGCGGGTaga                                                                                     <  1:595749/129‑1 (MQ=255)
                                                                                       gATCGGTGGGTAATTGGCCATGGGCCAGGTCAAGGCAGGATGAAAAACGAAAAAAGGAGCCAAGGCTCCTTTTGACTGCGATGACGTTGCTGGAACGCCATGCGGGACTGCGATGTGGAGCGGGTAGAGGGAATCGAACCCTCAACTaaa                                                               <  1:2075730/150‑1 (MQ=255)
                                                                                        aTCGGTGGGTAATTGGCCATGGGCCAGGTCAAGGCAGGATGAAAAACGAAAAAAGGAGCCAAGGCTCCTTTTGACTGCGATGACGTTGCTGGAACGCCATGCGGGACTGCGATGTGGAGCGGGTAGAGGGAATCGAACCCTCAACTAAAg                                                              >  2:1615204/1‑150 (MQ=255)
                                                                                        aTCGGTGGGTAATTGGCCATGGGCCAGGTCAAGGCAGGATGAAAAACGAAAAAAGGAGCCAAGGCTCCTTTTGACTGCGATGACGTTGCTGGAACGCCATGCGGGACTGCGATGTGGAGCGGGTAGAGGGAATCGAACCCTCAACTAAAg                                                              <  2:548990/150‑1 (MQ=255)
                                                                                        aTCGGTGGGTAATTGGCCATGGGCCAGGTCAAGGCAGGATGAAAAACGAAAAAAGGAGCCAAGGCTCCTTTTGACTGCGATGACGTTGCTGGAACGCCATGCGGGACTGCGATGTGGAGCGGGTAGAGGGAATCGAACCCTCAACTAAAg                                                              <  2:777099/150‑1 (MQ=255)
                                                                                           ggTGGGTAATTGGCCATGGGCCAGGTCAAGGCAGGATGAAAAACGAAAAAAGGAGCCAAGGCTCCTTTTGACTGCGATGACGTTGCTGGAACGCCATGCGGGACTGCGATGTGGAGCGGGTAGAGGGAATCGAACCCTCAACTAAAGCtt                                                           <  1:2005532/150‑1 (MQ=255)
                                                                                            gTGGGTAATTGGCCATGGGCCAGGTCAAGGCAGGATGAAAAACGAAAAAAGGAGCCAAGGCTCCTTTTGACTGCGATg                                                                                                                                  >  1:720380/1‑78 (MQ=255)
                                                                                            gTGGGTAATTGGCCATGGGCCAGGTCAAGGCAGGATGAAAAACGAAAAAAGGAGCCAAGGCTCCTTTTGACTGCGATg                                                                                                                                  <  2:720380/78‑1 (MQ=255)
                                                                                            gTGGGTAATTGGCCATGGGCCAGGTCAAGGCAGGATGAAAAACGAAAAAAGGAGCCAAGGCTCCTTTTGACTGCGATGACGTTGCTGGAACGCCATGCGGGACTGCGATGTGGAGCGGGTAGAGGGAATCGAACCCTCAACTAAAGCTTg                                                          >  2:1747187/1‑150 (MQ=255)
                                                                                             tGGGTAATTGGCCATGGGCCAGGTCAAGGCAGGATGAAAAACGAAAAAAGGAGCCAAGGCTCCTTTTGACTGCGATGACGt                                                                                                                              <  2:714567/81‑1 (MQ=255)
                                                                                             tGGGTAATTGGCCATGGGCCAGGTCAAGGCAGGATGAAAAACGAAAAAAGGAGCCAAGGCTCCTTTTGACTGCGATGACGt                                                                                                                              >  1:714567/1‑81 (MQ=255)
                                                                                             tGGGTAATTGGCCATGGGCCAGGTCAAGGCAGGATGAAAAACGAAAAAAGGAGCCAAGGCTCCTTTTGACTGCGATGACGTTGCTGGAACGCCATGCGGGACTGCGATGTGGAGCGGGTAGAGGGAATCGAACCCTCAACTAAAGCTTgg                                                         <  1:852120/150‑1 (MQ=255)
                                                                                             tGGGTAATTGGCCATGGGCCAGGTCAAGGCAGGATGAAAAACGAAAAAAGGAGCCAAGGCTCCTTTTGACTGCGATGACGTTGCTGGAACGCCATGCGGGACTGCGATGTGGAGCGGGTAGAGGGAATCGAACCCTCAACTAAAGCTTgg                                                         <  1:1615204/150‑1 (MQ=255)
                                                                                               ggTAATTGGCCATGGGCCAGGTCAAGGCAGGATGAAAAACGAAAAAAGGAGCCAAGGCTCCTTTTGACTGCGATGACGTTGCTGGAACGCCATGCGGGACTGCGATGTGGAGCGGGTAGAGGGAATCGAACCCTCAACTaaa                                                               <  2:9567/142‑1 (MQ=255)
                                                                                               ggTAATTGGCCATGGGCCAGGTCAAGGCAGGATGAAAAACGAAAAAAGGAGCCAAGGCTCCTTTTGACTGCGATGACGTTGCTGGAACGCCATGCGGGACTGCGATGTGGAGCGGGTAGAGGGAATCGAACCCTCAACTaaa                                                               >  1:9567/1‑142 (MQ=255)
                                                                                                 tAATTGGCCATGGGCCAGGTCAAGGCAGGATGAAAAACGAAAAAAGGAGCCAAGGCTCCTTTTGACTGCGATGACGTTGCTGGAACGCCATGCGGGACTGCGATGTGGAGCGGGTAGAGGGAATCGAACCCTCAACTAAAGCTTGGGAAg                                                     >  2:581805/1‑150 (MQ=255)
                                                                                                  aaTTGGCCATGGGCCAGGTCAAGGCAGGATGAAAAACGAAAAAAGGAGCCAAGGCTCCTTTTGACTGCGATGACGTTGCTGGAACGCCATGCGGGACTGCGATGTGGAGCGGGTAGAGGGAATCGAACCCTCAACTAAAGCTTGGGaa                                                      >  1:90803/1‑148 (MQ=255)
                                                                                                  aaTTGGCCATGGGCCAGGTCAAGGCAGGATGAAAAACGAAAAAAGGAGCCAAGGCTCCTTTTGACTGCGATGACGTTGCTGGAACGCCATGCGGGACTGCGATGTGGAGCGGGTAGAGGGAATCGAACCCTCAACTAAAGCTTGGGaa                                                      <  2:90803/148‑1 (MQ=255)
                                                                                                     tGGCCATGGGCCAGGTCAAGGCAGGATGAAAAACGAAAAAAGGAGCCAAGGCTCCTTTTGACTGCGATGACGTTGCTGGAACGCCATGCGGGACTGCTATGTGGAGCGGGTAGAGGGAATCGAACCCTCAACTAAAGCTTGGGAAGCttt                                                 <  2:1001883/150‑1 (MQ=255)
                                                                                                      ggCCATGGGCCAGGGCAAGGCAGGATGAAAAACGAAAAAAGGAGCCAAGGCTCCTTTTGACTGCGATGACGTTGCTGGAACGCCATGCGGGACTGCGATGTGGAGCGGGTAGAGGGAATCGAACCCTCAACTAAAGCTTGGGAAGCTTTc                                                <  1:1711084/150‑1 (MQ=255)
                                                                                                       gCCATGGGCCAGGTCAAGGCAGGATGAAAAACGAAAAAAGGAGCCAAGGCTCCTTTTGACTGCGATGACGTTGCTGGAACGCCATGCGGGACTGCGATGTGGAGCGGGTAGAGGGAATCGAACCCTCAACTAAAGCTTGGGaa                                                      <  2:2226204/143‑1 (MQ=255)
                                                                                                       gCCATGGGCCAGGTCAAGGCAGGATGAAAAACGAAAAAAGGAGCCAAGGCTCCTTTTGACTGCGATGACGTTGCTGGAACGCCATGCGGGACTGCGATGTGGAGCGGGTAGAGGGAATCGAACCCTCAACTAAAGCTTGGGaa                                                      >  1:2226204/1‑143 (MQ=255)
                                                                                                       gCCATGGGCCAGGTCAAGGCAGGATGAAAAACGAAAAAAGGAGCCAAGGCTCCTTTTGACTGCGATGACGTTGCTGGAACGCCATGCGGGACTGCGATGTGGAGCGGGTAGAGGGAATCGAACCCTCAACTAAAGCTTGGGAAGCTTTCg                                               >  2:2343267/1‑150 (MQ=255)
                                                                                                         cATGGGCCAGGTCAAGGCAGGATGAAAAACGAAAAAAGGAGCCAAGGCTc                                                                                                                                                 <  1:757316/50‑1 (MQ=255)
                                                                                                         cATGGGCCAGGTCAAGGCAGGATGAAAAACGAAAAAAGGAGCCAAGGCTc                                                                                                                                                 >  2:757316/1‑50 (MQ=255)
                                                                                                         cATGGGCCAGGTCAAGGCAGGATGAAAAACGAAAAAAGGAGCCAAGGCTCCTTTTGACTGCGATGACGTTGCTGGAACGCCATGCGGGACTGCGATGTGGAGCggg                                                                                         >  2:812550/1‑106 (MQ=255)
                                                                                                         cATGGGCCAGGTCAAGGCAGGATGAAAAACGAAAAAAGGAGCCAAGGCTCCTTTTGACTGCGATGACGTTGCTGGAACGCCATGCGGGACTGCGATGTGGAGCggg                                                                                         <  1:812550/106‑1 (MQ=255)
                                                                                                         cATGGGCCAGGTCAAGGCAGGATGAAAAACGAAAAAAGGAGCCAAGGCTCCTTTTGACTGCGATGACGTTGCTGGAACGCCATGCGGGACTGCGATGTGGAGCg                                                                                           <  2:175052/104‑1 (MQ=255)
                                                                                                         cATGGGCCAGGTCAAGGCAGGATGAAAAACGAAAAAAGGAGCCAAGGCTCCTTTTGACTGCGATGACGTTGCTGGAACGCCATGCGGGACTGCGATGTGGAGCg                                                                                           >  1:175052/1‑104 (MQ=255)
                                                                                                         cATGGGCCAGGTCAAGGCAGGATGAAAAACGAAAAAAGGAGCCAAGGCTCCTTTTGACTGCGATGACGTTGCTGGAACGCCATGCGGGACTGCGATGTGGAGCGGGTAGAGGGAATCGAACCCTCAACTAAAGCTTGGGAAGCTTTCGtt                                             >  1:1707917/1‑150 (MQ=255)
                                                                                                                    tCAAGGCAGGATGAAAAACGAAAAAAGGAGCCAAGGCTCCTTTTGACTGCGATGACGTTGCTGGAACGCCATGCGGGACTGCGATGTGGAGCGGGTAGAGGGAATCGAACCCTCAACTAAAGCTTGGGAAGCTTTCGTTTTGCCACTaaa                                  <  1:30376/150‑1 (MQ=255)
                                                                                                                      aaGGCAGGATGAAAAACGAAAAAAGGAGCCAAGGCTCCTTTTGACTGCGATGACGTTGCTGGAACGCCATGCGGGACTGCGATGTGGAGCGGGTAGAGGGAATCGAACCCTCAACTAAAGCTTGGGAAGCTTTCGTTTTGCCACTaaa                                  <  2:136137/148‑1 (MQ=255)
                                                                                                                      aaGGCAGGATGAAAAACGAAAAAAGGAGCCAAGGCTCCTTTTGACTGCGATGACGTTGCTGGAACGCCATGCGGGACTGCGATGTGGAGCGGGTAGAGGGAATCGAACCCTCAACTAAAGCTTGGGAAGCTTTCGTTTTGCCACTaaa                                  >  1:136137/1‑148 (MQ=255)
                                                                                                                       aGGCAGGATGAAAAACGAAAAAAGGAGCCAAGGCTCCTTTTGACTGCGATGACGTTGCTGGAACGCCATGCGGGACTGCGATGTGGAGCGGGTAGAGGGAATCGAACCCTCAACTAAAGCTTGGGAAGCTTTCGTTTTGCCACTAAACta                               >  1:1487847/1‑150 (MQ=255)
                                                                                                                         gCAGGATGAAAAACGAAAAAAGGAGCCAAGGCTCCTTTTGACTGCGATGACGTTGCTGGAACGCCATGCGGGACTGCGATGTGGAGCGGGTAGAGGGAATCGAACCCTCAACTAAAGCTTGGGAAGCTTTCGTTTTGCCACTAAACtata                             >  1:406685/1‑150 (MQ=255)
                                                                                                                         gCAGGATGAAAAACGAAAAAAGGAGCCAAGGCTCCTTTTGACTGCGATGACGTTGCTGGAACGCCATGCGGGACTGCGATGTGGAGCGGGTAGAGGGAATCGAACCCTCAACTAAAGCTTGGGAAGCTTTCGTTTTGCCACTAAACtata                             >  1:1540575/1‑150 (MQ=255)
                                                                                                                         gCAGGATGAAAAACGAAAAAAGGAGCCAAGGCTCCTTTTGACTGCGATGACGTTGCTGGAACGCCATGCGGGACTGCGATGTGGAGCGGGTAGAGGGAATCGAACCCTCAACTAAAGCTTGGGAAGCTTTCGTTTTGCCACTAAACtata                             >  1:2123316/1‑150 (MQ=255)
                                                                                                                          cAGGATGAAAAACGAAAAAAGGAGCCAAGGCTCCTTTTGACTGCGATGACGTTGCTGGAACGCCATGCGGGACTGCGATGTGGAGCGGGTAGAGGGAATCGAACCCTCAACTAAAGCTTGGGAAGCTTTCGTTTTGCCACTAAACTATAc                            >  1:769587/1‑150 (MQ=255)
                                                                                                                           aGGATGAAAAACGAAAAAAGGAGCCAAGGCTCCTTTTGACTGCGATGACGTTGCTGGAACGCCATGCGGGACTGCGATGTGGAGCGGGTAGATGGAATCGAACCCTCAACTAAAGCTTGGGAAGCTTTCGTTTTGCCACTAAACTATAcc                           <  1:104302/150‑1 (MQ=255)
                                                                                                                           aGGATGAAAAACGAAAAAAGGAGCCAAGGCTCCTTTTGACTGCGATGACGTTGCTGGAACGCCATGCGGGACTGCGATGTGGAGCGGGTAGAGGGAATCGAACCCTCaa                                                                    <  2:196722/109‑1 (MQ=255)
                                                                                                                           aGGATGAAAAACGAAAAAAGGAGCCAAGGCTCCTTTTGACTGCGATGACGTTGCTGGAACGCCATGCGGGACTGCGATGTGGAGCGGGTAGAGGGAATCGAACCCTCaa                                                                    >  1:196722/1‑109 (MQ=255)
                                                                                                                            ggATGAAAAACGAAAAAAGGAGCCAAGGCTCCTTTTGACTGCGATGACGTTGCTGGAACGCCATGCGGGACTGCGATGTGGAGCGGGTAGAGGGAATCGAACCCTCAACTAAAGCTTGGGAAGCTTTCGTTTTGCCACTAAACTATAccc                          <  1:1747187/150‑1 (MQ=255)
                                                                                                                                gAAAAACGAAAAAAGGAGCCAAGGCTCCTTTTGACTGCGATGACGTTGCTGGAACGCCATGCGGGACTGCGATGTGGa                                                                                              <  2:710028/78‑1 (MQ=255)
                                                                                                                                gAAAAACGAAAAAAGGAGCCAAGGCTCCTTTTGACTGCGATGACGTTGCTGGAACGCCATGCGGGACTGCGATGTGGa                                                                                              >  1:710028/1‑78 (MQ=255)
                                                                                                                                gAAAAACGAAAAAAGGAGCCAAGGCTCCTTTTGACTGCGATGACGTTGCTGGAACGCCATGCGGGACTGCGATGTGGAGCGGGTAGAGGGAATCGAACCCTCAACTAAAGCTTGGGAAGCTTTCGTTTTGCCACTAAACTATACCCGCAt                      >  1:1018662/1‑150 (MQ=255)
                                                                                                                                  aaaaCGAAAAAAGGAGCCAAGGCTCCTTTTGACTGCGATGACGTTGCTGGAACGCCATGCGGGACTGCGATGTGGAGCGGGTAGAGGGAATCGAACCCTCAACTaaa                                                               >  1:219171/1‑107 (MQ=255)
                                                                                                                                  aaaaCGAAAAAAGGAGCCAAGGCTCCTTTTGACTGCGATGACGTTGCTGGAACGCCATGCGGGACTGCGATGTGGAGCGGGTAGAGGGAATCGAACCCTCAACTaaa                                                               <  2:219171/107‑1 (MQ=255)
                                                                                                                                  aaaaCGAAAAAAGGAGCCAAGGCTCCTTTTGACTGCGATGACGTTGCTGGAACGCCATGCGGGACTGCGATGTGGAGCGGGTAGAGGGAATCGAACCCTCAACTAAAGCTTGGGAAGCTTTCGTTTTGCCGCTAAACTATAccc                          >  2:1585701/1‑144 (MQ=255)
                                                                                                                                  aaaaCGAAAAAAGGAGCCAAGGCTCCTTTTGACTGCGATGACGTTGCTGGAACGCCATGCGGGACTGCGATGTGGAGCGGGTAGAGGGAATCGAACCCTCAACTAAAGCTTGGGAAGCTTTCGTTTTGCCACTAAACTATAccc                          <  1:1585701/144‑1 (MQ=255)
                                                                                                                                     aCGAAAAAAGGAGCCAAGGCTCCTTTTGACTGCGATGACGTAGCTGGAACGCCATGCGGGACTGCGATGTGGAGCGGGTAGAGGGAATCGAACCCTCAACTAAAGCTTGGGAAGCTTTCGTTTTGCCACTAAACTATACCCGCATCTGat                 <  2:1651362/150‑2 (MQ=255)
                                                                                                                                          aaaaGGAGCCTTGGCTCCTTTTGACTGCGATGACGTTGCTGGAACGCCATGCGGGACTGCGATGTGGAGCGGGTAGAGGGAATCGAACCCTCAACTAAAGCTTGGGAAGCTTTCGTTTTGCCACTAAACTATACCCGCATCTGAGCTCGa            >  2:2228461/1‑150 (MQ=255)
                                                                                                                                          aaaaGGAGCCTTGGCTCCTTTTGACTGCGATGACGTTGCTGGAACGCCATGCGGGACTGCGATGTGGAGCGGGTAGAGGGAATCGAACCCTCAACTAAAGCTTGGGAAGCTTTCGTTTTGCCACTAAACTATACCCGCATCTGAGCTCGa            >  2:1458179/1‑150 (MQ=255)
                                                                                                                                          aaaaGGAGCCTTGGCTCCTTTTGACTGCGATGACGTTGCTGGAACGCCATGCGGGACTACGATGTGGAGCGGGTAGAGGGAATCGAACCCTCAACTAAAGCTTGGGAAGCTTTCGTTTTGCCACTAAACTATACCCGCATCTGAGCTCGa            >  1:2195305/1‑150 (MQ=255)
                                                                                                                                             aGGAGCCAAGGCTCCTTTTGACTGCGATGACGTTGCTGGAACGCCATGCGGGACTGCGATGTGGAGCGGGTAGAGGGAATCGAACCCTCAACTAAAGCTTGGGAAGCTTTCGTTTTGCCACTAAACTATACCCGCATCTGAGCTCGACtt         <  1:581805/150‑1 (MQ=255)
                                                                                                                                               gAGCCAAGGCTCCTTTTGACTGCGATGACGTTGCTGGAACGCCATGCGGGACTGCGATGTGGAGCGGGTAGAGGGAATCGAACCCTCAACTAAAGCTTGGGAAGCTTTCGTTTTGCCACTAAACTATACCCGCATCTGAGCTCGACtttt       <  2:1707917/150‑1 (MQ=255)
                                                                                                                                                  ccAAGGCTCCTTTTGACTGCGATGACGTTGCTGGAACGCCATGCGGGACTGCGATGTGGAGCGGGTAGAGGGAATCGAACCCTCAACTAAAGCTTGGGAAGCTTTCGTTTTGCCACTAAACTATACCCGCATCTGAGCTCGACTTTTTAc    <  1:1017788/150‑1 (MQ=255)
                                                                                                                                                    aaGGCTCCTTTTGACTGCGATGACGTTGCTGGAACGCCATGCGGGACTGCGATGTGGAGCGGGTAGAGGGAATCGAACCCTCAACTAAAGCTTGGGAAGCTTTCGTTTTGCCACTAAACTATACCCGCATCTGAGCTCGACTTTTTACCa  >  1:1301810/1‑150 (MQ=255)
                                                                                                                                                    aaGGCTCCTTTTGACTGCGATGACGTTGCTGGAACGCCATGCGGGACTGCGATGTGGAGCGGGTAGAGGGAATCGAACCCTCAACTAAAGCTTGGGAAGCTTTCGTTTTGCCACTAAACTATACCCGCATCTGAGCTCGACTTTTTACCa  <  1:1225002/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                    
CTTTCGAATAGTTTCACGAATCCAGTGAGCATATCCGCCACAAGCCGTTTCGCTTGCGCACGCAAAGGCCCGATGCGCATTGAAGCGGATCGGTGGGTAATTGGCCATGGGCCAGGTCAAGGCAGGATGAAAAACGAAAAAAGGAGCCAAGGCTCCTTTTGACTGCGATGACGTTGCTGGAACGCCATGCGGGACTGCGATGTGGAGCGGGTAGAGGGAATCGAACCCTCAACTAAAGCTTGGGAAGCTTTCGTTTTGCCACTAAACTATACCCGCATCTGAGCTCGACTTTTTACCA  >  NC_002947/1606301‑1606598

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