Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,642,334 T→C 6.9% intergenic (‑120/+32) PP_1442 ← / ← lon‑I hypothetical protein/DNA‑binding ATP‑dependent protease

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029471,642,3340TC 6.9% 305.4 / 8.9 102intergenic (‑120/+32)PP_1442/lon‑Ihypothetical protein/DNA‑binding ATP‑dependent protease
Reads supporting (aligned to +/- strand):  ref base T (48/47);  new base C (3/4);  total (51/51)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.49e-01

Only 100 of 102 total aligned reads displayed.

AGGGGAACGAGCAGACGAGGAGAGGTCATGGTGGCTCCGAAGGCTGGCGATGGCCTGAGGATACCCGATGATCGAGGAAATTTGTGCTGCCTGTGCCGGTCCTTTCGCGAGTAAGCCCGCGCCTGCAGGGGTAGCGTACCCGTGCAGGTGCGGGCTTGCCCGCGAAGAGGGCGGTGCAGGCTAGAACAGCACTTTCGCCACGTCGGCAAAGCGCTTGGCGAAATGCACTGTCAGGCCTTCGCGCAGGTAGTCCGGCAGTTCTTCGAAGTCGCCGCGGTTGGGCTCCGGCAGGATCA  >  NC_002947/1642186‑1642481
                                                                                                                                                     |                                                                                                                                                   
