Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0065841,563,2370TC55.9% ‑6.4 / 37.1 34V173A (GTT→GCT) LY180_07445transposase IS609
Reads supporting (aligned to +/- strand):  ref base T (8/7);  new base C (8/11);  total (16/18)
Fisher's exact test for biased strand distribution p-value = 7.30e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.54e-01
Rejected as consensus: E-value score below prediction cutoff.
Rejected as consensus: Frequency below/above cutoff threshold.

GGAAAACAGCCGTATTTTTCTGCCGAAACTGGGCTGGATGCGCTACCGGAATAGCCGTCAGGTCACGGGTGTTGTGAAAAATGTCACTGTCAGCCAGTCCTGCGGTAAGTGGTACATCAGTATTCAGACAGAAAGTGAAGTATCCACTCCGGTTCACCCTTCAGCATCAATGGTCGGGCTGGATGCTGGCGTGGCTAAACTCGCCACGCTGTCAGATGGCACAGTCTTTGAGCCTGTAAACAGTTTTCAGAAAAACCAGAAGACGCTGGCGAGACTTCAGCGACAGTTAAGCCGCAAGGTCAAATTCAGCAACAACTGGCAGAAGCAGAAACGCAAAATACAGCGACTGCATTCCTGTATCGCAAATATCCGCAGGGACTACCTTCACAAAGTCACAACGACCGTCAGCAAAAACCACGCAATGATAGTCATTGAGGATTTGAAG  >  CP006584/1563085‑1563529
                                                                                                                                                        |                                                                                                                                                                                                                                                                                                    
ggAAAACAGCCGTATTTTTCTGCCGAAACTGGGCTGGATGCGCTACCGGAACAGCCGTCAGGTCACGGGTGTTGTGAAAAATGTCACTGTCAGCCAGTCCTGCGGTAAGTGGTACATCAGTATTCAGACAGAAAGTGAAGTATCCACTCCGGCTCACCCTTCAGCATCAATGGTCGGGCTGGATGCTGGCGTGGCTAAACTCGCCACGCTGTCAGATGGCACAGTCTTTGAGCCTGTAAACAGTTTTCAGAAAAACCAGAAGACGCTGGCGAGACTTCAGCGACAGTTAAGCCGCAAg                                                                                                                                                     >  1:277571/1‑298 (MQ=11)
 gAAATCAGCCGTATTTTTCTGCCGGAGCGGGGCTGGATGCGCTACCGGAATAGCCGTCAGGTCACGGGTGTTGTGAAAAATGTCACTGTCAGCCAGTCCTGCGGTAAGTGGTACATCAGTATTCAGACAGAAAGTGAAGTATCGACTCCGGTTCACCCTTCAGCATCAATGGTCGGGCTGGATGCTGGCGTGGCTAAACTCGCCACGCTGTCAGATGGCACAGTCTTTGAGCCTGTAAACAGTTTTCAGAAAAACCAGAAGACGCTGGCGAGACTTCAGCGACAGTTAAGCCGCAAgg                                                                                                                                                    <  2:384079/298‑1 (MQ=9)
      cAGCCGTATTTTTCTGCCGAAACTGGGCTGGATGCGCTACCGGAATAGCCGTCAGGTCACGGGTGTTGTGAAAAATGTCACTGTCAGCCAGTCCTGCGGTAAGTGGTACATCAGTATTCAGACAGAAAGTGAAGTATCCACTCCGGTTCACCCTTCAGCATCAATGGTCGGGCTGGATGCTGGCGTGGCTAAACTCGCCACGCTGTCAGATGGCACAGTCTTTGAGCCTGTAAACAGTTTTCAGAAAAACCAGAAGACGCTGGCGAGACTTCAGCGACAGTTAAGCCGCAAGGTCaaa                                                                                                                                               <  1:339219/298‑1 (MQ=31)
        gCCGTATTTTTCTGCCGAAACTGGGCTGGATGCGCTACCGGAATAGCCGTCAGGTCACGGGTGTTGTGAAAAATGTCACTGTCAGCCAGTCCTGCGGTAAGTGGTACATCAGTATTCAGACAGAAAGTGAAGTATCCACTCCGGTTCACCCTTCAGCATCAATGGTCGGGCTGGATGCTGGCGTGGCTAAACTCGCCACGCTGTCAGATGGCACAGTCTTTGAGCCTGTAAACAGTTTTCAGAAAAACCAGAAGACGCTGGCGAGACTTCAGCGACAGTTAAGCCGCTAGGTCAAAtt                                                                                                                                             >  1:293429/1‑298 (MQ=14)
         ccGTATTTTTCTGCCGAAGCGGGGCTGTATGCGCTACCGGACCGGCCGTCAGGTCACGGGTGTTGTGAAAAATGTCACTGTCAGCCAGTCCTGCGGTAAGTGGTACATCAGTATTCAGACAGAAAGTGAAGTATCCACTCCGGCTCACCCTTCAGCATCAATGGTCGGGCTGGATGCTGGCGTGGCTAAACTCGCCACGCTGTCAGATGGCACAGTCTTTGAGCCTGTAAACAGTTTTCAGAAAAACCAGAAGACGCTGGCGAGACTTCAGCGACAGTTAAGCCGCAAGGTCAAATTc                                                                                                                                            <  1:184163/298‑1 (MQ=9)
            tATTTTTCTGCCGAAACTGGGCTGGATGCGCTACCGGAATAGCCGTCAGGTCACGGGTGTTGTGAAAAATGTCACTGTCAGCCAGTCCTGCGGTAAGTGGTACATCAGTATTCAGACAGAAAGTGAAGTATCCACTCCGGTTCACCCTTCAGCATCAATGGTCGGGCTGGATGCTGGCGTGGCTAAACTCGCCACGCTGTCAGATGGCACAGTCTTTGAGCCTGTAAACAGTTTTCAGAAAAACCAGAAGACGCTGGCGAGACTTCAGCGACAGTTAAGCCGCAAGGTCAAATTCAGc                                                                                                                                         >  1:112867/1‑298 (MQ=31)
                    tGCCGAAACTGGGCTGGATGCGCTACCGGAACAGCCGTCAGGTCACGGGTGTTGTGAAAAATGTCACTGTCAGCCAGTCCTGCGGTAAGTGGTACATCAGTATTCAGACAGAAAGTGAAGTATCCACTCCGGCTCACCCTTCAGCATCAATGGTCGGGCTGGATGCTGGCGTGGCTAAACTCGCCACGCTGTCAGATGGCACAGTCTTTGAGCCTGTAAACAGTTTTCAGAAAAACCAGAAGACGCTGGCGAGACTTCAGCGACAGTTAAGCCGCAAGGTCAAATTCAGCAACAACTg                                                                                                                                 <  1:328443/298‑1 (MQ=14)
                     gCCGAAACTGGGCTGGATGCGCGACCGGAACAGCCGTCAGGTCACGGGTGTTGTGAAAAATGTCACTGTCAGCCAGTCCTGCGGTAAGTGGTACATCAGTATTCAGACAGAAAGTGAAGTATCCACTCCGGCTCACCCTTCAGCATCAATGGTCGGGCTGGATGCTGGCGTGGCTAAACTCGCCACGCTGTCAGATGGCACAGTCTTTGAGCCTGTAAACAGTTTTCTGAAAAACCAGAAGACGCTGGCGAGACTTCAGCGACAGTTAAGCCGCAAGGTCAAATTCAGCAACAACTgg                                                                                                                                <  2:350932/298‑1 (MQ=14)
                           aCTGGGCTGGATGCGCTACCGGAATAGCCGTCAGGTCACGGGTGTTGTGAAAAATGTCACTGTCAGCCAGTCCTGCGGTAAGTGGTACATCAGTATTCAGACAGAAAGTGAAGTATCCACTCCGGTTCACCCTTCAGCATCAATGGTCGGGCTGGATGCTGGCGTGGCTAAACTCGCCACGCTGTCAGATGGCACAGTCTTTGAGCCTGTAAACAGTTTTCAGAAAAACCAGAAGACGCTGGCGAGACTTCAGCGACAGTTAAGCCGCAAGGTCAAATTCAGCAACAACTGgcagaag                                                                                                                          >  1:211270/1‑298 (MQ=31)
                                gCTGGATGGGCTACCGGAATAGCCTTCAGGTTACGGGTGTTGTGAAAAATGTCACTGTCAGCCAGTCCTGCGGTAAGTGGTACATCAGTATTCAGACAGAAAGTGAAGTATCCACTCCGGTTCACCCTTCAGCATCGATGGTCGGGCTGGATGCTGGCGTGGCTAAACTCGCCACGCTGTCAGATGGCACAGTCTTTGAGCCTGTAAACAGTTTTCAGAAAAACCAGAAGACGCTGGCGAGACTTCAGCGACAGTTAAGCCGCAAGGTCAAATTCAGCAACAACTGgcagaagcagaa                                                                                                                     <  2:402709/298‑1 (MQ=14)
                                        cgcTACCGGTATAGCCGTCAGGTCCAGGGTGTTGTGAAAAATGTCACTGTCAGCCAGTCCTGCGGTAAGTGGTACATCAGTATTCAGACAGAAAGTGAAGTATCCACTCCGGTTCACCCTTCGGCAGCAATGGTCGGGCTGGATGCTGGCGTGGCTAAACTCGCCACGCTGTCAGATGGCACAGTCTTTGAGCCTGTAAACAGTTTTCAGAAAAACCAGAAGACGCTGGCGAGACTTCAGCGACAGTTAAGCCGCAAGGTCAAATTCAGCAACAACTGGCAGAAGCAGAAACGCaaaa                                                                                                             <  2:154241/298‑1 (MQ=9)
                                          cTACCGGAACAGCCGTCAGGTCACGGGTGTTGTGAAAAATGTCACTGTCAGCCAGTCCTGCGGTAAGTGGTACATCAGTATTCAGACAGAAAGTGAAGTATCCACTCCGGCTCACCCTTCAGCATCAATGGTCGGGCTGGATGCTGGCGTGGCTAAACTCGCCACGCTGTCAGATGGCACAGTCTTTGAGCCTGTAAACAGTTTTCAGAAAAACCAGAAGACGCTGGCGAGACTTCAGCGACAGTTAAGCCGCAAGGTCAAATTCAGCAACAACTGGCAGAAGCAGAAACGCAAAATa                                                                                                           >  1:169588/1‑298 (MQ=14)
                                               ggAATAGCCGTCAGGTCACGGGTGTTGTGAAAAATGTCACTGTCAGCCAGTCCTGCGGTAAGTGGTACATCAGTATTCAGACAGAAAGTGAAGTATCCACTCCGGTTCACCCTTCAGCATCAATGGTCGGGCTGGATGCTGGCGTGGCTAAACTCGCCACGCTGTCAGATGGCACAGTCTTTGAGCCTGTAAACAGTTTTAAGAAAAACCAGAAGACGCTGGCGAGACTTCAGCGACAGTTAAGCCGCAAGGGCAAATTCAGCAAAAACTGGCAGAAGCAGAAACGCAACATACAGcc                                                                                                      >  2:332644/1‑297 (MQ=14)
                                                gAACAGCCGTCAGGTCACGGGTGTTGTGAAAAATGTCACTGTCAGCCAGTCCTGCGGTAAGTGGTACATCAGTATTCAGACAGAAAGTGAAGTATCCACTCCGGCTCACCCTTCAGCATCAATGGTCGGGCTGGATGCTGGCGTGGCTAAACTCGCCACGCTGTCAGATGGCACAGTCTTTGAGCCTGTAAACAGTTTTCAGAAAAACCAGAAGACGCTGGCGAGACTTCAGCGACAGTTAAGCCGCAAGGTCAAATTCAGCAACAACTGGCAGAAGCAGAAACGCAAAATACAGCGa                                                                                                     <  2:131481/298‑1 (MQ=14)
                                                  aTAGCCGTCAGGTCACGGGTGTTGTGAAAAATGTCACTGTCAGCCAGTCCTGCGGTAAGTGGTACATCAGTATTCAGACAGAAAGTGAAGTATCCACTCCGGTTCACCCTTCAGCATCAATGGTCGGGCTGGATGCTGGCGTGGCTAAACTCGCCACGCTGTCAGATGGCACAGTCTTTGAGCCTGTAAACAGTTTTCAGAAAAACCAGAAGACGCTGGCGAGACTTCAGCGACAGTTAAGCCGCAAGGTCAAATTCAGCAACAACTGGCAGAAGCAGAAACGCAAAATACAGCGACt                                                                                                   <  1:309577/298‑1 (MQ=31)
                                                   tAGCCGTCAGGTCACGGGTGTTGTGAAAAATGTCACTGTCAGCCAGTCCTGCGGTAAGTGGTACATCAGTATTCAGACAGAAAGTGAAGTATCCACTCCGGTTCACCCTTCAGCATCAATGGTCGGGCTGGATGCTGGCGTGGCTAAACTCGCCACGCTGTCAGATGGCACAGTCTTTGAGCCTGTAAACAGTTTTCAGAAAAACCAGAAGACGCTGGCGAGACTTCAGCGACAGTTAAGCCGCAAGGTCAAATTCAGCAACAACTGGCAGAAGCAGAAACGCAAAATACAGCGACTg                                                                                                  >  1:315958/1‑298 (MQ=31)
                                                      ccgtggtgTCACGGGTGTTGGAAAAAATGTCACTGTCAGCCAGTCCTGCGGTAAGTGGTACATCAGTATTCAGAAAGAAAGTGAAGTATCCACTACGGCTCACCCTTCAGCATCAATGGTCGGGCTGGATGCTGGCGTGGCTAAACTCGCCACGCTGTCAGATGGCACAGTCTTTGAGCCTGTAAACAGTTTTCAGAAAAACCAGAAGACGCTGGCGAGACTTCAGCGACAGTTAAGCCGCAAGGTCAAATTCAGCAACAACTGGCAGAAGCAGAAACGCAAAATACAGCGACTGCAt                                                                                               <  2:369744/291‑1 (MQ=9)
                                                            ggTCACGGGTGTTGTGAAAAATGTCACTGTCAGCCAGTCCTGCGGTAAGTGGTACATCAGTATTCAGACAGAAAGTGAAGTATCCACTCCGGCTCACCCTTCAGCATCAATGGTCGGGCTGGATGCTGGCGTGGCTAAACTCGCCACGCTGTCAGATGGCACAGTCTTTGAGCCTGTAAACAGTTTTCAGAAAAACCAGAAGACGCTGGCGAGACTTCAGCGACAGTTAAGCCGCAAGGTCAAATTCAGCAACAACTGGCAGAAGCAGAAACGCAAAATACAGCGACTGCATTCCTGt                                                                                         >  1:64937/1‑298 (MQ=17)
                                                                  gggTGTTGTGAAAAATGTCACTGTCAGCCAGTCCTGCGGTAAGTGGTACATCAGTATTCAGACAGAAAGTGAAGTATCCACTCCGGTTCACCCTTCAGCATCAATGGTCGGGCTGGATGCTGGCGTGGCTAAACTCGCCACGCTGTCAGATGGCACAGTCTTTGAGCCTGTAAACAGTTTTCAGAAAAACCAGAAGACGCTGGGGAGACTTCAGCGACAGTTAAGCCGCAAGGTCAAATTCAGCAACAACTGGCAGAAGCAGAAACGCAAAATACAGCGACTGCATTcc                                                                                            >  2:112210/1‑289 (MQ=14)
                                                                  gggTGTTGTGAAAAATGTCACTGTCAGCCAGTCCTGCGGTAAGTGGTACATCAGTATTCAGACAGAAAGTGAAGTATCCACTCCGGTTCACCCTTCAGCATCAATGGTCGGGCTGGATGCTGGCGTGGCTAAACTCGCCACGCTGTCAGATGGCACAGTCTTTGAGCCTGTAAACAGTTTTCAGAAAAACCAGAAGACGCTGGGGAGACTTCAGCGACAGTTAAGCCGCAAGGTCAAATTCAGCAACAACTGGCAGAAGCAGAAACGCAAAATACAGCGACTGCATTcc                                                                                            <  1:112210/289‑1 (MQ=14)
                                                                   ggTGTTGTGAAAAATGTCACTGTCAGCCAGTCCTGCGGTAAGTGGTACATCAGTATTCAGACAGAAAGTGAAGTATCCACTCCGGCTCACCCTTCAGCATCAATGGTCGGGCTGGATGCTGGCGTGGCTAAACTCGCCACGCTGTCAGATGGCACAGTCTTTGAGCCTGTAAACAGTTTTCAGAAAAACCAGAAGACGCTGGCGAGACTTCAGCGACAGTTAAGCCGCAAGGTCAAATTCAGCAACAACTGGCAGAAGCAGAAACGCAAAATAAAGCGACGGCATTCCttgaacgaaa                                                                                  >  2:210771/1‑289 (MQ=9)
                                                                    gtgtTGTGAAAAATGTCACTGTCAGCCAGTCCTGCGGTAAGTGGTACATCAGTATTCAGACAGAAAGTGAAGTATCCACTCCGGCTCACCCTTCAGCATCAATGGTCGGGCTGGATGCTGGCGTGGCTAAACTCGCCACGCTGTCAGATGGCACAGTCTTTGAGCCTGTAAACAGTTTTCAGAAAAACCAGAAGACGCTGGCGAGACTTCAGCGACAGTTAAGCCGCAAGGTCAAATTCAGCAACAACTGGCAGAAGCAGAAACGCAAAATACAGCGACTGCATTCCTGTATCGCaaa                                                                                 <  1:210771/298‑1 (MQ=17)
                                                                          taaaaaatttacatGTCAGCCAGTCCTGGGATAAGTGGTACATCAGTATTCAGACAGAAAGGGAAGTATCCACTCCGGCTCACCCTTCAGCATCAATGGTCGGGCTGGATGCTTGGGTGGCTAAACTCGCCACGCTGTCAGATGGCACAGTCTTTGAGCCTGTAAACAGTTTTCAGAAAAACCAGAAGACGCTGGCGAGACTTCAGCGACAGTTAAGCCGCAAGGTCAAATTCAGCAACAACTGGCAGAAGCAGAAACGCAAAATACAGCGACTGCATTCCTGTATCGCAAATATCCg                                                                           <  2:242609/285‑1 (MQ=11)
                                                                           gAAAAATGTCACTGTCAGCCAGTCCTGCGGTAAGTGGTACATCAGTATTCAGACAGAAAGTGAAGTATCCACTCCGGCTCACCCTTCAGCATCAATGGTCGGGCTGGATGCTGGCGTGGCTAAACTCGCCACGCTGTCAGATGGCACAGTCTTTGAGCCTGTAAACAGTTTTCAGAAAAACCAGAAGACGCTGGCGAGACTTCAGCGACAGTTAAGCCGCAAGGTCAAATTCAGCAACAACTGGCAGAAGCAGAAACGCAAAATACAGCGACTGCATTCCTGTATCGCAAATATCCGc                                                                          >  1:242609/1‑298 (MQ=17)
                                                                                       tGTCAGCCAGTCCTGCGGTAAGTGGTACATCAGTATTCAGACAGAAAGTGAAGTATCCACTCCGGCTCACCCTTCAGCATCAATGGTCGGGCTGGATGCTGGCGTGGCTAAACTCGCCACGCTGTCAGATGGCACAGTCTTTGAGCCTGTAAACAGTTTTCAGAAAAACCAGAAGACGCTGGCGAGACTTCAGCGACAGTTAAGCCGCAAGGTCAAATTCAGCAACAACTGGCAGAAGCAGAAACGCAAAATACAGCGACTGCATTCCTGTATCGCAAATATCCGCAGGGACTACCtt                                                              <  1:79419/298‑1 (MQ=17)
                                                                                                  ccTGCGGTAAGTGGTACATCAGTATTCAGACAGAAAGTGAAGTATCCACTCCGGCTCACCCTTCAGCATCAATGGTCGGGCTGGATGCTGGCGTGGCTAAACTCGCCACGCTGTCAGATGGCACAGTCTTTGAGCCTGTAAACAGTTTTCAGAAAAACCAGAAGACGCTGGCGAGACTTCAGCGACAGTTAAGCCGCAAGGTCAAATTCAGCAACAACTGGCAGAAGCAGAAACGAAAAATACAGCGACTGCCTTCCTGTATCGCAAATATCCCCAGGCACTACCTCCACaaatcccc                                                   >  2:102948/1‑293 (MQ=12)
                                                                                                                         aTTCAGAGAGAAAGTGAATGATCCACTCCGGCTCACCCTTCCGCATCAATGGTCGGGCTGGATGCTGGCGTGGCTAAACTCGCCACGCGGTCAGATGGCACAGTCTGTGAGCCTGTAAACAGTTTTCAGAAAAACCAGAAGACGCTGGCGAGACTTCAGCGACAGTTAAGCCGCAAGGTCAAATTCAGCAACAACTGGCAGAAGCAGAAACGCAAAATACAGCGACTGCATTCCTGTATCGCAAATATCCGCAGGGACTACCTTCACAAAGTCACAACGGCCGTCAGCAAAAACCACg                            <  2:409845/298‑1 (MQ=255)
                                                                                                                          ttCAGACAGAAAGTGAAGTATCCACTCCGGTTCACCCTTCAGCATCAATGGTCGGGCTGGATGCTGGCGTGGCTAAACTCGCCACGCTGTCAGATGGCACAGTCTTTGAGCCTGTAAACAGTTTTCAGAAAAACCAGAAGACGCTGGCGAGACTTCAGCGACAGTTAAGCCGCAAGGTCAAATTCAGCAACAACTGGCAGAAGCAGAAACGCAAAATACAGc                                                                                                       >  1:313268/1‑222 (MQ=31)
                                                                                                                          ttCAGACAGAAAGGGAAGTATCCACTCCGGTTCACCCTTCAGCATCAATGGTCGGGCTGGATGCTGGCGTGGCTAAACTCGCCACGCTGTCAGATGGCACAGTCTTTGAGCCTGTAAACAGTTTTCAGAAAAACCAGAAGACGCTGGCGAGACTTCAGCGACAGTTAAGCCGCAAGGTCAAATTCAGCAACAACTGGCAGAAGCAGAAACGCAAAATACAGc                                                                                                       <  2:313268/222‑1 (MQ=14)
                                                                                                                               acagaAAGTGAAGTATCCACTCCGGCTCACCCTTCAGCATCAATGGTCGGGCTGGATGCTGGCGTGGCTAAACTCGCCACGCTGTCAGATGGCACAGTCTTTGAGCCTGTAAACAGTTTTCAGAAAAACCAGAAGACGCTGGCGAGACTTCAGCGACAGTTAAGCCGCAAGGTCAAATTCAGCAACAACTGGCAGAAGCAGAAACGCAAAATACAGCGACTGCATTCCTGTATCGCAAATATCCGCAGGGACTACCTTCACAAAGTCACAACGGCCGTCAGCAAAAACCACGCAATGa                      <  1:320189/298‑1 (MQ=17)
                                                                                                                                    aaGTGAAGTATCCACTCCGGCTCACCCTTCAGCATCAATGGTCGGGCTGGATGCTGGCGTGGCTAAACTCGCCACGCTGTCAGATGGCACAGTCTTTGAGCCTGTAAACAGTTTTCAGAAAAACCAGAAGACGCTGGCGAGACTTCAGCGACAGTTAAGCCGCAAGGTCAAATTCAGCAACAACTGGCAGAAGCAGAAACGCAAAATACAGCGACTGCATTCCTGTATCGCAAATATCCGCAGGGACTACCTTCACAAAGTCACAACGGCCGTCAGCAAAAACCACGCAATGATAGTc                 >  1:95410/1‑298 (MQ=17)
                                                                                                                                         aaGTATCCACTCCGGCTCACCCTTCAGCATCAATGGTCGGGCTGGATGCTGGCGTGGCTAAACTCGCCACGCTGTCAGATGGCACAGTCTTTGAGCCTGTAAACAGTTTTCAGAAAAACCAGAAGACGCTGGCGAGACTTCAGCGACAGTTAAGCCGCAAGGTCAAATTCAGCAACAACTGGCAGAAGCAGAAACGCAAAATACAGCGACTGCATTCCTGTATCGCAAATATCCGCAGGGACTACCTTCACAAAGTCACAACGGCCGTCAGCAAAAACCACGCAATGATAGTCATTGa            <  1:330133/298‑1 (MQ=17)
                                                                                                                                                cACTCCGGTTCACCCTTCAGCATCAATGGTCGGGCTGGATGCTGGCGTGGCTAAACTCGCCACGCTGTCAGATGGCACAGTCTTTGAGCCTGTAAACAGTTTTCAGAAAAACCAGAAGACGCTGGCGAGACTTCAGCGACAGTTAAGCCGCAAGGTCAAATTCAGCAACAACTGGCAGAAGCAGAAACGCAAAATACAGCGACTGCATTCCTGTATCGCAAATATCCGCAGGGACTACCTTCACAAAGTCACAACGGCCGTCAGCAAAAACCACGCAATGATAGTCATTGAGGATTTg     >  1:281868/1‑298 (MQ=17)
                                                                                                                                                   tCCGGCTCACCCTTCAGCATCAATGGTCGGGCTGGATGCTGGCGTGGCTAAACTCGCCACGCTGTCAGATGGCACAGTCTTTGAGCCTGTAAACAGTTTTCAGAAAAACCAGAAGACGCTGGCGAGACTTCAGCGACAGTTAAGCCGCAAGGTCAAATTCAGCAACAACTGGCAGAAGCAGAAACGCAAAATACAGCGACTGCATTCCTGTATCGCAAATATCCGCAGGGACTACCTTCACAAAGTCACAACGGCCGTCAGCAAAAACCACGCAATGATAGTCATTGAGGATTTGAAg  >  1:306414/1‑298 (MQ=17)
                                                                                                                                                        |                                                                                                                                                                                                                                                                                                    
GGAAAACAGCCGTATTTTTCTGCCGAAACTGGGCTGGATGCGCTACCGGAATAGCCGTCAGGTCACGGGTGTTGTGAAAAATGTCACTGTCAGCCAGTCCTGCGGTAAGTGGTACATCAGTATTCAGACAGAAAGTGAAGTATCCACTCCGGTTCACCCTTCAGCATCAATGGTCGGGCTGGATGCTGGCGTGGCTAAACTCGCCACGCTGTCAGATGGCACAGTCTTTGAGCCTGTAAACAGTTTTCAGAAAAACCAGAAGACGCTGGCGAGACTTCAGCGACAGTTAAGCCGCAAGGTCAAATTCAGCAACAACTGGCAGAAGCAGAAACGCAAAATACAGCGACTGCATTCCTGTATCGCAAATATCCGCAGGGACTACCTTCACAAAGTCACAACGACCGTCAGCAAAAACCACGCAATGATAGTCATTGAGGATTTGAAG  >  CP006584/1563085‑1563529

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