Predicted mutation
evidence seq id position mutation annotation gene description
RA CP006584 1,964,699 C→T V180I (GTC→ATC)  LY180_09345 ← transposase IS609

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0065841,964,6990CT93.6% 78.4 / ‑5.5 48V180I (GTC→ATC) LY180_09345transposase IS609
Reads supporting (aligned to +/- strand):  ref base C (3/0);  new base T (23/21);  total (26/22)
Fisher's exact test for biased strand distribution p-value = 2.42e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.88e-01

TTTGACATGTTGCTGACCTTCAAATCCTCAATGACTATCATTGCGTGGTTTTTGCTGACGGTCGTTGTGACTTTATGAAGGTAGTCCCTGCGGATATTTGCGATACGGGAATGCAGTCGCTGTATTTTGCGTTTCTGCTTCTGCCAGTTGTTGCTGAATTTGACCTTGCGGCTTAACTGTCGCTGAAGTCTCGCCAGCGTCTTCTGGTTTTTCTGAAAACTGTTTACAGGCCCAAAGACTGTGCCATCTGACAGCGTGGCGAGTTTAGCCACGCCAGCATCCAGCCCGACCATTGATGCTGAAGGGTGAACCGGAGTGGATACTTCACTTTCTGTCTGAATACTGATGTACCACTTACCGCAGGACTGGCTGACAGTGACATTTTTCACAACACCCGTGACCTGACGGCTATTCCGGTAGCGCATCCAGCCCAGTTTCGGCAGAAAAATACGGCTGTTTTCCTGATCGAGCTTTACGCCCTGCGGGTAGCGGAATGCATCATTCTGTCCCCGCTTTTTGAATCGGGGAAAAGCAGCCCGCTTCCGGAAGAAGTTTTTGTAAGCCCGCTCAAGGTCTTTCAGTG  >  CP006584/1964410‑1964992
                                                                                                                                                                                                                                                                                                 |                                                                                                                                                                                                                                                                                                     