aGGGGAACGAGCAGACGAGGAGAGGTCATGGTGGCTCCGAAGGCTGGCGATGGCCTGAGGATACCCGATGATCGAGGAAATTTGTGCTGCCTGTGCCGGTCCTTTCGCGAGTAAGCCCGCGCCTGCAGGGGTAGCGTACCCGTGCAGGTg                                                                                                                                                    >  2:1145731/1‑150 (MQ=255)
aGGGGAACGAGCAGACGAGGAGAGGTCATGGTGGCTCCGAAGGCTGGCGATGGCCTGAGGATACCCGATGATCGAGGAAATTTGTGCTGCCTGTGCCGGTCCTTTCGCGAGTAAGCCCGCGCCTGCAGGGGTAGCGTACCCGTGCAGGTg                                                                                                                                                    <  1:1160853/150‑1 (MQ=255)
  gggAACGAGCAGACGAGGAGAGGTCATGGTGGCTCCGAAGGCTGGCGATGGCCTGAGGATACCCGATGATCGAGGAAATTTGTGCTGCCTGTGCCGGTCCTTTCGCGAGTAAGCCCGCGCCTGCAGGGGTAGCGTACCCGTGCAGGTGCg                                                                                                                                                  <  2:693746/150‑1 (MQ=255)
   ggAACGAGGAGGCGAGGAGAGGTCATGGTGGCTCCGGAGGCCGGCGATGGCCTGAGGATACCCGATGATCGAGGAAATTTGTTCTGCCTGTGCCGGTTCGTTCGCGAGTAAGCCCGCGCCTGCAAGGGTAGCGTACCCGTGCAGGTGCgg                                                                                                                                                 <  1:1531135/150‑1 (MQ=255)
   ggAACGAGCAGACGAGGAGAGGTCATGGTGGCTCCGAAGGCTGGCGATGGCCTGAGGATACCCGATGATCGAGGAAATTTGTGCTGCCTGTGCCGGTCCTTTCGCGAGTAAGCCCGCGCCTGCAGGGGTAGCGTACCCGTGCAGGTGCgg                                                                                                                                                 <  2:1625358/150‑1 (MQ=255)
    gAACGAGCAGACGAGGAGAGGTCATGGTGGCTCCGAAGGCTGGCGATGGCCTGAGGATACCCGATGATCGAGGAAATTTGTGCTGCCTGTGCCGGTCCTTTCGCGAGTAAGCCCGCGCCTGCAGGGGTAGCGTACCCGTGCAGGTGCggg                                                                                                                                                >  1:320406/1‑150 (MQ=255)
     aaCGAGCAGACGAGGAGAGGTCATGGTGGCTCCGAAGGCTGGCGATGGCCTGAGGATACCCGATGATCGAGGAAATTTGTGCTGCCTGTGCCGGTCCTTTCGCGAGTAAGCCCGCGCCTGCAGGGGTAGCGTACCCGTGCAGGTGCGGGc                                                                                                                                               <  1:1145731/150‑1 (MQ=255)
       cGAGCAGACGAGGAGAGGTCATGGTGGCTCCGAAGGCTGGCGATGGCCTGAGGATACCCGATGATCGAGGCAATTCGTGCTGCCTGTGCCGGTCCTTTCGCGAGTAAGCCCGCGCCTGCAGGGGTAGCGTACCCGTGCAGGTGCGGGCtt                                                                                                                                             <  1:199302/150‑1 (MQ=255)
         aGCAGACGAGGAGAGGTCATGGTGGCTCCGAAGGCTGGCGATGGCCTGAGGATACCCGATGATCGAGGAAATTTGTGCTGCCTGTGCCGGTCCTTTCGCGGGCAAGCCCGCACCTGCACGGGTACGCTACCCCTGCAGGCGCGGGCttac                                                                                                                                           >  2:1446631/1‑148 (MQ=25)
          gCAGACGAGGAGAGGTCATGGTGGCTCCGAAGGCTGGCGATGGCCTGAGGATACCCGATGATCGAGGAAATTTGTGCTGCCCGTGCCGGTCCTTTCGCGAGTAAGCCCGCGCCTGCAGGGGTAGCGTACCCGTGCAGGTGCGGGCTTGcc                                                                                                                                          >  2:783911/1‑150 (MQ=255)
             gACGAGGAGAGGTCATGGTGGCTCCGAAGGCTGGCGATGGCCTGAGGATACCCGATGATCGAGGAAATTTGTGCTGCCTGTGCCGGTCCTTTCGCGAGTAAGCCCGCGCCTGCAGGGGTAGCGTACCCGTGCAGGTGCGGGCTTGCCagc                                                                                                                                       >  2:343058/1‑150 (MQ=255)
              aCGAGGAGAGGTCATGGTGGCTCCGAAGGCTGGCGATGGCCTGAGGATACCCGATGATCGAGGAAATTTGTGCTGCCTGTGCCGGTCCTTTCGCGAGTAAGCCCGCGCCTGCAGGGGTGGCGTACCCGTGCAGATGCGGGCTTGcgcgcg                                                                                                                                      >  1:657296/1‑150 (MQ=255)
                gaggagAGGTCATGGTGGCTCCGAAGGCTGGCGATGGCCTGAGGATACCCGATGATCGAGGAAATTTGTGCTGCCTGTGCCGGTCCTTTCGCGAGTAAGCCCGCGCCTGTAGGGGTAGCGTACCCGTGCAGGTGCGGGCTTGCCCGCGaa                                                                                                                                    >  2:1248198/1‑150 (MQ=255)
                 aggagAGGTCATGGTGGCTCCGAAGGCTGGCGATGGCCTGAGGATACCCGATGATCGAGGAAATTTGTGCTGCCTGTGCCGGTCCTTTCGCGGGCAAGCCCGCACCTGCACGGGTACGCTACCCCTGCAGGCGCGGGCTTACTCGCGaaa                                                                                                                                   <  1:1446631/150‑2 (MQ=18)
                 aggagAGGTCATGGTGGCTCCGAAGGCTGGCGATGGCCTGAGGATACCCGATGATCGAGGAAATTTGTGCTGCCTGTGCCGGTCCTTTCGCGAGTAAGCCCGCACCTGCACGGGTACGCTACCCCTGCAGGCGCGGGCTTACTCGCGaaa                                                                                                                                   <  2:587545/150‑2 (MQ=25)
                     gagGTCATGGTGGCTCCGAAGGCTGGCGATGGCCTGAGGATACCCGATGATCGAGGAAATTTGTGCTGCTTGTGCCGGTCCTTTCGCGAGTAAGCCCGCGCCTGCAGGGGTAGCGTACCCGTGCAGGTGCGGGCTTGCCCGCGAAGAggg                                                                                                                               <  1:27603/150‑1 (MQ=255)
                       ggTCATGGTGGCTCCGAAGGCTGGCGATGGCCTGAGGATACCCGATGATCGAGGAAATTTGTGCTGCCTGTGCCGGTCCTTTCGCGAGTAAGCCCGCGCCTGCAGGGGTAGCGTACCCGTGCAGGTGCGGGCTTGCCCGCGAAGAGGGCg                                                                                                                             <  2:1447385/150‑1 (MQ=255)
                        gTCATGGTGGCTCCGAAGGCTGGCGATGGCCTGAGGATACCCGATGATCGAGGAAATTTGTGCTGCCTGTGCCGGTCCTTTCGCGAGTAAGCCCGCGCCTGCAGGGGTAGCGTACCCGTGCAGGTGCGGGCTTGCCCGCGAAGAggg                                                                                                                               >  1:108848/1‑147 (MQ=255)
                        gTCATGGTGGCTCCGAAGGCTGGCGATGGCCTGAGGATACCCGATGATCGAGGAAATTTGTGCTGCCTGTGCCGGTCCTTTCGCGAGTAAGCCCGCGCCTGCAGGGGTAGCGTACCCGTGCAGGTGCGGGCTTGCCCGCGAAGAGGGCgg                                                                                                                            >  2:342216/1‑150 (MQ=255)
                         tCATGGTGGCTCCGAAGGCTGGCGATGGCCTGAGGATACCCGATGATCGAGGAAATTTGTGCTGCCTGTGCCGGTCCTTTCGCGAGTAAGCCCGCGCCTGCAGGGGTAGCGTACCCGTGCAGGTGCGGGCTTGCCCGCGAAGAGGGCGGt                                                                                                                           >  1:762788/1‑150 (MQ=255)
                         tCATGGTGGCTCCGAAGGCTGGCGATGGCCTGAGGATACCCGATGATCGAGGAAATTTGTGCTGCCTGTGCCGGTCCTTTCGCGAGTAAGCCCGCGCCTGCAGGGGTAGCGTACCCGTGCAGGTGCGGGCTTGCCCGCGAAGAGGGCGGt                                                                                                                           >  2:1552882/1‑150 (MQ=255)
                         tCATGGTGGCTCCGAAGGCTGGCGATGGCCGGAGGATACCCGATGATCGAGGAAATTTGTGCTGCCTGTGCCGGTCCTTTCGCGGGCAAGCCCGCACCTGCACGGGTACGCTACCCCTGCAGGCGCgg                                                                                                                                                 <  2:1118443/128‑1 (MQ=21)
                         tCATGGTGGCTCCGAAGGCTGGCGATGGCCAGAGGATACCCGATGATCGAGGAAATTTGTGCTGCCTGTGCCGGTCCTTTCGCGGGCAAGCCCGCACCTGCACGGGTACGCTACCCCTGCAGGCGCgg                                                                                                                                                 >  1:1118443/1‑128 (MQ=21)
                           aTGGTGGCTCCGAAGGCTGGCGATGGCCTGAGGATACCCGATGATCGAGGAAATTTGTGCTGCCTGTGCCGGTCCTTTCGCGAGTAAGCCCGCGCCTGCAGGGGTAGCGTACCCGTGCAGGTGCGGGCTTGCCCGCGAAGAGGGCGGTGc                                                                                                                         >  1:665519/1‑150 (MQ=255)
                              gtggCTCCGAAGGCTGGCGATGGCCTGAGGATACCCGATGATCGAGGAAATTTGTGCTGCCTGTGCCGGTCCTTTCGCGAGTAAGCCCGCGCCTGCAGGGGTAGCGTACCCGTGCAGGTGCGGGCTTGCCCGCGAAGAGGGCGGTGCAgg                                                                                                                      <  1:1552882/150‑1 (MQ=255)
                                   tCCGAAGGCTGGCGATGGCCTGAGGATACCCGATGATCGAGGAAATTTGTGCTGCCTGTGCCGGTCCTTTCGCGAGTAAGCCCGCGCCTGCAGGGGTAGCGTACCCGTGCAGGTGCGGGCTTGCCCGCGAAGAGGGCGGTGCAGGCTAGa                                                                                                                 >  1:26278/1‑150 (MQ=255)
                                    ccGAAGGCTGGCGATGGCCTGAGGATACCCGATGATCGAGGAAATTTGTGCTGCCTGTGCCGGTCCTTTCGCGAGTAAGCCCGCGCCTGCAGGGGTAGCGTACCCGTGCAGGTGCGGGCTTGCCCGCGAAGAGGGCGGTGCAGGCTAGaa                                                                                                                >  2:235967/1‑150 (MQ=255)
                                        aGGCTGGCGATGGCCTGAGGATACCCGATGATCGAGGAAATTTGTGCTGCCTGTGCCGGTCCTTTCGCGAGTAAGCCCGCGCCTGCAGGGGTAGCGTACCCGTGCAGGTGCGGGCTTGCCCGCGAAGAGGGCGGTGCAGGCTAGAACAGc                                                                                                            <  2:762788/150‑1 (MQ=255)
                                                           gATACCCGATGATCGAGGAAATTTGTGCTGCCTGTGCCGGTCCTTTCGCGAGTAAGCCCGCGCCTGCAGGGGTAGCGTACCCGTGCAGGTGCGGGCTTGCCCGCGAAGAGGGCGGTGCAGGCTAGAACAGCACTTTCGCCACGTCgg                                                                                            >  1:1520450/1‑147 (MQ=255)
                                                           gATACCCGATGATCGAGGAAATTTGTGCTGCCTGTGCCGGTCCTTTCGCGAGTAAGCCCGCGCCTGCAGGGGTAGCGTACCCGTGCAGGTGCGGGCTTGCCCGCGAAGAGGGCGGTGCAGGCTAGAACAGCACTTTCGCCACGTCgg                                                                                            <  2:1520450/147‑1 (MQ=255)
                                                               cccGATGATCGAGGAAATTTGTGCTGCCTGTGCCGGTCCTTTCGCGAGTAAGCCCGCGCCTGCAGGGGTAGCGTACCCGTGCAGGTGCGGGATTGCCCGCGAAGAGGGCGGTGCAGGCTAGAACAGCACTTTCGCCACGTCGGCAAAgcg                                                                                     >  2:219883/1‑150 (MQ=255)
                                                                ccGATGATCGAGGAAATTTGTGCTGCCTGTGCCGGTCCTTTCGCGAGTAAGCCCGCGCCTGCAGGGGTAGCGTACCCGTGCAGGTGCGGGCTTGCCCGCGAAGAGGGCGGTGCAGGCTAGaa                                                                                                                >  1:179872/1‑122 (MQ=255)
                                                                ccGATGATCGAGGAAATTTGTGCTGCCTGTGCCGGTCCTTTCGCGAGTAAGCCCGCGCCTGCAGGGGTAGCGTACCCGTGCAGGTGCGGGCTTGCCCGCGAAGAGGGCGGTGCAGGCTAGaa                                                                                                                <  2:179872/122‑1 (MQ=255)
                                                                ccGATGATCGAGGAAATTTGTGCTGCCTGTGCCGGTCCTTTCGCGAGTAAGCCCGCGCCTGCAGGGGTAGCGTACCCGTGCAGGTGCGGGCTTGCCCGCGAAGAGGGCGGTGCAGGCTAGAACAGCa                                                                                                           <  2:570128/127‑1 (MQ=255)
                                                                ccGATGATCGAGGAAATTTGTGCTGCCTGTGCCGGTCCTTTCGCGAGTAAGCCCGCGCCTGCAGGGGTAGCGTACCCGTGCAGGTGCGGGCTTGCCCGCGAAGAGGGCGGTGCAGGCTAGAACAGCa                                                                                                           >  1:570128/1‑127 (MQ=255)
                                                                ccGATGATCGAGGAAATTTGTGCTGCCTGTGCCGGTCCTTTCGCGAGTAAGCCCGCGCCTGCAGGGGTAGCGTACCCGTGCAGGTGCGGGCTTGCCCGCGAAGAGGGCGGTGCAGGCTAGAACAGCACTTTCGCCACGTCGGCAAAgcgc                                                                                    >  2:559401/1‑150 (MQ=255)
                                                                ccGATGATCGAGGAAATTTGTGCTGCCTGTGCCGGTCATTTCGCGCGCACGAAAGCGCCTGCAGGGGTAGCGTACCCGTGCAGGTGCGGGCTTGCCCGCGAAGAGGGCGGTGCAGGCTAGAACAGCACTTTCGCCACGTCGGCAAAgcgc                                                                                    >  1:1104661/1‑150 (MQ=255)
                                                                  gatgatCGAGGAAATTTGTGCTGCCTGTGCCGGTCCTTTCGCGAGTAAGCCCGCGCCTGCAGGGGTAGCGTACAAGTGCAGGTGCGGGCTTGAGGGCTTCGAGGGCGGTGCAGGCTAGAACAGCACTTTCGCCACGTCGGCAAAGCGCtt                                                                                  >  2:387976/1‑150 (MQ=255)
                                                                     gatCGAGGAAATTTGTGCTGCCTGTGCCGGTCCTTTCGCGAGTAAGCCCGCGCCTGCAGGGGTAGCGTACCCGTGCAGGTGCGGCCTTGCCCGCGAAGAGGGCGGTGCAGGCTAGAACAGCACTTTCGCCACGTCGGCAAAGCGCTTGGc                                                                               <  2:300710/150‑1 (MQ=255)
                                                                       tCGAGGAAATTTGTGCTGCCTGTGCCGGTCCTTTCGCGCGTAAGCCCGCGCCTGCAGGGGTAGCGTACCCGTGCAGGTGCGGGCTTGCCCGCGAAGAGGGCGGTGCAGGCTAGAACAGCACTTTCGCCACGTCGGCAAAGCGCTTGGCGa                                                                             <  1:343058/150‑1 (MQ=255)
                                                                       tCGAGGAAATTTGTGCTGCCTGTGCCGGTCCTTTCGCGAGTAAGCCCGCGCCTGCAGGGGTAGCGTACCCGTGCAGGTGCGGGCTTGCCCGCGAAGAGGGCGGTGCAGGCTAGAACAGCACTTTCGCCACGTCGGCAAAGCGCTTGGCGa                                                                             >  1:1510420/1‑150 (MQ=255)
                                                                          aGGAAATTTGTGCTGCCTGTGCCGGTCCTTTCGCGAGTAAGCCCGCGCCTGCAGGGGTAGCGTACCCGTGCAGGTGCGGGCTTGCCCGCGAAGAGGGCGGTGCAGGCTAGAACAGCACTTTCGCCACGTCGGCAAAGCGCTTGGCGAAAt                                                                          >  1:695737/1‑150 (MQ=255)
                                                                          aGGAAATTTGTGCTGCCTGTGCCGGTCCTTTCGCGAGTAAGCCCGCGCCTGCAGGGGTAGCGTACCCGTGCAGGTGCGGGCTTGCCCGCGAAGAGGGCGGTGCAGGCTAGAACAGCACTTTCGCCACGTCGGCAAAGCGCTTGGCGAAAt                                                                          <  1:1204417/150‑1 (MQ=255)
                                                                                tttGTGCTGCCTGTGCCGGTCCTTTCGCGAGTAAGCCCGCGCCTGCAGGGGTAGCGTACCCGTGCAGGTGCGGGCTTGCCCGCGAAGAGGGCGGTGCAGGCTAGAACAGCACTTTCGCCACGTCGGCAAAGCGCTTGGCGAAATGCACTg                                                                    >  1:941560/1‑150 (MQ=255)
                                                                                    tgctgcCTGTGCCGGTCCTTTCGCGAGTAAGCCCGCGCCTGCAGGGGTAGCGTACCCGTGCAGGTGCGGGCTTGCCCGCGAAGAGGGCGGTGCAGGCTAGAACAGCACTTTCGCCACGTCGGCAAAGCGCTTg                                                                                 >  2:598672/1‑133 (MQ=255)
                                                                                    tgctgcCTGTGCCGGTCCTTTCGCGAGTAAGCCCGCGCCTGCAGGGGTAGCGTACCCGTGCAGGTGCGGGCTTGCCCGCGAAGAGGGCGGTGCAGGCTAGAACAGCACTTTCGCCACGTCGGCAAAGCGCTTg                                                                                 <  1:598672/133‑1 (MQ=255)
                                                                                    tgctgcCTGTGCCGGTCCTTTCGCGAGTAAGCCCGCGCCTGCAGGGGTAGCGTACCCGTGCAGGTGCGGGCTTGCCCGCGAAGAGGGCGGTGCAGGCTAGAACAGCACTTTCGCCACGTCGGCAAAGCGCTTGGCGAAATGCACTGTCAg                                                                >  2:105515/1‑150 (MQ=255)
                                                                                     gctgcCTGTGCCGGTCCTTTCGCGAGTAAGCCCGCGCCTGCAGGGGTAGCGTACCCGTGCAGGTGCGGGCTTGCCCGCGAAGAGGGCGGTGCAGGCTAGAACAGCACTTTCGCCACGTCGGCAAAGCGCTTGGCGAAATGCACTGTCAgg                                                               >  1:378013/1‑150 (MQ=255)
                                                                                      ctgcCTGTGCCGGTCCTTTCGCGAGTAAGCCCGCGCCTGCAGGGGTAGCGTACCCGTGCAGGTGCGGGCTTGCCCGCGAAGAGGGCg                                                                                                                             <  2:660841/87‑1 (MQ=255)
                                                                                      ctgcCTGTGCCGGTCCTTTCGCGAGTAAGCCCGCGCCTGCAGGGGTAGCGTACCCGTGCAGGTGCGGGCTTGCCCGCGAAGAGGGCg                                                                                                                             >  1:660841/1‑87 (MQ=255)
                                                                                      ctgcCTGTGCCGGTCCTTTCGCGAGTAAGCCCGCGCCTGCAGGGGTAGCGTACCCGTGCAGGTGCGGGCTTGCCCGCGAAGAGGGCGGTGCAGGCTAGAACAGCACTTTCGcc                                                                                                   <  1:1324304/113‑1 (MQ=255)
                                                                                      ctgcCTGTGCCGGTCCTTTCGCGAGTAAGCCCGCGCCTGCAGGGGTAGCGTACCCGTGCAGGTGCGGGCTTGCCCGCGAAGAGGGCGGTGCAGGCTAGAACAGCACTTTCGcc                                                                                                   >  2:1324304/1‑113 (MQ=255)
                                                                                       tgcCTGTGCCGGTCCTTTCGCGAGTAAGCCCGCGCCTGCAGGGGTAGCGTACCCGTGCAGGTGCGGGCTTGCCCGCGAAGAGGGCgg                                                                                                                            >  1:1320303/1‑87 (MQ=255)
                                                                                       tgcCTGTGCCGGTCCTTTCGCGAGTAAGCCCGCGCCTGCAGGGGTAGCGTACCCGTGCAGGTGCGGGCTTGCCCGCGAAGAGGGCgg                                                                                                                            <  2:1320303/87‑1 (MQ=255)
                                                                                               gatctCCTTTCGCGAGTAAGCCCGCGCCTGCAGGGGTAGCGTACCCGTGCAGGTGCGGGCTTGCCCGCGAAGAGGGCGGTGCAGGCTAGAACAGCACTTTCGCCACGTCGGCAAAGCGCTTGGCGAAATGCACTGTCAGGCCTTCGCGCa                                                     <  2:1265711/146‑1 (MQ=255)
                                                                                                cGGTCCTTTCGCGAGTAAGCCCGCGCCTGCAGGGGTAGCGTACCCGTGCAGGTGCGGGCTTGCCCGCGAAGAGGGCGGTGCAGGCTAGAACAGCACTTTCGCCACGTCGGCAAAGCGCTTGGCGAAATGCACTGTCAGGCCTTCGCGCAg                                                    >  2:364118/1‑150 (MQ=255)
                                                                                                cGGTCCTTTCGCGAGTAAGCCCGCGCCTGCAGGGGTAGCGTACCCGTGCAGGTGCGGGCTTGCCCGCGAAGAGGGCGGTGCAGGCTAGAACAGCACTTTCGCCACGTCGGCAAAGCGCTTGGCGAAATGCACTGTCAGGCCTTCGCGCAg                                                    >  2:1486052/1‑150 (MQ=255)
                                                                                                cGGTCCTTTCGCGAGTAAGCCCGCGCCTGCAGGGGTAGCGTACCCGTGCAGGTGCGGGCTTGCCCGCGAAGAGGGCGGTGCAGGCTAGAACAGCACTTTCGCCACGTCGGCAAAGCGCTTGGCGAAATGCACTGTCAGGCCTTCGCGCAg                                                    <  1:559401/150‑1 (MQ=255)
                                                                                                 ggTCCTTTCGCGAGTAAGCCCGCGCCTGCAGGGGTAGCGTACCCGTGCAGGTGCGGGCTTGCCCGCGAAGAGGGCGGTGCAGGCTAGAACAGCACTTTCGCCACGTCGGCAAAGCGCTTGGCGAAATGCACTGTCAGGCCTTCGCGCAgg                                                   <  2:1104661/150‑1 (MQ=255)
                                                                                                     cTTTCGCGAGTAAGCCCGCGCCTGCAGGGGTAGCGTACATGTGAAGGTGCGGGCTTGCCCGCGAAGAGGGCGGTGCAGGCTAGAACAGCACTTTCGCCACGTCGGCAAAGCGCTTGGCGAAATGCACTGTCAGGCCTTCGCGCAgagatc                                               >  1:1265711/1‑145 (MQ=255)
                                                                                                       ttCGCGAGTAAGCCCGCGCCTGCAGGGGTAGCGTACCCGTGCAGGTGCGGGCTTGCCCGCGAAGAGGGCGGTGCAGGCTAGAACa                                                                                                              >  1:1012332/1‑85 (MQ=255)
                                                                                                       ttCGCGAGTAAGCCCGCGCCTGCAGGGGTAGCGTACCCGTGCAGGTGCGGGCTTGCCCGCGAAGAGGGCGGTGCAGGCTAGAACa                                                                                                              <  2:1012332/85‑1 (MQ=255)
                                                                                                       ttCGCGAGTAAGCCCGCGCCTGCAGGGGTAGCGTACCCGTGCAGGTGCGGGCTTGCCCGCGAAGAGGGCGGTGCAGGCTAGAACAGCACTTTCGCCACGTCGGCAAAGCGCTTGGCGAAATGCACTGTCAGGCCTTCGAGAAGGTAGTcc                                             <  2:1510420/150‑1 (MQ=255)
                                                                                                            gAGTAAGCCCGCGCCTGCAGGGGTAGCGTACCCGTGCAGGTGCGGGCTTGCCCGCGAAGAGGGCGGTGCAGGCTAGAACAGCACTTTCGCCACGTCGGCAAAGCGCTTGGc                                                                               >  1:625484/1‑111 (MQ=255)
                                                                                                            gAGTAAGCCCGCGCCTGCAGGGGTAGCGTACCCGTGCAGGTGCGGGCTTGCCCGCGAAGAGGGCGGTGCAGGCTAGAACAGCACTTTCGCCACGTCGGCAAAGCGCTTGGc                                                                               <  2:625484/111‑1 (MQ=255)
                                                                                                               tAAGCCCGCGCCTGCAGGGGTAGCGTACCCGTGCAGGTGCGGGCTTGCCCGCGAAGAGGGCGGTGCAGGCTAGAACAGCACTTTCGCCACGTCGGCAAAGCGCTTGGCGAAATGCACTGTCAGGCCTTCGCGCAGGTAGTCCGGCAGttc                                     >  1:804923/1‑150 (MQ=255)
                                                                                                                  gCCCGCGCCTGCAGGGGTAGCGTACCCGTGCAGGTGCGGGCTTGCCCGCGAAGAGGGCGGTGCAGGCTAGAACAGCACTTTCGCCACGTCGGCAAAGCGCTTGGCGAAATGCACTGTCAGGCCTTCGCGCAGGTAGTCCGGCAGttcttc                                  <  1:235967/150‑1 (MQ=255)
                                                                                                                   cccGCGCCTGCAGGGGTAGCGTACCCGTGCAGGTGCGGGCTTGCCCGCGAAGAGGGCGGTGCAGGCTAGAACAGCACTTTCGCCACGTCAGCAAAGCGCTTGGCGAAATGCACTGTCAGGCCTTCGCGCAGGTAGTCCGGCAGTTCTTCg                                 >  2:575004/1‑150 (MQ=255)
                                                                                                                    ccGCGCCTGCAGGGGTAGCGTACCCGTGCAGGTGCGGGCTTGCCCGCGAAGAGGGCGGTGCAGGCTAGAACAGCACTTTCGCCACGTCGGCAAAGCGCTTGGCGAAATGCACTGTCAGGCCTTCGCGCCGGTAGTCCGGCAGTTCTTCGa                                >  2:495111/1‑150 (MQ=255)
                                                                                                                    ccGCGCCTGCAGGGGTAGCGTACCCGTGCAGGTGCGGGCTTGCCCGCGAAGAGGGCGGTGCAGGCTAGAACAGCACTTTCGCCACGTCGGCAAAGCGCTTGGCGAAATGCACTGTCAGGCCTTCGCGCAGGTAGTCCGGCAGTTCTTCGa                                >  2:495272/1‑150 (MQ=255)
                                                                                                                    ccGCACCTGCACGGGTACGCTACCCCTGCAGGCGCGGGCTTACTCGCGAAGAGGGCGGTGCAGGCTAGAACAGCACTTTCGCCACGTCGGCAAAGCGCTTGGCGAAATGCACTGt                                                                   >  2:727625/1‑115 (MQ=18)
                                                                                                                    ccGCACCTGCACGGGTACGCTACCCCTGCAGGCGCGGGCTTACTCGCGAAGAGGGCGGTGCAGGCTAGAACAGCACTTTCGCCACGTCGGCAAAGCGCTTGGCGAAATGCACTGt                                                                   <  1:727625/115‑1 (MQ=18)
                                                                                                                     cgcgCCTGCAGGGGTAGCGTACCCGTGCAGGTGCGGGCTTGCCCGCGAAGAGGGCGGTGCAGGCTAGAACAGCACTTTCGCCACGTCGGCAAAGCGCTTGGCGAAATGCACTGTCAGGCCTTCGCGCAGGTAGTCCGGCAGTTCTTCGaa                               <  2:804923/150‑1 (MQ=255)
                                                                                                                           tGCTGGTGTAGCGTACCCGTGCAGGTGCTGGCTTCCCCGCGAAGAGGGCGGTGCAGGCTAGAACAGCCCTTTCGCCACGTCGGCAAAGCGCTTGGCGAAATGCTCTGTCAGGCCTTCGCGCAGGTAGTCCGGCAGTTCTTCGAAGTcgcc                         <  2:657296/150‑1 (MQ=255)
                                                                                                                           tGCAGGGGTAGCGTACCCGTGCAGGTGCGGGCTTGCCCGCGAAGAGGGCGGTGCAGGCTAGAACAGCACTTTCGCCACGTCGGCAAAGCGCTTGGCGAAATGCACTGTCAGGCCTTCGCGCAGGTAGTCCGGCAGTTCTTCGAAGTcgcc                         >  2:745462/1‑150 (MQ=255)
                                                                                                                              aGGGGTAGCGTACCCGTGCAGGTGCGGGCTTGCCCGCGAAGAGGGCGGTGCAGGCTAGAACAGCACTTTCGCCACGTCGGCAAAGCGCTTGGCGAAATGCACTGTc                                                                  >  2:120715/1‑106 (MQ=255)
                                                                                                                              aGGGGTAGCGTACCCGTGCAGGTGCGGGCTTGCCCGCGAAGAGGGCGGTGCAGGCTAGAACAGCACTTTCGCCACGTCGGCAAAGCGCTTGGCGAAATGCACTGTc                                                                  <  1:120715/106‑1 (MQ=255)
                                                                                                                                 ggTAGCGTACCCGTGCAGGTGCGGGCTTGCCCGCGAAGAGGGCGGTGCAGGCTAGAACAGCACTTTCGCCACGTCGGCAAAGCGCTTGGCGAAATGCACTGTCAGGCCTTCGCGCAGGTAGTCCGGCAGTTCTTCGAAGTCGCCGCGGtt                   <  2:264445/150‑1 (MQ=255)
                                                                                                                                  gTAGCGTACCCGTGCAGGTGCGGGCTTGCCCGCGAAGAGGGCGGTGCAGGCTAGAACAGCACTTTCGCCACGTCGGCAACGCGCTTGGCGAAATGCACTGTCAGGCCTTCGCGCAGGTAGTCCGGCAGTTCTTCGAAGTCGCCGCGGTTg                  <  1:1266188/150‑1 (MQ=255)
                                                                                                                                  gTAGCGTACCCGTGCAGGTGCGGGCTTGCCCGCGAAGAGGGCGGTGCAGGCTAGAACAGCACTTTCGCCACGTCGGCAAAGCGCTTGGCGAAATGCACTGTCAGGCCTTCGCGCAGGTAGTCCGGCAGTTCTTCGAAGTCGCCGCGGTTg                  >  2:443827/1‑150 (MQ=255)
                                                                                                                                  gTAGCGTACCCGTGCAGGTGCGGGCTTGCCCGCGAAGAGGGCGGTGCAGGCTAGAACAGCACTTTCGCCACGTCGGCAAAGCGCTTGGCGAAATGCACTGTCAGGCCTTCGCGCAGGTAGTCCGGCAGTTCTTCGAAGTCGCCGCGGTTg                  <  2:1165704/150‑1 (MQ=255)
                                                                                                                                  gTAGCGTACCCGTGCAGGTGCGGGCTTGCCCGCGAAGAGGGCGGTGCAGGCTAGAACAGCACTTTCGCCACGTCGGCAAAGCGCTTGGCGAAATGCACTGTCAGGCCTTCGCGCAGGTAGTCCGGCAGTTCTTCGAAGTCGCCGCGGTTg                  <  1:105515/150‑1 (MQ=255)
                                                                                                                                    aGCGTACCCGTGCAGGTGCGGGCTTGCCCGCGAAGAGGGCGGTGCAGGCTAGAACAGCACTTTCGCCACGTCGGCAAAGCGCTTGGCGAAATGCACTGTCAGGCCTTCGCGCAGGTAGTCCGGCAGTTCTTCGAAGTCGCCGCGGTTggg                <  1:575004/150‑1 (MQ=255)
                                                                                                                                     gCGTACCCGTGCAGGTGCGGGCTTGCCCGCGAAGAGGGCGGTGCAGGCTAGAACAGCACTTTCGCCACGTCGGCAAAGCGCTTGGCGAAATGCACTGTCAGGCCTTCGCGCAGGTAGTCCGGCAGTTCTTCGAAGTCGCCGCGGTTGGGc               <  1:1502240/150‑1 (MQ=255)
                                                                                                                                      cGTACCCGTGCAGGTGCGGGCTTGCCCGCGAAGAGGGCGGTGCAGGCTAGAACAGCACTTTCGCCACGTCGGCAAAGCGCTTGGCGAAATGCACTGTCAGGCCTTCGCGCAGGTAGTCCGGCAGTTCTTCGaa                               <  2:1488197/133‑1 (MQ=255)
                                                                                                                                      cGTACCCGTGCAGGTGCGGGCTTGCCCGCGAAGAGGGCGGTGCAGGCTAGAACAGCACTTTCGCCACGTCGGCAAAGCGCTTGGCGAAATGCACTGTCAGGCCTTCGCGCAGGTAGTCAGGCAGTTCTTCGaa                               >  1:1488197/1‑133 (MQ=255)
                                                                                                                                       gTACCCGTGCAGGTGCGGGCTTGCCCGCGAAGAGGGCGGTGCAGGCTAGAACAGCACTTTCGCCACGTCGGCAAAGCGCTTGGCGAAATGCACTGTCAGGCCTTCGCGCAGGTAGTCCGGCAGttctt                                   <  1:103379/128‑1 (MQ=255)
                                                                                                                                       gTACCCGTGCAGGTGCGGGCTTGCCCGCGAAGAGGGCGGTGCAGGCTAGAACAGCACTTTCGCCACGTCGGCAAAGCGCTTGGCGAAATGCACTGTCAGGCCTTCGCGCAGGTAGTCCGGCAGttctt                                   >  2:103379/1‑128 (MQ=255)
                                                                                                                                       gTACCCGTGCAGGTGCGGGCTTGCCCGCGAAGAGGGCGGTGCAGGCTAGAACAGCACTTTCGCCACGTCGGCAAAGCGCTTGGCGAAATGCACTGTCAGGCCTTCGCGCAGGTAGTCCGGCAGTTCTTCGAAGTCGCCGCGGTTGGGCTc             <  2:695737/150‑1 (MQ=255)
                                                                                                                                         tACCCGTGCAGGTGCGGGCTTGCCCGCGAAGAGGGCGGTGCAGGCTAGAACAGCACTTTCGCCACGTCGGCAAAGCGCTTGGCGAAATGCACTGTCAGGCCTTCGCGCAGGTAGTCCGGCAGTTCTTCGAAGTCGCCGCGGTTGGGCTcc            >  2:1604950/1‑150 (MQ=255)
                                                                                                                                         tACCCGTGCAGGTGCGGGCTTGCCCGCGAAGAGGGCGGTGCAGGCTAGAACAGCACTTTCGCCACGTCGGCAAAGCGCTTGGCGAAATGCACTGTCAGGCCTTCGCGCAGGTAGTCCGGCAGTTCTTCGAAGTCGCCGCGGTTGGGCTcc            >  2:1406977/1‑150 (MQ=255)
                                                                                                                                         tACCCGTGCAGGTGCGGGCTTGCCCGCGAAGAGGGCGGTGCAGGCTAGAACAGCACTTTCGCCACGTCGGCAAAGCGCTTGGCGAAATGAACTGTCAGGCCTTCGCGCAGGTAATCCGGCAGTTCTTCGAAGTCGCCGCCGTTGGGCTcc            <  2:1089945/150‑1 (MQ=255)
                                                                                                                                          aCCCGTGCAGGTGCGGGCTTGCCCGCGAAGAGGGCGGTGCAGGCTAGAACAGCACTTTCGCCACGTCGGCAAAGCGCTTGGCGAAATGCACTGTCAGGCCTTCGCGCAg                                                    >  2:1613417/1‑109 (MQ=255)
                                                                                                                                          aCCCGTGCAGGTGCGGGCTTGCCCGCGAAGAGGGCGGTGCAGGCTAGAACAGCACTTTCGCCACGTCGGCAAAGCGCTTGGCGAAATGCACTGTCAGGCCTTCGCGCAg                                                    <  1:1613417/109‑1 (MQ=255)
                                                                                                                                            ccGTGCAGGTGCGGGCTTGCCCGCGAAGAGGGCGGTGCAGGCTAGAACAGCACTTTCGCCACGTCGGCAAAGCGCTTGGCGAAATGCATTGTCAGGCCTTCGCGCAGGTAGTCCGGCAGTTCTTCGAAGTCGCCGCGGTTGGGCTCCGGc         <  1:1248198/150‑1 (MQ=255)
                                                                                                                                               tGCAGGTGCGGGCTTGCCCGCGAAGAGGGCGGTGCAGGCTAGAACAGCACTTTCGCCACGTCGGCAAAGCGCTTGGCGAAATGCACTGTCa                                                                 <  1:905052/91‑1 (MQ=255)
                                                                                                                                               tGCAGGTGCGGGCTTGCCCGCGAAGAGGGCGGTGCAGGCTAGAACAGCACTTTCGCCACGTCGGCAAAGCGCTTGGCGAAATGCACTGTCa                                                                 >  2:905052/1‑91 (MQ=255)
                                                                                                                                                 cAGGTGCGGGCTTGCCCGCGAAGAGGGCGGTGCAGGCTAGAACAGCACTTTCGCCACGTCGGCAAAGCGCTTGGCGAAATGCACTGTCAGGCCTTCGCGCAGGTAGTCCGGCAGTTCTTCGAAGTCGCCGCGGTTGGGCTCCGGCAGGAt    <  2:665519/150‑1 (MQ=255)
                                                                                                                                                 cAGGTGCGGGCTTGCCCGCGAAGAGGGCGGTGCAGGCTAGAACAGCACTTTCGCCACGTCGGCAAAGCGCTTGGCGAAATGCACTGTCAGGCCTTCGCGCAGGTAGTCCGGCAGTTCTTCGAAGTCGCCGCGGTTGGGCTCCGGCAGGAt    <  1:342216/150‑1 (MQ=255)
                                                                                                                                                  aGGTGCGGGCTTGCCCGCGAAGAGGGCGGTCAGGCTAGAACAGCACTTTCGCCACGTCGGCAAAGCGCTTGGCGAAATGCACTGTCAGGCCTTCGCGCAGGTAGTCCGGCAGTTCTTCGAAGTCGCCGCGGTTGGGCTCCGGCAGGATCa  >  1:657249/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                   
AGGGGAACGAGCAGACGAGGAGAGGTCATGGTGGCTCCGAAGGCTGGCGATGGCCTGAGGATACCCGATGATCGAGGAAATTTGTGCTGCCTGTGCCGGTCCTTTCGCGAGTAAGCCCGCGCCTGCAGGGGTAGCGTACCCGTGCAGGTGCGGGCTTGCCCGCGAAGAGGGCGGTGCAGGCTAGAACAGCACTTTCGCCACGTCGGCAAAGCGCTTGGCGAAATGCACTGTCAGGCCTTCGCGCAGGTAGTCCGGCAGTTCTTCGAAGTCGCCGCGGTTGGGCTCCGGCAGGATCA  >  NC_002947/1642186‑1642481

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