Predicted mutation
evidence seq id position mutation annotation gene description
RA CP006584 1,563,642 G→A R308H (CGC→CAC)  LY180_07445 → transposase IS609

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0065841,563,6420GA92.7% 59.3 / 0.8 42R308H (CGC→CAC) LY180_07445transposase IS609
Reads supporting (aligned to +/- strand):  ref base G (3/0);  new base A (22/16);  total (26/16)
Fisher's exact test for biased strand distribution p-value = 2.68e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.80e-01

TGGCGAGACTTCAGCGACAGTTAAGCCGCAAGGTCAAATTCAGCAACAACTGGCAGAAGCAGAAACGCAAAATACAGCGACTGCATTCCTGTATCGCAAATATCCGCAGGGACTACCTTCACAAAGTCACAACGACCGTCAGCAAAAACCACGCAATGATAGTCATTGAGGATTTGAAGGTCAGCAACATGTCAAAGTCAGCAGCGGGTACGGTCAGCCAGCCGGGGCGCAATGTCCGGGCAAAATCAGGTTTAAACCGTTCGATACTGGATCAGGGCTGGTATGAAATGCGCCGCCAGCTTGAGTACAAGCAGCTCTGGCGTGGCGGTCAGGTGCTTGCTGTTCCGCCAGCGTACACAAGCCAGCGTTGCGCGTGCTGTGGTCATACAGCGAAAGAAAATCGCCTGTCACAAAGTAAATTCAGATGCCAGGTATGTGGATATACAGCGAACGCCGATGTAAATGGCGCTCGTAACATTTTAGCGGCGGGGCACGCCGTT‑CTTGCCTGTG‑GAG‑AGATGGTGCAGTCAGGCCGCCCGTTGAAGCAGGAACCCACCGAAATGATTCAGGCGACAGCCTGAACGTAGCA  >  CP006584/1563351‑1563936
                                                                                                                                                                                                                                                                                                   |                                                                                                                                                                                                                                                                                                         
tGGCGAGACTTCAGCGACAGTTAAGCCGCAAGGTCAAATTCAGCAACAACTGGCAGAAGCAGAAACGCAAAATACAGCGACTGCATTCCTGTATCGCAAATATCCGCAGGGACTACCTTCACAAAGTCACAACGGCCGTCAGCAAAAACCACGCAATGATAGTCATTGAGGATTTGAAGGTCAGCAACATGTCAAAGTCAGCAGCGGGTACGGTCAGCCAGCCGGGGCGCAATGTCCGGGCAAAATCAGGTTTAAACCGTTCGATACTGGATCAGGGCTGGTATGAAATGCACCGCCa                                                                                                                                                                                                                                                                                                     >  1:759402/1‑298 (MQ=17)
tGGCGAGACTTCAGCGACAGTTAAGCCGCAAGGTCAAATTCAGCAACAACTGGCAGAAGCAGAAACGCAAAATACAGCGACTGCATTCCTGTATCGCAAATATCCGCAGGGACTACCTTCACAAAGTCACAACGGCCGTCAGCAAAAACCACGCAATGATAGTCATTGAGGATTTGAAGGTCAGCAACATGTCAAAGTCAGCAGCGGGTACGGTCAGCCAGCCGGGGCGCAATGTCCGGGCAAAATCAGGTTTAAACCGTTCGATACTGGATCAGGGCTGGTATGAAATGCACCGCCa                                                                                                                                                                                                                                                                                                     <  1:704403/298‑1 (MQ=17)
    gagaCTTCAGCGACAGTTAAGCCGCAAGGTCAAATTCAGCAACAACTGGCAGAAGCAGAAACGCAAAATACAGCGACTGCATTCCTGTATCGCAAATATCCGCAGGGACTACCTTCACAAAGTCACAACGGCCGTCAGCAAAAACCACGCAATGATAGTCATTGAGGATTTGAAGGTCAGCAACATGTCAAAGTCAGCAGCGGGTACGGTCAGCCAGCCGGGGCGCAATGTCCGGGCAAAATCAGGTTTAAACCGTTCGAAACGGGATCAGGGCTGGTATGAAATGCACCGCCCGCtt                                                                                                                                                                                                                                                                                                 >  2:285723/1‑298 (MQ=12)
      gaCTTCAGCGACAGTTAAGCCGCAAGGTCCAATTCAGCAACAACTGGCCGAAGCAGAAACGCAAAATACAGCGACTGCATTCCTGTATCGCAAATATCCGCAGGGACTACCTTCACAAAGTCACAACGGCCGTCAGCAAAAACCACGCACTGATACTCATTGAGGATTTGCAGGTCAGCAACAAGCCAACGTCATCAGCGGGTACGGTCACCCAGCCGGGGCCCCACTTCCGGGCAAAACCAGGTTTCAACCCTCCGATACTGGATCAGCGCTGGTATAAAATCCACGGCCAGCCTGa                                                                                                                                                                                                                                                                                               >  2:422937/1‑298 (MQ=11)
             gCGACAGTTAAGCCGCAAGGTCAAATTCAGCAACAACTGGCAGAAGCAGAACCGCAAAATACAGCGACTGCATTCCTGTATCGCAACTATCCGCAGGGACTACCTTCACAAAGTCACAACCGCCGTCAGCAAAAACCACGCCATGATAGTCATTGAGGATTTGAAGGTCAGCAACAAGCCAAAGCCAGCAGCGGGGACGGTCAGCCAGGCGGGGCGCAATGTCCGGCAAAACACCGCATTTAACCGGTCCATCCTGGCTCAGGGCCGGTCAGAAATGCACCGCACACTTGAGttaaag                                                                                                                                                                                                                                                                                        >  2:514513/1‑293 (MQ=11)
                    tgaagacaaaaGGTCAAAATCAACAACAAGTAGCAGAAGGAGAAACGAAAAATACAGCGACTGCATTCCTGTATCGAAAAAATCCGCAGGGACTACCTTCGCAAAGTCACAACGGCCGTCAGCAAAAACCACGCAATGATAGTCATGGAGGATTTGAAGGTCAGCAACATGTCAAAGTCAGCAGCGGGTACGGTCAGCCAGCCGGGGCGCAATGTCCGGGCAAAATCAGGTTTAAACCGTTCGATACTGGATCAGGGCTGGTATGAAATGCACCGCCAGCTTGAGTACAAGCAGctct                                                                                                                                                                                                                                                                                 <  2:570092/289‑1 (MQ=11)
                          cacAAGGCAAAAGTAAGCAACAACTGGAGGAAGCAGAAACGCAAAATACAGCCACTGCATTCCTGTATCGGAAATGTCCGCAGGGACTACCTTCACAAAGTCACAACGGCCGTCAGCAAAAACCACGCAAGGATAGTCATTGAGGATTTGAAGGTCAGCAACATGTCAAAGTCAGCAGCGGGTACGGTCAGCCAGCCGGGGCGCAATGTCCGGGCAAAAGCAGGTTTAAACCGTTCGATACTGGATCAGGGCTGGTAGGAAATGCACCGCCAGCTTGAGTACAAGCAGCTCtggcgtg                                                                                                                                                                                                                                                                           <  2:568413/296‑1 (MQ=11)
                           gCAAGGTCAAATTCAGCAACAACTGGCAGAAGCAGAAACGCAAAATACAGCGACTGCATTCCTGTATCGCAAATATCCGCAGGGACTACCTTCACAAAGTCACAACGGCCGTCAGCAAAAACCACGCAATGATAGTCATTGAGGATTTGAAGGTCAGCAACATGTCAAAGTCAGCAGCGGGTACGGTCAGCCAGCCGGGGCGCAATGTCCGGGCAAAATCAGGTTTAAACCGTTCGATACTGGATCAGGGCTGGTATGAAATGCACCGCCAGCTTGAGTACAAGCAGCTCtggcgtgg                                                                                                                                                                                                                                                                          >  1:809944/1‑298 (MQ=17)
                              aGGTCAAATTCAGCAACAACTGGCAGAAGCAGAAACGCAAAATACAGCGACTGCATTCCTGTATCGCAAATATCCGCAGGGACTACCTTCACAAAGTCACAACGGCCGTCAGCAAAAACCACGCAATGATAGTCATTGAGGATTTGAAGGTCAGCAACATGTCAAAGTCAGCAGCGGGTACGGTCAGCCCGCCGGGGCGCACTGTCCGGGCAAACTCAGGTTTAAACCGTTCGATACTGGACCAGGGCTGGTACGAATTGCACCGCCAGCTGGATTACAAGCAGCTCTGGCGTGCCgg                                                                                                                                                                                                                                                                       >  2:790028/1‑298 (MQ=11)
                              aGGTCAAATTCAGCAACAACTGGCAGAAGCAGAAACGCAAAATACAGCGACTGCATTCCTGTATCGCAAATATCCGCAGGGACTACCTTCACAAAGTCACAACGGCCGTCAGCAAAAACCACGCAATGATAGTCATTGAGGATTCGAAGGTCAGCAACATGTCAAAGTCAGCAGCGGGTACGGTCAGCCAGCCGGGGCGCAATGTCCGGGCAAAATCAGGTTTAAACCGTTCGATACTGGATCAGGGCTGGTATGAAATGCACCGCCAGCTTGAGTACAAGCAGCTCTGGCGTGGCgg                                                                                                                                                                                                                                                                       >  1:348719/1‑298 (MQ=14)
                                      ttccaggaagaatGGCAGAAGCAGAAACGAAAAAAACAGCGACTGCATTGCTGTATCGCAAATATCCGCAGGGACTACCTTCACAAAGTCACAACGGCCGTCAGCAAAAACCACGCAATGATAGTCATTGAGGATTTGAAGGTCAGCAACATGTCAAAGTCAGCAGCGGGTACGGTCAGCCAGCCGGGGCGCAATGTCCGGGCAAAATCAGGTTTAAACCGTTCGATACTGGATCAGGGCTGGTATGAAATGCACCGCCAGCTTGAGTACAAGCAGCTCTGGCGTGGCGGTCAGGtgc                                                                                                                                                                                                                                                               <  2:286351/286‑1 (MQ=255)
                                         aGCAACAACTGGCAGAAGCAGAAACGCAAAATACAGCGACTGCATTCCTGTATCGCAAATATCCGCAGGGACTACCTTCACAAAGTCACAACGGCCGTCAGCAAAAACCACGCAATGATAGTCATTGAGGATTTGAAGGTCAGCAACCTGTCAAAGTCAGCAGCGGGTACGGTCAGCCAGCCGGGGCGCAATGTCCGGGCAAAATCAGGTTTAAACCGTCGGATACTGGATCAGGGCTGGGATGAAATGCCCCGCCAGCTTGAGTACAAGCAGCCCTGGCCGGGCGGGCAGGtgcctt                                                                                                                                                                                                                                                            >  2:310525/1‑295 (MQ=11)
                                         aGCAACAACTGGCAGAAGCAGAAACGCAAAATACAGCGACTGCATTCCTGTATCGCAAATATCCGCAGGGACTACCTTCACAAAGTCACAACGGCCGTCAGCAAAAACCACGCAATGATAGTCATTGAGGATTTGAAGGTCAGCAACATGTCAAAGTCAGCAGCGGGTACGGTCAGCCAGCCGGGGCGCAATGTCCGGGCAAAATCAGGTTTAAACCGTTCGATACTGGATCAGGGCTGGTATGAAATGCACCGCCAGCTTGAGTACAAGCAGCTCTGGCGTGGCGGTCAGGtgcttg                                                                                                                                                                                                                                                            >  1:726499/1‑298 (MQ=14)
                                                  tGGCAGAAGCAGAAACGCAAAATACAGCGACTGCATTCCTGTATCGCAAATATCCGCAGGGACTACCTTCACAAAGTCACAACGGCCGTCAGCAAAAACCACGCAATGATAGTCATTGAGGATTTGAAGGTCAGCAACATGTCAAAGTCAGCAGCGGGTACGGTCAGCCAGCCGGGGCGCAATGTCCGGGCAAAATCAGGTTTAAACCGTTCGATACTGGATCAGGGCTGGTATGAAATGCACCGCCAGCTTGAGTACAAGCAGCTCTGGCGTGGCGGTCAGGTGCTTGCTGTTCCGc                                                                                                                                                                                                                                                   <  1:652800/298‑1 (MQ=14)
                                                  tGGCAGAAGCAGAAACGCAAAATACAGCGACTGCATTCCTGTATCGCAAATATCCGCAGGGACTACCTTCACAAAGTCACAACGGCCGTCAGCAAAAACCACGCAATGATAGTCATTGAGGATTTGAAGGTCAGCAACATGTCAAAGTCAGCAGCGGGTACGGTCAGCCAGCCGGGGCGCAATGTCCGGGCAAAATCAGGTTTAAACCGTTCGATACTGGATCAGGGCTGGTATGAAATGCACCGCCAGCTTGAGTACAAGCAGCTCTGGCGTGGCGGTCAGGTGCTTGCTGTTCCGc                                                                                                                                                                                                                                                   >  2:308311/1‑298 (MQ=14)
                                                   ggCAGAAGCAGAAACGCAAAATACAGCGACTGCATTCCTGTATCGCAAATATCCGCAGGGACTACCTTCACAAAGTCACAACGGCCGTCAGCAAAAACCACGCAATGATAGTCATTGAGGATTTGAAGGTCAGCAACATGTCAAAGTCAGCAGCGGGTACGGTCAGCCAGCCGGGGCGCAATGTCCGGGCAAAATCAGGTTTAAACCGTTCGATACTGGATCAGGGCTGGTATGAAATGCACCGCCAGCTTGAGTACAAGCAGCTCTGGCGTGGCGGTCAGGTGCTTGCTGTTCCGcc                                                                                                                                                                                                                                                  >  1:191628/1‑298 (MQ=14)
                                                      agaagcagaaACGCAAAATACAGCGACTGCATTCCTGTATCGCAAATATCCGCAGGGACTACCTTCACAAAGTCACAACGGCCGTCAGCAAAAACCACGCAATGATAGTCATTGAGGATTTGAAGGTCAGCAACATGTCAAAGTCAGCAGCGGGTACGGTCAGCCAGCCGGGGCGCAATGTCCGGGCAAAATCAGGTTTAAACCGTTCGATACTGGATCAGGGCTGGTATGAAATGCACCGCCAGCTTGAGTACAAGCAGCTCTGGCGTGGCGGTCAGGTGCTTGCTGTTCCGCCAGc                                                                                                                                                                                                                                               >  1:332138/1‑298 (MQ=14)
                                                          gcagaaACGCAAAATACAGCGACTGCATTCCTGTATCGCAAATATCCGCAGGGACTACCTTCACAAAGTCACAACGGCCGTCAGCAAAAACCACGCAATGATAGTCATTGAGGATTTGAAGGTCAGCAACATGTCAAAGTCAGCAGCGGGTACGGTCAGCCAGCCGGGGCGCAATGTCCGGGCAAAATCAGGTTTAAACCGTTCGATACTGGATCAGGGCTGGTATGAAATGCACCGCCAGCTTGAGTACAAGCAGCTCTGGCGTGGCGGTCAGGTGCTTGCTGTTCCGCCAGCGTat                                                                                                                                                                                                                                           >  1:409070/1‑297 (MQ=14)
                                                               aaCGCAAAATACAGCGACTGCATTCCTGTATCGCAAATATCCGCAGGGACTACCTTCACAAAGTCACAACGGCCGTCAGCAAAAACCACGCAATGATAGTCATTGAGGATTTGAAGGTCAGCAACATGTCAAAGTCAGCAGCGGGTACGGTCAGCCAGCCGGGGCGCAATGTCCGGGCAAAATCAGGTTTAAACCGTTCGATACTGGATCAGGGCTGGTATGAAATGCACCGCCAGCTTGAGTACAAGCAGCTCTGGCGTGGCGGTCAGGTGCTTGCTGTTCCGCCAGCGTATACAAg                                                                                                                                                                                                                                      >  1:341019/1‑298 (MQ=14)
                                                                 cGCAAAAAACAGCGACTGCATTCCTGTATCGCAAATATCCGCAGGGACTACCTTCACAAAGTCACAACGGCCGTCAGCAAAAACCACGCAATGATAGTCATTGAGGATTTGAAGGTCAGCAACATGTCAAAGTCAGCAGCGGGTACGGTCAGCCAGCCGGGGCGCAATGTCCGGGCAAAATCAGGTTTAAACCGTTCGATACTGGATCAGGGCTGGTATGAAATGCACCGCCAGCTTGAGTACAAGCAGCTCTGGCGTGGCGGTCAGGTGCTTGCTGTTCCGCCAGCGTATACAAGcc                                                                                                                                                                                                                                    <  2:603514/298‑1 (MQ=14)
                                                                   cAAAATACAGCGACTGCATTCCTGTATCGCAAATATCCGCAGGGACTACCTTCACAAAGTCACAACGGCCGTCAGCAAAAACCACGCAATGATAGTCATTGAGGATTTGAAGGTCAGCAACATGTCAAAGTCAGCAGCGGGTACGGTCAGCCAGCCGGGGCGCAATGTCCGGGCAAAATCAGGTTTAAACCGTTCGATACTGGATCAGGGCTGGTATGAAATGCACCGCCAGCTTGAGTACAAGCAGCTCTGGCGTGGCGGTCAGGTGCTTGCTGTTCCGCCAGCGTATACAAGCCAg                                                                                                                                                                                                                                  <  1:308311/298‑1 (MQ=14)
                                                                           aGCGACTGCATTCCTGTATCGCAAATATCCGCAGGGACTACCTTCACAAAGTCACAACGGCCGTCAGCAAAAACCACGCAATGATAGTCATTGAGGATTTGAAGGTCAGCAACATGTCAAAGTCAGCAGCGGGTACGGTCAGCCAGCCGGGGCGCAATGTCCGGGCAAAATCAGGTTTAAACCGTTCGATACTGGATCAGGGCTGGTATGAAATGCACCGCCAGCTTGAGTACAAGCAGCTCTGGCGTGGCGGTCAGGtgctt                                                                                                                                                                                                                                                             <  1:112979/263‑1 (MQ=14)
                                                                           aGCGACTGCATTCCTGTATCGCAAATATCCGCAGGGACTACCTTCACAAAGTCACAACGGCCGTCAGCAAAAACCACGCAATGATAGTCATTGAGGATTTGAAGGTCAGCAACATGTCAAAGTCAGCAGCGGGTACGGTCAGCCAGCCGGGGCGCAATGTCCGGGCAAAATCAGGTTTAAACCGTTCGATACTGGATCAGGGCTGGTATGAAATGCACCGCCAGCTTGAGTACAAGCAGCTCTGGCGTGGCGGTCAGGtgctt                                                                                                                                                                                                                                                             >  2:112979/1‑263 (MQ=14)
                                                                               aCTGCATTCCTGTATCGCAAATATCCGCAGGGACTACCTTCACAAAGTCACAACGGCCGTCAGCAAAAACCACGCAATGATAGTCATTGAGGATTTGAAGGTCAGCAACATGTCAAAGTCAGCAGCGGGTACGGTCAGCCAGCCGGGGCGCAATGTCCGGGCAAAATCAGGTTTAAACCGTTCGATACTGGATCAGGGCTGGTATGAAATGCACCGCCAGCTTGAGTACAAGCAGCTCTGGCGTGGCGGTCAGGTGCTTGCTGTTCCGCCAGCGTATACAAGCCAGCGTTGCGCGtac                                                                                                                                                                                                                      >  1:702604/1‑296 (MQ=14)
                                                                                   cATTCCTGTATCGCAAATATCCGCAGGGACTACCTTCACAAAGTCACAACGGCCGTCAGCAAAAACCACGCAATGATAGTCATTGAGGATTTGAAGGTCAGCAACATGTCAAAGTCAGCAGCGGGTACGGTCAGCCAGCCGGGGCGCAATGTCCGGGCAAAATCAGGTTTAAACCGTTCGATACTGGATCAGGGCTGGTATGAAATGCACCGCCAGCTTGAGTACAAGCAGCTCTGGCGTGGCGGTCAGGTGCTTGCTGTTCCGCCAGCGTATACAAGCCAGCGTTGCGCGTACtgtg                                                                                                                                                                                                                  <  1:39728/298‑1 (MQ=14)
                                                                                     ttCCTGTATCGCAAATACCCGCAGGGACTACCTTCAAAAAGTCACAACGGCCGTCAGCAAAAACCACGCAATGATAGTCATTGAGGATTTGAAGGTCAGCAACATGTCAAAGTCAGCAGCGGGTACGGTCAGCCAGCCGGGGCGCAATGTCCGGGCAAAATCAGGTTTAAACCGTTCGATACTGGATCAGGGCTGGTATGAAATGCACCGCCAGCTTGAGTACAAGCAGCTCTGGCGTGGCGGTCAGGTGCTTGCTGTTCCGCCAGCGTATACAAGCCAGCGTTGCGCGTACTGTGGt                                                                                                                                                                                                                <  1:155431/298‑1 (MQ=9)
                                                                                                                                                                                                                                               gaaaaaTCAGGTTTAAACCGTTCGCAACTGGATCAGGGCTGGTATGAAATGCACCGCCAGCTTGAGTACAAGCAGCTCTGGCGTGGCGGTCAGGGGCTTGCTGTTCCGCCAGCGTATACAAGCCAGCGTTGCGCGTACTGTGGTCATACAGCGAAAGAGAACCGCCTGTCACAAAGTAAATTCAGATGCCAGGTATGTGGATATACAGCGAACGCCGATGTAAATGGCGCTCGTAACATTTTAGCGGCGGGGCACGCCGTT‑CTTGCCTGTG‑GAG‑AGATGGTGCAGTCAGGCCGCCCGt                                                   <  2:723137/296‑1 (MQ=9)
                                                                                                                                                                                                                                               gCAAAATCAGGTTTAAACCGTTCGATACTGGATCAGGGCTGGTATGAAATGCACCGCCAGCTTGAGTACAAGCAGCTCTGGCGTGGCGGTCAGGTGCTTGCTGTTCCGCCAGCGTATACAAGCCAGCGTTGCGCGTACTGTGGTCATACAGCGAAAGAGAACCGCCTGTCACAAAGTAAATTCAGATGCCAGGTATGTGGATATACAGCGAACGCCGATGTAAATGGCGCTCGTAACATTTTAGCGGCGGGGCACGCCGTT‑CTTGCCTGTG‑GAG‑AGATGGTGCAGTCAGGCCGCCCGt                                                   >  1:328821/1‑298 (MQ=9)
                                                                                                                                                                                                                                                cAAAATCAGGTTTAAACCGTTCGATACTGGATCAGGGCTGGTATGAAATGCACCGCCAGCTTGAGTACAAGCAGCTCTGGCGTGGCGGTCAGGTGCTTGCTGTTCCGCCAGCGTATACAAGCCAGCGTTGCGCGTACTGTGGTCATACAGCGAAAGAGAACCGCCTGTCACAAAGTAAATTCAGATGCCAGGTATGTGGATATACAGCGAACGCCGATGTAAATGGCGCTCGTAACATTTTAGCGGCGGGGCACGCCGTT‑CTTGCCTGTG‑GAG‑AGATGGTGCAGTCAGGCCGCCCGtt                                                  >  1:624421/1‑298 (MQ=9)
                                                                                                                                                                                                                                                 aaaaTCAGGTTTAAACCGTTCGATACTGGATCAGGGCTGGTATGAAATGCGCCGCCAGCTTGAGTACAAGCAGCTCTGGCGTGGCGGTCAGGTGCTTGCTGTTCCGCCAGCGTACACAAGCCAGCGTTGCGCGTGCTGTGGTCATACAGCGAAAGAAAATCGCCTGTCACAAAGTAAATTCAGATGCCAGGTATGTGGATATACAGCGAACGCCGATGTAAATGGCGCTCGTAACATTTTAGCGGCGGGGCACGCCGTT‑CTTGCCTGTG‑GAG‑AGATGGTGCAGTCAGGCCGCCCGTTg                                                 >  1:693931/1‑298 (MQ=31)
                                                                                                                                                                                                                                                 aaaaTCAGGTTTAAACCGTTCGATACTGGATCAGGGCTGGTATGAAATGCACCGCCAGCTTGAGTACAAGCAGCTCTGGCGTGGCGGTCAGGTGCTTGCTGTTCCGCCAGCGTATACAAGCCAGCGTTGCGCGTACTGTGGTCATACAGCGAAAGAGAACCGCCTGTCACAAAGTAAATTCAGATGCCAGGTATGTGGATATACAGCGAACGCCGATGTAAATGGCGCTCGTAACATTTTAGCGGCGGGGCACGCCGTT‑CTTGCCTGTG‑GAG‑AGATGGTGCAGTCAGGCCGCCCGTTg                                                 >  1:752932/1‑298 (MQ=9)
                                                                                                                                                                                                                                                   aaTCAGGTTTAAACCGTTCGATACTGGATCAGGGCTGGTATGAAATGCGCCGCCAGCTTGAGTACAAGCAGCTCTGGCGTGGCGGTCAGGTGCTTGCTGTTCCGCCAGCGTACACAAGCCAGCGTTGCGCGTGCTGTGGTCATACAGCGCAAGAAAATCGCCTGTCACAAAGTAAATTACGATGCCAGGTATGTGGATATACAGCGAACGCCGCTGTAAATGGCACTCGTAACATTCTAGCGGCGGGGCCCGCCGTC‑CTT‑CCTGTG‑GGGTGGATGTTGCAGTCCGgcctatcgatgac                                               >  2:246767/1‑288 (MQ=1)
                                                                                                                                                                                                                                                     caaGGTTTAAACGGTTTGATAATGGATCGGGGCTGGAATGAAAGGTACCGCCAGCTTGGGGACAAGCAGCTCTGGCGGGGCGGTCAGGTGCTTGCTGTTCCGCCAGGGTATACAAGCCAGCGTTGGGCGTACTGTGGTCATACAGCGAAAGAGAACCGCCTGTCAAAAAGTAAATTCAGATGCCAGGTATGTGGATATACAGCGAACGCCGATGTAAATGGCGCTCGTAACATTTTAGCGGCGGGGCACGCCGTT‑CTTGCCTGTG‑GAG‑AGATGGTGCAGTCAGGCCGCCCGTTGAAGc                                             <  2:810642/296‑1 (MQ=11)
                                                                                                                                                                                                                                                        ggTTTAAACCGTTCGATACTGGATCAGGGCTGGTATGAAATGCACCGCCAGCTTGAGTACAAGCAGCTCTGGCGTGGCGGTCAGGTGCTTGCTGTTCCGCCAGCGTATACAAGCCAGCGTTGCGCGTACTGTGGTCATACAGCGAAAGAGAACCGCCTGTCACAAAGTAAATTCAGATGCCAGGTATGTGGATATACAGCGAACGCCGATGTAAATGGCGCTCGTAACATTTTAGCGGCGGGGCACGCCGTT‑CTTGCCTGTG‑GAG‑AGATGGTGCAGTCAGGCCGCCCGTTGAAGCAgg                                          >  1:9360/1‑298 (MQ=9)
                                                                                                                                                                                                                                                                   ttCGATACTGGATCAGGGCTGGTATGAAATGCACCGCCAGCTTGAGTACAAGCAGCTCTGGCGTGGCGGTCAGGTGCTTGCTGTTCCGCCAGCGTATACAAGCCAGCGTTGCGCGTACTGTGGTCATACAGCGAAAGAGAACCGCCTGTCACAAAGTAAATTCAGATGCCAGGTATGTGGATATACAGCGAACGCCGATGTAACTGGCGCTCGTAACATTTTAGCGGCGGGGCACCCCCTTCCTCCCCTGTGCGAG‑AG‑TGGTCCGGCCAGCCCGCCCGTTG‑AGCAGGAAccccaccaa                               >  2:797576/1‑292 (MQ=11)
                                                                                                                                                                                                                                                                     cGATACTGGATCAGGGCTGGTATGAAATGCACCGCCAGCTTGAGTACAAGCAGCTCTGGCGTGGCGGTCAGGTGCTTGCTGTTCCGCCAGCGTATACAAGCCAGCGTTGCGCGTACTGTGGTCATACAGCGAAAGAGAACCGCCTGTCACAAAGTAAATTCAGATGCCAGGTATGTGGATATACAGCGAACGCCGATGTAAATGGCGCTCGTAACATTTTAGCGGCGGGGCACGCCGTT‑CTTGCCTGTG‑GAG‑AGATGGTGCAGTCAGGCCGCCCGTTGAAGCAGGAACCCACCGAAAt                             <  2:332138/298‑1 (MQ=9)
                                                                                                                                                                                                                                                                                         tATGAAATGCACCGCCAGCTTGAGTACAAGCAGCTCTGGCGTGGCGGTCAGGTGCTTGCTGTTCCGCCAGCGTATACAAGCCAGCGTTGCGCGTACTGTGGTCATACAGCGAAAGAGAACCGCCTGTCACAAAGTAAATTCAGATGCCAGGTATGTGGATATACAGCGAACGCCGATGTAAATGGCGCTCGTAACATTTTAGCGGCGGGGCACGCCGTT‑CTTGCCTGTG‑GAG‑AGATGGTGCAGTCAGGCCGCCCGTTGAAGCAGGAACCCACCGAAATGATTCAGGCGACAGCCTGaa         >  1:633976/1‑298 (MQ=9)
                                                                                                                                                                                                                                                                                          aTGAAATGCGCCGCCAGCTTGAGTACAAGCAGCTCTGGCGTGGCGGTCAGGTGCTTGCTGTTCCGCCAGCGTACACAAGCCAGCGTTGCGCGTGCTGTGGTCATACAGCGAAAGAAAATCGCCTGTCACAAAGTAAATTCAGATGCCAGGTAGGTGGATATACAGCGAACGCCGATGTAAATGGCGCTCGTACCATTTTAGCGGCCGGGCACGCCGGT‑CTCTCCTGTG‑CAG‑AGATGGTGCAGTCAGGCCGCCCGTTGCAGCAGTAACCCACCCACATAGTTAAGGCGACAGCCTGaaa        >  2:22220/1‑297 (MQ=11)
                                                                                                                                                                                                                                                                                            gAAATGCACCGCCAGCTTGAGTACAAGCAGCTCTGGCGTGGCGGTCAGGTGCTTGCTGTTCCGCCAGCGTATACAAGCCAGCGTTGCGCGTACTGTGGTCATACAGCGAAAGAGAACCGCCTGTCACAAAGTAAATTCAGATGCCAGGTATGTGGATATACAGCGAACGCCGATGTAAATGGCGCTCGTAACATTTTAGCGGCGGGGCACGCCGTT‑CTTGCCTGTG‑GAG‑AGATGGTGCAGTCAGGCCGCCCGTTGAAGCAGGAACCCACCGAAATGATTCAGGCGACAGCCTGAACGt      <  1:322359/298‑1 (MQ=9)
                                                                                                                                                                                                                                                                                             aaaTGCACCGCCAGCTTGAGTACAAGCAGCTCTGGCGTGGCGGTCAGGTGCTTGCTGTTCCGCCAGCGTATACAAGCCAGCGTTGCGCGTACTGTGGTCATACAGCGAAAGAGAACCGCCTGTCACAAAGTAAATTCAGATGCCAGGTATGTGGATATACAGCGAACGCCGATGTAAATGGCGCTCGTAACATTTTAGCGGCGGGGCACGCCGTT‑CTTGCCTGTG‑GAG‑AGATGGTGCAGTCAGGCCGCCCGTTGAAGCAGGAACCCACCGAAATGATTCAGGCGACAGCCTGAACGTa     >  1:215718/1‑298 (MQ=9)
                                                                                                                                                                                                                                                                                              aaTGCACCGCCAGCTTGAGTACAAGCAGCTCTGGCGTGGCGGTCAGGTGCTTGCTGTTCCGCCAGCGTATACAAGCCAGCGTTGCGCGTACTGTGGTCATACAGCGAAAGAGAACCGCCTGTCACAAAGTAAATTCAGATGCCAGGTATGTGGATATACAGCGAACGCCGATGTAAATGGCGCTCGTAACATTTTAGCGGCGGGGCACGCCGTT‑CTTGCCTGTG‑GAG‑AGATGGTGCAGTCAGGCCGCCCGTTGAAGCAGGAACCCACCGAAATGATTCAGGCGACAGCCTGAACGTAg    <  1:454285/298‑1 (MQ=12)
                                                                                                                                                                                                                                                                                                tGCACCGCCAGCTTGAGTACAAGCAGCTCTGGCGTGGCGGTCAGGTGCTTGCTGTTCCGCCAGCGTATACAAGCCAGCGTTGCGCGTACTGTGGTCATACAGCGAAAGAGAACCGCCTGTCACAAAGTAAATTCAGATGCCAGGTATGTGGATATACAGCGAACGCCGATGTAAATGGCGCTCGTAACATTTTAGCGGCGGGGCACGCCGTT‑CTTGCCTGTG‑GAG‑AGATGGTGCAGTCAGGCCGCCCGTTGAAGCAGGAACCCACCGAAATGATTCAGGCGACAGCCTGAACGTAGCa  <  1:750889/298‑1 (MQ=12)
                                                                                                                                                                                                                                                                                                   |                                                                                                                                                                                                                                                                                                         
TGGCGAGACTTCAGCGACAGTTAAGCCGCAAGGTCAAATTCAGCAACAACTGGCAGAAGCAGAAACGCAAAATACAGCGACTGCATTCCTGTATCGCAAATATCCGCAGGGACTACCTTCACAAAGTCACAACGACCGTCAGCAAAAACCACGCAATGATAGTCATTGAGGATTTGAAGGTCAGCAACATGTCAAAGTCAGCAGCGGGTACGGTCAGCCAGCCGGGGCGCAATGTCCGGGCAAAATCAGGTTTAAACCGTTCGATACTGGATCAGGGCTGGTATGAAATGCGCCGCCAGCTTGAGTACAAGCAGCTCTGGCGTGGCGGTCAGGTGCTTGCTGTTCCGCCAGCGTACACAAGCCAGCGTTGCGCGTGCTGTGGTCATACAGCGAAAGAAAATCGCCTGTCACAAAGTAAATTCAGATGCCAGGTATGTGGATATACAGCGAACGCCGATGTAAATGGCGCTCGTAACATTTTAGCGGCGGGGCACGCCGTT‑CTTGCCTGTG‑GAG‑AGATGGTGCAGTCAGGCCGCCCGTTGAAGCAGGAACCCACCGAAATGATTCAGGCGACAGCCTGAACGTAGCA  >  CP006584/1563351‑1563936

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