Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,642,306 G→A 5.9% intergenic (‑92/+60) 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,3060GA 5.9% 351.2 / 4.6 121intergenic (‑92/+60)PP_1442/lon‑Ihypothetical protein/DNA‑binding ATP‑dependent protease
Reads supporting (aligned to +/- strand):  ref base G (60/52);  new base A (6/2);  total (67/54)
Fisher's exact test for biased strand distribution p-value = 2.93e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 127 total aligned reads displayed.

GCGCAAGCAGAAAGCAGGGCGAAGCTCAGGGGAACGAGCAGACGAGGAGAGGTCATGGTGGCTCCGAAGGCTGGCGATGGCCTGAGGATACCCGATGATCGAGGAAATTTGTGCTGCCTGTGCCGGTCCTTTCGCGAGTAAGCCCG‑‑CGCCTGCAGGGGTAGCGTACCCGTGCAGGTGCGGGCTTGCCCGCGAAGAGGGCGGTGCAGGCTAGAACAGCACTTTCGCCACGTCGGCAAAGCGCTTGGCGAAATGCACTGTCAGGCCTTCGCGCAGGTAGTCCGGCAGTTCTTCGAAGTC  >  NC_002947/1642159‑1642455
                                                                                                                                                     |                                                                                                                                                      
gcgcAAGCAGAAAGCAGGGCGAAGCTCAGGGGAACGAGCAGACGAGGAGAGGTCATGGTGGCTCCGAAGGCTGGCGATGGCCTGAGGATACCCGATGATCGAGGAAATTTGTGCTGCCTGTGCCGGTCCTTTCGCGAGTAAGCCCG‑‑CGcc                                                                                                                                                      >  1:1461064/1‑150 (MQ=255)
    aaGCAGAAAGCAGGGCGAAGCTCAGGGGAACGAGCAGACGAGGAGAGGTCATGGTGGCTCCGAAGGCTGGCGATGGCCTGAGGATACCCGATGATCGAGGAAATTTGTGCTGCCTGTGCCGGTCCTTTCGCGAGTAAGCCCG‑‑CACCTGCa                                                                                                                                                  >  1:2132613/1‑150 (MQ=255)
     aGCAGAAAGCAGGGCGAAGCGCAGGGGAACGAGCAGACGAGGAGAGGTCATGGTGGCTCCGAAGGCTGGCGATGGCCTGAGGATACCCGATGATCGAGGAAATTTGTGCTGCCTGTGCCGGTCCTTTCGCGAGTAAGCCCG‑‑CGCCTGCAg                                                                                                                                                 <  1:1602357/150‑1 (MQ=255)
      gCAGAAAGCAGGGCGAAGCTCAGGGGAACGAGCAGACGAGGAGAGGTCATGGTGGCTCCGAAGGCTGGCGATGGCCTGAGGATACCCGATGATCGAGGAAATTTGTGCTGCCTGTGCCGGTCCTTTCGCGAGTAAGCCCG‑‑CGCCTGCAgg                                                                                                                                                >  2:742145/1‑150 (MQ=255)
        aGAAAGCAGGGCGAAGCTCAGGGGAACGAGCAGACGAGGAGAGGTCATGGTGGCTCCGAAGGCTGGCGATGGCCTGAGGATACCCGATGATCGAGGAAATTTGTGCTGCCTGTGCCGGTCCTTTCGCGAGTAAGCCCG‑‑CACCTGCAcggg                                                                                                                                              >  2:1714498/1‑150 (MQ=255)
          aaaGCAGGGCGAAGCTCAGGGGAACGAGCAGACGAGGAGAGGTCATGGTGGCTCCGAAGGCTGGCGATGGCCTGAGGATACCCGATGATCGAGGAAATTTGTGCTGCCTGTGCCGGTCCTTTCGCGAGTAAGCCCG‑‑CGCCTGCAGGGGTa                                                                                                                                            <  2:914211/150‑1 (MQ=255)
            aGCAGGGCGAAGCTCAGGGGAACGGGCACACGGGGAGATGTCACGGTGGCTCCGCCGGCATCGGAGGGTTATGTGATACCCGATGATCGAGGAAATTTGTGCTGCCTGTGCCGGTCCTTTCGCGAGTAAGCCCG‑‑CGCCTGCAGGGGTAGc                                                                                                                                          <  1:1224490/150‑1 (MQ=255)
                       gCTCAGGGGAACGAGCAGACGAGGAGAGGTCATGGTGGCTCCGAAGGCTGGCGATGGCCTGAGGATACCCGATGATCGAGGAAATTTGTGCTGCCTGTGCTGTTCCTGTCGCGACTAAGCCCG‑‑CGCCTGCAGGGGTAGCGTACCCGTGCa                                                                                                                              >  2:663425/1‑150 (MQ=255)
                         tCAGGGGAACGAGCAGACGAGGAGAGGTCATGGTGGCTCGGAAGGCTGGCGATGGCCTGAGGATACCAGATGATAGAGGAAATTTGTGCTGCCTGTGCCGGTCCTTTCGCGAGTAAGCACG‑‑CGCCTGGAGGGGTAGCGTAACCGTGCAgg                                                                                                                            >  1:909862/1‑150 (MQ=255)
                          cAGGGGAACGAGCAGACGAGGAGAGGTCATGGTGGCTCCGAAGGCTGGCGATGGCCTGAGGATACCCGATGATCGAGGAAATTTGTGCTGCCTGTGCCGGTCCTTTCGCGAGTAAGCCCG‑‑CGCCTGCAGGGGTAGCGTACCCGTGCAGGt                                                                                                                           <  1:655412/150‑1 (MQ=255)
                                aaCGAGCAGACGAGGAGAGGTCGTGGTGGCTCCGAAGGCTGGCGATGGCCTGAGGATACCCGATGATCGAGGAAATTTGTGCTGCCTGTGCCGGTCCTTTCGCGAGTAAGCCCG‑‑CGCCTGCAGGGGTAGCGTACCCGTGCAGGTGCGGGc                                                                                                                     <  2:1198033/150‑1 (MQ=255)
                                     gCAGACGAGGAGAGGTCATGGTGGCTCCGAAGGCTGGCGATGGCCTGAGGATACCCGATGATCGAGGAAATTTGTGCTGCCTGTGCCGGTCCTTTCGCGAGTAAGCCCG‑‑CACCTGCACGGGTACGCTACCCCTGCAGGCGCGGGCttact                                                                                                                <  2:2132613/150‑4 (MQ=37)
                                        gACGAGGAGAGGTCATGGTGGCTCCGAAGGCTGGCGATGGCCTGAGGATACCCGATGATCGAGGAAATTTGTGCTGCCTGTGCCGGTCCTTTCGCGAGTAAGCCCG‑‑CGCCTGCAGGGGTAGCGTACCCGTGCAGGTGCGGGCTTGCCcgc                                                                                                             >  2:1771799/1‑150 (MQ=255)
                                           gaggagAGGTCATGGTGGCTCCGAAGGCTGGCGATGGCCTGAGGATACCCGATGATCGAGGAAATTTGTGCTGCCTGTGCCGGTCCTTTCGCGAGTAAGCCCG‑‑CGCCTGCAGGGGTAGCGTACCCGTGCAGGTGCGGGCTTGCCCGCGaa                                                                                                          <  1:628349/150‑1 (MQ=255)
                                           gaggagAGGTCATGGTGGCTCCGAAGGCTGGCGATGGCCTGAGGATACCCGATGATCGAGGAAATTTGTGCTGCCTGTGCCGGTCCTTTCGCGAGTAAGCCCG‑‑CGCCTGCAGGGGTAGCGTACCCGTGCAGGTGCGGGCTTGCCCGCGaa                                                                                                          <  1:1771799/150‑1 (MQ=255)
                                            aggagAGGTCATGGTGGCTCCGAAGGCTGGCGATGGCCTGAGGATACCCGATGATCGAGGAAATTTGTGCTGCCTGTGCCGGTCCTTTCGCGAGTAAGCCCG‑‑CGCCTGCAGGGGTAGCGTACCCGTGCAGGTGCGGGCTTGCCCGCGAag                                                                                                         >  2:2315112/1‑150 (MQ=255)
                                                gagGTCATGGTGGCTCCGAAGGCTGGCGATGGCCTGAGGATACCCGATGATCGAGGAAATTTGTGCTGCCTGTGCCGGTCCTTTCGCGAGTAAGCCCG‑‑CGCCTGCAGGGGTAGCGTACCCGTGCAGGTGCGGGCTTGCCCGCGAAGAggg                                                                                                     <  1:201583/150‑1 (MQ=255)
                                                gagGTCATGGTGGCTCCGAAGGCTGGCGATGGCCTGAGGATACCCGATGATCGAGGAAATTTGTGCTGCCTGTGCCGGTCCTTTCGCGAGTAAGCCCG‑‑CGCCTGCAGGGGTAGCGTACCCGTGCAGGTGCGGGCTTGCCCGCGAAGAggg                                                                                                     >  2:2301628/1‑150 (MQ=255)
                                                 agGTCATGGTGGCTCCGAAGGCTGGCGATGGCCTGAGGATACCCGATGATCGAGGAAATTTGTGCTGCCTGTGCCGGTCCTTTCGCGAGTAAGCCCG‑‑CGCCTGCAGGGGTAGCGTACCCGTGCAGGTGCGGGCTTGCCCGCGAAGAGGGc                                                                                                    >  1:1449027/1‑150 (MQ=255)
                                                  ggTCATGGTGGCTCCGAAGGCTGGCGATGGCATGAGGATACCCGATGATCAAGGAAATTTGTGCTGCCTGTGCCGGTCCTTTCGCGAGTAAGCCCG‑‑CGCCTGCAGGGGTAGCGTACCCGTGTAGGTGCGGGCTTGCCCGCGAAGAGGGCg                                                                                                   <  2:1250098/150‑1 (MQ=255)
                                                   gTCATGGTGGCTCCGAAGGCTGGCGATGGCCTGAGGATACCCGATGATCGAGGAAATTTGTGCTGCCTGTGCCGGTCCTTTCGCGAGTAAGCCCG‑‑CGCCTGCAGGGGTAGCGTACCCGTGCAGGTGCGGGCTTGCCCGCGAAGAGGGCgg                                                                                                  >  2:941910/1‑150 (MQ=255)
                                                   gTCATGGTGGCTCCGAAGGCTGGCGATGGCCTGAGGATACCCGATGATCGAGGAAATTTGTGCTGCCTGTGCCGGTCCTTTCGCGAGTAAGCCCG‑‑CGCCTGCAGGGGTAGCGTACCCGTGCAGGTGCCGGCTTGCCCGCGAAGAGGGCgg                                                                                                  <  1:742145/150‑1 (MQ=255)
                                                    tCATGGTGGCTCCGAAGGCTGGCGATGGCCTGAGGATACCCGATGATCGAGGAAATTTGTGCTGCCTGTGCCGGTCCTTTCGCGAGTAAGCCCG‑‑CGCCTGCAGGGGTAGCGTACCCGTGCAGGTGCGGGCTTGCCCGCGAAGAGGGCGGt                                                                                                 <  1:2106202/150‑1 (MQ=255)
                                                        ggtggCTCCGAAGGCTGGCGATGGCCTGAGGATACCCGATGATCGAGGAAATTTGTGCTGCCTGTGCCGGTCCTTTCGCGAGTAAGCCCG‑‑CGCCTGCAGGGGTAGCGTACCCGTGCAGGTGCGGGCTTGCCCGCGAAGAGGGCGGTGcgt                                                                                             <  1:2301628/150‑3 (MQ=255)
                                                        ggtggCTCCGAAGGCTGGCGATGGCCTGAGGATACCCGATGATCGAGGAAATTTGTGCTGCCTGTGCCGGTCCTTTCGCGAGTAAGCCCG‑‑CGCCTGCAGGGGTAGCGTACCCGTGCAGGTGCGGGCTTGCCCCCGAAGCCAGCGCTGCAg                                                                                             <  2:1602585/150‑1 (MQ=255)
                                                         gtggCTCCGAAGGCTGGCGATGGCCTGAGGATACCCGATGATCGAGGAAATTTGTGCTGCCTGTGCCGGTCCTTTCGCGAGTAAGCCCG‑‑CGCCTGCAGGGGTAGCGTACCCGTGCAGGTGCGGGCTTGCCCGCGAAGAGGGCGGTGCAgg                                                                                            >  1:2358548/1‑150 (MQ=255)
                                                         gtggCTCCGAAGGCTGGCGATGGCCTGAGGATACCCGATGATCGAGGAAATTTGTGCTGCCTGTGCCGGTCCTTTCGCGAGTAAGCCCG‑‑CGCCTGCAGGGGTAGCGTACCCGTGCAGGTGCGGGCTTGCCCGCGAAGAGGGCGGTGCAgg                                                                                            <  2:1449027/150‑1 (MQ=255)
                                                         gtggCTCCGAAGGCTGGCGATGGCCTGAGGATACCCGATGATCGAGGAAATTTGTGCTGCCTGTGCCGGTCCTTTCGCGAGTAAGCCCG‑‑CGCCTGCAGGGGTAGCGTACCCGTGCAGGTGCGGGCTTGCCCGCGAAGAGGGCGGTGCAgg                                                                                            >  1:477740/1‑150 (MQ=255)
                                                          tggCTCCGAAGGCTGGCGATGGCCTGAGGATACCCGATGATCGAGGAAATTTGTGCTGCCTGTGCCGGTCCTTTCGCGAGTAAGCCCG‑‑CGCCTGCAGGGGTAGCGTACCCGTGCa                                                                                                                              <  2:899793/115‑1 (MQ=255)
                                                          tggCTCCGAAGGCTGGCGATGGCCTGAGGATACCCGATGATCGAGGAAATTTGTGCTGCCTGTGCCGGTCCTTTCGCGAGTAAGCCCG‑‑CGCCTGCAGGGGTAGCGTACCCGTGCa                                                                                                                              >  1:899793/1‑115 (MQ=255)
                                                          tggCTCCGAAGGCTGGCGATGGCCTGAGGATACCCGATGATCGAGGAAATTTGTGCTGCCTGTGCCGGTCCTTTCGCGAGTAAGCCCG‑‑CGCCTGCAGGGGTAGCGTACCCGTGCAGGTGCGGGCTTGCCCGCGAAGAGGGCGATGCAGGc                                                                                           >  2:1224232/1‑150 (MQ=255)
                                                            gCTCCGAAGGCTGGCGATGGCCTGAGGATACCCGATGATCGAGGAAATTTGTGCTGCCTGTGCCGGTCCTTTCGCGAGTAAGCCCG‑‑CGCCTGCAGGGGTAGCGTACCCGTGCAGGTGCGGGCTTGCCCGCGAAGAGGGCGGTGCAGGCTa                                                                                         <  1:1711910/150‑1 (MQ=255)
                                                              tCCGAAGGCTGGCGATGGCCTGAGGATACGAGCTGACCCTTGAAAGTTGTGCTGCCTGTGACGGTCCTTTAGCACACGTCTGAACTCCAGTCACGATAGCGGACCCGTGCAGGTGCGGGCTTGCCCGCGAAGAGGCCGGGGCAGGCTAGa                                                                                       >  1:1126167/1‑150 (MQ=255)
                                                              tCCGAAGGCTGGCGATGGCCTGAGGATACCCGATGATCGAGGAAATTTGTGCTGCCTGTGCCGGTCCTTTCGCGAGTAAGCCCG‑‑CGCCTGCAGGGGTAGCGTACCCGTGCAGGTGCGGGCTTGCCCGCGACGAGGGCGGTGCAGGCTAGa                                                                                       >  1:119845/1‑150 (MQ=255)
                                                              tCCGAAGGCTGGCGATGGCCTGAGGATACCCGATGATCGAGGAAATTTGTGCTGCCTGTGCCGGTCCTTTCGCGAGTAAGCCCG‑‑CGCCTGCAGGGGTAGCGTACCCGTGCAGGTGCGGGCTTGCCCGCGAAGAGGGCGGTGCAGGCTAGa                                                                                       >  1:166137/1‑150 (MQ=255)
                                                               ccGAAGGCTGGCGATGGCCTGAGGATACCCGATGATCGAGGAAATTTGTGCTGCCTGTGCCGGTCCTTTCGCGAGTAAGCCCG‑‑CGCCTGCAGGGGTAGCGTACCCGTGCAGGTGCGGGCTTGCCCGCGAAGAGGGCGGTGCAGGCTAGaa                                                                                      <  1:941910/150‑1 (MQ=255)
                                                                 gAAGGCTGGCGATGGCCTGAGGATACCCGATGATCGAGGAGATTTGTGCTGCCTGTGCCGGTCCTTTCGCGAGTAAGCCCG‑‑CACCTGCACGGGTACGCTACCCCTGCAGGCGCgg                                                                                                                       >  2:839109/1‑115 (MQ=25)
                                                                 gAAGGCTGGCGATGGCCTGAGGATACCCGATGATCGAGGAGATTTGTGCTGCCTGTGCCGGTCCTTTCGCGAGTAAGCCCG‑‑CACCTGCACGGGTACGCTACCCCTGCAGGCGCgg                                                                                                                       <  1:839109/115‑1 (MQ=25)
                                                                 gAAGGCTGGCGATGGCCTGAGGATACCCGATGATCGAGGAAATTTGTGCTGCCTGTGCCGGTCCTTTCGCGAGTAAGCCCG‑‑CGCCTGCAGGGGTAGCGTACCCGTGCAGGTGCGGGCTTGCCCGCGAAGAGGGCGGTGCAGGCTAGAACa                                                                                    <  2:1461064/150‑1 (MQ=255)
                                                                       tGGCGATGGCCTGAGGATACCCGATGATCGAGGAAATTTGTGCTGCCTGTGCCGGTCCTTTCGCGAGTAAGCCCG‑‑CGCCTGCAGGGGTAGCGTACCCGTGCAGGTGCGGGCTTGCCCGCGAAGAGGGCGGTGCAGGCTAGAACAGCACtt                                                                              >  2:1413578/1‑150 (MQ=255)
                                                                       tGGCGATGGCCTGAGGATACCCGATGATCGAGGAAATTTGTGCTGCCTGTGCCGGTCCTTTCGCGAGTAAGCCCG‑‑CGCCTGCAGGGGTAGCGTACCCGTGCAGGTGCGGGCTTGCCCGCGAAGAGGGCGGTGCAGGCTAGAACAGCACtt                                                                              >  2:49837/1‑150 (MQ=255)
                                                                           gATGGCCTGAGGATACCCGATGATCGAGGAAATTTGTGCTGCCTGTGCCGGTCCTTTCGCGAGTAAGCCCG‑‑CGCCTGCAGGGGTAGCGTACCCGTGCAGGTGCGGGCTTGCCCGCGAAGAGGGCGGTGCAGGCTAGAACAGCACTTTCGc                                                                          >  2:504292/1‑150 (MQ=255)
                                                                                 ccgcgcaTACCCGATGACAGAGGCCCTCTGTGCACCTTGTGCCGGTCCTCTCGCGAGTAAGCCCG‑‑CGCCTGCAGGGGTAGCGTACCAGGGCAGGTGCGGGCTTGCCCGCGAAGAGGGCGGTGCAGTCTAGAACAGCACTTTCGCCACGTc                                                                    >  2:1207428/7‑150 (MQ=255)
                                                                                 cTGAGGATACCCGATGATCGAGGAAATTTGTGCTGCCTGTGCCGGTCCTTTCGCGAGTAAGCCCG‑‑CGCCTGCAGGGGAAAACGCGCCGTGCAGGTGCGGGCTTGCCCGCGAAGAGGGCGGTGCAGGCTAGAACAGCACTTTCGCCACGTc                                                                    >  1:1970932/1‑150 (MQ=255)
                                                                                   gAGGATACCCGATGATCGAGGAAATTTGTGCTGCCTGTGCCGGTCCTTTCGCGAGTAAGCCCG‑‑CGCCTGccg                                                                                                                                                 <  2:482227/72‑3 (MQ=255)
                                                                                   gAGGATACCCGATGATCGAGGAAATTTGTGCTGCCTGTGCCGGTCCTTTCGCGAGTAAGCCCG‑‑CGCCTGCAg                                                                                                                                                 >  1:482227/1‑72 (MQ=255)
                                                                                         aCCCGATGACAGAGGAAAATTGTGCTGCCTGTGCCGGTCCTTTCGCGAGTAAGCCCG‑‑CGCCTGCAGGGGTAGCGTACCCGTGCAGGTGCGGGCTTGCCCGCGAAGAGGGCGGTGCAGGCTAGAACAGCACTTTCGCCACGTCGGCAAAgc                                                            >  2:1919849/1‑150 (MQ=255)
                                                                                           ccGATGATCGAGGAAATTTGTGCTGCCTGTGCCGGTCCTTTCGCGAGTAAGCCCG‑‑CGCCTGCAGGGGTAGCGTACCCGTGCAGGTGCGGGCTTGCCCGCGAAGAGGGCGGTGCAGGCTAGAACAGCACTTTCGCCACGTCGGCAAAgcgc                                                          >  1:2327994/1‑150 (MQ=255)
                                                                                           ccGATGATCGAGGAAATTTGTGCTGCCTGTGCCGGTCCTTTCGCGAGTAAGCCCG‑‑CGCCTGCAGGGGTAGCGTACCCGTGCAGGTGCGGGCTTGCCCGCGAAGAGGGCGGTGCAGGCTAGAACAGCACTTTCGCCACGTCGGCAAAgcgc                                                          >  1:1064267/1‑150 (MQ=255)
                                                                                           ccGATGATCGAGGAAATTTGTGCTGCCTGTGCCGGTCCTTTCGCGAGTAAGCCCG‑‑CGCCTGCAGGGGTAGCGTACCCGTGCAGGTGCGGGCTTGCCCGCGAAGAGGGCGGTGCAGGCTAGAACAGCACTTTCGCCACGTCGGCAAAgcgc                                                          >  1:454340/1‑150 (MQ=255)
                                                                                           ccGATGATCGAGGAAATTTGTGCTGCCTGTGCCGGTCCTTTCGCGAGTAAGCCCG‑‑CGCCTGCAGGGGTAGCGTACCCGTGCAGGTGCGGGCTTGCCCGCGAAGAGGGCGGTGCAGGCTAGAACAGCACTTTCGCCACGTCGGCAAAgcgc                                                          >  2:1675230/1‑150 (MQ=255)
                                                                                            cGATGATCGAGGAAATTTGTGCTGCCTGTGCCGGTCCTTTCGCGAGTAAGCCCG‑‑CGCCTGCAGGGGTAGCGTACCCGTGCAGGTGCGGGCTTGCCCGCGAAGAGGGCGGTGCAGGc                                                                                           >  2:2287160/1‑116 (MQ=255)
                                                                                            cGATGATCGAGGAAATTTGTGCTGCCTGTGCCGGTCCTTTCGCGAGTAAGCCCG‑‑CGCCTGCAGGGGTAGCGTACCCGTGCAGGTGCGGGCTTGCCCGCGAAGAGGGCGGTGCAGGc                                                                                           <  1:2287160/116‑1 (MQ=255)
                                                                                             gatgatCGAGGAAATTTGTGCTGCCTGTGCCGGTCCTTTCGCGAGTAAGCCCG‑‑CGCCTGCAGGGGTAGCGTACCCGTGCAGGTGCGGGCTTGCCCGCGAAGAGGGCGGTGCAGGCTAGAACAGCACTTTCGCCACGTCGGCAAAGCGCtt                                                        >  2:48620/1‑150 (MQ=255)
                                                                                             gatgatCGAGGAAATTTGTGCTGCCTGTGCCGGTCCTTTCGCGAGTAAGCCCG‑‑CGCCTGCAGGGGTAGCGTACCCGTGCAGGTGCGGGCTTGCCCGCGAAGAGGGCGGTGCAGGCTAGAACAGCACTTTCGCCACGTCGGCAAAGCGCtt                                                        >  1:482623/1‑150 (MQ=255)
                                                                                              atgatCGAGGAAATTTGTGCTGCCTGTGCCGGTCCTTTCGCGAGTAAGCCCG‑‑CGCCTGCAGGGGTAGCGTACCCGTGCAGGTGCGGGCTTGCCCGCGAAGAGGGCGGTGCAGGCTAGAACAGCACTTTCGCCACGTCGGCAAAGCGCTTg                                                       <  2:1118426/150‑1 (MQ=255)
                                                                                              atgatCGAGGAAATTTGTGCTGCCTGTGCCGGTCCTTTCGCGAGTAAGCCCG‑‑CGCCTGCAGGGGTAGCGTACCCGTGCAGGTGCGGGCTTGCCCGCGAAGAGGGCGGTGCAGGCTAGAACAGCACTTTCGCCACGTCGGCAAAGCGCTTg                                                       <  1:325617/150‑1 (MQ=255)
                                                                                                gatCGAGGAAATTTGTGCTGCCTGTGCCGGTCCTTTCGCGAGTAAGCCCG‑‑CGCCTGCAGGGGTAGCGTACCCGTGCAGGTGCGGGCTTGCCCGCGAAGAGGGCGGTGCAGGCTAGAACAGCACTTTCGCCACGTCGGCAAAGCGCTTGGc                                                     <  2:454340/150‑1 (MQ=255)
                                                                                                gatCGAGGAAATTTGTGCTGCCTGTGCCGGTCCATTCGCGAGTAAGCCCG‑‑CGCCTGCAGGGGGAGCGTACGCGTGCAGGTGCGGGCTTGCCCGCGAAGAGGGCGGTGCAGGCTAGAACAGCACTTTCGCCACGTCGGCAAAGCGCTTGGc                                                     >  1:649303/1‑150 (MQ=255)
                                                                                                  tCGAGGAAATTTGTGCTGCCTGTGCCGGTCCTTTCGCGAGTAAGCCCG‑‑CGCCTGCAGGGGTAGCGTACCCGTGCAGGTGCGGGCTTGCCCGCGAAGAGGGCGGTGCAGGCTAGAACAGCACTTTCGCCACGTCGGCAAAACGCTTGGCgc                                                   >  1:2202984/1‑149 (MQ=255)
                                                                                                   cGAGGAAATTTGTGCTGCCTGTGCCGGTCCTTTCGCGAGTAAGCCCG‑‑CGCCTGCAGGGGTAGCGTACCCGTGCAGGTGCGGGCTTGCCCGCGAAGAGGGCGGTGCAGGCTAGAACAGCACTTTCGCCACGTCGGCAAAGCGCTTGGCGaa                                                  <  2:2327994/150‑1 (MQ=255)
                                                                                                    gAGGAAATTTGTGCTGCCTGTGCCGGTCCTTTCGCGAGTAAGCCCG‑‑CGCCTGCAGGGGTAGCGTACCCGTGCAGGTGCGGGCTTGCCCGCGAAGAGGGCGGTGCAGGCTAGAACAGCACTTTCGCCACGTCGGCAAAGCGCTTGGCGaaa                                                 <  2:1126167/150‑1 (MQ=255)
                                                                                                    gAGGAAATTTGTGCTGCCTGTGCCGGTCCTTTCGCGAGTAAGCCCG‑‑CGCCTGCAGGGGTAGCGTACCCGTGCAGGTGCGGGCTTGCCCGCGAAGAGGGCGGTGCAGGCTAGAACAGCACTTTCGCCACGTCGGCAAAGCGCTTGGCGaaa                                                 <  2:482623/150‑1 (MQ=255)
                                                                                                    gAGGAAATTTGTGCTGCCTGTGCCGGTCCTTTCGCGAGTAAGCCCG‑‑CGCCTGCAGGGGTAGCGTACCCGTGCAGGTGCGGGCTTGCCCGCGAAGAGGGCGGTGCAGGCTAGAACAGCACTTTCGCCACGTCGGCAAAGCGCTTGGCGaaa                                                 >  2:611895/1‑150 (MQ=255)
                                                                                                     aGGAAATTTGTGCTGCCTGTGCCGGTCCTTTCGCGAGTAAGCCCG‑‑CGCCTGCAGGGGTAGCGTACCCGTGCAGGTGCGGGCTTGCCCGCGAAGAGGGCGGTGCAGGCTa                                                                                         <  2:2160253/109‑1 (MQ=255)
                                                                                                       gAAATTTGTGCTGCCTGTGCCGGTCCTTTCGCGAGTAAGCCCG‑‑CGCCTGCAGGGGTAGCGTACCCGTGCAGGTGCGGGCTTGCCCGCGAAGAGGGCGGTGCAGGCTAGAACAGCACTTTCGCCACGTCGGCAAAGCGCTTGGCGAAATGc                                              >  2:2419546/1‑150 (MQ=255)
                                                                                                       gAAATTTGTGCTGCCTGTGCCGGTCCTTTCGCGAGTAAGCCCG‑‑CGCCTGCAGGGGTAGCGTACCCGTGCAGGTGCGGGCTTGCCCGCGAAGAGGGCGGTGCAGGCTAGAACAGCACTTTCGCCACGTCAGCAAAGCGCTTGGCGAAATGc                                              <  1:2083454/150‑1 (MQ=255)
                                                                                                       gAAATTTGTGCTGCCTGTGCCGATCCTTTCGCGAGTAAGCCCG‑‑CGCCTGCAGGGGTAGCGTACCCGGGCAGGTGCGGGCTTGCCCGCGAAGAGGGCGGTGCAGGCTAGAACAGCACTTTCGCCACGTCGGCAAAGCGCTTGGCGAAATGc                                              >  2:868802/1‑150 (MQ=255)
                                                                                                         aaTTTGTGCTGCCTGTGCCGGTCCTTTCGCGAGTAAGCCCG‑‑CGCCTGCAGGGGTAGCGTACCCGTGCAGGTGCGGGCTTGCCCGCGAAGAGGGCGGTGCAGGCTAGAACAGCACTTTCGCCACGTCGGCAAAGCGCTTGGCGAAATGCAc                                            <  1:2419546/150‑1 (MQ=255)
                                                                                                             tgtgCTGCCTGTGCCGGTCCTTTCGCGAGTAAGCCCG‑‑CGCCTGCAGGGGTAGCGTACCCGTGCAGGTGCGGGCTTGCCCGCGAAGAGGGCGGTGCAGGCTAGAACAGCACTTTCGCCACGTCGGCAAAGCGCTTGGCGAAATGCACTGTc                                        <  2:1434268/150‑1 (MQ=255)
                                                                                                                    ccTGTGCCGGTCCTTTCGCGAGTAAGCCCG‑‑CGCCTGCAGGGGTAGCGTACCCGTGCAGGTGCGGGCTTGCCCGCGAAGAGGGCGGTGCAGGCTAGAAAAGCACTTTCGCCACGTCGGCGAAGCGCTTGGCGAAATGCACTGTCAGGCCtt                                 >  2:1245657/1‑150 (MQ=255)
                                                                                                                     cTGTGCCGGTCCTTTCGCGAGTAAGCCCG‑‑CGCCTGCAGGGGTAGCGTACCCGTGCAGGTGCGGGCTTGCCCGCGAAGAGGGCGGTGCAGGCTAGAACAGCACTTTc                                                                            <  2:2202521/106‑1 (MQ=255)
                                                                                                                     cTGTGCCGGTCCTTTCGCGAGTAAGCCCG‑‑CGCCTGCAGGGGTAGCGTACCCGTGCAGGTGCGGGCTTGCCCGCGAAGAGGGCGGTGCAGGCTAGAACAGCACTTTc                                                                            >  1:2202521/1‑106 (MQ=255)
                                                                                                                     cTGTGCCGGTCCTTTCGCGAGTAAGCCCG‑‑AGCCTGCAGGGGTAGCGTACCCGTGCAGGTGCGGGCTTGCCCGCGAAGAGGGCGGTGCAGGCTAGAACAGCACTTTCGCCACGTCGGCAAAGCGCTTGGCGAAATGCACTGCCAGGCCTTc                                >  1:312809/1‑150 (MQ=255)
                                                                                                                          ccGGTCCTTTCGCGAGTAAGCCCG‑‑CGCCTGCAGGGGTAGCGTACCCGTGCAGGTGCGGGCTTGCCCGCGAAGAGGGCGGTGCAGGCTAGAACa                                                                                    <  1:1387966/93‑1 (MQ=255)
                                                                                                                          ccGGTCCTTTCGCGAGTAAGCCCG‑‑CGCCTGCAGGGGTAGCGTACCCGTGCAGGTGCGGGCTTGCCCGCGAAGAGGGCGGTGCAGGCTAGAACa                                                                                    >  2:1387966/1‑93 (MQ=255)
                                                                                                                          ccGGTCCTTTCGCGAGTAAGCCCG‑‑CGCCTGCAGGGGTAGCGTACCCGTGCAGGTGCGGGCTTGCCCGCGAAGAGGGCGGTGCAGGCTAGAACAGCACTTTCGCCACGTCGGCAAAGCGCTTGGCGAAATGCACTGTCAGGCCTTCGCGCa                           <  2:939122/150‑1 (MQ=255)
                                                                                                                            ggTCCTTTCGCGAGTAAGCCCG‑‑CGCCTGCAGGGGTAGCGTACCCGTGCAGGTGCGGGCTTGCCCGCGAAGAGGGCGGTGCAGGCTAGAACAGCACTTTCGCCACGTCGGCAAAGCGCTTGGCGAAATGCACTGTCAGGCCTTCGCGCAg                          <  1:1990987/149‑1 (MQ=255)
                                                                                                                            ggTCCTTTCGCGAGTAAGCCCG‑‑CGCCTGCAGGGGTAGCGTACCCGTGCAGGTGCGGGCTTGCCCGCGAAGAGGGCGGTGCAGGCTAGAACAGCACTTTCGCCACGTCGGCAAAGCGCTTGGCGAAATGCACTGTCAGGCCTTCGCGCAg                          >  2:1990987/1‑149 (MQ=255)
                                                                                                                              tCCTTTCGCGAGTAAGCCCG‑‑CGCCTGCAGGGGTAGCGTACCCGTGCAGGTGCGGGCTTGCCCGCGAAGAGGGCGGTGCAGGCTAGAACa                                                                                    <  2:769600/89‑1 (MQ=255)
                                                                                                                              tCCTTTCGCGAGTAAGCCCG‑‑CGCCTGCAGGGGTAGCGTACCCGTGCAGGTGCGGGCTTGCCCGCGAAGAGGGCGGTGCAGGCTAGAACa                                                                                    >  1:769600/1‑89 (MQ=255)
                                                                                                                                 tttCGCGAGTAAGCCCG‑‑CGCCTGCAGGGGTAGCGTACCCGTGCAGGTGCGGGCTTGCCCGCGAAGAGGGAGGTGCAGGCTAGAACAGCACTTTCGCCACGTCGGCAAAGCGCTTGGCGAAATGCACTGTCAGGCCTTCGCGCAGGTAGTc                    >  2:750626/1‑150 (MQ=255)
                                                                                                                                   tCGCGAGTAAGCCCG‑‑CGCCTGCAGGGGTAGCGTACCCGTGCAGGTGCGGGCTTGCCCGCGAAGAGGGCGGTGCa                                                                                              >  1:611999/1‑74 (MQ=255)
                                                                                                                                   tCGCGAGTAAGCCCG‑‑CGCCTGCAGGGGTAGCGTACCCGTGCAGGTGCGGGCTTGCCCGCGAAGAGGGCGGTGCa                                                                                              <  2:611999/74‑1 (MQ=255)
                                                                                                                                    cgcgAATAAGCCCG‑‑CGCCTGCAGGGGTAGCGTACCCGTGCAGGTGCGGGCTTGCCCGCGAAGAGGGCGGTGCGGGCTAGAACAGCACTTTCGCCACGTCGGCAAAGCGCTTGGCGAAATGCACTGTCAGGCCTTCGCGCAGGTAGTCCgg                 <  2:1751992/150‑1 (MQ=255)
                                                                                                                                     gcgAGTAAGCCCG‑‑CGCCTGCAGGGGTAGCGTACCCGTGCAGGTGCGGGCTTGCCCGCGAAGAGGGCGGTGCAGGCTAGAACAGCACTTTCGCCACGTCGGCAAAGCGCTTGGCGAAATGCACTGTCAGGCCTTCGCGCAGGTAGTCCGGc                >  2:1660039/1‑150 (MQ=255)
                                                                                                                                      cgAGTAAGCCCG‑‑CGCCTGCAGGGGTAGCGTACCCGTGCAGGTGCGGGCTTGCCCGCGAAGAGGGCGGTGCAGGCTAGAACAGCACTTTCGCCCAGTTTAACTAGCGCTTGGCGAAATGCACTGTCAGGCCTTCGCGCAGGTAGTCCGGCa               >  2:1243729/1‑150 (MQ=255)
                                                                                                                                             gCCCG‑‑CGCCTGCAGGGGTAGCGTACCCGTGCAGGTGCGGGCTTGCCCGCGAAGAGGGCGGTGCAGGCTAGAACAGCACTTTCGCCACGTCGGCAAAGCGCTTGGCGAAATGCAc                                            <  2:784725/114‑1 (MQ=255)
                                                                                                                                             gCCCG‑‑CGCCTGCAGGGGTAGCGTACCCGTGCAGGTGCGGGCTTGCCCGCGAAGAGGGCGGTGCAGGCTAGAACAGCACTTTCGCCACGTCGGCAAAGCGCTTGGCGAAATGCAc                                            >  1:784725/1‑114 (MQ=255)
                                                                                                                                              cccG‑‑CGCCTGCAGGGGTAGCGTACCCGTGCAGGTGCGGGCTTGCCCGCGAAGAGGGCGGTGCAGGCTAGAACAGCACTTTCGCCACGTCGGCAAAGCGCTTGGCGAAATGCa                                             <  2:2009850/112‑1 (MQ=255)
                                                                                                                                              cccG‑‑CGCCTGCAGGGGTAGCGTACCCGTGCAGGTGCGGGCTTGCCCGCGAAGAGGGCGGTGCAGGCTAGAACAGCACTTTCGCCACGTCGGCAAAGCGCTTGGCGAAATGCACTGTCAGGCCTTCGCGCAGGTAGTCCGGCAGTTCTTCg       >  2:497578/1‑150 (MQ=255)
                                                                                                                                              cccG‑‑CGCCTGCAGGGGTAGCGTACCCGTGCAGGTGCGGGCTTGCCCGCGAAGAGGGCGGTGCAGGCTAGAACAGCACTTTCGCCACGTCGGCAAAGCGCTTGGCGAAATGCACTGTCAGGCCTTCGCGCAGGTAGTCCGGCAGTTCTTCg       >  1:936061/1‑150 (MQ=255)
                                                                                                                                               ccG‑‑CGCCTGCAGGGGTAGCGTACCCGTGCAGGTGCGGGCTTGCCCGCGAAGAGGGCGGTGCAGGCTAGAACAGCACTTTCGCCAc                                                                       >  1:1760091/1‑85 (MQ=255)
                                                                                                                                               ccG‑‑CGCCTGCAGGGGTAGCGTACCCGTGCAGGTGCGGGCTTGCCCGCGAAGAGGGCGGTGCAGGCTAGAACAGCACTTTCGCCAc                                                                       <  2:1760091/85‑1 (MQ=255)
                                                                                                                                               ccG‑‑CACCTGCACGGGTACGCTACCCCTGCAGGCGCGGGCTTACTCGCGAAGAGGGCGGTGCAGGCTAGAACAGCACTTTCGCCACGTCGGCAAAGCGCTTGGCGAAATGCACTGTCAGGCCTTCGCGCAGGTAGTCCGGCAGTTCTTCGa      >  2:881660/1‑150 (MQ=25)
                                                                                                                                               ccG‑‑CACCTGCACGGGTACGCTACCCCTGCAGGCGCGGGCTTACTCGCGAAGAGGGCGGTGCAGGCTAGAACAGCACTTTCGCCACGTCGGCAAAGCGCTTGGCGAAATGCACTGTCAGGCCTTCGCGCAGGTAGTCCGGCAGTTCTTCGa      >  1:2191365/1‑150 (MQ=25)
                                                                                                                                                cg‑‑cgCCTGCAGGGGTAGCGTACCCGTGCAGGTGCGGGCTTGCCCGCGAAGAGGGCGGTGCAGGCTAGAACAGCACTTTCGCCACGTCGGCAAAGCGCTTGGCGAAATGCACTGTCAGGCCTTCGCGCAGGTAGTCCGGCAGTTCTTCGaa     <  1:967548/150‑1 (MQ=255)
                                                                                                                                                cg‑‑cgCCTGCAGGGGTAGCGTACCCGTGCAGGTGCGGGCTTGCCCGCGAAGAGGGCGGTGCAGGCTAGAACAGCACTTTCGCCACGTCGGCAAAGCGCTTGGCGAAATGCACTGTCAGGCCTTCGCGCAGGTAGTCCGGCAGTTCTTCGaa     <  1:49837/150‑1 (MQ=255)
                                                                                                                                                 g‑‑cgcCTGCAGGGGTAGCGTACCCGTGCAGGTGCGGGCTTGCCCGCGAAGAGGGCGGTGCAGGCTAGAACAGCACTTTCGCCACGTCGGCAAAGCGCTTGGCGAAATGCACTGTCAGGCCTTCGCGCAGGTAGTCCGGCAGTTCTTCGaa     <  1:2001520/149‑1 (MQ=255)
                                                                                                                                                 g‑‑cgcCTGCAGGGGTAGCGTACCCGTGCAGGTGCGGGCTTGCCCGCGAAGAGGGCGGTGCAGGCTAGAACAGCACTTTCGCCACGTCGGCAAAGCGCTTGGCGAAATGCACTGTCAGGCCTTCGCGCAGGTAGTCCGGCAGTTCTTCGaa     >  2:2001520/1‑149 (MQ=255)
                                                                                                                                                     |                                                                                                                                                      
GCGCAAGCAGAAAGCAGGGCGAAGCTCAGGGGAACGAGCAGACGAGGAGAGGTCATGGTGGCTCCGAAGGCTGGCGATGGCCTGAGGATACCCGATGATCGAGGAAATTTGTGCTGCCTGTGCCGGTCCTTTCGCGAGTAAGCCCG‑‑CGCCTGCAGGGGTAGCGTACCCGTGCAGGTGCGGGCTTGCCCGCGAAGAGGGCGGTGCAGGCTAGAACAGCACTTTCGCCACGTCGGCAAAGCGCTTGGCGAAATGCACTGTCAGGCCTTCGCGCAGGTAGTCCGGCAGTTCTTCGAAGTC  >  NC_002947/1642159‑1642455

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