Predicted mutation
evidence seq id position mutation annotation gene description
RA CP006584 2,506,026 T→C intergenic (‑1383/‑130) LY180_11895 ← / → LY180_11905 hypothetical protein/phosphate acetyltransferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0065842,506,0260TC97.9% 160.7 / ‑3.2 47intergenic (‑1383/‑130)LY180_11895/LY180_11905hypothetical protein/phosphate acetyltransferase
Reads supporting (aligned to +/- strand):  ref base T (0/0);  major base C (24/22);  minor base . (1/0);  total (25/22)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.43e-01

CCGCTATGTTGAAGACAACTACGCGACGAAAGAAGACGCGAAGCGCGCAATGGACGTTTACTGCCACCGCCTGGCGAAATACATCGGTGCCTACACTGCGCTGATGGATGGTCGTCTGGACGCTGTTGTATTCACTGGTGGTATCGGTGAAAATGCCGCGATGGTTCGTGAACTGTCTCTGGGCAAACTGGGCGTGCTGGGCTTTGAAGTTGATCATGAACGCAACCTGGCTGCACGTTTCGGCAAATCTGGTTTCATCAACAAAGAAGGTACCCGTCCTGCGGTGGTTATCCCAACCAACGAAGAACTGGTTATCGCGCAAGACGCGAGCCGCCTGACTGCCTGATTTCACACCGCCAGCTCAGCTGGCGGTGCTGTTTTGTAACCCGCCAAATCGGCGGTAACGAAAGAGGATAAACCGTGTCCCGTATTATTATGCTGATCCCTACCGGAACCAGCGTCGGTCTGACCAGCGTCAGCCTTGGCGTGATCCGTGCAATGGAACGCAAAGGCGTTCGTCTGAGCGTTTTCAAACCTATCGCTCAGCCGCGTACCGGTGGCGATGCGCCCGATCAGA  >  CP006584/2505736‑2506312
                                                                                                                                                                                                                                                                                                   |                                                                                                                                                                                                                                                                                               
ccGCTATGTTGAAGACAACTACGCGACGAAAGAAGACGCGAAGCGCGCAATGGACGTTTACTGCCACCGCCTGGCGAAATGCATCGGTGCCTACACTGCGCTGATGGATGGTCGTCTGGACGCTGTTGTATTCACTGGTGGTATCGGTGAAAATGCCGCGATGGTTCGTGAACTGTCTCTGGGCAAACTGGGCGTGCTGGGCTTTGAAGTTGATCATGAACGCAACCTGGCTGCACGTTTCGGCAAATCTGGTTTCATCAACAAAGAAGGTACCCGTCCTGCGGTGGTTACCccaacca                                                                                                                                                                                                                                                                                         <  1:517136/299‑1 (MQ=255)
 cGCTATGTTGAAGACAACTACGCGACGAAAGAAGACGCGAAGCGCGCAATGGACGTTTACTGCCACCGCCTGGCGAAATGCATCGGTGCCTACACTGCGCTGATGGATGGTCGTCTGGACGCTGTTGTATTCACTGGTGGTATCGGTGAAAATGCCGCGATGGTTCGTGAACTGTCTCTGGGCAAACTGGGCGTGCTGGGCTTTGAAGTTGATCATGAACGCAACCTGGCTGCACGTTTCGGCAAATCTGGTTTCATCAACAAAGAAGGTACCCGTCCTGCGGTGGTTACCccaaccaa                                                                                                                                                                                                                                                                                       >  1:33481/1‑299 (MQ=255)
   cTAGGTTGAATACAACTACGCCACGGAAGAAGACGCGAAGCGCGGAATGGCAGTTTACTGCCACCGCCTGGCGAAATGCAGCGGTGCCTACACTGGGCTGATGGGTGGTCGTCTGGGCGCTGTTGTATTCACTGGTGGTATCGGTGAAAATGCCGCGATGGTTCGTGAACTGTCTCTGGGCAAACTGGGCGTGCTGGGCTTTGAAGTTGATCATGAACGCAACCTGGCTGCACGTTTCGGCAAATCTGGTTTCATCAACAAAGAAGGTACCCGTCCTGCGGTGGTTACCCCAACCAACg                                                                                                                                                                                                                                                                                     <  2:611304/299‑1 (MQ=255)
               cAACTACGCGACGAAAGAAGACGCGAAGCGCGCAATGGACGTTTACTGCCACCGCCTGGCGAAATGCATCGGTGCCTACACTGCGCTGATGGATGGTCGTCTGGACGCTGTTGTATTCACTGGTGGTATCGGTGAAAATGCCGCGATGGTTCGTGAACTGTCTCTGGGCCAACTGGGCGTGCTGGGCTTTGAAGTTGATCATGAACGCAACCTGGCTGCACGTTTCGGCAACTCTGGTTCCATCAACAAAGAAGGTACCCGTCCTGCGGTGGTGACCCCAACCAACGAAGAACAGGttc                                                                                                                                                                                                                                                                         >  2:225687/1‑298 (MQ=255)
                aaatACGGGACGAAAGAGGCAGGGAGGGGCGAATGGAGGTTTGACTGCGCCCGCTTGGCGAAGTGGATGCGGGCCGCCACTGCGCGGATGGATGGTCGTCTGGACGCTGTTGTATTCACTGGTGGTATCGGTGAAAATGCCGCGATGGTTCGTGAACTGTCTCTGGGCAAACTGGGCGTGCTGGGCTTTGAAGTTGATCATGAACGCAACCTGGCTGCACGTTTCGGCAAATCTGGTTTCATCAACAAAGAAGGTACCCGTCCTGCGGTGGTTACCCCAACCAACGAAGAACTGGTTAt                                                                                                                                                                                                                                                                        <  2:293145/296‑1 (MQ=255)
                           gaAAGAAGACGCGAAGCGCGCAATGGACGTTTACTGCCACCGCCTGGCGAAATGCATCGGTGCCTACACTGCGCTGATGGATGGTCGTCTGGACGCTGTTGTATTCACTGGTGGTATCGGTGAAAATGCCGCGATGGTTCGTGAACTGTCTCTGGGCAAACTGGGCGTGCTGGGCTTTGAAGTTCATGATGAACGCAACCTGGCTGCCCGTTTCGGCAAACCGGGTTTCATCAAAAAAGCAGGTCCCCGTCCACGGGGTGTTCCCCCCACCAACGAAGAACGAGTTACCGCCAAAGaag                                                                                                                                                                                                                                                             >  2:520447/1‑297 (MQ=255)
                                             cgcAATGGACGTTTACTGCCACCGCCTGGCGAAATGCATCGGTGCCTACACTGCGCTGATGGATGGTCGTCTGGACGCTGTTGTATCCACTGGTGGTATCGGTGAAAATGCCGCGCTGGCTCGTGAACTGTCTCTGGTCAAACTGGGCGTGCTGGGCTTTGAAGTTGATCATGAACGCAACCTGGCTGCCCGTTTCGGCAAATCTGGTTTCATCACCAAAGAAGGTACCCGCCCGTCGGTGGTTACCCAACAACACGAACAAATGGTTACCGCGCAAAACCCGAgcccgccgccgaccc                                                                                                                                                                                                                                           >  2:447204/1‑287 (MQ=255)
                                                 aTGGACGTTTACTGCCACCGCCTGGCGAAATGCATCGGTGCCTACACTGCGCTGATGGATGGTCGTCTGGACGCTGTTGTATTCACTGGTGGTATCGGTGAAAATGCCGCGATGGTTCGTGAACTGTCTCTGGGCAAACTGGGCGTGCTGGGCTTTGAAGTTGATCATGAACGCAACCTGGCTGCACGTTTCGGCAAATCTGGTTTCATCAACAAAGAAGGTACCCGTCCTGCGGTGGTTACCCCAACCAACGAAGAACTGGTTATCGCGCAAGACGCGAGCCGCCTGACTGCCTGAtt                                                                                                                                                                                                                                       <  1:431727/299‑1 (MQ=255)
                                                         ttACTGCCACCGCCTGGCGAAATGCATCGGTGCCTACCCTGCGCTGATGGATGGTCGTCTGGACGCTGTTGTATTCACTGGTGGTATCGGTGAAAATGCCGCGATGGTTCGTGAACTGTCTCTGGGCAAACTGGGCGTGCTGGGCTTTGAAGTTGATCATGAACGCAACCTGGCTGCACGTTTCGGCAAATCTGGTTTCATCAACAAAGAAGGTACCCGTCCTGCGGTGGTTAcccc                                                                                                                                                                                                                                                                                              <  2:404076/237‑1 (MQ=255)
                                                         ttACTGCCACCGCCTGGCGAAATGCATCGGTGCCTACACTGCGCTGATGGATGGTCGTCTGGACGCTGTTGTATTCACTGGTGGTATCGGTGAAAATGCCGCGATGGTTCGTGAACTGTCTCTGGGCAAACTGGGCGTGCTGGGCTTTGAAGTTGATCATGAACGCAACCTGGCTGCACGTTTCGGCAAATCTGGTTTCATCAACAAAGAAGGTACCCGTCCTGCGGTGGTTAcccc                                                                                                                                                                                                                                                                                              >  1:404076/1‑237 (MQ=255)
                                                                        tGGCGAAATGCATCGGTGCCTACACTGCGCTGATGGATGGTCGTCTGGACGCTGTTGTATTCACTGGTGGTATCGGTGAAAATGCCGCGATGGTTCGTGAACTGTCTCTGGGCAAACTGGGCGTGCTGGGCTTTGAAGTTGATCATGAACGCAACCTGGCTGCACGTTTCGGCAAATCTGGTTTCATCAACAAAGAAGGTACCCGTCCTGCGGTGGTTACCCCAACCAACGAAGAACTGGTTATCGCGCAAGACGCGAGCCGCCTGACTGCCTGATTTCACACCGCCAGCTCAGCTGGc                                                                                                                                                                                                                 >  1:642667/1‑299 (MQ=255)
                                                                                           cTACACTGCGCTGATGGATGGTCGTCTGGACGCTGTTGTATTCACTGGTGGTATCGGTGAAAATGCCGCGATGGTTCGTGAACTGTCTCTGGGCAAACTGGGCGTGCTGGGCTTTGAAGTTGATCATGAACGCAACCTGGCTGCACGTTTCGGCAAATCTGGTTTCATCAACAAAGAAGGTACCCGTCCTGCGGTGGTTACCCCAACCAACGAAGAACTGGTTATCGCGCAAGACGCGAGCCGCCTGACTGCCTGATTTCACACCGCCAGCTCAGCTGGCGGTGCTGTTTTGTAACCCg                                                                                                                                                                                              <  1:447204/299‑1 (MQ=255)
                                                                                               acTGCGCTGATGGATGGTCGTCTGGACGCTGTTGTATTCACTGGTGGTATCGGTGAAAATGCCGCGATGGTTCGTGAACTGTCTCTGGGCAAACTGGGCGTGCTGGGCTTTGAAGTTGATCATGAACGCAACCTGGCTGCACGTTTCGGCAAATCTGGTTTCATCAACAAAGAAGGTACCCGTCCTGCGGTGGTTACCCCAACCAACGAAGAACTGGTTATCGCGCAAGACGCGAGCCGCCTGACTGCCTGATTTCACACCGCCAGCTCAGCTGGCGGGGCTGTTTTGTAACCCGCCaa                                                                                                                                                                                          >  1:394394/1‑299 (MQ=255)
                                                                                                     cTGATGGATGGTCGTCTGGACGCTGTTGTATTCACTGGTGGTATCGGTGAAAATGCCGCGATGGTTCGTGAACTGTCTCTGGGCAAACTGGGCGTGCTGGGCTTTGAAGTTGATCATGAACGCAACCTGGCTGCACGTTTCGGCAAATCTGGTTTCATCAACAAAGAAGGTACCCGTCCTGCGGTGGTTACCCCAACCAACGAAGAACTGGTTATCGCGCAAGACGCGAGCCGCCTGACTGCCTGATTTCACACCGCCAGCTCAGCTGGCGGTGCTGTTTTGTAAccagatcggaagag                                                                                                                                                                                    >  1:81711/1‑287 (MQ=255)
                                                                                                               gtcgtcTGGACGCTGTTGTATTCACTGGTGGTATCGGTGAAAATGCCGCGATGGTTCGTGAACTGTCTCTGGGCAAACTGGGCGTGCTGGGCTTTGAAGTTGATCATGAACGCAACCTGGCTGCACGTTTCGGCAAATCTGGTTTCATCAACAAAGAAGGTACCCGTCCTGCGGTGGTTACCCCAACCAACGAAGAACTGGTTATCGCGCAAGACGCGAGCCGCCTGACTGCCTGATTTCACACCGCCAGCTCAGCTGGCGGTGCTGTTTTGTAACCCGCCAAATCGGCGGTAACGaaa                                                                                                                                                                          >  1:575889/1‑299 (MQ=255)
                                                                                                                    cTGGACGCTGTTGTATTCACTGGTGGTATCGGTGAAAATGCCGCGATGGTTCGTGAACTGTCTCTGGGCAAACTGGGCGTGCTGGGCTTTGAAGTTGATCATGAACGCAACCTGGCTGCACGTTTCGGCAAATCTGGTTTCATCAACAAAGAAGGTACCCGTCCTGCGGTGGTTACCCCAACCAACGAAGAACTGGTTATCGCGCAAGACGCGAGCCGCCTGACTGCCTGATTTCACACCGCCAGCTCAGCTGGCGGTGCTGTTTTGTAACCCGCCAAATCGGCGGTAACGAAAGAGGa                                                                                                                                                                     <  1:14057/299‑1 (MQ=255)
                                                                                                                       gACGCTGTTGTATTCACTGGTGGTATCGGTGAAAATGCCGCGATGGTTCGTGAACTGTCTCTGGGCAAACTGGGCGTGCTGGGCTTTGAAGTTGATCATGAACGCAACCTGGCTGCACGTTTCGGCCAATCTGGTTTCATCAACAACGAAGGTACCAGTCCTGCGGTGGTTACCCCAACCAACGAAGAACTGGTTATCGCGCAAGACGCGAGCCGCCTGACTGCCTGATTTCACACCGCCAGCTCAGCTGGCGGTGCTGTTTTGTAACCCGCCAAATCGGCGGTAACGAAAGAGGATaa                                                                                                                                                                  >  2:255944/1‑299 (MQ=255)
                                                                                                                         tgtggttttATTTATTGGTGTGAACGGGGAAAATAGCGCGTGTGTTCGTGAGCTGTCTCTGGGCAAACTGGGGGTGCTGGGCTTTGAAGTTGATCATGAACGCAACCTGGCGGCACGTTTCGGCAAATCTGGTTTCATCAACAAAGAAGGTACCCGTCCTGCGGTGGTTACCCCAACCAACGAAGAACTGGTTATCGCGCAAGACGCGAGCCGCCTGACTGCCTGATTTCACACCGCCAGCTCAGCTGGCGGTGCTGTTTTGTAACCCGCCAAATCGGCGGTAACGAAAGAGGATAAAc                                                                                                                                                                <  2:33481/295‑1 (MQ=255)
                                                                                                                            tgttgtATTCACTGGTGGTATCGGTGAAAATGCCGCGATGGTTCGTGAACTGTCTCTGGGCAAACTGGGCGTGCTGGGCTTTGAAGTTGATCATGAACGCAACCTGGCTGCACGTTTCGGCAAATCTGGTTTCATCAACAAAGAAGGTACCCGTCCTGCGGTGGTTACCCCAACCAACGAAGAACTGGTTATCGCGCAAGACGCGAGCCGCCTGACTGCCTGATTTCACACCGCCAGCTCAGCTGGCGGTGCTGTTTTGTAACCCGCCAAATCGGCGGTAACGAAAGAGGATAAACCgt                                                                                                                                                             <  1:372790/299‑1 (MQ=255)
                                                                                                                             gttgtATTCACTGGTGGTATCGGTGAAAATGCCGCGATGGTTCGTGAACTGTCTCTGGGCAAACTGGGCGTGCTGGGCTTTGAAGTTGATCATGAACGCAACCTGGCTGCACGTTTCGGCAAATCTGGTTTCATCAACAAAGAAGGTCCCCGTCCTGCGGAGGTTACCCCAACCAACGCAGAAACGGTTCTCGCGCACGACGCGACCCGCCTCAATGCCGGATTTCACACCCCCCGCTCACACGTCGGGGCTGTTTTGTAACCCCCCAAACCCGCAGTACCGGAcgagaaggaacagtg                                                                                                                                                            >  2:548521/1‑288 (MQ=255)
                                                                                                                             gttgtATTCACTGGTGGTATCGGTGAAAATGCCGCGATGGTTCGTGAACTGTCTCTGGGCAAACTGGGCGTGCTGGGCTTTGAAGTTGATCATGAACGCAACCTGGCTGCACGTTTCGGCAAATCTGGTTTCATCAACAAAGAAGGTACCCGTCCTGCGGTGGTTACCCCAACCAACGAAGAACTGGTTATCGCGCAAGACGCGAGCCGCCTGACTGCCTGATTTCACACCGCCAGCTCAGCTGGCGGTGCTGTTTTGTAACCCGCCAAATCGGCGGTAACGAAAGAGGATAAACCgtg                                                                                                                                                            <  1:542587/299‑1 (MQ=255)
                                                                                                                                 tATTCACTGGTGGTATCGGTGAAAATGCCGCGATGGTTCGTGAACTGTCTCTGGGCAAACTGGGCGTGCTGGGCTTTGAAGTTGATCATGAACGCAACCTGGCTGCACGTTTCGGCAAATCTGGTTTCATCAACAAAGAAGGTACCCGTCCTGCGGTGGTTACCCCAACCAACGAAGAACTGGTTATCGCGCAAGACGCGAGCCGCCTGACTGCCTGATTTCACACCGCCAGCTCAGCTGGCGGTGCTGTTTTGTAACCCGCCAAATCGGCGGTAACGAAAGAGGATAAACCGTGTccc                                                                                                                                                        >  1:250698/1‑299 (MQ=255)
                                                                                                                                 tATTCACTGGTGGTATCGGTGAAAATGCCGCGATGGTTCGTGAACTGTCTCTGGGCAAACTGGGCGTGCTGGGCTTTGAAGTTGATCATGAACGCAACCTGGCTGCACGTTTCGGCAAATCTGGTTTCATCAACAAAGAAGGTACCCGTCCTGCGGTGGTTACCCCAACCAACGAAGAACTGGTTATCGCGCAAGACGCGAGCCGCCTGACTGCCTGATTTCACACCGCCAGCTCAGCTGGCGGTGCTGTTTTGTAACCCGCCAAATCGGCGGTAACGAAAGAGGATAAACCGTGTccc                                                                                                                                                        >  1:611412/1‑299 (MQ=255)
                                                                                                                                  aTTCACTGGTGGTATCGGTGAAAATGCCGCGATGGTTCGTGAACTGTCTCTGGGCAAACTGGGCGTGCTGGGCTTTGAAGTTGATCATGAACGCAACCTGGCTGCACGTTTCGGCAAATCTGGTTTCATCAACAAAGAAGGTACCCGTCCTGCGGTGGTTACCCCAACCAACGAAGAACTGGTTATCGCGCAAGACGCGAGCCGCCTGACTGCCTGATTTCACACCGCCAGCTCAGCTGGCGGTGCTGTTTTGTAACCCGCCAAATCGGCGGTAACGAAAGAGGATAAACCGTGTCCCg                                                                                                                                                       <  1:451743/299‑1 (MQ=255)
                                                                                                                                       cTGGTGGTATCGGTGAAAATGCCGCGATGGTTCGTGAACTGTCTCTGGGCAAACTGGGCGTGCTGGGCTTTGAAGTTGATCATGAACGCAACCTGGCTGCACGTTTCGGCAAATCTGGTTTCATCAACAAAGAAGGTACCCGTCCTGCGGTGGTTACCCCAACCAACGAAGAACTGGTTATCGCGCAAGACGCGAGCCGCCTGACTGCCTGATTTCACACCGCCAGCTCAGCTGGCGGTGCTGTTTTGTAACCCGCCAAATCGGCGGTAACGAAAGAGGATAAACCGTGTCCCGtatta                                                                                                                                                  >  1:136150/1‑299 (MQ=255)
                                                                                                                                       cTGGTGGTATCGGTGAAAATGCCGCGATGGTTCGTGAACTGTCTCTGGGCAAACTGGGCGTGCTGGGCTTTGAAGTTGATCATGAACGCAACCTGGCTGCACGTTTCGGCAAATCTGGTTTCATCAACAAAGAAGGTACCCGTCCTGCGGTGGTTACCCCAACCAACGAAGAACTGGTTATCGCGCAAGACGCGAGCCGCCTGACTGCCTGATTTCACACCGCCAGCTCAGCTGGCGGTGCTGTTTTGTAACCCGCCAAATCGGCGGTAACGAAAGAGGATAAACCGTGTCCCGtatta                                                                                                                                                  <  1:25365/299‑1 (MQ=255)
                                                                                                                                        tggtggTATCGGTGAAAATGCCGCGATGGTTCGTGAACTGTCTCTGGGCAAACTGGGCGTGCTGGGCTTTGAAGTTGATCATGAACGCAACCTGGCTGCACGTTTCGGCAAATCTGGTTTCATCAACAAAGAAGGTACCCGTCCTGCGGTGGTTACCCCAACCAACGAAGAACTGGTTATCGCGCAAGACGCGAGCCGCCTGACTGCCTGATTTCACACCGCCAGCTCAGCTGGCGGTGCTGTTTTGTAACCCGCCAAATCGGCGGTAACGAAAGAGGATAAACCGTGTCCCGtattat                                                                                                                                                 <  1:221272/299‑1 (MQ=255)
                                                                                                                                          tgggtATCGGTGGAAATGCCGCGATGGTTCATGAACTGTATCTGGGCAAACTGGGCGTGCTGGGCTTTGAAGTTGATCATGAACGCAACCTGGCTGCACGTTTCGGCAAATCTGGTTTCATCAACAAAGAAGGTACCCGTCCTGCGGTGGTTACCCCAACCAACGAAGAACTGGTTATCGCGCAAGACGCGAGCCGCCTGACTGCCTGATTTCACACCGCCAGCTCAGCTGGCGGTGCTGTTTTGTAACCCGCCAAATCGGCGGTAACGAAAGAGGATAAACCGTGTCCCGtattatta                                                                                                                                               <  2:136150/297‑1 (MQ=255)
                                                                                                                                            ggtATCGGTGAAAATGCCGCGATGGTTCGTGAACTGTCTCTGGGCAAACTGGGCGTGCTGGGCTTTGAAGTTGATCATGAACGCAACCTGGCTGCACGTTTCGGCAAATCTGGTTTCATCAACAAAGAAGGTACCCGTCCTGCGGTGGTTACCCCAACCAACGAAGAACTGGTTATCGCGCAAGACGCGAGCCGCCTGACTGCCTGATTTCACACCGCCAGCTCAGCTGGCGGTGCTGTTTTGTAACCCGCCAAATCGGCGGTAACGAAAGAGGATAAACCGTGTCCCGTATTATTATg                                                                                                                                             <  1:191705/299‑1 (MQ=255)
                                                                                                                                                       aaaTGCCTCGATGGTTCGTGAACTGTCTCTGGGCAAATTGGGCGTGCTGTGCTTTGAAGTTGATCATGATCGCAACCTTGCTGCACGTTTCGGCAAATCTGGTTTCATCAACAAAGAAGGTACCCGTCCTGCGGTGGTTACCCCAACCAACGAAGAACTGGTTATCGCGCAAGACGCGAGCCGCCTGACTGCCTGATTTCACACCGCCAGCTCAGCTGGCGGTGCTGTTTTGTAACCCGCCAAATCGTCGGTAACGAAAGAGGATAAACCGTGTCCCGTATTATTATGCTGATCCCTAc                                                                                                                                  <  1:27959/299‑1 (MQ=255)
                                                                                                                                                          tGCCGCGATGGTTCGTGAACTGTCTCTGGGCAAACTGGGCGTGCTGGGCTTTGAAGTTGATCATGAACGCAACCTGGCTGCACGTTTCGGCAAATCTGGTTTCATCAACAAAGAAGGTACCCGTCCTGCGGTGGTTACCCCAACCAACGAAGAACTGGTTATCGCGCAAGACGCGAGCCGCCTGACTGCCTGATTTCACACCGCCAGCTCAGCTGGCGGTGCTGTTTTGTAACCCGCCAAATCGGCGGTAACGAAAGAGGATAAACCGTGTCCCGTATTATTATGCTGATCCCTACCgg                                                                                                                               >  1:640351/1‑299 (MQ=255)
                                                                                                                                                               cgATGGTTCGTGAACTGTCTCTGGGCAAACTGGGCGTGCTGGGCTTTGAAGTTGATCATGAACGCAACCTGGCTGCACGTTTCGGCAAATCTGGTTTCATCAACAAAGAAGGTACCCGTCCTGCGGTGGTTACCCCAACCAACGAAGAACTGGTTATCGCGCAAGACGCGAGCCGCCTGACTGCCTGATTTCACACCGCCAGCTCAGCTGGCGGTGCTGTTTTGTAACCCGCCAAATCGGCGGTAACGAAAGAGGATAAACCGTGTCCCGTATTATTATGCTGAtagatcggaagagca                                                                                                                          >  1:626320/1‑285 (MQ=255)
                                                                                                                                                                        gTGAACTGTCTCTGGGCAAACTGGGCGTGCTGGGCTTTGAAGTTGATCATGAACGCAACCTGGCTGCACGTTTCGGCAAATCTGGTTTCATCAACAAAGAAGGTACCCGTCCTGCGGTGGTTACCCCAACCAACGAAGAACTGGTTATCGCGCAAGACGCGAGCCGCCTGACTGCCTGATTTCACACCGCCAGCTCAGCTGGCGGTGCTGTTTTGTAACCCGCCAAATCGGCGGTAACGAAAGAGGATAAACCGTGTCCCGTATTATTATGCTGATCCCTACCGGAACCAGCGTCGGTc                                                                                                                 >  1:484018/1‑299 (MQ=255)
                                                                                                                                                                                    tGGGCAAACTGGGGGTGTGGGGCTTTGAAGTTGATCATGAACGCAACCTGGCTTCACGTTTCGGCAAATCTGGTTTCATCAACAAAGAGGGTACCCGTCCTGCGGTGGTTACCCCAACCAACGAAGAACTGGTTATCGCGCAAGACGCGAGCCGCCTGACTGCCTGATTTCACACCGCCAGCTCAGCTGGCGGTGCTGTTTTGTAACCCGCCAAATCGGCGGTAACGAAAGAGGATAAACCGTGTCCCGTATTATTATGCTGATCCCTACCGGAACCAGCGTCGGTCTGACCAGCGTCa                                                                                                     <  2:575889/299‑1 (MQ=255)
                                                                                                                                                                                        cAAACTGGGCGTGCTGGGCTTTGAAGTTGATCATGAACGCAACCTGGCTGCACGTTTCGGCAAATCTGGTTTCATCAACAAAGAAGGTACCCGTCCTGCGGTGGTTACCCCAACCAACGAAGAACTGGTTATCGCGCAAGACGCGAGCCGCCTGACTGCCTGATTTCACACCGCCAGCTCAGCTGGCGGTGCTGTTTTGTAACCCGCCAAATCGGCGGTAACGAAAGAGGATAAACCGTGTCCCGTATTATTATGCTGATCCCTACCGGAACCAGCGTCGGTCTGACCAGCGTCAGCCt                                                                                                 >  1:143163/1‑299 (MQ=255)
                                                                                                                                                                                          caCTGGGCGTGCTGGGCTTTGAAGTTGATCATGAACGCAACCTGGCTGCACGTTTCGGCAAATCTGGTTTCATCAACAAAGAAGGTACCCGTCCAGCGGTGGTTACCCCAACCAACGAAGAACTGGTTATCGCGCAAGACGCGAGCCGCCTGACTGCCTGATTTCACACCGCCAGCTCAGCTCGCGGTGCTGTTTTGTAACCCGCCAAATCGGCGGTAACGAAAGAGGATAAACCGTGTCCCGTATTCTTATGCTGATCCCTACGGGAACCAGCGACGGTCTGCCCAGCGTCAGCCtta                                                                                               >  2:137667/2‑298 (MQ=255)
                                                                                                                                                                                                                          gAACGCAACCTGGCTGCACGTTTCGGCAAATCTGGTTTCATCAACAAAGAAGGTACCCGTCCTGCGGTGGTTACCCCAACCAACGAAGAACTGGTTATCGCGCAAGACGCGAGCCGCCTGACTGCCTGATTTCACACCGCCAGCTCAGCTGGCGGTGCTGTTTTGTAACCCGCCAAATCGGCGGTAACGAAAGAGGATAAACCGTGTCCCGTATTATTATGCTGATCCCTACCGGAACCAGCGTCGGTCTGACCAGCGTCAGCCTTGGCGTGATCCGTGCAATGGAACGCAAAgcgggt                                                               >  1:585574/1‑294 (MQ=255)
                                                                                                                                                                                                                             cGCAACCTGGCTGCACGTTTCGGCAAATCTGGTTTCATCAACAAAGAAGGTACCCGTCCTGCGGTGGTTACCCCAACCAACGAAGCACTGGTTATCGCGCAAGACGCGAGCCGCCTGACTGCCTGATTTCACACCGCCAGCTCAGCTGGCGGTGCTGTTTTGTAACCCGCCAAATCGGCGGTAACGAAAGAGGATAAACCGTGTCCCGTATTATTATGCTGACCCCTCCCGGAACCAGCGTCGGTCCGACCAGCGTCAGCCTTGGCCTTACCCGCGCCATGGACCGCAAAGCCGTTCGt                                                            >  2:348865/1‑299 (MQ=255)
                                                                                                                                                                                                                               cAACCTGGCTGCACGTTTCGGCAAATCTGGTTTCATCAACAAAGAAGGTACCCGTCCTGCGGTGGTTACCCCAACCAACGAAGAACTGGTTATCGCGCAAGACGCGAGCCGCCTGACTGCCTGATTTCACACCGCCAGCTCAGCTGGCGGTGCTGTTTTGTAACCCGCCAAATCGGCGGTAACGAAAGAGGATAAACCGTGTCCCGTATTATTATGCTGATCCCTACCGGAACCAGCGTCGGTCTGACCAGCGTCAGCCTTGGCGTGATCCGTGCAATGGAACGCAAAGGCGTTCGTCt                                                          <  1:92782/299‑1 (MQ=255)
                                                                                                                                                                                                                                           aCGTTTCGGCAAATCTGGTTTCATCAACAAAGAAGGTACCCGTACTGCGGTGGTTACCCCAACCAACGAAGAACTGGTTATCGCGCAAGACGCGAGCCGCCTGACTGCCTGATTTCACACCGCCAGCTCAGCTGGCGGTGCTGTTTTGTAACCCGCCAAATCGGCGGTAACGAAAGAGGATAAACCGTGTCCCGTATTATTATGCTGATCCCTACCGGAACCAGCGTCGGTCTGACCAGCGTCAGCCTTGGCGTGATCCGTGCAATGGAACGCAAAGGCGTTCGTCTGAGCGTTTTCaa                                              <  1:330653/299‑1 (MQ=255)
                                                                                                                                                                                                                                              tttCGGCAAATCTGGTTTCATCAACAAAGAAGGTACCCGTCCTGCGGTGGTTACCCCAACCAACGAAGAACTGGTTATCGCGCAAGACGCGAGCCGCCTGACTGCCTGATTTCACACCGCCAGCTCAGCTGGCGGTGCTGTTTTGTAACCCGCCAAATCGGCGGTAACGAAAGAGGATAAACCGTGTCCCGTATTATTATGCTGATCCCTACCGGAACCAGCGTCGGTCTGACCAGCGTCAGCCTTGGCGTGATCCGTGCAATGGAACGCAAAGGCGTTCGTCTGAGCGTTTTCAAAcc                                           <  1:520447/299‑1 (MQ=255)
                                                                                                                                                                                                                                                 cGGCAAATCTGGTTTCATCAACAAAGAAGGTACACGTCCTGCGGTGGTTACCCCAACCAACGAAGAAATGGTTATCGCGCAAGACGCGAGCCGACTGACTGCCTGATTTCACACCGCCAGCTCAGCTGGCGGTGCTGGTTTGTAACACGCCAAATCGGCGGTAACGCAAGAGGATAAACCGTGTACCGTATTATTAGGCTTATGCATACCGGCACCAGACAGAGTCGTACCAGCGGCCACATTGGCGTGGTCCGTGACGTCCAACCACAAGGGGTTCTTTCTAGCGTTTTCAAActtct                                        >  2:621664/1‑295 (MQ=255)
                                                                                                                                                                                                                                                         cTGGTTTCATCAACAAAGAAGGTACCCGTCCTGCGGTGGTTACCCCAACCAACGAAGAACTGGTTATCGCGCAAGACGCGAGCCGCCTGACTGCCTGATTTCACACCGCCAGCTCAGCTGGCGGTGCTGTTTTGTAACCCGCCAAATCGGCGGTAACGAAAGAGGATAAACCGTGTCCCGTATTATTATGCTGATCCCTACCGGAACCAGCGTCGGTCTGACCAGCGTCAGCCTTGGCGTGATCCGTGCAATGGAACGCAAAGGCGTTCGTCTGAGCGTTTTCAAACCTATCGCTCAGc                                >  1:366690/1‑299 (MQ=255)
                                                                                                                                                                                                                                                         cTGGTTTCATCAACAAAGAAGGTACCCGTCCTGCGGTGGTTACCCCAACCAACGAAGAACTGGTTATCGCGCAAGACGCGAGCCGCCTGACTGCCTGATTTCACACCGCCAGATCAGCTGGCGGTGCTGTTTTGTAACCCGACAAATCCGAGGTAACGACAGAGGATAAACCGTGTCCCGTATTATTATGCTGATCCCTACCGGAACCAGCGTCGGTCTGACCAGCGTCAGCCTTGGCGTGACCCGTGCCATGGACCGCAACGGCGCTCGTCGCAGCGTTTTCACAcccacagcccacg                                >  2:116006/1‑288 (MQ=255)
                                                                                                                                                                                                                                                               tcatcaACAAAGAAGGTACCCGTCCTGCGGTGGTTACCCCAACCAACGAAGAACTGGTTATCGCGCAAGACGCGAGCCGCCTGACTGCCTGATTTCACACCGCCAGCTCAGCTGGCGGTGCTGTTTTGTAACCCGCCAAATCGGCGGTAACGAAAGAGGATAAACCGTGTCCCGTATTATTATGCTGATCCCTACCGGAACCAGCGTCGGTCTGACCAGCGTCAGCCTTGGCGTGATCCGTGCAATGGAACGCAAAGGCGTTCGTCTGAGCGTTTTCAAACCTATCGCTCAGCCGCGTa                          >  1:291062/1‑299 (MQ=255)
                                                                                                                                                                                                                                                                tgtaaaaaaaGAAGGCACCGTGCCGGCGGTGTGTACCCCAACACGAGAAGAACTGGTTATCGCGCAAGGCGGGAGCCGCCTGTCTGCTTGATTTCAGACCGCCAGCTCAGCTGGCGGTGCTGTGTTGTAACCCGCCAAATCGGCGGTAACGAAAGAGGATAAACCGTGTCCCGGATTATTATGCTGATCCCTACCGGAACCAGCGTCGGTCTGACCAGCGTCAGCCTTGGCGTGATCCGTGCAATGGAACGCAAAGGCGTTCGTCTGAGCGTTTTCAAACCTATCGCTCAGCCGCGTAc                         <  2:40979/292‑1 (MQ=255)
                                                                                                                                                                                                                                                                    accaaAGAAGGTACCCGTCCTGCGGTGGTTACCCCAACCAACGAAGAACTGGTTATCGCGCAAGACGCGAGCCGCCTGACTGCCTGATTTCACACCGCCAGCTCAGCTGGCGGTGCTGTTTTGTAACCCGCCAAATCGGCGGTAACGAAAGAGGATAAACCGTGTCCCGTATTATTATGCTGATCCCTACCGGAACCAGCGTCGGTCTGACCAGCGTCAGCCTTGGCGTGATCCGTGCAATGGAACGCAAAGGCGTTCGTCTGAGCGTTTTCAAACCTATCGCTCAGCCGCGTACCGGt                     <  1:734/297‑1 (MQ=255)
                                                                                                                                                                                                                                                                                      ccTGCGGTGGTTACCCCAACCAACGAAGAACTGGTTATCGCGCAAGACGCGAGCCGCCTGACTGCCTGATTTCACACCGCCAGCTCAGCTGGCGGTGCTGTTTTGTAACCCGCCAAATCGGCGGTAACGAAAGAGGATAAACCGTGTCCCGTATTATTATGCTGATCCCTACCGGAACCAGCGTCGGTCTGACCAGCGTCAGCCTTGGCGTGATCCGTGCAATGGAACGCAAAGGCGTTCGTCTGAGCGTTTTCAAACCTATCGCTCAGCCGCGTACCGGTGGCGATGCGCCCGATCAg   <  1:255944/299‑1 (MQ=255)
                                                                                                                                                                                                                                                                                       cTGCGGTGGTTACCCCAACCAACGAAGAACTGGTTATCGCGCAAGACGCGAGCCGCCTGACTGCCTGATTTCACACCGCCAGCTCAGCTGGCGGTGCTGTTTTGTAACCCGCCAAATCGGCGGTAACGAAAGAGGATAAACCGTGTCCCGTATTATTATGCTGATCCCTACCGGAACCAGCGTCGGTCTGACCAGCGTCAGCCTTGGCGTGATCCGTGCAATGGAACGCAAAGGCGTTCGTCTGAGCGTTTTCAAACCTATCACTCAGCCGCGTACCGGTGGCGATGCGCCCGATCAGa  >  1:560956/1‑299 (MQ=255)
                                                                                                                                                                                                                                                                                                   |                                                                                                                                                                                                                                                                                               
CCGCTATGTTGAAGACAACTACGCGACGAAAGAAGACGCGAAGCGCGCAATGGACGTTTACTGCCACCGCCTGGCGAAATACATCGGTGCCTACACTGCGCTGATGGATGGTCGTCTGGACGCTGTTGTATTCACTGGTGGTATCGGTGAAAATGCCGCGATGGTTCGTGAACTGTCTCTGGGCAAACTGGGCGTGCTGGGCTTTGAAGTTGATCATGAACGCAACCTGGCTGCACGTTTCGGCAAATCTGGTTTCATCAACAAAGAAGGTACCCGTCCTGCGGTGGTTATCCCAACCAACGAAGAACTGGTTATCGCGCAAGACGCGAGCCGCCTGACTGCCTGATTTCACACCGCCAGCTCAGCTGGCGGTGCTGTTTTGTAACCCGCCAAATCGGCGGTAACGAAAGAGGATAAACCGTGTCCCGTATTATTATGCTGATCCCTACCGGAACCAGCGTCGGTCTGACCAGCGTCAGCCTTGGCGTGATCCGTGCAATGGAACGCAAAGGCGTTCGTCTGAGCGTTTTCAAACCTATCGCTCAGCCGCGTACCGGTGGCGATGCGCCCGATCAGA  >  CP006584/2505736‑2506312

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 9 ≤ ATCG/ATCG < 21 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG < 38 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: