Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,642,334 T→C 5.2% 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 5.2% 408.6 / 7.6 134intergenic (‑120/+32)PP_1442/lon‑Ihypothetical protein/DNA‑binding ATP‑dependent protease
Reads supporting (aligned to +/- strand):  ref base T (64/63);  new base C (4/3);  total (68/66)
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.24e-01

Only 100 of 135 total aligned reads displayed.

CAGGGGAACGAGCAGACGAGGAGAGGTCATGGTGGCTCCGAAGGCTGGCGATGGCCTGAGGATACCCGATGATCGAGGAAATTTGTGCTGCCTGTGCCGGTCCTTTCGCGAGTAAGCCCG‑‑CGCCTGCAGGGGTAGCGTACCCGTGCAGGTGCGGGCTTGCCCGCGAAGAGGGCGGTGCAGGCTAGAACAGCACTTTCGCCACGTCGGCAAAGCGCTTGGCGAAATGCACTGTCAGGCCTTCGCGCAGGTAGTCCGGCAGTTCTTCGAAGTCGCCGCGGTTGGGCTCCGGCAGGATCA  >  NC_002947/1642185‑1642481
                                                                                                                                                        |                                                                                                                                                   
cAGGGGAACGAGCAGACGAGGAGAGGTCATGGTGGCTCCGAAGGCTGGCGATGGCCTGAGGATACCCGATGATCGAGGAAATTTGTGCTGCCTGTGCCGGTCCTTTCGCGAGTAAGCCCG‑‑CGCCTGCAGGGGTAGCGTACCCGTGCAGGt                                                                                                                                                     <  1:655412/150‑1 (MQ=255)
   gggAACGAGCAGACGAGGAGAGGTCATGGTGGCTCCGAAGGCTGGCGATGGCCTGAGGATACCCGATGATCGAGGAAATTTGTGCTGCCTGTGCCGGTCCTTTCGCGAGTAAGCCCG‑‑CGCCTGCAGGGGTAGCGTACCCGTGCAGGTGCg                                                                                                                                                  <  2:533256/150‑1 (MQ=255)
    ggAACGAGCAGACGAGGAGAGGTCATGGTGGCTCCGAAGGCTGGCGATGGCCTGAGGATACCCGATGATCGAGGAAATTTGTGCTGCCTGTGCCGGTCCTTTCGCGAGTAAGCCCG‑‑CGCCTGCAGGGGTAGCGTACCCGTGCAGGTGCgg                                                                                                                                                 <  2:1233996/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)
                       agGTCATGGTGGCTCCGAAGGCTGGCGATGGCCTGAGGATACCCGATGATCGAGGAAATTTGTGCTGCCTGTGCCGGTCCTTTCGCGAGTAAGCCCG‑‑CGCCTGCAGGGGTAGCGTACCCGTGCAGGTGCGGGCTTGCCCGCGAAGAGGGc                                                                                                                              >  1:1449027/1‑150 (MQ=255)
                         gTCATGGTGGCTCCGAAGGCTGGCGATGGCCTGAGGATACCCGATGATCGAGGAAATTTGTGCTGCCTGTGCCGGTCCTTTCGCGAGTAAGCCCG‑‑CGCCTGCAGGGGTAGCGTACCCGTGCAGGTGCGGGCTTGCCCGCGAAGAGGGCgg                                                                                                                            >  2:941910/1‑150 (MQ=255)
                         gTCATGGTGGCTCCGAAGGCTGGCGATGGCCTGAGGATACCCGATGATCGAGGAAATTTGTGCTGCCTGTGCCGGTCCTTTCGCGAGTAAGCCCG‑‑CGCCTGCAGGGGTAGCGTACCCGTGCAGGTGCCGGCTTGCCCGCGAAGAGGGCgg                                                                                                                            <  1:742145/150‑1 (MQ=255)
                          tCATGGTGGCTCCGAAGGCTGGCGATGGCCTGAGGATACTCGATGATCGAGGAAATTTGTGCTGCCTGTGCCGGTCCTTTCGCGAGTAAGCCCG‑‑CGCCTGCAGGGGTAGCGTACCCGTGGAGGTGCGGGCTTGCCCGCGAAGAGGGCGGt                                                                                                                           <  2:1406949/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:477740/1‑150 (MQ=255)
                               gtggCTCCGAAGGCTGGCGATGGCCTGAGGATACCCGATGATCGAGGAAATTTGTGCTGCCTGTGCCGGTCCTTTCGCGAGTAAGCCCG‑‑CGCCTGCAGGGGTAGCGTACCCGTGCAGGTGCGGGCTTGCCCGCGAAGAGGGCGGTGCAgg                                                                                                                      <  2:1449027/150‑1 (MQ=255)
                               gtggCTCCGAAGGCTGGCGATGGCCTGAGGATACCCGATGATCGAGGAAATTTGTGCTGCCTGTGCCGGTCCTTTCGCGAGTAAGCCCG‑‑CGCCTGCAGGGGTAGCGTACCCGTGCAGGTGCGGGCTTGCCCGCGAAGAGGGCGGTGCAgg                                                                                                                      >  1:2358548/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:1889325/150‑1 (MQ=255)
                                     ccGAAGGCTGGCGATGGCCTGAGGATACCCGATGATCGAGGAAATTTGTGCTGCCTGTGCCGGTCCTTTCGCGAGTAAGCCCG‑‑CGCCTGCAGGGGTAGCGTACCCGTGCAGGTGCGGGCTTATGCGCGAAGAGGGCGGTGCAGGCTAGaa                                                                                                                <  1:1243203/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:49837/1‑150 (MQ=255)
                                             tGGCGATGGCCTGAGGATACCCGATGATCGAGGAAATTTGTGCTGCCTGTGCCGGTCCTTTCGCGAGTAAGCCCG‑‑CGCCTGCAGGGGTAGCGTACCCGTGCAGGTGCGGGCTTGCCCGCGAAGAGGGCGGTGCAGGCTAGAACAGCACtt                                                                                                        >  2:1413578/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)
                                                               aCCCGATGACAGAGGAAAATTGTGCTGCCTGTGCCGGTCCTTTCGCGAGTAAGCCCG‑‑CGCCTGCAGGGGTAGCGTACCCGTGCAGGTGCGGGCTTGCCCGCGAAGAGGGCGGTGCAGGCTAGAACAGCACTTTCGCCACGTCGGCAAAgc                                                                                      >  2:1919849/1‑150 (MQ=255)
                                                                 ccGATGATCGAGGAAATTTGTGCTGCCTGTGCCGGTCCTTTCGCGAGTAAGCCCG‑‑CGCCTGCAGGGGTAGCGTACCCGTGCAGGTGCGGGCTTGCCCGCGAAGAGGGCGGTGCAGGCTAGAACAGCACTTTCGCCACGTCGGCAAAgcgc                                                                                    >  2:1675230/1‑150 (MQ=255)
                                                                 ccGATGATCGAGGAAATTTGTGCTGCCTGTGCCGGTCCTTTCGCGAGTAAGCCCG‑‑CGCCTGCAGGGGTAGCGTACCCGTGCAGGTGCGGGCTTGCCCGCGAAGAGGGCGGTGCAGGCTAGAACAGCACTTTCGCCACGTCGGCAAAgcgc                                                                                    >  1:939122/1‑150 (MQ=255)
                                                                 ccGATGATCGAGGAAATTTGTGCTGCCTGTGCCGGTCCTTTCGCGAGTAAGCCCG‑‑CGCCTGCAGGGGTAGCGTACCCGTGCAGGTGCGGGCTTGCCCGCGAAGAGGGCGGTGCAGGCTAGAACAGCACTTTCGCCACGTCGGCAAAgcgc                                                                                    >  1:1064267/1‑150 (MQ=255)
                                                                 ccGATGATCGAGGAAATTTGTGCTGCCTGTGCCGGTCCTTTCGCGAGTAAGCCCG‑‑CGCCTGCAGGGGTAGCGTACCCGTGCAGGTGCGGGCTTGCCCGCGAAGAGGGCGGTGCAGGCTAGAACAGCACTTTCGCCACGTCGGCAAAgcgc                                                                                    >  1:454340/1‑150 (MQ=255)
                                                                 ccGATGATCGAGGAAATTTGTGCTGCCTGTGCCGGTCCTTTCGCGAGTAAGCCCG‑‑CGCCTGCAGGGGTAGCGTACCCGTGCAGGTGCGGGCTTGCCCGCGAAGAGGGCGGTGCAGGCTAGAACAGCACTTTCGCCACGTCGGCAAAgcgc                                                                                    >  1:2327994/1‑150 (MQ=255)
                                                                  cGATGATCGAGGAAATTTGTGCTGCCTGTGCCGGTCCTTTCGCGAGTAAGCCCG‑‑CGCCTGCAGGGGTAGCGTACCCGTGCAGGTGCGGGCTTGCCCGCGAAGAGGGCGGTGCAGGc                                                                                                                     <  1:2287160/116‑1 (MQ=255)
                                                                  cGATGATCGAGGAAATTTGTGCTGCCTGTGCCGGTCCTTTCGCGAGTAAGCCCG‑‑CGCCTGCAGGGGTAGCGTACCCGTGCAGGTGCGGGCTTGCCCGCGAAGAGGGCGGTGCAGGc                                                                                                                     >  2:2287160/1‑116 (MQ=255)
                                                                   gatgatCGAGGAAATTTGTGCTGCCTGTGCCGGTCCTTTCGCGAGTAAGCCCG‑‑CGCCTGCAGGGGTAGCGTACCCGTGCAGGTGCGGGCTTGCCCGCGAAGAGGGCGGTGCAGGCTAGAACAGCACTTTCGCCACGTCGGCAAAGCGCtt                                                                                  >  2:48620/1‑150 (MQ=255)
                                                                      gatCGAGGAAATTTGTGCTGCCTGTGCCGGTCCTTTCGCGAGTAAGCCCG‑‑CGCCTGCAGGGGTAGCGTACCCGTGCAGGTGCGGGCTTGCCCGCGAAGAGGGCGGTGCAGGCTAGAACAGCACTTTCGCCACGTCGGCAAAGCGCTTGGc                                                                               <  2:454340/150‑1 (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)
                                                                           aGGAAATTTGTGCTGCCTGTGCCGGTCCTTTCGCGAGTAAGCCCG‑‑CGCCTGCAGGGGTAGCGTACCCGTGCAGGTGCGGGCTTGCCCGCGAAGAGGGCGGTGCAGGCTa                                                                                                                   >  1:2160253/1‑109 (MQ=255)
                                                                             gAAATTTGTGCTGCCTGTGCCGGTCCTTTCGCGAGTAAGCCCG‑‑CGCCTGCAGGGGTAGCGTACCCGTGCAGGTGCGGGCTTGCCCGCGAAGAGGGCGGTGCAGGCTAGAACAGCACTTTCGCCACGTCGGCAAAGCGCTTGGCGAAATGc                                                                        >  2:2419546/1‑150 (MQ=255)
                                                                             gAAATTTGTGCTGCCTGTGCCGATCCTTTCGCGAGTAAGCCCG‑‑CGCCTGCAGGGGTAGCGTACCCGGGCAGGTGCGGGCTTGCCCGCGAAGAGGGCGGTGCAGGCTAGAACAGCACTTTCGCCACGTCGGCAAAGCGCTTGGCGAAATGc                                                                        >  2:868802/1‑150 (MQ=255)
                                                                               aaTTTGTGCTGCCTGTGCCGGTCCTTTCGCGAGTAAGCCCG‑‑CGCCTGCAGGGGTAGCGTACCCGTGCAGGTGCGGGCTTGCCCGCGAAGAGGGCGGTGCAGGCTAGAACAGCACTTTCGCCACGTCGGCAAAGCGCTTGGCGAAATGCAc                                                                      <  1:2419546/150‑1 (MQ=255)
                                                                                        tgcCTGTGCCGGTCCTTTCGCGAGTAAGCCCG‑‑CGCCTGCAGGGGTAGCGTACCCGTGCAGGTGCGGGCTTGCCCGCGAAGAGGGCGGTGCAGGCTAGAACAGCACTTTCGCCACGTCGGCAAAGCGCTTGGCGAAATGCACTGTCAGGcc                                                             >  1:2419057/1‑150 (MQ=255)
                                                                                        tgcCTGTGCCGGTCCTTTCGCGAGTAAGCCCG‑‑CGCCTGCAGGGGTAGCGTACCCGTGCAGGTGCGGGCTTGCCCGCGAAGAGGGCGGTGCAGGCTAGAACAGCACTTTCGCCACGTCGGCAAAGCGCTTGGCGAAATGCACTGTCAGGcc                                                             >  2:483334/1‑150 (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                                                                                                              >  2:1387966/1‑93 (MQ=255)
                                                                                                ccGGTCCTTTCGCGAGTAAGCCCG‑‑CGCCTGCAGGGGTAGCGTACCCGTGCAGGTGCGGGCTTGCCCGCGAAGAGGGCGGTGCAGGCTAGAACa                                                                                                              <  1:1387966/93‑1 (MQ=255)
                                                                                                ccGGTCCTTTCGCGAGTAAGCCCG‑‑CGCCTGCAGGGGTAGCGTACCCGTGCAGGTGCGGGCTTGCCCGCGAAGAGGGCGGTGCAGGCTAGAACAGCACTTTCGCCACGTCGGCAAAGCGCTTGGCGAAATGCACTGTCAGGCCTTCGCGCa                                                     <  2:939122/150‑1 (MQ=255)
                                                                                                  ggTCCTTTCGCGAGTAAGCCCG‑‑CGCCTGCAGGGGTAGCGTACCCGTGCAGGTGCGGGCTTGCCCGCGAAGAGGGCGGTGCAGGCTAGAACAGCACTTTCGCCACGTCGGCAAAGCGCTTGGCGAAATGCACTGTCAGGCCTTCGCGCAg                                                    >  2:1990987/1‑149 (MQ=255)
                                                                                                    tCCTTTCGCGAGTAAGCCCG‑‑CGCCTGCAGGGGTAGCGTACCCGTGCAGGTGCGGGCTTGCCCGCGAAGAGGGCGGTGCAGGCTAGAACa                                                                                                              >  1:769600/1‑89 (MQ=255)
                                                                                                    tCCTTTCGCGAGTAAGCCCG‑‑CGCCTGCAGGGGTAGCGTACCCGTGCAGGTGCGGGCTTGCCCGCGAAGAGGGCGGTGCAGGCTAGAACa                                                                                                              <  2:769600/89‑1 (MQ=255)
                                                                                                       tttCGCGAGTAAGCCCG‑‑CGCCTGCAGGGGTAGCGTACCCGTGCAGGTGCGGGCTTGCCCGCGAAGAGGGCGGTGCAGGCTAGAACAg                                                                                                             >  2:835831/1‑87 (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)
                                                                                                            cgAGTAAGCCCG‑‑CGCCTGCAGGGGTAGCGTACCCGTGCAGGTGCGGGCTTGCCCGCGAAGAGGGCGGTGCAGGCTAGAACAGCACTTTCGCCACGTCGGCAAAGCGCTTGGCGAAATGCACTGTCAGGCCTTCGCGCAGGTAGTCCGGCa                                         <  1:611895/150‑1 (MQ=255)
                                                                                                                   gCCCG‑‑CGCCTGCAGGGGTAGCGTACCCGTGCAGGTGCGGGCTTGCCCGCGAAGAGGGCGGTGCAGGCTAGAACAGCACTTTCGCCACGTCGGCAAAGCGCTTGGCGAAATGCAc                                                                      >  1:784725/1‑114 (MQ=255)
                                                                                                                    cccG‑‑CGCCTGCAGGGGTAGCGTACCCGTGCAGGTGCGGGCTTGCCCGCGAAGAGGGCGGTGCAGGCTAGAACAGCACTTTCGCCACGTCGGCAAAGCGCTTGGCGAAATGCa                                                                       >  1:2009850/1‑112 (MQ=255)
                                                                                                                    cccG‑‑CGCCTGCAGGGGTAGCGTACCCGTGCAGGTGCGGGCTTGCCCGCGAAGAGGGCGGTGCAGGCTAGAACAGCACTTTCGCCACGTCGGCAAAGCGCTTGGCGAAATGCa                                                                       <  2:2009850/112‑1 (MQ=255)
                                                                                                                    cccG‑‑CGCCTGCAGGGGTAGCGTACCCGTGCAGGTGCGGGCTTGCCCGCGAAGAGGGCGGTGCAGGCTAGAACAGCACTTTCGCCACGTCGGCAAAGCGCTTGGCGAAATGCACTGTCAGGCCTTCGCGCAGGTAGTCCGGCAGTTCTTCg                                 >  1:936061/1‑150 (MQ=255)
                                                                                                                      cg‑‑cgCCTGCAGGGGTAGCGTACCCGTGCAGGTGCGGGCTTGCCCGCGAAGAGGGCGGTGCAGGCTAGAACAGCACTTTCGCCACGTCGGCAAAGCGCTTGGCGAAATGCACTGTCAGGCCTTCGCGCAGGTAGTCCGGCAGTTCTTCGaa                               <  1:49837/150‑1 (MQ=255)
                                                                                                                      cg‑‑cgCCTGCAGGGGTAGCGTACCCGTGCAGGTGCGGGCTTGCCCGCGAAGAGGGCGGTGCAGGCTAGAACAGCACTTTCGCCACGTCGGCAAAGCGCTTGGCGAAATGCACTGTCAGGCCTTCGCGCAGGTAGTCCGGCAGTTCTTCGaa                               <  1:967548/150‑1 (MQ=255)
                                                                                                                      cg‑‑caCCTGCACGGGTACGCTACCCCTGCAGGCGCGGGCTTACTCGCGAAGAGGGCGGTGCAGGCTAGAACAGCACTTTCGCCACGTCGGCAAAGCGCTTGGCGAAATGCACTGTCAGGCCTTCGCGCAGGTAGTCCGGCAGTTCTTCGaa                               >  1:454951/1‑150 (MQ=25)
                                                                                                                       g‑‑cgcCTGCAGGGGTAGCGTACCCGTGCAGGTGCGGGCTTGCCCGCGAAGAGGGCGGTGCAGGCTAGAACAGCACTTTCGCCACGTCGGCAAAGCGCTTGGCGAAATGCACTGTCAGGCCTTCGCGCAGGTAGTCCGGCAGTTCTTCGaa                               <  1:2001520/149‑1 (MQ=255)
                                                                                                                           gcCTGCAGGGGTAGCGTACCCGTGCAGGTGCGGGCTTGCCCGCGAAGAGGGCGGTGCAGGCTAGAACAGCACTTTCGCCACGTCGGCAAAGCGCTTGGCGAAATGCACTGTCAGGCCTTCGCGCAGGTAGTCCGGCAGTTCTTCGAAGTc                            >  2:1847804/1‑150 (MQ=255)
                                                                                                                                cAGGGGTAGCGTACCCGTGCAGGTGCGGGCTTGCCCGCGAAGAGGGCGGTGCAGGCTAGAACAGCACTTTCGCCACGTCGGCAAAGCGCTTGGCGAAATGCACTGTCAGGCCTTCGCGCAGGTAGTCCGGCAGtt                                      >  2:488916/1‑135 (MQ=255)
                                                                                                                                cAGGGGTAGCGTACCCGTGCAGGTGCGGGCTTGCCCGCGAAGAGGGCGGTGCAGGCTAGAACAGCACTTTCGCCACGTCGGCAAAGCGCTTGGCGAAATGCACTGTCAGGCCTTCGCGCAGGTAGTCCGGCAGtt                                      <  1:488916/135‑1 (MQ=255)
                                                                                                                                    ggTAGCGTACCCGTGCAGGTGCGGGCTTGCCCGCGAAGAGGGCGGTGCAGGCTAGAACAGCa                                                                                                           <  1:1440376/62‑1 (MQ=255)
                                                                                                                                    ggTAGCGTACCCGTGCAGGTGCGGGCTTGCCCGCGAAGAGGGCGGTGCAGGCTAGAACAGCa                                                                                                           >  2:1440376/1‑62 (MQ=255)
                                                                                                                                     gTAGCGTACCCGTGCAGGTGCGGGCTTGCCCGCGAAGAGGGCGGTGCAGGCTAGAACAGCACTTTCGCCACGTCGGCAAAGCGCTTGGCGAAATGCACTGTCAGGCCTTCGCGCAGGTAGTCCGGCAGTTCTTCGAAGTCGCCGCGGTTg                  <  1:48620/150‑1 (MQ=255)
                                                                                                                                     gTAGCGTACCCGTGCAGGTGCGGGCTTGCCCGCGAAGAGGGCGGTGCAGGCTAGAACAGCACTTTCGCCACGTCGGCAAAGCGCTTGGCGAAATGCACTGTCAGGCCTTCGCGCAGGTAATCCTACAGTTCTTCGAAGTCGCCGCGGTTg                  <  1:663425/150‑1 (MQ=255)
                                                                                                                                     gTAGCGTACAAGTGCAGGTGAGGGATTGCCCGCGAAGAGGGCTGTGCAGGCTAGAACAGCACTTTCGCCACGTCGGCAAAGCGCTTGGCGAAATGCACTGTCAGGCCTTCGCGCAGGTAGTCCGGCAGTTCTTCGAAGTCGCCGCGGTTg                  <  2:2358548/150‑1 (MQ=255)
                                                                                                                                      tAGCGTACCAGTGCAGGTGCGGGCTTGCCCGCGAAGAGGGCGGTGCAGGCTAGAACAGCACTTTCGCCACGTCGGCAAAGCGCTTGGCGAAATGCACTGTCAGGCCTTCGCGCAGGTAGTACGGCAGTTCTTCGAAGTCGCCGCGGTTgg                 >  2:1089838/1‑150 (MQ=255)
                                                                                                                                       aGCGTACCCGTGCAGGTGCGGGCTTGCCCGAGAAGAGGGCGGTGCAGGCTAGAACAGCACTTTCGCCACGTCGGCAAAGCGCTTGGCGAAATGCACTGTCAGGCCTTCGCGCAGCTAGTCCGGCAGTTCTTCGAAGTCGCCGCGGTTggg                <  2:1277319/150‑1 (MQ=255)
                                                                                                                                          gTACCCGTGCAGGTGCGGGCTTGCCCGCGAAGAGGGCGGTGCAGGCTAGAACAGCACTTTCGCCACGTCGGCAAAGCGCTTGGCGAAATGCACTGTCAGGCCTTCGCGCAGGTAGTCCGGTTCTTCTTCGAAGTCGCCGCGGTTGGGCTc             <  1:1243729/150‑1 (MQ=255)
                                                                                                                                          gTACCCGTGCAGGTGCGGGCTTGCCCGCGAAGAGGGCGGTGCAGGATAGAACAGCACTTTCGCCACGTCGGCAAAGCGCTTGGCGAAATGCACTGTCAGGCCTTCGCGCAg                                                    <  2:1526036/111‑1 (MQ=255)
                                                                                                                                          gTACCCGTGCAGGTGCGGGCTTGCCCGCGAAGAGGGCGGTGCAGGATAGAACAGCACTTTCGCCACGTCGGCAAAGCGCTTGGCGAAATGCACTGTCAGGCCTTCGCGCAg                                                    >  1:1526036/1‑111 (MQ=255)
                                                                                                                                               ccGTGCAGGTGCGGGCTTGCCCGCGAAGAGGGCGGTGCAGGCTAGAACAGCACTTTCGCCACGTCGGCAAAGCGCTTGGCGAAATGCACTGTCAGGCCTTCGCGCAGGTAGTCCGGCAGAACATCGAAGTCGCCGCGGTTGGGCTCCGGc         >  2:2134388/1‑150 (MQ=255)
                                                                                                                                                cGTGCAGGTGCGGGCTTGCCCGCGAAGAGGGCGGTGCAGGCTAGAACAGCACTTTCGCCACGTCGGCAAAGCGCTTGGCGAAATGCACTGTCAGGCCTTCGCGCAGGTAGTCCGGCAGTTCTTCGCAGTCGCCGCGGTTTGGCTCCGGCa        <  2:1970932/150‑1 (MQ=255)
                                                                                                                                                cGTGCAGGTGCGGGCTTGCCCGCGAAGAGGGCGGTGCAGGCTAGAACAGCACTTTCGCCACGTCGGCAAAGCGCTTGGCGAAATGCACTGTCAGGCCTTCGCGCAGGTAGTCCGGCAGTTCTTCGAAGTCGCCGCGGTTGGGCTCCGGCa        <  1:1675230/150‑1 (MQ=255)
                                                                                                                                                 gTGCAGGTGCGGGCTTGCCCGCGAAGAGGGCGGTGCAGGCTAGAACAGCACTTTCGCCACGTCGGCAAAGCGCTTGGCGAAATGCACTGTCAGGCCTTCGCGCAGGTAGTCCGGCAGTTCTTCGAAGTCGCCGCGGTTGGGCTCCGGCAg       >  2:1872847/1‑150 (MQ=255)
                                                                                                                                                 gTGCAGGTGCGGGCTTGCCCGCGAAGAGGGCGGTGCAGGCTAGAACAGCACTTTCGCCACGTCGGCAAAGCGCTTGGCGAAATGCACTGTCAGGCCTTCGCGCAGGTAGTCCGGCAGTTCTTCGAAGTCGCCGCGGTTGGGCTCCGGCAg       >  2:1310499/1‑150 (MQ=255)
                                                                                                                                                 gTGCAGGTGCGGGCTTGCCCGCGAAGAGGGCGGTGCAGGCTAGAACAGCACTTTCGCCACGTCGGCAAAGCGCTTGGCGAAATGCACTGTCAGGCCTTCGCGCAGGTAGTCCGGCAGTTCTTCGAAGTCGCCGCGGTTGGGCTCCGGCAg       >  2:101354/1‑150 (MQ=255)
                                                                                                                                                      ggTGCGGGCTTGCCCGCGAAGAGGGCGGTGCAGGCTAGAACAGCACTTTCGCCACGTCGGCAAAGCGCTTGGCGAAATGCACTGTCAGGCCTTCGCGCAGGTAGTCCGGCAGTTCTTCGAAGTCGCCGCGGTTGGGCTCCGGCAGGATCa  <  2:477740/150‑1 (MQ=255)
                                                                                                                                                      ggTGCGGGCTTGCCCGCGAAGAGGGCGGTGCAGGCTAGAACAGCACTTTCGCCACGTCGGCAAAGCGCTTGGCGAAATGCACTGTCAGGCCTTCGCGCAGGTAGTCCGGCAGTTCTTCGAAGTCGCCGCGGTTGGGCTCCGGCAGGATCa  <  2:1677494/150‑1 (MQ=255)
                                                                                                                                                        |                                                                                                                                                   
CAGGGGAACGAGCAGACGAGGAGAGGTCATGGTGGCTCCGAAGGCTGGCGATGGCCTGAGGATACCCGATGATCGAGGAAATTTGTGCTGCCTGTGCCGGTCCTTTCGCGAGTAAGCCCG‑‑CGCCTGCAGGGGTAGCGTACCCGTGCAGGTGCGGGCTTGCCCGCGAAGAGGGCGGTGCAGGCTAGAACAGCACTTTCGCCACGTCGGCAAAGCGCTTGGCGAAATGCACTGTCAGGCCTTCGCGCAGGTAGTCCGGCAGTTCTTCGAAGTCGCCGCGGTTGGGCTCCGGCAGGATCA  >  NC_002947/1642185‑1642481

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: