Predicted mutation
evidence seq id position mutation annotation gene description
RA CP006584 1,964,559 G→A N226N (AAC→AAT)  LY180_09345 ← transposase IS609

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0065841,964,5590GA95.0% 66.1 / ‑3.5 42N226N (AAC→AAT) LY180_09345transposase IS609
Reads supporting (aligned to +/- strand):  ref base G (2/0);  new base A (19/19);  total (22/20)
Fisher's exact test for biased strand distribution p-value = 4.88e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.97e-01

GTACTCAAGCTGGCGGCGCATTTCATACCAGCCCTGATCCAGTATCGAACGGTTTAAACCTGATTTTGCCCGGACATTGCGCCCCGGCTGGCTGACCGTACCCGCTGCTGACTTTGACATGTTGCTGACCTTCAAATCCTCAATGACTATCATTGCGTGGTTTTTGCTGACGGTCGTTGTGACTTTATGAAGGTAGTCCCTGCGGATATTTGCGATACGGGAATGCAGTCGCTGTATTTTGCGTTTCTGCTTCTGCCAGTTGTTGCTGAATTTGACCTTGCGGCTTAACTGTCGCTGAAGTCTCGCCAGCGTCTTCTGGTTTTTCTGAAAACTGTTTACAGGCCCAAAGACTGTGCCATCTGACAGCGTGGCGAGTTTAGCCACGCCAGCATCCAGCCCGACCATTGATGCTGAAGGGTGAACCGGAGTGGATACTTCACTTTCTGTCTGAATACTGATGTACCACTTACCGCAGGACTGGCTGACAGTGACATTTTTCACAACACCCGTGACCTGACGGCTATTCCGGTAGCGCATCCAGCCCAGTTTCGGCAGA  >  CP006584/1964298‑1964853
                                                                                                                                                                                                                                                                     |                                                                                                                                                                                                                                                                                                      
gTACTCAAGCTGGCGGCGCATTTCATACCAGCCCTGATCCAGTATCGAACGGTTTAAACCTGATTTTGCCCGGACATTGCGCCCCGGCTGGCTGACCGTACCCGCTGCTGACTTTGACATGTTGCTGACCTTCAAATCCTCAATGACTATCATTGCGTGGTTTTTGCTGACGGTCGTTGTGACTTTATGAAGGTAGTCCCTGCGGATATTTGCGATACGGGAATGCAGTCGCTGTATTTTGCGTTTCTGCTTCTGCCAGTTATTGCTGAATTTGACCTTGCGGCTTAACTGTCGCTGa                                                                                                                                                                                                                                                                    <  1:470984/298‑1 (MQ=14)
           tccgatctaTTTCATACCAGCCCTGTTCCAGTAGCGGACGGTTTAAACCTGATTTTGCCCGGACATTGCGCCCCGGCTGGCTGACCGTACCCGCTGCTGACTTTGACATGTTGCTGACCTTCAAATCCTCAATGACTATCATTGCGGGGTTTTTGCTGACGGTCGTTGTGACTTTATGAAGGTAGTCCCTGCGGATATTTGCGATACGGGAATGCAGTCGCTGTATTTTGCGTTTCTGCTTCTGCCAGTTATTGCTGAATTTGACCTTGCGGCTTAACTGTCGCTGAAGTCTCGCCAg                                                                                                                                                                                                                                                         <  2:236668/290‑1 (MQ=9)
               gcgcATTGCGTACCAGCCCTGATCCAGTATCGGACGGTTTAAACCTGAGTTTGCCCGGACATTGCGCCCCGGCTGGCTGGCCGTACCCGCTGCTGACTTTGACATGTTGCTGACCTTCAAATCCTCAATGACTATCATTGCGTGGTTTTTGCTGACGGTCGTTGTGACGTTATGAAGGTAGTCCCTGCGGATATTTGCGATACGGGAATGCAGTCGCTGTATTTTGCGTTTCTGCTTCTGCCAGTTATTGCTGAATTTGACCTTGCGGCTTAACTGTCGCTGAAGTCTCGCCAGCGtc                                                                                                                                                                                                                                                     <  2:574407/298‑1 (MQ=9)
                    tttCATACCAGCCCTGATCCAGTATCGAACGGTTTAAACCTGATTTTGCCCGGACATTGCGCCCCGGCTGGCTGACCGTACCCGCTGCTGACTTTGACATGTTGCTGACCTTCAAATCCTCAATGACTATCATTGCGTGGTTTTTGCTGACGGTCGTTGTGACTTTATGAAGGTAGTCCCTGCGGATATTTGCGATACGGGAATGCAGTCGCTGTATTTTGCGTTTCTGCTTCTGCCAGTTATTGCTGAATTTGACCTTGCGGCTTAACTGTCGCTGAAGTCTCGCCAGcagatcgga                                                                                                                                                                                                                                                >  1:236668/1‑290 (MQ=14)
                            cAGCCCTGATCCAGTATCGAACGGTTTAAACCTGATTTTGCCCGGACATTGCGCCCCGGCTGGCTGACCGTACCCGCTGCTGACTTTGACATGTTGCTGACCTTCAAATCCTCAATGACTATCATTGCGTGGTTTTTGCTGACGGTCGTTGTGACTTTATGAAGGTAGTCCCTGCGGATATTTGCGATACGGGAATGCAGTCGCTGTATTTTGCGTTTCTGCTTCTGCCAGTTATTGCTGAATTTGACCTTGCGGCTTAACTGTCGCTGAAGTCTCGCCAGCGTCTTCTGGTTTTTCt                                                                                                                                                                                                                                        <  1:70141/298‑1 (MQ=14)
                                ccTGATCCAGTATCGAACGGTTTAAACCTGATTTTGCCCGGACATTGCGCCCCGGCTGGCTGACCGTACCCGCTGCTGACTTTGACATGTTGCTGACCTTCAAATCCTCAATGACTATCATTGCGTGGTTTTTGCTGACGGTCGTTGTGACTTTATGAAGGTAGTCCCTGCGGATATTTGCGATACGGGAATGCAGTCGCTGTATTTTGCGTTTCTGCTTCTGCCAGTTATTGCTGAATTTGACCTTGCGGCTTAACTGTCGCTGAAGTCTCGCCAGCGTCTTCTGGTTTTTCTGaaa                                                                                                                                                                                                                                    >  1:410190/1‑298 (MQ=14)
                                  agctCCAATAGCGAACGTTGTAAACCTGAGTTTGCACTGACATTGCGCCCCGGCTGGCTCACCGTACCCGCTGCTGCAGTTGACATGTTGCTGACCTGCAAATCATCAATGACTATCATTGCGTGGTTTTTGCTGACGGTCGTTGTGACTTTATGAAGGTAGTCCCTGCGGATATTTTCGATACGGGAATGCAGTCGCTGTATTTTGCGTTTCTGCTTCTGCCAGTTATTGCTGAATTTGACCTTGCTGCTTAACTGTCGCTGAAGTCTCGCCAGCGTCTTCTGGTTTTTCTGAAAAc                                                                                                                                                                                                                                  <  2:476883/295‑1 (MQ=11)
                                        aGTATCGAACGGTTTAAACCTGATTTTGCCCGGACATTGCGCCCCGGCTGGCTGACCGTACCCGCTGCTGACTTTGACATGTTGCTGACCTTCAAATCCTCAATGACTATCATTGCGTGGTTTTTGCTGACGGTCGTTGTGACTTTATGAAGGTAGTCCCTGCGGATATTTGCGATACGGGAATGCAGTCGCTGTATTTTGCGTTTCTGCTTCTGCCAGTTATTGCTGAATTTGACCTTGCGGCTTAACTGTCGCTGAAGTCTCGCCAGCGTCTTCTGGTTTTTCTGAAAACTGTTTa                                                                                                                                                                                                                            >  1:605073/1‑298 (MQ=14)
                                         gTATCGAACGGTTTAAACCTGATTTTGCCCGGACATTGCGCCCCGGCTGGCTGACCGTACCCGCTGCTGACTTTGACATGTCGCTGACCTTCAAATCCTCAATGACTATCATTGCGTGGTTTTTGCTGACGGTCGTTGTGACTTTATTAAGGTAGTCCCTGCGGATATTTGCGCTACGGGAATGCAGTCGCCGGATTTTGCTTTTCCGCTTCTGCCAGTTATTTCTGAATTTGACCTTGCGTGTTACATTTAGCTGCATTCTCCCAGGCGTCTTCTGGTTTTTCTCTAAACTTTTTAc                                                                                                                                                                                                                           >  2:364856/1‑298 (MQ=11)
                                                     ttAAACCTGATTTTGCCCGGACATTGCGCCCCGGCTGGCTGACCGTACCCGCTGCTGACTTTGACATGTTGCTGACCTTCAAATCCTCAATGACTATCATTGCGGGGTTTTTGCTGACGGTCGTTGTGACTTTGTGAAGGTAGTCCCTGCGGATATTCGCGATACGGGAATGCAGTCGCTGTCTTTTGCGTTTCTGCTTCTGCCAGTTGTTGCTGACTTTGTCCTTGCGGCTTAACTGTCGTGGAAGTCTGGCCAACGTCTTCTGGTTTTTTTGCAAACTGTTTACAGGccccacgag                                                                                                                                                                                                               >  2:580174/1‑292 (MQ=11)
                                                         aCCTGATTTTGCCCGGACATTGCGCCCCGGCTGGCTGACCGTACCCGCTGCTGACTTTGACATGTTGCTGACCTTCAAATCCTCAATGACTATCATTGCGTGGTTTTTGCTGACGGTCGTTGTGACTTTATGAAGGTAGTCCCTGCGGATATTTGCGATACGGGAATGCAGTCGCTGTATTTTGCGTTTCTGCTTCTGCCAGTTATTGCTGAATTTGACCTTGCGGCTTAACTGTCGCTGAAGTCTCGCCAGCGTCTTCTGGTTTTTCTGAAAACTGTTTACAGGCCCAAAGACtgtg                                                                                                                                                                                                           <  1:396658/298‑1 (MQ=17)
                                                                   gCCCGGACATTGCGCCCCGGCTGGCTGACCGTACCCGCTGCTGACTTTGACATGTTGCTGACCTTCAAATCCTCAATGACTATCATTGCGTGGTTTTTGCTGACGGTCGTTGTGACTTTATGAAGGTAGTCCCTGCGGATATTTGCGATACGGGAATGCAGTCGCTGTATTTTGCGTTTCTGCTTCTGCCAGTTATTGCTGAATTTGACCTTGCGGCTTAACTGTCGCTGAAGTCTCGCCAGCGTCTTCTGGTTTTTCTGAAAACTGTTTACAGGCCCAAAGACTGTGCCATCTGACa                                                                                                                                                                                                 <  1:326433/298‑1 (MQ=17)
                                                                      cGGACATTGCGCCCCGGCTGGCTGACCGTACCCGCTGCTGACTTTGACATGTTGCTGACCTTCAAATCCTCAATGACTATCATTGCGTGGTTTTTGCTGACGGTCGTTGTGACTTTATGAAGGTAGTCCCTGCGGATATTTGCGATACGGGAATGCAGTCGCTGTATTTTGCGTTTCTGCTTCTGCCAGTTATTGCTGAATTTGACCTTGCGGCTTAACTGTCGCTGAAGTCTCGCCAGCGGCCTCTGGTTTTTCTGAAAACCTTTTACAGGCACAAAGACTGTGCCATTTGCAAGCg                                                                                                                                                                                              >  2:632377/1‑298 (MQ=255)
                                                                          cATTGCGCCCCGGCTGGCTGACCGTACCCGCTGCTGACTTTGACATGTTGCTGACCTTCAAATCCTCAATGACTATCATTGCGTGGTTTTTGCTGACGGTCGTTGTGACTTTGTGAAGGTAGTCCCTGCGGATATTTGCGATACGGGAATGCAGTCGCTGTATTTTGCGTTTCTGCTTCTGCCAGTTGTTGCTGAATTTGACCTTGCGGCTTAACTGTCGCTGAAGTCTCGCCAGCGTCTTCTGGTTTTTCTGAAAACTGTTTACAGGCCCAAAGACTGTGCCCTCTgtaacgcgtgt                                                                                                                                                                                          >  2:229716/1‑288 (MQ=9)
                                                                               cgcCCCGGCTGGCTGACCGTACCCGCTGCTGACTTTGACATGTTGCTGACCTTCAAATCCTCAATGACTATCATTGCGTGGTTTTTGCTGACGGTCGTTGTGACTTTATGAAGGTAGTCCCTGCGGATATTTGCGATACGGGAATGCAGTCGCTGTATTTTGCGTTTCTGCTTCTGCCAGTTATTGCTGAATTTGACCTTGCGGCTTAACTGTCGCTGAAGTCTCGCCAGCGTCTTCTGGTTTTTCTGAAAACTGTTTACAGGCCCAAAGACTGTGCCATCTGACAGCGTGGCGAGtt                                                                                                                                                                                     <  1:573883/298‑1 (MQ=17)
                                                                                 ccccGGCTGGCTGACCGTACCCGCTGCTGACTTTGACATGTTGCTGACCTTCAAATCCTCAATGACTATCATTGCGTGGTTTTTGCTGACGGTCGTTGTGACTTTATGAAGGTAGTCCCTGCGGATATTTGCGATACGGGAATGCAGTCGCTGTATTTTGCGTTTCTCCTTCTGCCAGTTCTTGCTGACTTTGACCTTGCGCCTTAACTGTTGCTGAAGTCTCGCCAGCGTCTTCTGGTCTTTCTGACAACTGCTTACAGGCCCAAACACTTTGCCATCCGACAGCGAGGACAGTTTa                                                                                                                                                                                   >  2:8345/1‑298 (MQ=11)
                                                                                                                tttGACATGTTGCTGACCTTCAAAT‑CTCAATGACTATCATTGCGTGGTTTTTGCTGACGGTCGTTGTGACTTTATGAAGGTAGTCCCTGCGGATATTTGCGATACGGGAATGCAGTCGCTGTATTTTGCGTTTCTGCTTCTGCCAGTTATTGCTGAATTTGACCTTGCGGCTTAACTGTCGCTGAAGTCTCGCCAGCGTCTTCTGGTTTTTCTGAAAACTGTTTACAGGCCCAAAGACTGTGCCATCTGACAGCGTGGCGAGTTTAGCCACGCCAGCATCCAGCCCGATCATTGATGc                                                                                                                                                   >  1:53287/1‑298 (MQ=17)
                                                                                                                       tGTTGCTGACCTTCAAATCCTCAATGACTATCATTGCGTGGTTTTTGCTGACGGTCGTTGTGACTTTATGAAGGTAGTCCCTGCGGATATTTGCGATACGGGAATGCAGTCGCTGTATTTTGCGTTTCTGCTTCTGCCAGTTATTGCTGAATTTGACCTTGCGGCTTAACTGTCGCTGAAGTCTCGCCAGCGTCTTCTGGTTTTTCTGAAAACTGTTTACAGGCCCAAAGACTGTGCCATCTGACAGCGTGGCGAGTTTAGCCACGCCAGCATCCAGCCCGATCATTGATGCTGAAgg                                                                                                                                             >  1:320636/1‑298 (MQ=17)
                                                                                                                        ggtgctgcGCTTCAATGCCTCTGTGACGTTCATGTCGTGGTTTTTGGGTGGGGGCGTTGGTACTTTAGGAAGCTGTGCCCTGCGGATAGTTGCGATATGGGAATGCAGGCGCGTTATGTTGCGTTTCTTCTTCTGCCAGTTAGTGCTGAATTTGACCTAGCGGCTTAACTGTCGCTGAAGTATCGCCAGCGTCTTCTGGTTTTTCTGAAAACTGTTTACAGGCCCAAAGACTGTGCCATCTGACAGCGTGGCGAGTTTAGCCACGCCAGCATCCAGCCCGATCATTGATGCTGAAggg                                                                                                                                            <  2:234438/296‑1 (MQ=11)
                                                                                                                         ttGCTGACCTTCAAATCCTCAATGACTATCATTGCGTGGTTTTTGCTGACGGTCGTTGTGACTTTATGAAGGTAGTCCCTGCGGATATTTGCGATACGGGAATGCAGTCGCTGTATTTTGCGTTTCTGCTTCTGCCAGTTATTGCTGAATTTGACCTTGCGGCTTAACTGTCGCTGAAGTCTCGCCAGCGTCTTCTGGTTTTTCTGAAAACTGTTTACAGGCCCAAAGACTGTGCCATCTGACAGCGTGGCGAGTTTAGCCACGCCAGCATCCAGCCCGATCATTGATGCTGAAGGGt                                                                                                                                           <  1:193094/298‑1 (MQ=17)
                                                                                                                               aCCTTCAAATCCTCAATGACTATCATTGCGTGGTTTTTGCTGACGGTCGTTGTGACTTTATGAAGGTAGTCCCTGCGGATATTTGCGATACGGGAATGCAGTCGCTGTATTTTGCGTTTCTGCTTCTGCCAGTTATTGCTGAATTTGACCTTGCGGCTTAACTGTCGCTGAAGTCTCGCCAGCGTCTTCTGGTTTTTCTGAAAACTGTTTACAGGCCCAAACACTGTTCCATCTGACCGCCTGGCGAGTTTAGCCCCGCCAGCAACCACCCCGACCTTTTATGCACCAGGGCGAACCg                                                                                                                                     >  2:524990/1‑298 (MQ=11)
                                                                                                                                           tCAATGACTATCATTGCGTGGTTTTTGCTGACGGTCGTTGTGACTTTATGAAGGTAGTCCCTGCGGATATTTGCGATACGGGAATGCAGTCGCTGTATTTTGCGTTTCTGCTTCTGCCAGTTATTGCTGAATTTGACCTTGCGGCTTAACTGTCGCTGAAGTCTCGCCAGCGTCTTCTGGTTTTTCTGAAAACTGTTTACAGGCCCAAAGACTGTGCCATCTGACAGCGTGGCGAGTTTAGCCACGCCAGCATCCAGCCCGATCATTGATGCTGAAGGGTGAACCGGAGTTGATACtt                                                                                                                         >  1:82024/1‑298 (MQ=17)
                                                                                                                                              aTGACTATCATTGCGTGGTTTTTGCTGACGGTCGTTGTGACTTTATGAAGGTAGTCCCTGCGGATATTTGCGATACGGGAATGCAGTCGCTGTATTTTGCGTTTCTGCTTCTGCCAGTTATTGCTGAATTTGACCTTGCGGCTTAACTGTCGCTGAAGTCTCGCCAGCGTCTTCTGGTTTTTCTGAAAACTGTTTACAGGCCCAAAGACTGTGCCATCTGACAGCGTGGCGAGTTTAGCCACGCCAGCATCCAGCCCGATCATTGATGCTGAAGGGTGAACCGGAGTTGCTACTTCAc                                                                                                                      >  1:153183/1‑298 (MQ=17)
                                                                                                                                                    ttcttttggttgTTTGTGCTGCCGGGCGTGTTTACTTTATGAATGGAGTCCCTGGGTGTATTTGCGTTACGGGAATGCGGTGTCTTTATTTTTCGTTTCTGCTTCTTCCAGTTATTGCTGAATTTGACCTTGCGGCTTAACTGTCGCTGAAGTCTCGCCAGCGTCTTCTGGTTTTTCTGAAAACTGTTTACAGGCCCAAAGACTGTGCCATCTGACAGCGTGGCGAGTTTAGCCACGCCAGCATCCAGCCCGATCATTGATGCTGAAGGGTGAACCGGAGTTGATACTTCACTTtctg                                                                                                                <  2:380992/287‑1 (MQ=11)
                                                                                                                                                    aTCATTGCGTGGTTTTTGCTGACGGTCGTTGTGACTTTATGAAGGTAGTCCCTGCGGATATTTGCGATACGGGAATGCAGTCGCTGTATTTTGCGTTTCTGCTTCTGCCAGTTATTGCTGAATTTGACCTTGCGGCTTAACTGTCGCTGAAGTCTCGCCAGCGTCTTCTGGTTTTTCTGAAAACTGTTTACAGGCCCAAAGACTGTGCCATCTGACAGCGTGGCGAGTTTAGCCACGCCAGCATCCAGCCCGATCATTGATGCTGAAGGGTGAACCGGAGTTGATACTTCACTTtctg                                                                                                                >  1:406986/1‑298 (MQ=17)
                                                                                                                                                      cATTGCGTGGTTTTTGCTGACGGTCGTTGTGACTTTATGAAGGTAGTCCCTGCGGATATTTGCGATACGGGAATGCAGTCGCTGTATTTTGCGTTTCTGCTTCTGCCAGTTATTGCTGAATTTGACCTTGCGGCTTAACTGTCGCTGAAGTCTCGCCAGCGTCTTCTGGTTTTTCTGAAAACTGTTTACAGGCCCAAAGACTGTGCCATCTGACAGCGTGGCGAGTTTAGCCACGCCAGCATCCAGCCCGATCATTGATGCTGAAGGGTGAACCGGAGTTGATACTTCACTTtctgtc                                                                                                              <  1:632377/298‑1 (MQ=17)
                                                                                                                                                                         aCGGTCGTTGTGACTTTATGAAGGTAGTCCCTGCGGATATTTGCGATACGGGAATGCAGTCGCTGTATTTTGCGTTTCTGCTTCTGCCAGTTATTGCTGAATTTGACCTTGCGGCTTAACTGTCGCTGAAGTCTCGCCAGCGTCTTCTGGTTTTTCTGAAAACTGTTTACAGGCCCAAAGACTGTGCCATCTGACAGCGTGGCGAGTTTAGCCACGCCAGCATCCAGCCCGATCATTGATGCTGAAGGGTGAACCGGAGTTGATACTTCACTTTCTGTCTGAATACTGATGTACCACt                                                                                           >  1:218217/1‑298 (MQ=17)
                                                                                                                                                                         aCGGTCGTTGTGACTTTATGAAGGTAGTCCCTGCGGATATTTGCGATACGGGAATGCAGTCGCTGTATTTTGCGTTTCTGCTTCTGCCAGTTATTGCTGAATTTGACCTTGCGGCTTAACTGTCGCTGAAGTCTCGCCAGCGTCTTCTGGTTTTTCTGAAAACTGTTTACAGGCCCAAAGACTGTGCCATCTGACAGCGTGGCGAGTTTAGCCACGCCAGCATCCAGCCCGATCATTGATGCTGAAGGGTGAACCGGAGTTGATACTTCACTTTCTGTCTGAATACTGATGTACCACt                                                                                           >  1:493556/1‑298 (MQ=17)
                                                                                                                                                                                ttGTGACTTTATGAAGGTAGTCCCTGCGGATATTTGCGATACGGGAATGCAGTCGCTGTATTTTGCGTTTCTGCTTCTGCCAGTTATTGCTGAATTTGACCTTGCGGCTTAACTGTCGCTGAAGTCTCGCCAGCGTCTTCTGGTTTTTCTGAAAACTGTTTACACGCCCAAAGACTGTGCCATCTGAAAGCGTGGCGAGGTTAGCAACGCCAGAATCCAGCCCGACCACTGATGCTGAAGGGTGACCCGGAGTGGATAATTCACTCTTTGTCTCCATCTGGATCTACAACTTTCCGcc                                                                                    >  2:238668/1‑297 (MQ=11)
                                                                                                                                                                                      gttgtttgtggTAGTCCCTCCGGAGTTTTGGGATCCGGGTGTGCAGGCGCTGTATTTTGGGTTTTTGCTTCTGCCAGTTAGTGCTGAATGTGTCCTTGCGGCTTAACTGTCGCTGAAGTCTCGCCATCGTCTTCTGGTTTTTCTGAAAACTGTTTACAGGCCCAAAGACTGTGCCATCTGACAGCGTGGCGAGTTTAGCCACGCCAGCATCCAGCCCGATCATTGATGCTGAAGGGTGAACCGGAGTTGATACTTCACTTTCTGTCTGAATACTGCTGTACCACTTACCGCAGGACTg                                                                              <  2:381656/289‑1 (MQ=11)
                                                                                                                                                                                         tgtGAAGGTAGGTCCTTTGGGTATTTGCGGTACGGGAATGGAGGTGCTGTGTTGTGCGTTTCTGCTTCTGCCACTTATTGCGTAATTTGACCTTGCGGCTTAACTGTCGCTGAAGTCTCGCCAGCGTCTTCTGGTTTTTCTGAAAACTGTTTACAGGCCCAAAGACTGTGCCATCTGACAGCGTGGCGAGTTTAGCCACGCCAGCATCCAGCCCGATCATTGATGCTGAAGGGTGAACCGGAGTTGATACTTCACTTTCTGTCTGAATACTGATGTACCACTTACCGCAGGACTGGCt                                                                           <  2:469195/296‑1 (MQ=11)
                                                                                                                                                                                                 tAGTCCCTTGTGATATTTGCGAGACGGGAATGCAGTCGTCTTTGTTTGAGTTTCTGCGTCTGCCAGTTATTGCTGACTTTGACCTTGCGGCTTAACTGTCGCTGCAGTCTCGCCAGCGTCTTCTGGTTTTTCTGAGAACTGTTTACAGGCCCCAAGACTGTGCCATCTGACAGCGTTGCGAGTTTAGCCACGCCAGCATCCAGCCCGATCCTTGATGCTGAAGGGTGAACCGGAGTTGATACTTCACTTTCTGTCTGAATACTGATGTACCACTTACCGCAGGACTGGCTGACCGTGa                                                                   <  2:529028/298‑1 (MQ=11)
                                                                                                                                                                                                             atatTTGCGATACGGGAATGCAGTCGCTGTATTTTGCGTTTCTGCTTCTGCCAGTTATTGCTGAATTTGACCTTGCGGCTTAACTGTCGCTGAAGTCTCGCCAGCGTCTTCTGGTTTTTCTGAAAACTGTTTACAGGCCCAAAGACTGTGCCATCTGACAGCGTGGCGAGTTTAGCCACGCCAGCATCCAGCCCGATCATTGATGCTGAAGGGTGAACCGGAGTTGATACTTCACTTTCTGTCTGAATACTGATGTACCACTTACCGCAGGACTGGCTGACCGTGACATTTTTCacaa                                                       >  1:475751/1‑298 (MQ=17)
                                                                                                                                                                                                                 ttGCGATACGGGAATGCAGTCGCTGTATTTTGCGTTTCTGCTTCTGCCAGTTATTGCTGAATTTGACCTTGCGGCTTAACTGTCGCTGAAGTCTCGCCAGCGTCTTCTGGTTTTTCTGAAAACTGTTTACAGGCCCAAAGACTGTGCCATCTGACAGCGTGGCGAGTTTAGCCACGCCAGCATCCAGCCCGATCATTGATGCTGAAGGGTGAACCGGAGTTGATACTTCACTTTCTGTCTGAATACTGATGTACCACTTACCGCAGGACTGGCTGACCGTGACATTTTTCACAACAcc                                                   >  1:633572/1‑298 (MQ=17)
                                                                                                                                                                                                                      aTACGGGAATGCAGTCGCTGTATTTTGCGTTTCTGCTTCTGCCAGTTATTGCTGAATTTGACCTTGCGGCTTAACTGTCGCTGAAGTCTCGCCAGCGTCTTCTGGTTTTTCTGAAAACTGTTTACAGGCCCAAAGACTGTGCCATCTGACAGCGTGGCGAGTTTAGCCACGCCAGCATCCAGCCCGATCATTGATGCTGAAGGGTGAACCGGAGTTGATACTTCACTTTCTGTCTGAATACTGATGTACCACTTACCGCAGGACTGGCTGACCGTGACATTTTTCACAACACCCGTGa                                              >  1:477332/1‑298 (MQ=17)
                                                                                                                                                                                                                                gcgtacgttGTATTTCGGATTTCTGTTTCTGCCAGTTATTGCTGCATTTGGCCTTACGGCTTAGCGTTCGCTGCAGTCTCGCCTGCTGCTTCTTGTTTTTCTGCAAACTGTTTACAGGCCCAAAGACTGTGCCATCTGACAGCGTGGCGAGTTTAGCCACGCCAGCATCCAGCCCGATCATTGATGCTGAAGGGTGAACCGGAGTTGATACTTCACTTTCTGTCTGAATACTGATGTACCACTTACCGCAGGACTGGCTGACCGTGACATTTTTCACAACACCCGTGACCTGACGGCt                                    <  2:258011/290‑1 (MQ=11)
                                                                                                                                                                                                                                gCAGTCGCTGTATTTTGCGTTTCTGCTTCTGCCAGTTATTGCTGAATTTGACCTTGCGGCTTAACTGTCGCTGAAGTCTCGCCAGCGTCTTCTGGTTTTTCTGAAAACTGTTTACAGGCCCAAAGACTGTGCCATCTGACAGCGTGGCGAGTTTAGCCACGCAAGCATCCAGCCCGATCATTGATGCTGAAGGGTGAACCGGAGTTGATACTTCACTTTCTGTCTGAATACTGATGTACCACTTACCGCAGGACTGGCTGACCGTGACATTTTTCACAACACCCGTGACCTGACGGCt                                    <  1:141362/298‑1 (MQ=9)
                                                                                                                                                                                                                                   gTCTCTGTATTTTGCGTTTCTGATTCGGCCAGTTATTGCGGGATTTGACCTTGGGGCTTAATGGTCGGTGAAGTCTCGCCAGCGTCTTCTGGTTTTTCTGAAAACTGTTTACAGGCCCAAAGACTGTGCCATCTGACAGCGTGGCGAGGTTAGCCACGCCAGCATCCAGCCCGATCATTGATGCTAAAGGGTGAACCGGAGTTGATACTTCACTTTCTGTCTGAATACTGATGTACCACTTACCGCAGGACTGGCTGACCGTGACATTTTTCACAACACCCGTGACCTGACGGCTAtt                                 <  2:43871/298‑1 (MQ=11)
                                                                                                                                                                                                                                   gTCGCTGTATTTTGCGTTTCTGCTTCTGCCAGTTATTGCTGAATTTGACCTTGCGGCTTAACTGTCGCTGAAGTCTCGCCAGCGTCTTCTGGTTTTTCTGAAAACTGTTTACAGGCCCAAAGACTGTGCCATCTGACCGCGTGGCGAGTTTAGCCACCCCAGCATCCAGCCCGACAATTCATCCTGAAGGGTGAACCGGAGTTGATACCTCACTTTCTGGCTGAATACCTAAGTACCACTTACCGCAGGACTGGCCTACCGTGCCATTTTTTGACACACCCGTCCCCTCACggttctt                                 >  2:328908/1‑293 (MQ=11)
                                                                                                                                                                                                                                          tATTTTGCGTTTCTGCTTCTGGCAGTTTTGGCTGATTTTGACCTGGCGGCTTAACTGTCGCTGAAGTCTCGCCATCGTCTTCGGGTTTTTCTGAAAGCTTTTTACAGGCCCCAAGACTGTGCCATCTGACAGCGTGGCGAGTTGGGCCACGCCAGCATACAGCCCGATCATTGATGCTGAAGGGTGAACCGGAGTTGATACTTCACTTTCTGTCTGAGTACTGATGTACCACTTACCGCAGGACTGGCTGACCGTGACATTTTTCACAACACCCGTGACCTGACGGCTATTCCGGTAg                          <  2:243625/298‑1 (MQ=11)
                                                                                                                                                                                                                                                               ccAGTTATTGCTGAATTTGACCTTGCGGCTTAACTGTCGCTGAAGTCTCGCCAGCGTCTTCTGGTTTTTCTGAAAACTGTTTACAGGCCCAAAGACTGTGCCATCTGACAGCGTGGCGAGTTTAGCCACGCCAGCATCCAGCCCGATCATTGATGCTGAAGGGTGAACCGGAGTTGATACTTCACTTTCTGTCTGAATACTGATGTACCACTTACCGCAGGACTGGCTGACCGTGACATTTTTCACAACACCCGTGACCTGACGGCTATTCCGGTAGCGCATCCAGCCCAGTTTCGGc     <  1:238668/298‑1 (MQ=14)
                                                                                                                                                                                                                                                                cAGTTATTGCTGAATTTGACCTTGCGGCTTAACTGTCGCTGAAGTCTCGCCAGCGTCTTCTGGTTTTTCTGAAAACTGTTTACAGGCCCAAAGACTGTGCCATCTGACAGCGTGGCGAGTTTAGCCACGCCAGCATCGAGCCCGATCATTGATGCTGAAGGGGGAACCGGAGTTGATACTTCACTTTCTGTCTGAATACTGATGTAACACTTACCGCAGGACTGGCTGACCGTGACATTTTTCACAACACCCGTAACCTGACGGCTATTCCGGTATCGCCTCCAGCCCAGTTTCGGCa    >  1:68034/1‑298 (MQ=11)
                                                                                                                                                                                                                                                                  gttattGCTGAATTTGACCTTGCGGCTTAACTGTCGCTGAAGTCTCGCCAGCGTCTTCTGGTTTTTCTGAAAACTGTTTACAGGCCCAAAGACTGTGCCATCTGACAGCGTGGCGAGTTTAGCCACGCCAGCATCCAGCCCGATCATTGATGCTGAAGGGTGAACCGGAGTTGATACTTCACTTTCTGTCTGAATACTGATGTACCACTTACCGCAGGACTGGCTGACCGTGACATTTTTCACAACACCCGTGACCTGACGGCTATTCCGGTAGCGCATCCAGCCCAGTTTCGGCAGa  >  1:487036/1‑298 (MQ=14)
                                                                                                                                                                                                                                                                     |                                                                                                                                                                                                                                                                                                      
GTACTCAAGCTGGCGGCGCATTTCATACCAGCCCTGATCCAGTATCGAACGGTTTAAACCTGATTTTGCCCGGACATTGCGCCCCGGCTGGCTGACCGTACCCGCTGCTGACTTTGACATGTTGCTGACCTTCAAATCCTCAATGACTATCATTGCGTGGTTTTTGCTGACGGTCGTTGTGACTTTATGAAGGTAGTCCCTGCGGATATTTGCGATACGGGAATGCAGTCGCTGTATTTTGCGTTTCTGCTTCTGCCAGTTGTTGCTGAATTTGACCTTGCGGCTTAACTGTCGCTGAAGTCTCGCCAGCGTCTTCTGGTTTTTCTGAAAACTGTTTACAGGCCCAAAGACTGTGCCATCTGACAGCGTGGCGAGTTTAGCCACGCCAGCATCCAGCCCGACCATTGATGCTGAAGGGTGAACCGGAGTGGATACTTCACTTTCTGTCTGAATACTGATGTACCACTTACCGCAGGACTGGCTGACAGTGACATTTTTCACAACACCCGTGACCTGACGGCTATTCCGGTAGCGCATCCAGCCCAGTTTCGGCAGA  >  CP006584/1964298‑1964853

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