Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,606,450 A→T 12.2% 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,4500AT12.2% 176.4 / 12.5 64intergenic (‑158/+53)PP_1407/PP_t64sulfate transporter/tRNA‑Gly
Reads supporting (aligned to +/- strand):  ref base A (25/31);  new base T (3/5);  total (28/36)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CTTTCGAATAGTTTCACGAATCCAGTGAGCATATCCGCCACAAGCCGTTTCGCTTGCGCACGCAAAGGCCCGATGCGCATTGAAGCGGATCGGTGGGTAATTGGCCATGGGCCAGGTCAAGGCAGGATGAAAAACGAAAAAAGGAGCCAAGGCTCCTTTTGACTGCGATGACGTTGCTGGAACGCCATGCGGGACTGCGATGTGGAGCGGGTAGAGGGAATCGAACCCTCAACTAAAGCTTGGGAAGCTTTCGTTTTGCCACTAAACTATACCCGCATCTGAGCTCGACTTTTTACCA  >  NC_002947/1606301‑1606598
                                                                                                                                                     |                                                                                                                                                    
cTTTCGAATAGTTTCACGAATCCAGTGAGCATATCCGCCACAAGCCGTTTCGCTTGCGCACGCAAAGGCCCGATGCGCATTGAAGCGGATCGGTGGGTAATTGGCCATGGGCCAGGTCAAGGCAGGATGAAAAACGAAAAAAGGAGCCaa                                                                                                                                                      >  1:262752/1‑150 (MQ=255)
  ttCGTATTGTTTCACGAATCCAGTGAGCATATCCGCCACAAGCCGTTTCGCTTGCGCACGCAAAGGCCCGACGCGCATTGAAGCGGATCGGTGGGTAATTGGCCATGGGCCAGGTCAAGGCAGGATGAAAAACGAAAAAAGGAGCCAAgg                                                                                                                                                    <  1:1763106/150‑1 (MQ=255)
  ttCGAATAGCTTCACGAATCCAGTGAGCATATCCGCCACAAGCCGTTTCGCTTGCGCACGCAAAGGCCCGATGCGCATTGAAGCTGATCGGTGGGTAATTGCCCATGTGCCATGTCAAGGCAGGATGAAAAACGAAAAAAGGAGCCAAgg                                                                                                                                                    <  2:363976/150‑1 (MQ=255)
      aaTAGTTTCACGAATCCAGTGAGCATATCCGCCACAAGCCGTTTCGCTTGCGCACGCAAAGGCCCGATGCGCATTGAAGCGGATCGGTGGGTAATTGGCCATGGGCCAGGTCAAGGCAGGATGAAAAACGAAAAAAGGAGCCAAGGCTcc                                                                                                                                                >  1:397022/1‑150 (MQ=255)
       aTAGTTTCACGAATCCAGTGAGCATATCCGCCACAAGCCGTTTCGCTTGCGCACGCAAAGGCCCGATGCGCATTGAAGCGGATCGGTGGGTAATTGGCCATGGGCCAGGTCAAGGCAGGATGAAAAACGAAAAAAGGAGCCAAGGCTCCt                                                                                                                                               >  1:1901338/1‑150 (MQ=255)
         agtcccACCAATCCAGTCAGCACATCCGCCACAAGCCGTTTCGCTTGCGCACGCAAAGGCCCGATGCGCATTGAAGCGGATCGGTGGGTAATTGGCCATGGGCCAGGTCAAGGCAGGATGAAAAACGAAAAAAGGAGCCAAGGCTCCttt                                                                                                                                             <  2:1545264/145‑1 (MQ=255)
          gTTTCCCGAATTCAGTGAGCATCTCCCCCACATGCCGTTTTGCTTGCGCCCGCAAAGGCCCGATGCGCATTGAAGCGGATCGGTGGGTAATTGGCCATGGGCCAGGTCAAGGCAGGATGAAAAACGAAAAAAGGAGCCTTGGCTCCtttt                                                                                                                                            <  1:150934/150‑1 (MQ=255)
             tCACGAATCCAGTGAGCATATCCGCCACAAGCCGTTTCGCTTGCGCACGCAAAGGCCCGATGCGCATTGAAGCGGATCGGTGGGTAATTGGCCATGGGCCAGGTCAAGGCAGGATGAAAAACGAAAAAAGGAGCCAAGGCTCCTTTTGAc                                                                                                                                         >  1:267578/1‑150 (MQ=255)
               aCGAATCCAGTGAGCATATCCGCCACAAGCCGTTTCGCTTGCGCACGCAAAGGCCCGATGCGCATTGAAGCGGATCGGTGGGTAATTGGCCATGGGCCAGGTCAAGGCAGGATGAAAAACGAAAAAAGGAGCCAAGGCTCCTTTTGACTg                                                                                                                                       >  1:1413182/1‑150 (MQ=255)
                 gAATCCAGTGAGCATATCCGCCACAAGCCGTTTCGCTTGCGCACGCAAAGGCCCGATGCGCATTGAAGCGGATCGGTGGGTAATTGGCCATGGGCCAGGTCAAGGCAGGATGAAAAACGAAAAAAGGAGCCTTGGCTCCttttttcgttt                                                                                                                                     <  1:147154/150‑8 (MQ=255)
                   gtacAGTGAGCATATCCGCCACAAGCCGTTTCGCGTGCGCAGGCAAAGGCCCGATGCGCATTGAAGCGGATCGGTGGGTAATTGGCCATGGGCCAGGTCAAGGCAGGATGAAAAACGAAAAAAGGAGCCAAGGCTCCTTTTGACTGCGAt                                                                                                                                   >  2:829123/4‑150 (MQ=255)
                   aTCCAGTGAGCATATCCGCCACAAGCCGTTTCGCTTGCGCACGCAAAGGCCCGATGCGCATTGAAGCGGATCGGTGGGTAATTGGCCATGGGCCAGGTCAAGGCAGGATGAAAAACGAAAAAAGGAGCCAAGGCTCCTTTTGACTGCGAt                                                                                                                                   >  1:1486408/1‑150 (MQ=255)
                     ccAGCGAGCATACCCGCGACAGCCCGGTGCGCGTGCGCACGCAAAGGCCCGATGCGCATTGAAGGGGATCGGTGGGTAATTGGCCATGGGCCAGGTCAAGGCAGGATGAAAAACGAAAAAAGGAGCCAAGGCTCCTTTTGACTGCGATGa                                                                                                                                 <  1:1912193/150‑1 (MQ=255)
                       aGTGAGCATATCCGCCACAAGCCGTTTCGCTTGCGCACGCAAAGGCCCGATGCGCATTGAAGCGGATCGGTGGGTAATTGGCCATGGGCCAGGTCAAGGCAGGATGAAAAACGAAAAAAGGAGCCAAGGCTCCTTTTGACTGCGATGACg                                                                                                                               >  1:669072/1‑150 (MQ=255)
                       aGTGAGCATATCCGCCACAAGCCGTTTCGCTTGCGCACGCAAAGGCCCGATGCGCATTGAAGCGGATCGGTGGGTAATTGGCCATGGGCCAGGTCAAGGCAGGATGAAAAACGAAAAAAGGAGCCAAGGCTCCTTTTGACTGCGATGACg                                                                                                                               >  1:248887/1‑150 (MQ=255)
                        gTGAGCATATCCGCCACAAGCCGTTTCGCTTGCGCACGCAAAGGCCCGATGCGCATTGAAGCGGATCGGTGGGTAATTGGCCATGGGCCAGGTCAAGGCAGGATGAAAAACGAAAAAAGGAGCCAAGGCTCCTTTTGACTGCGATGACGt                                                                                                                              >  2:1804907/1‑150 (MQ=255)
                           aGCATATCCGCCACAAGCCGTTTCGCTTGCGCACGCAAAGGCCCGATGCGCATTGAAGCGGATCGGTGGGTAATTGGCCATGGGCCAGGTCAAGGCAGGATGAAAAACGAAAAAAGGAGCCAAGGCTCCTTTTGACTGCGATGACGTTGc                                                                                                                           >  1:1828532/1‑150 (MQ=255)
                                atCCGCCACAAGCCGTTTCGCTTGCGCACGCAAAGGCCCGATGCGCATTGAAGCGGATCGGTGGGTAATTGGCCATGGGCCAGGTCAAGGCAGGATGAAAAACGAAAAAAGGAGCCAAGGCCCCTTTTGACTGCGATGACGTTGCTGGaa                                                                                                                      >  1:277596/1‑150 (MQ=255)
                                 tCCGCCACAAGCCGTTTCGCTTGCGCACGCAAAGGCCCGATGCGCATTGAAGCGGATCGGTGGGTAATTGGCCATGGGCCAGGTCAAGGCAGGATGAAAAACGAAAAAAGGAGCCAAGGCTCCTTTTGACTGCGATGACGTTGCTGGAAc                                                                                                                     >  2:110837/1‑150 (MQ=255)
                                   cGCCACAAGCCGTTTCGCTTGCGCACGCAAAGGCCCGATGCGCATTGAAGCGGATCGGTGGGTAATTGGCCATGGGCCAGGTCAAGGCAGGATGAAAAACGAAAAAAGGAGCCAAGGCTCCTTTTGACTGCGATGACGTTGCTGGAACGc                                                                                                                   >  2:1167223/1‑150 (MQ=255)
                                   cGCCACAAGCCGTTTCGCTTGCGCACGCAAAGGCCCGATGCGCATTGAAGCGGATCGGTGGGTAATTGGCCATGGGCCAGGTCAAGGCAGGATGAAAAACGAAAAAAGGAGCCAAGGCTCCTTTTGACTGCGATGACGTTGCTGGAACGc                                                                                                                   <  2:1901338/150‑1 (MQ=255)
                                        caAGCCGTTTCGCTTGCGCACGCAAAGCCCCGATGCGCATTGAAGCGGATCGGTGGGTAATTGGCCATGGGCCAGGTCCAGGCAGGATGAAACGGTATGCAGCTTTCTATTGCTGCTTTTGACTGCGATGACGTTGCTGGAACGCCATGc                                                                                                              <  1:1470123/150‑1 (MQ=255)
                                               tttCGCTTGCGCACGCAAAGGCCCGATGCGCATTGAAGCGGATCGGTGGGTAATTGGCCATGGGCCAGGTCAAGGCAGGATGAAAAACGAAAAAAGGAGCCAAGGCTCCTTTTGACTGCGATGACGTTGCTGGAACGCCATGCGGGACTg                                                                                                       <  2:1413182/150‑1 (MQ=255)
                                               tttCGCTTGCGCACGCAAAGGCCCGATGCGCATTGAAGCGGATCGGTGGGTAATTGGCCATGGGCCAGGTAAAGGCAGGATGAAAAACGAAAAAAGGAGCCAAGGCTCCTTTTGAATGCGATGACGTTGCTGGAACGCCATGCGGGACTg                                                                                                       >  2:792954/1‑150 (MQ=255)
                                                ttCGCTTGCGCACGCAAAGGCCCGATGCGCATTGAAGCGGATCGGTGGGTAATTGGCCATGGGCCAGGTCAAGGCAGGATGAAAAACGAAAAAAGGAGCCAAGGCTCCTTTTGACTGCGATGACGTTGCTGGAACGCCATGCGGGACTGc                                                                                                      <  2:397022/150‑1 (MQ=255)
                                                   gCTTGCGCACGCAAAGGCCCGATGCGCATTGAAGCGGATCGGTGGGTAATTGGCCATGGGCCAGGTCAAGGCAGGATGAAAAACGAAAAAAGGAGCCAAGCCACCTTTTGACTGCGATGACGTTGCTGGAACGCCATGCGGGACTGCGAt                                                                                                   <  1:902763/150‑1 (MQ=255)
                                                         gcacgcaAAGGCCCGATGCGCATTGAAGCGGATCGGTGGGTAATTGGCCATGGGCCAGGTCAAGGCAGGATGAAAAACGAAAAAAGGAGCCAAGGCTCCTTTTGACTGCGATGACGTTGc                                                                                                                           >  2:1916339/1‑120 (MQ=255)
                                                         gcacgcaAAGGCCCGATGCGCATTGAAGCGGATCGGTGGGTAATTGGCCATGGGCCAGGTCAAGGCAGGATGAAAAACGAAAAAAGGAGCCAAGGCTCCTTTTGACTGCGATGACGTTGc                                                                                                                           <  1:1916339/120‑1 (MQ=255)
                                                            cgcaAAGGCCCGATGCGCATTGAAGCGGATCGGTGGGTAATTGGCCATGGGCCAGGTCAAGGCAGGATGAAAAACGAAAAAAGGAGCCAAGGCTCCTTTTGACTGCGATGACGTTGCTGGAACGCCATGCGGGACTGCGATGTGGAGCgg                                                                                          <  1:1850625/150‑1 (MQ=255)
                                                            cgcaAAGGCCCGATGCGCATTGAAGCGGATCGGTGGGTAATTGGCCATGGGCCAGGTCAAGGCAGGATGAAAAACGAAAAAAGGAGCCAAGGCTCCTTTTGACTGCGATGACGTTGCTGGAACGCCATGCGGGACTGCGATGTGGAGCgg                                                                                          <  2:669072/150‑1 (MQ=255)
                                                             gcaAAGGCCCGATGCGCATTGAAGCGGATCGGTGGGTAATTGGCCATGGGCCAGGTCAAGGCAGGATGAAAAACGAAAAAAGGAGCCAAGGCTCCTTTTGACTGCGATGACGTTGCTGGAACGCCATGCGGGACTGCGATGTGGAGCggg                                                                                         <  2:1828532/150‑1 (MQ=255)
                                                             gcaAAGGCCCGATGCGCATTGAAGCGGATCGGTGGGTAATTGGCCATGGGCCAGGTCAAGGCAGGATGAAAAACGAAAAAAGGAGCCAAGGCTCCTTTTGACTGCGATGACGTTGCTGGAACGCCATGCGGGACTGCGATGTGGAGCggg                                                                                         >  2:389228/1‑150 (MQ=255)
                                                                 aGGCCCGATGCGCATTGAAGCGGATCGGTGGGTAATTGGCCATGGGCCAGGTCAAGGCAGGATGAAAAACGAAAAAAGGAGCCAAGGCTCCTTTTGACTGCGATGACGTTGCTGGAACGCCATGCGGGACTGCGATGTGGAGCGGGTaga                                                                                     <  2:248887/150‑1 (MQ=255)
                                                                  gtcccGATGCGCATTGAAGCGGAGCGGTGGGTAATTGGGGGGGGGGGGGGGGAGGGCAGGATGAAAAACGAAAAAAGGAGCCAAGGCTCCTTTTGACTGCGATGACGTTGCTGGAACGCCATGCGGGACTGCGATGTGGAGCGGGTagag                                                                                    <  1:461104/148‑1 (MQ=255)
                                                                  ggaccGATGCGCATTGCAGCGGACCGGTGGGTAATTGGCCATGGGCCAGGTCAAGGCAGGATAAAAAAAAAAAAAAGGAGCCAAGGCTCCTTTTGACTGCGATGACGTTGCTGGAACGCCATGCGGGACTGCGATGTGGAGCGGGTagag                                                                                    <  2:956983/147‑1 (MQ=255)
                                                                  ggCCCGATGCGCATTGAAGCGGATCGGTGGGTAATTGGCCATGGGCCAGGTCAAGGCAGGATGAAAAACGAAAAAAGGAGCCAAGGCTCCTTTTGACTGCGATGACGTTGCTGGAACGCCATGCGGGACTGCGATGTGGAGCGGGTagag                                                                                    <  2:579916/150‑1 (MQ=255)
                                                                   gCCCGATGCGCATTGAAGCGGATCGGTGGGTAATTGGCCATGGGCCAGGTCAAGGCAGGATGAAAAACGAAAAAAGGAGCCAAGGCTCCTTTTGACTGCGATGACGTTGCTGGAACGCCATGCGGGACTGCGATGTGGAGCGGGTAGAgg                                                                                   <  1:1593947/150‑1 (MQ=255)
                                                                     ccGATGCGCATTGAAGCGGATCGGTGGGTAATTGGCCATGGGCCAGGTCAAGGCAGGATGAAAAACGAAAAAAGGAGCCAAGGCTCCTTTTGACTGCGATGACGTTGCTGGAACGCCATGCGGGACTGCGATGTGGAGCGGGTAGAGGGa                                                                                 <  1:389228/150‑1 (MQ=255)
                                                                        aTGCGCATTGAAGCGGATCGGTGGGTAATTGGCCATGGGCCAGGTCAAGGCAGGATGAAAAACGAAAAAAGGAGCCAAGGCTCCTTTTGACTGCGATGACGTTGCTGGAACGCCATGCGGGACTGCGATGTGGAGCGGGTAGAGGGAATc                                                                              <  2:28582/150‑1 (MQ=255)
                                                                        aTGCGCATTGAAGCGGATCGGTGGGTAATTGGCCATGGGCCAGGTCAAGGCAGGATGAAAAACGAAAAAAGGAGCCAAGGCTCCTTTTGACTGCGATGACGTTGCTGGAACGCCATGCGGGACTGCGATGTGGAGCGGGTAGAGGGAATc                                                                              <  1:1381308/150‑1 (MQ=255)
                                                                            gcATTGAAGCGGATCGGTGGGTAATTGGCCATGGGCCAGGTCAAGGCAGGATGAAAAACGAAAAAAGGAGCCAAGGCTCCTTTTGACTGCGATGACGTTGCTGGAACGCCATGCGGGACTGCGATGTGGAGCGGGTAGAGGGAATCGAAc                                                                          <  1:829123/150‑1 (MQ=255)
                                                                              aTTGAAGCGGATCGGTGGGTAATTGGCCATGGGCCAGGTCAAGGCAGGATGAAAAACGAAAAAAGGAGTCAAGGCTCCTTTTGACTGCGATGACGTTGCTGGAACGCCATGCGGGACTGCGATGTGGAGc                                                                                            >  1:893259/1‑130 (MQ=255)
                                                                              aTTGAAGCGGATCGGTGGGTAATTGGCCATGGGCCAGGTCAAGGCAGGATGAAAAACGAAAAAAGGAGTCAAGGCTCCTTTTGACTGCGATGACGTTGCTGGAACGCCATGCGGGACTGCGATGTGGAGc                                                                                            <  2:893259/130‑1 (MQ=255)
                                                                                       gATCGGTGGGTAATTGGCCATGGGCCAGGTCAAGGCAGGATGAAAAACGAAAAAAGGAGCCAAGGCTCCTTTTGACTGCGATGACGTTGCTGGAACGCCATGCGGGACTGCGATGTGGAGCGGGTAGAGGGAATCGAACCCTCAACTaaa                                                               >  1:1742664/1‑150 (MQ=255)
                                                                                            gTGGGTAATTGGCCATGGGCCAGGTCAAGGCAGGATGAAAAACGAAAAAAGGATCCAAGGCTCCTTTTGACTGCGATGACGTTGCTGGAACGCCATGCGGGACTGCGATGTGGAGCGGGTAGAGGGAATCGAACCCTCAACTAAAGCTTg                                                          >  1:1429127/1‑150 (MQ=255)
                                                                                             tGGGTAATTGGCCATGGGCCAGGTCAAGGCAGGATGAAAAACGAAAAAAGGAGCCTTGGCTCCttttt                                                                                                                                           <  2:1896541/68‑2 (MQ=38)
                                                                                             tGGGTAATTGGCCATGGGCCAGGTCAAGGCAGGATGAAAAACGAAAAAAGGAGCCTTGGCTCCttttt                                                                                                                                           >  1:1896541/1‑67 (MQ=38)
                                                                                               ggTAATTGGCCATGGGCCAGGTCAAGGCAGGATGAAAAACGAAAAAAGGAGCCAAGGCTCCTTTTGACTGCGATGACGTTGCTGGAACGCCATGCGGACCGCCGATGTGGAGCGGATATAGGGAATCGAACCCTCAACTAAAGCTTGGGa                                                       >  2:990735/1‑150 (MQ=255)
                                                                                                      ggCCATGGGTCAGGTCAAAGCAGGATGAAAAACGAAAAAAGGAGCCAAGGCTCCTTTTGACTGCGATGACGTTGCTGGAACGCCATGCGGGACTGCGATGTGGAGCGGGTAGAGGGAATCGAACCCTCAACTAAAGCTTGGGAAGCTTTc                                                <  2:504340/150‑1 (MQ=255)
                                                                                                           tGGGCCAGGTCAAGGCAGGATGAAAAACGAAAAAAGGAGCCAAGGCTCCTTTTGACTGCGATGACGTTGCTGGAACGCCATGCGGGACTGCGATGTGGAGCGGGTAGAGGGAATCGAACCCTCAACTAAAGCTTGGGAAGCTTTCGtttt                                           <  1:932229/150‑1 (MQ=255)
                                                                                                             ggCCAGGTCAAGGCAGGATGAAAAACGAAAAAAGGAGCCAAGGCTCCTTTTGACTGCGATGACGTTGCTGGAACGCCATGCGGGACTGCGATGTGGAGCGGGTAGAGGGAATCGAACCCTCAACTAAAGCTTGGGAAGCTTTCGTTTTGc                                         <  2:267578/150‑1 (MQ=255)
                                                                                                                    tCAAGGCAGGATGAAAAACGAAAAAAGGAGCCAAGGCTCCTTTTGACTGCGATGACGTTGCTGGAACGCCATGCGGGACTGCGATGTGGAGCGGGTAGAGGGAATCGAACCCTCAACTAAAGCTTGGGAAGCTTTCGTTTTGCCACTaaa                                  <  2:1486408/150‑1 (MQ=255)
                                                                                                                      aaGGCAGGATGAAAAACGAAAAAAGGAGCCAAGGCTCCTTTTGACTGCGATGACGTTGCTGGAACGCCATGCGGGACTGCGa                                                                                                    <  1:1012037/82‑1 (MQ=255)
                                                                                                                      aaGGCAGGATGAAAAACGAAAAAAGGAGCCAAGGCTCCTTTTGACTGCGATGACGTTGCTGGAACGCCATGCGGGACTGCGa                                                                                                    >  2:1012037/1‑82 (MQ=255)
                                                                                                                      aaGGCAGGATGAAAAACAAACCAAGATGTTAAGGTTTGTCTCTACTGCGATGACGTGGCTTTAACGCCATGCTGGACTGCGAGGGGGAGCGGGTAGAGGGAATCGAACCCTCAACTAAAGCTTGGGAAACTTTCGTTGTGCCACTAAACt                                >  2:960775/1‑150 (MQ=255)
                                                                                                                                       gtcaaaaGGAGCCTTGGCTCCTTTTGACTGCGATGACGTTGCTGGAACGCCATGCGGGACTGCGATGTGGAGCGGGTAGAGGGAATCGAACCCTCAACTAAAGCTTGGGAAGCTTTCGTTTTGCCACTAAACTATACCCGCATCTGAGCt               >  2:672457/4‑150 (MQ=255)
                                                                                                                                        aaaaaaGGAGCCAAGGCTCCTTTTGACTGCGATGACGTTGCTGGAACGCCATGCGGGACTGCGATGTGGAGCGGGTAGAGGGAATCGAACCCTCAACTAAAGCTTGGGAAGCTTTCGTTTTGCCACTAAACTATACCCGCATCTGAGCTc              >  1:1424147/1‑150 (MQ=255)
                                                                                                                                          aaaaGGAGCCTTGGCTCCTTTTGACTGCGATGACGTTGCTGGAACGCCATGCGGGACTGCGATGTGGAGCGGGTagag                                                                                    <  2:1190835/78‑1 (MQ=39)
                                                                                                                                          aaaaGGAGCCTTGGCTCCTTTTGACTGCGATGACGTTGCTGGAACGCCATGCGGGACTGCGATGTGGAGCGGGTagag                                                                                    >  1:1190835/1‑78 (MQ=39)
                                                                                                                                           aaaGGAGCCAAGGCTCCTTTTGACTGCGATGACGTTGATTTAACGCCATGCGGTACTGCGATGTGGAGCGGGTAGAGGGAATCGAACCCTCAACTAAAGCTTGGGAAGCTTTCGTTTTGCCACTAAACTATACCCGCATCTGAGCTCGAc           <  1:1374269/150‑1 (MQ=255)
                                                                                                                                            aaGGAGCCAAGGCTCCTTTTGACTGCGATGACGTTGCTGGAACGCCATGCGGGACTGCGATGTGGAGCGGGTAGAGGGAATCGaa                                                                           >  1:1365652/1‑85 (MQ=255)
                                                                                                                                            aaGGAGCCAAGGCTCCTTTTGACTGCGATGACGTTGCTGGAACGCCATGCGGGACTGCGATGTGGAGCGGGTAGAGGGAATCGaa                                                                           <  2:1365652/85‑1 (MQ=255)
                                                                                                                                             aGGAGCCAAGGCTCCTTTTGACTGCGATGACGTTGCTGGAACGCCATGCGGGACTGCGATGTGGAGCGGGTAGAGGGAATCGAACCCTCAACTAAAGCTTGGGAAGCTTTCGTTTTGCCACTAAACTATACCCGCATCTGAGCTCGACtt         >  2:533200/1‑150 (MQ=255)
                                                                                                                                                 gCCAAGGCTCCTTTTGACTGCGATGACGTTGCTGGAACGCCATGCGGGACTGCGATGTGGAGCGGGTAGAGGTAATCGAACCCTCAACTAAAGCTTGGGAAGCTTTCGTTTTGCCACTAACCTATACCCGCATCTGAGCTCGACTTTTTa     <  2:596708/150‑1 (MQ=255)
                                                                                                                                                   cAAGGCTCCTTTTGACTGCGATGACGTTGCTGGAACGCCATGCGGGACTGCGATGTGGAGCGGGTAGAGGGAATCGAACCCTCAACTAAAGCTTGGGAAGCTTTCGTTTTGCCACTAAACTATACCCGCATCTGAGCTCGACTTTTTAcc   <  1:1167223/150‑1 (MQ=255)
                                                                                                                                                    aaGGCTCCTTTTGACTGCGATGACGTTGCTGGAACGCCATGCGGGACTGCGATGTGGAGCGGGTAGAGGGAATCGAACCCTCAACTAAAGCTTGGAAAAGTATCGTTTTGCCACTAAACTATACCCGAATCTGAGCTCGACTTTTTACCa  >  1:1551937/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                    
CTTTCGAATAGTTTCACGAATCCAGTGAGCATATCCGCCACAAGCCGTTTCGCTTGCGCACGCAAAGGCCCGATGCGCATTGAAGCGGATCGGTGGGTAATTGGCCATGGGCCAGGTCAAGGCAGGATGAAAAACGAAAAAAGGAGCCAAGGCTCCTTTTGACTGCGATGACGTTGCTGGAACGCCATGCGGGACTGCGATGTGGAGCGGGTAGAGGGAATCGAACCCTCAACTAAAGCTTGGGAAGCTTTCGTTTTGCCACTAAACTATACCCGCATCTGAGCTCGACTTTTTACCA  >  NC_002947/1606301‑1606598

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