tttGACATGTTGCTGACCTTCAAATCTCAATGACTATCATTGCGTGGTTTTTGCTGACGGTCGTTGTGACTTTATGAAGGTAGTCCCTGCGGATATTTGCGATACGGGAATGCAGTCGCTGTATTTTGCGTTTCTGCTTCTGCCAGTTATTGCTGAATTTGACCTTGCGGCTTAACTGTCGCTGAAGTCTCGCCAGCGTCTTCTGGTTTTTCTGAAAACTGTTTACAGGCCCAAAGACTGTGCCATCTGACAGCGTGGCGAGTTTAGCCACGCCAGCATCCAGCCCGATCATTGATGc                                                                                                                                                                                                                                                                                              >  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:493556/1‑298 (MQ=17)
                                                         aCGGTCGTTGTGACTTTATGAAGGTAGTCCCTGCGGATATTTGCGATACGGGAATGCAGTCGCTGTATTTTGCGTTTCTGCTTCTGCCAGTTATTGCTGAATTTGACCTTGCGGCTTAACTGTCGCTGAAGTCTCGCCAGCGTCTTCTGGTTTTTCTGAAAACTGTTTACAGGCCCAAAGACTGTGCCATCTGACAGCGTGGCGAGTTTAGCCACGCCAGCATCCAGCCCGATCATTGATGCTGAAGGGTGAACCGGAGTTGATACTTCACTTTCTGTCTGAATACTGATGTACCACt                                                                                                                                                                                                                                      >  1:218217/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)
                                                                                                                                                              tttGACCTTGCGGCTTAACTGTCGCTGAAGTCTCGCCAGCGTCTTCTGGTTTTTCTGAAACCTGTTTACAGGCCCAAAGACTGTGCCATCTGACAGCGTGGCGAGTTTAGCCACGCCAGAATCCAGCCCGATCATTGATGCTGAAGGGGGAACCGGAGTTGATAATTCACTTTCTGTCTGAATACTGATTTACCACTTACCGCAGGACTGGCTGACCGGCACATTTTTGACAACACCCGTGACCTGACGGCTCTTCCGGTAGCCGAGCCAGCCCATTTTGTGCAGCACAATACGGCtt                                                                                                                                 >  2:310163/1‑297 (MQ=11)
                                                                                                                                                                        cGGCTTAACTGTCGCTGAAGTCTCGCCAGCGTCTTCTGGTTTTTCTGAAAACTGTTTACAGGCCCAAAGACTGTGCCATCTGACAGCGTGGCGAGTTTAGCCACGCCAGCATCCAGCCCGATCATTGATGCTGAAGGGTGAACCGGAGTTGATACTTCACTTTCTGTCTGAATACTGATGTACCACTTACCGCAGGACTGGCTGACCGTGACATTTTTCACAACACCCGTGACCTGACGGCTATTCCGGTAGCGCATCCAGCCCAGTTTCGGCAGAAAAATACGGCTGTTTTCCTGAt                                                                                                                       <  1:555863/298‑1 (MQ=14)
                                                                                                                                                                            ttAACTGTCGCTGAAGTCTCGCCAGCGTCTTCTGGTTTTTCTGAAAACTGTTTACAGGCCCAAAGACTGTGCCATCTGACAGCGTGGCGAGTTTAGCCACGCCAGCATCCAGCCCGATCATTGATGCTGAAGGGTGAACCGGAGTTGATACTTCACTTTCTGTCTGAATACTGATGTACCACTTACCGCAGGACTGGCTGACCGTGACATTTTTCACAACACCCGTGACCTGACGGCTATTCCGGTAGCGCATCCAGCCCAGTTTCGGCAGAAACATACGGCTGTTTTctggtccgga                                                                                                                   >  2:214979/1‑289 (MQ=9)
                                                                                                                                                                            ttAACTGTCGCTAAAGTCTCGCCAGCGTCTTCTGGTTTTTCTGAAAACTGTTTACAGGCCCAAAGACTGTGCCATCTGACAGCGTGGCGAGTTTAGCCACGCCAGCATCCAGCCCGATCATTGATGCTGAAGGGTGAACCGGAGTTGATACTTCACTTTCTGTCTGAATACTGATGTACCACTTACCGCAGGACTGGCTGACCGTGACATTTTTCACAACACCCGTGACCTGACGGCTATTCCGGTAGCGCATCCAGCCCAGTTTCGGCAGAAAAATACGGCTGTTTTCCTGATCGAg                                                                                                                   >  1:173123/1‑298 (MQ=14)
                                                                                                                                                                             tAACTGTCGCTGAAGTCTCGCCAGCGTCTTCTGGTTTTTCTGAAAACTGTTTACAGGCCCAAAGACTGTGCCATCTGACAGCGTGGCGAGTTTAGCCACGCCAGCATCCAGCCCGATCATTGATGCTGAAGGGTGAACCGGAGTTGATACTTCACTTTCTGTCTGAATACTGATGTACCACTTACCGCAGGACTGGCTGACCGTGACATTTTTCACAACACCCGTGACCTGACGGCTATTCCGGTAGCGCATCCAGCCCAGTTTCGGCAGAAAAATACGGCTGTTTTCCTGATCGAGc                                                                                                                  >  1:93185/1‑298 (MQ=14)
                                                                                                                                                                                      cTGAAGTCTCGCCAGCGTCTTCTGGTTTTTCTGAAAACTGTTTACAGGCCCAAAGACTGTGCCATCTGACAGCGTGGCGAGTTTAGCCACGCCAGCATCCAGCCCGATCATTGATGCTGAAGGGTGAACCGGAGTTGATACTTCACTTTCTGTCTGAATACTGATGTACCACTTACCGCAGGACTGGCTGACCGTGACATTTTTCACAACACCCGTGACCTGACGGCTATTCCGGTAGCGCATCCAGCCCAGTTTCGGCAGAAAAATACGGCTGTTTTCCTGATCGAGCTTTACGccc                                                                                                         <  1:214979/298‑1 (MQ=17)
                                                                                                                                                                                        gAAGTCTCGCCAGCGTCTTCTGGTTTTTCTGAAAACTGTTTACAGGCCCAAAGACTGTGCCATCTGACAGCGTGGCGAGTTTAGCCACGCCAGCATCCAGCCCGATCATTGATGCTGAAGGGTGAACCGGAGTTGATACTTCACTTTCTGTCTGAATACTGATGTACCACTTACCGCAGGACTGGCTGACCGTGACATTTTTCACAACACCCGTGACCTGTCCGCTATTCCGGGAGCGCATCCACCCCAGTTTCGGCAGAAAAAACCGGATGTTTTCCTGACCGAGCTTTACccccct                                                                                                       >  2:88745/1‑296 (MQ=11)
                                                                                                                                                                                            tctcGCCAGCGTCTTCTGGTTTTTCTGAAAACTGTTTACAGGCCCAAAGACTGTGCCATCTGACAGCGTGGCGAGTTTAGCCACGCCAGCATCCAGCCCGATCATTGATGCTGAAGGGTGAACCGGAGTTGATACTTCACTTTCTGTCTGAATACTGATGTACCACTTACCGCAGGACTGGCTGACCGTGACATTTTTCACAACACCCGTGACCTGACGGCTATTCCGGTAGCGCATCCAGCCCAGTTTCGGCAGAAAAATACGGCTGTTTTCCTGATCGAGCTTTACGCCCTGCggg                                                                                                   >  1:642355/1‑298 (MQ=17)
                                                                                                                                                                                                     cGTCTTCTGGTTTTTCTGAAAACTGTTTACAGGCCCAAAGACTGTGCCATCTGACAGCGTGGCGAGTTTAGCCACGCCAGCATCCAGCCCGATCATTGATGCTGAAGGGTGAACCGGAGTTGATACTTCACTTTCTGTCTGAATACTGATGTACCACTTACCGCAGGACTGGCTGACCGTGACATTTTTCACAACACCCGTGACCTGACGGCTATTCCGGTAGCGCATCCGGCCCAGTTTCGGCAGAAAAATACGGCTGTTTTCCTGATCGAGCTTTACGCCCTGCGGGTAGCGGAAt                                                                                          <  1:484959/298‑1 (MQ=17)
                                                                                                                                                                                                                    cTGAAAACTGTTTACAGGCCCAAAGACTGTGCCATCTGACAGCGTGGCGAGTTTAGCCACGCCAGCATCCAGCCCGATCATTGATGCTGAAGGGTGAACCGGAGTTGATACTTCACTTTCTGTCTGAATACTGATGTACCACTTACCGCAGGACTGGCTGACCGTGACATTTTTCACAACACCCGTGACCTGACGGCTATTCCGGTAGCGCATCCAGCCCAGTTTCGGCAGAAAAATACGGCTGTTTTCCTGATCGAGCTTTACGCCCTGCGGGTAGCGGAATGCATCATTCTGTccc                                                                           <  1:8345/298‑1 (MQ=17)
                                                                                                                                                                                                                          acggTTTACAGGCCCAAAGACTGTGCCATCTGACAGCGTGGCGAGTTTAGCCACGCCAGCATCCAGCCCGATCATTGATGCTGAAGGGTGAACCGGAGTTGATACTTCACTTTCTGTCTGAATACTGATGTACCACTTACCGCAGGACTGGCTGACCGTGACATTTTTCACAACACCCGTGACCTGACGGCTATTCCGGTAGCGCATCCAGCCCAGTTTCGGCAGAAAAATACGGCTGTTTTCCTGAt                                                                                                                       <  2:445174/245‑1 (MQ=14)
                                                                                                                                                                                                                          aCTGTTTACAGGCCCAAAGACTGTGCCATCTGACAGCGTGGCGAGTTTAGCCACGCCAGCATCCAGCCCGATCATTGATGCTGAAGGGTGAACCGGAGTTGATACTTCACTTTCTGTCTGAATACTGATGTACCACTTACCGCAGGACTGGCTGACCGTGACATTTTTCACAACACCCGTGACCTGACGGCTATTCCGGTAGCGCATCCAGCCCAGTTTCGGCAGAAAAATACGGCTGTTTTCCTGAt                                                                                                                       >  1:445174/1‑248 (MQ=14)
                                                                                                                                                                                                                                            gACTGTGCCATCTGACAGCGTGGCGAGTTTAGCCACGCCAGCATCCAGCCCGATCATTGATGCTGAAGGGTGAACCGGAGTTGATACTTCACTTTCTGTCTGAATACTGATGTACCACTTACCGCAGGACTGGCTGACCGTGACATTTTTCACAACACCCGTGACCTGACGGCTATTCCGGTAGCGCATCCAGCCCAGTTTCGGCAGAAAAATACGGCTGTTTTCCTGATCGAGCTTTACGCCCTGCGGGTAGCGGAATGCATCATTCTGTCCCCGCTTTTTGAATCGGGGAAAagca                                                   <  1:326598/298‑1 (MQ=14)
                                                                                                                                                                                                                                             aCTGTGCCATCTGACAGCGTGGCGAGTTTAGCCACGCCAGCATCCAGCCCGATCATTGATGCTGAAGGGTGAACCGGAGTTGATACTTCACTTTCTGTCTGAATACTGATGTACCACTTACCGCAGGACTGGCTGACCGTGACATTTTTCAAAACACCCGTGACCTGACGGCTATTCCGGTAGCGCATCCAGCCCAGTTTCGACAGAAAAATACGGCTGCTTTCCTGATCGAGCTTTACGCCCTGCCGGTCGCGGGTAGCATCCTCCTTTCCCCGCCTTTTGAATCGGGGAAaagaag                                                  >  2:281108/1‑298 (MQ=11)
                                                                                                                                                                                                                                              cTGTGCCATCTGACAGCGTGGCGAGTTTAGCCACGCCAGCATCCAGCCCGATCATTGATGCTGAAGGGTGAACCGGAGTTGATACTTCACTTTCTGTCTGAATACTGATGTACCACTTACCGCAGGACTGGCTGACCGTGACATTTTTCACAACACCCGTGACCTGACGGCTATTCCGGTAGCGCATCCAGCCCAGTTTCGGCAGAAAAATACGGCTGTTTTCCTGATCGAGCTTTACGCCCTGCGGGTAGCGGAATGCATCATTCTGTCCCCGCTTTTTGAATCGGGGAAAagcagc                                                 <  1:524990/298‑1 (MQ=14)
                                                                                                                                                                                                                                                             gCGTGGCGAGTTTAGCCACGCCAGCATCCAGCCCGATCATTGATGCTGAAGGGTGAACCGGAGTTGATACTTCACTTTCTGTCTGAATACTGATGTACCACTTACCGCAGGACTGGCTGACCGTGACATTTTTCACAACACCCGTGACCTGACGGCTATTCCGGTAGCGCATCCAGCCCAGTTTCGGCAGAAAAATACGGCTGTTTTCCTGATCGAGCTTTACGCCCTGCGGGTAGCGGAATGCATCATTCTGTCCCCGCTTTTTGAATCGGGGAAAAGCAGCCCGCTTCCGgaagaa                                  <  1:364856/298‑1 (MQ=14)
                                                                                                                                                                                                                                                                  cggattttAGCCATGGCAACGTCTGGCCCGTGCATTGATGCTGAAGGATGAACCGCAGTTGATACTTCACTTTCTTTCTGAATTATTATGTACCCCTTACCGCAGGACTGGCTGACCGTGACATTTTTCACAACACCCGTGACCTGACGGCTATTCCGGTAGCGCATCCAGCACAGTTTCGGCAGAAAAATACGGCTGTTTTCCTGATCGAGCTTTACGCCCTGCGGGTAGCGGAATGCATCATTCTGTCCCCGCTTTTTGAATCGGGGAAAAGCAGCCCGCTTCCGGAAGAAGtttt                             <  2:320636/293‑1 (MQ=11)
                                                                                                                                                                                                                                                                       tgttgcccccgcAGCATCCAGCCAGGTCATTGATGTTGAAGGGTGAACCGGAGTTTATAATTCACTTTCTGTCTGAATACTGATGTACCACTTACCGCAGGACTGGCTAACCGTGACATTTTTCACAACACCCGTGACCTGACGGCTATTCCGGTAGCGCATCCAGCCCAGTTTCGGCAGAAAAATACGGCTGTTTTCCTGATCGAGCTTTACGCCCTGCGGGTAGCGGAATGCATCATTCTGTCCCCGCTTTTTGAATCGGGGAAAAGCAGCCCGCTTCCGGAAGAAGTTTTTGTaa                        <  2:410190/287‑1 (MQ=11)
                                                                                                                                                                                                                                                                                gCCAGCATCCAGCCCGATCATTGATGCTGAAGGGTGAACCGGAGTTGATACTTCACTTTCTGTCTGAATACTGATGTACCACTTACCGCAGGACTGGCTGACCGTGACATTTTTCACAACACCCGTGACCTGACGGCTATTCCGGTAGCGCATCCAGCCCAGCTTCGGCAGAAAAATACGGCTGTTTTCCTGATCGAGCTTTACGCCCTGCGGGTAGCGGAATGCATCATTCTGTCCCAGCTTTTTGAATTGGGGAAAGGCAGCCCGCTTCCGCGAGAATGTTTTGGAAGCCCGCCCa               >  2:320063/1‑298 (MQ=11)
                                                                                                                                                                                                                                                                                          aGCCCGATCATTGATGCTGAAGGGTGAACCGGAGTTGATACTTCACTTTCTGTCTGAATACTGATGTACCACTTACCGCAGGACTGGCTGACCGTGACATTTTTCACAACACCCGTGACCTGACGGCTATTCAGGTAGCGCATCCAGCCCAGCTTCGGCAGAAAAACACGGCTGTTTTCCTGATCGAGCTGTACGCCCTGCGGGTAGAGGAATGCATTATTCTGTCCCCGCTTTTTGAATCGGGGAAAAGCAGCCCGCTTCCAGAAGAAGGTTGTGTAAGCCCGCTCAAGGTttttga     >  2:442767/1‑296 (MQ=11)
                                                                                                                                                                                                                                                                                             ccGATCATTGATGCTGAAGGGTGAACCGGAGTTGATACTTCACTTTCTGTCTGAATACTGATGTACCACTTACCGCAGGACTGGCTGACCGTGACATTTTTCACAACACCCGTGACCTGACGGCTATTCCGGTAGCGCATCCAGCCCAGTTTCGGCAGAAAAATACGGCTGTTTTCCTGATCGAGCTTTACGCCCTGCGGGTAGCGGAATGCATCATTCTGTCCCCGCTTTTTGAATCGGGGAAAAGCAGCCCGCTTCCGGAAGAAGTTTTTGTAAGCCCGCTCAAGGTCTTTCAGTg  >  1:335266/1‑298 (MQ=14)
                                                                                                                                                                                                                                                                                                 |                                                                                                                                                                                                                                                                                                     
TTTGACATGTTGCTGACCTTCAAATCCTCAATGACTATCATTGCGTGGTTTTTGCTGACGGTCGTTGTGACTTTATGAAGGTAGTCCCTGCGGATATTTGCGATACGGGAATGCAGTCGCTGTATTTTGCGTTTCTGCTTCTGCCAGTTGTTGCTGAATTTGACCTTGCGGCTTAACTGTCGCTGAAGTCTCGCCAGCGTCTTCTGGTTTTTCTGAAAACTGTTTACAGGCCCAAAGACTGTGCCATCTGACAGCGTGGCGAGTTTAGCCACGCCAGCATCCAGCCCGACCATTGATGCTGAAGGGTGAACCGGAGTGGATACTTCACTTTCTGTCTGAATACTGATGTACCACTTACCGCAGGACTGGCTGACAGTGACATTTTTCACAACACCCGTGACCTGACGGCTATTCCGGTAGCGCATCCAGCCCAGTTTCGGCAGAAAAATACGGCTGTTTTCCTGATCGAGCTTTACGCCCTGCGGGTAGCGGAATGCATCATTCTGTCCCCGCTTTTTGAATCGGGGAAAAGCAGCCCGCTTCCGGAAGAAGTTTTTGTAAGCCCGCTCAAGGTCTTTCAGTG  >  CP006584/1964410‑1964992

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