Predicted mutation
evidence seq id position mutation annotation gene description
RA CP006584 1,588,241 A→C E403D (GAA→GAC LY180_07560 → type IV secretion protein Rhs

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0065841,588,2410AC95.7% 120.6 / NA 47E403D (GAA→GACLY180_07560type IV secretion protein Rhs
Reads supporting (aligned to +/- strand):  ref base A (2/0);  new base C (27/18);  total (29/18)
Fisher's exact test for biased strand distribution p-value = 5.17e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GGGGCATCCGCAGGCGAACCTGAACCGGGAATGGCAGGTGGTGGCGAGCGATCTTCACGGCGAACAGCCGCAGGCGGTACCGGGACGCAGGGGTTCAGGCACCACGCTGGATAACCATTTTGCGGTAATACCGGCAGACAGAACATGGCGACCACAGCCGTTGCTGAAACCGCTGGTGGACGGCCCGCAGAGCGCCGTCGTGACGGGACCGGCAGGCGAGGAAATCTTCTGTGATGAACATGGCCGCGTGCGGGTGAAATTTAACTGGGACCGCTATAACCCGTCAAACCAGGAAAGTTCATGCTGGATCCGTGTGGCACAGGCGTGGGCAGGCACCGGCTTTGGCAACCTGGCGATACCGCGAGTGGGTCAGGAGGTGATTGTGGACTTCCTCAACGGCGACCCGGACCAGCCGATTATTATGGGGCGTACCTACCACCATGAAAACCGCACACCCGGCAGCCTGCCGGGGACGAAGACGCAGATGACGATCCGCTCGAAAACCTATAAGGGCAGCGGGTTTAATGAACTGAAGTTTGACGATGCGACAGGGAAAGAACAGGTCTACATCCACGCGCAGAAGAACAT  >  CP006584/1587947‑1588534
                                                                                                                                                                                                                                                                                                      |                                                                                                                                                                                                                                                                                                      
ggggGATCCGGAGCGGGACCTGAAGCGGGAAGGGCAGGGGGGGGCGGGCGATGTTCACGGCGAACAGCCGCAGGCGGTAACGGGACGCAGGGGTGCGGGCACCACGCGGGATAACCATTTTGCGGTAATGCCGGGAGACAGAGCATGGCGACCACAGCCGTTGCTGAAACCGCTGGTGGACGGCGCGCAGAGCGCCGTCGTGACGGGACCGGCAGGCGAGGAAATCTTCTGTGATGAACATGGCCGCGTGCGGGTGAAATTTAACTGGGACCGCTATAACCCGTCAAACCAGGACAGtt                                                                                                                                                                                                                                                                                                    <  2:158378/299‑1 (MQ=11)
   gccgaggcgggggAACGTGAACGGGGGAGGGGGGGTGGGGGCGAGCGGTCTTCAGGGCGAACAGCGGCAGGGGGTACCGGGACGCAGGGGTTCAGGCAGCACGCTGGATAACCAGTTTGCGGTAATACCGGCAGACAGAACATGGCGACCACAGCCGTTGCTGAAACCGCTGGTGGACGGCCCGCAGAGCGCCGTCGTGACGGGACCGGCAGGCGAGGAAATCTTCTGTGATGAACATGGCCGCGTGCGGGTGAAATTTAACTGGGACCGCTATAACCCGTCAAACCAGGACAGTTCAt                                                                                                                                                                                                                                                                                                 <  2:356331/287‑1 (MQ=11)
    cATCCGCAGGCGAACCTGAACCGGGAATGGCAGGTGGTGGCGAGCGATCTTCACGGCGAACAGCCGCAGGCGGTACCGGGACGCAGGGGTTCAGGCACCACGCTGGATAACCATTTTGCGGTAAGACCGGCAGACAGAACATGGCGACCACAGACGTCGCTGAAACCGCTGCTCGACGGCCCGCGGAGAGCCGGCGTGGCGGGACCGGCAGGCGAGGACAACTCGTGTGATCAACACAGCCGAGTACGGGGGGAATTTAAATGGTGCCACTAGAACCAGTCCACCCACGCAAggcagct                                                                                                                                                                                                                                                                                                >  2:137485/1‑293 (MQ=11)
        cGCAGGCGAACCTGAACCGGGAATGGCAGGTGGTGGCGAGCGATCTTCACGGCGAACAGCCGCAGGCGGTACCGGGACGCAGGGGTTCAGGCACCACGCTGGATAACCATTTTGCGGTAATACCGGCAGACAGAACATGGCGACCACAGCCGTTGCTGAAACCGCTGGTGGACGGCCCGCAGAGCGCCGTCGTGACGGGACCGGCAGGCGAGGAAATCTTCTGTGATGAACATGGCCGCGTGCGGGTGAAATTTAACTGGGACCGCTATAACCCGTCAAACCAGGACAGTTCATGCTgg                                                                                                                                                                                                                                                                                            >  1:20457/1‑299 (MQ=25)
                aaCCTGAACCGGGAATGGCAGGTGGTGGCGAGCGATCTTCACGGCGAACAGCCGCAGGCGGTACCGGGACGCAGGGGTTCAGGCACCACGCTGGATAACCATTTTGCGGTAATACCGGCAGACAGAACATGGCGACCACAGCCGTTGCTGAAACCGCTGGTGGACGGCCGGCAGAGCGCCGTCGTGACGGGACCGGCAGGCGAGGAAATCTTCTGGGATGAACATGGCCGCGTGGGGGTGAAATTTAACTGGGACCGCATTAACCCGTCAAACCAGGACAGTTCATGCTGGATCCgtgt                                                                                                                                                                                                                                                                                    >  1:300453/1‑299 (MQ=255)
                               tGGCAGGTGGTGGCGAGCGATCTTCACGGCGAACAGCCGCAGGCGGTACCGGGACGCAGGGGTTCAGGCACCACGCTGGATAACCATTTTGCGGTAATACCGGCAGACAGAACATGGCGACCACAGCCGTTGCTGAAACCGCTGGTGGACGGCCCGCAGAGCGCCGTCGTGACGGGACCGGCAGGCGAGGAAATCTTCTGTGATGAACATGGCCGCGTGCGGGTGAAATTTAACTGGGACCGCTATAACCCGTCAAAACAGGACAGTTCAGGCGGGGTCCGTGTGGCCAAGGCGTGggc                                                                                                                                                                                                                                                                     >  1:86293/1‑299 (MQ=255)
                                  cAGGTGGTGGCGAGCGATCTTCACGGCGAACAGCCGCGGGCGGTACCGGGACGCAGGGGTTCAGGCACCACGCTGGATAACCATTTTGCGGTAATACCGGCAGACAGAACATGGCGACCACAGCCGTTGCTGAAACCGCTGGTGGACGGCCCGCAGAGCGCCGTCGTGACGGGACCGGCAGGCGAGGAAATCTTCTGTGATGAACATGGCCGCGTGCGGGTGAAATTTAACTGGGAACGCTATAACCCGTCAAACCAGGACAGTTCATGCTGGATCCAGGTGGCACAGGCGTGggcagg                                                                                                                                                                                                                                                                  >  1:622202/1‑299 (MQ=255)
                                  cAGGTGGTGGCGAGCGATCTTCACGGCGAACAGCCGCAGGCGGTACCGGGACGCAGGGGTCCAGGCACCACGCTGGATACCCATTTTGCGGTAATACCGGCAGACAGAACATGGCGCCCACAGCCGTTGCTGAACCCGCTGGTGGACGGCGCGCAGAGCGCCGTCGTGACGCACCCTGCAGGCGAGGAAACAGTCAGTGCTGAACATGGCCGCGTGCGGGTGCCCTCTAACTGAGAGCGCGATAACCCGCCAACCCAAGACAGATCACGCTGGTAACGTGTGGcacgcgccagcgcaga                                                                                                                                                                                                                                                                  >  2:95156/1‑286 (MQ=11)
                                                 gaTCTTCACGGCGAACAGCCGCAGGCGGTACCGGGACGCAGGGGTTCAGGCACCACGCTGGATAACCATTTTGCGGTAATACCGGCAGACAGAACATGGCGACCACAGCCGTTGCTGAAACCGCTGGTGGACGGCCCGCAGAGCGCCGTCGTGACGGGACCGGCAGGCGAGGAAATCTTCTGTGATGAACATGGCCGCGTGCGGGTGAAATTTAACTGGGACCGCTATAACCCGTCAAACCAGGACAGTTCATGCTGGATCCGTGTGGCACAGGCGTGGGCAGGCACCGGCTTTGGCaa                                                                                                                                                                                                                                                   >  1:108884/1‑299 (MQ=37)
                                                              aaCAGCCGCAGGCGGTACCGGGACGCAGGGGTTCAGGCACCACGCTGGATAACCATTTTGCGGTAATACCGGCAGACAGAACATGGCGACCACAGCCGTTGCTGAAACCGCTGGTGGACGGCCCGCAGAGCGCCGTCGTGACGGGACCGGCAGGCGAGGAAATCTTCTGTGATGAACATGGCCGCGTGCGGGTGAAATTTAACTGGGACCGCTATAACCCGTCAAACCAGGACAGTTCATGCTGGATCCGTGTGGCACAGGCGTGGGCAGGCACCGGCTTTGGCAACCTGGCGATACcg                                                                                                                                                                                                                                      >  1:546636/1‑299 (MQ=37)
                                                                             tACCGGGACGCAGGGGTTCAGGCATCACGCTGGATAACCATTTGGGGGTAATTCCGGCAGACAGAACATGGCGACCACCGCCGGGGCTGAAACCGCTGGTGGACGGCCCGCAGAGCGCCGTCGTGACGGGACCGGCAGGCGAGGAAATCTTCTGGGATGAACATGGCCGCGTGCGGGTGCAATGTAACTGGGACCGCTATAACCCGTCAAACCAGGACAGTTCATGCTGGATCCGTGTGGCACAGGCGTGGGCAGGCACCGGCTTTGGCAACCTGGCGATACCGCGAGTGGGTCaggat                                                                                                                                                                                                                       <  1:185327/299‑2 (MQ=11)
                                                                                  ggACGCAGGGGTTCAGGCACCACGCTGGATAACCATTTTGCGGTAATACCGGCAGACAGAACATGGCGACCACAGCCGTTGCTGAAACCGCTGGTGGACGGCCCGCAGAGCGCCGTCGTGACGGGACCGGCCGGCGAGGAAATCTTCTGTGATGAACATGGCCGCGTGCGGGTGAAATGTAACTCGGACCGCGATAACCCGTCAAACCAGGACAGTCCACGCTGGATCCGTGTGGCACCGGCGTGGGCCGGCCCCGGCTTTTGCAAACAGGCGCAACCGCGCCGGGGCCAGGAGGGGag                                                                                                                                                                                                                  >  2:260439/1‑298 (MQ=11)
                                                                                   gACGCAGGGGTTCAGGCACCACGCTGGATAACCATTTTGCGGTAATACCGGCAGACAGAACATGGCGACCACAGCCGTTGCTGAAACCGCTGGTGGACGGCCCGCAGAGCGCCGTCGTGACGGGACCGGCAGGCGAGGAAATCTTCTGTGATGAACATGGCCGCGTGCGGGTGAAATTTAACTGGGACCGCTATAACCCGTCAAACCAGGACAGTTCATGCTGGATCCGTGTGGCACAGGCGTGGGCAGGCACCGGCTTTGGCAACCTGGCGATAACGCGCGTGGGTAAGGCGGTGCtt                                                                                                                                                                                                                 >  1:644859/1‑299 (MQ=255)
                                                                                    aCGCAGGGGTTCAGGCACCACGCTGGATAACCATTTTGCGGTAATACCGGCAGACAGAACATGGCGACCACAGCCGTTGCTGAAACCGCTGGTGGACGGCCCGCAGAGCGCCGTCGTGACGGGACCGGCAGGCGAGGAAATCTTCTGTGATGAACATGGCCGCGTGCGGGTGAAATTTAACTGGGACCGCTATAACCCGTCAAACCAGGACAGTTCATGCTGGATCCGTGTGGCACAGGCGTGGGCAGGCACCGGCTTTGGCAACCTGGCGATACCGCGAGTGGGTCAGGAGGTGATtg                                                                                                                                                                                                                <  1:131626/299‑1 (MQ=39)
                                                                                                                      tttGCGGTAATACCGGCAGACAGAACATGGCGACCACAGCCGTTGCTGAAACCGCTGGTGGACGGCCCGCAGAGCGCCGTCGTGACGGGACCGGCAGGCGAGGAAATCTTCTGTGATGAACATGGCCGCGTGCGGGTGAAATTTAACTGGGACCGCTATAACCCGTCAAACCAGGACAGTTCATGCTGGATCCGTGTGGCACAGGCGTGGGCAGGCACCGGCTTTGGCAACCTGGCGATACCGCGAGTGGGTCAGGAGGTGATTGTGGACTTCCTCAACGGCGACCCGGACCAGCCGat                                                                                                                                                                              >  1:241730/1‑299 (MQ=39)
                                                                                                                                      cagacagaACATGGCGACCACAGCCGTTGCTGAAACCGCTGGTGGACGGCCCGCAGAGCGCCGTCGTGACGGGACCGGCAGGCGAGGAAATCTTCTGTGATGCACATGGCCGCGTGCGGGTGAAATTTAACTGGGACCGCTATAACCCGTCAAACCAGGACAGTTCATGCTGGATCCGTGTGGCACAg                                                                                                                                                                                                                                                                             >  2:124543/1‑188 (MQ=25)
                                                                                                                                      cagacagaACATGGCGACCACAGCCGTTGCTGAAACCGCTGGTGGACGGCCCGCAGAGCGCCGTCGTGACGGGACCGGCAGGCGAGGAAATCTTCTGTGATGAACATGGCCGCGTGCGGGTGAAATTTAACTGGGACCGCTATAACCCGTCAAACCAGGACAGTTCATGCTGGATCCGTGTGGCACAg                                                                                                                                                                                                                                                                             <  1:124543/188‑1 (MQ=25)
                                                                                                                                      cagacagaACATGGCGACCACAGCCGTTGCTGAAACCGCTGGTGGACGGCCCGCAGAGCGCCGTCGTGACGGGACCGGCAGGCGAGGAAATCTTCTGTGATGAACATGGCCGCGTGCGGGTGAAATTTAACTGGGACCGCTATAACCCGTCAAACCAGGACAGTTCATGCTGGATCCGTGTGGCACAGGCGTGGGCAGGCACCGGCTTTGGCAACCTGGCGATACCGCGAGTGGGTCAGGAGGTGATTGTGGACTTCCTCAACGGCGACCCGGACCAGCCGATTATTATGGGGCGtacc                                                                                                                                                              >  1:355431/1‑299 (MQ=39)
                                                                                                                                          cagaACATGGCGACCACAGCCGTTGCTGAAACCGCTGGTGGACGGCCCGCAGAGCGCCGTCGTGACGGGACCGGCAGGCGAGGAAATCTTCTGTGATGAACATGGCCGCGTGCGGGTGAAATTTAACTGGGACCGCTATAACCCGTCAAACCAGGACAGTTCATGCTGGATCCGTGTGGCACAGGCGTGGGCAGGCACCGGCTTTGGCAACCTGGCGATACCGCGAGTGGGTCAGGAGGTGATTGTGGACTTCCTCAACGGCGACCCGGACCAGCCGATTATTATGGGGCGtacctacc                                                                                                                                                          >  1:74083/1‑299 (MQ=38)
                                                                                                                                             aaCATGGCGACCACAGCCGTTGCTGAAACCGCTGGTGGACGGCCCGCAGAGCGCCGTCGTGACGGGACCGGCAGGCGAGGAAATCTTCTGTGATGAACATGGCCGCGTGCGGGTGAAATTTAACTGGGACCGCTATAACCCGTCAAACCAGGACAGTTCATGCTGGATCCGTGTGGCACAGGCGTGGGCAGGCACCGGCTTTGGCAACCTGGCGATACCGCGAGTGGGTCAGGAGGTGATTGTGGACTTCCTCAACGGCGACCCGGACCAGCCGATTATTATGGGGCGTACCCaccacc                                                                                                                                                       >  1:516558/1‑299 (MQ=37)
                                                                                                                                               cATGGCGACCACAGCCGTTGCTGAAACCGCTGGTGGACGGCCCGCAGAGCGCCGTCGTGACGGGACCGGCAGGCGAGGAAATCTTCTGTGATGAACATGGCCGCGTGCGGGTGAAATTTAACTGGGACCGCTATAACCCGTCAAACCAGGACAGTTCATGCTGGATCCGTGTGGCACAGGCGTGGGCAGGCACCGGCTTTGGCAACCTGGCGATACCGCGAGTGGGTCAGGAGGTGATTGTGGACTTCCTCAACGGCGACCCGGACCAGCCGATTATTATGGGGCGTACCTACCACCAt                                                                                                                                                     <  1:564791/299‑1 (MQ=37)
                                                                                                                                                        cacaGCCGTTGCTGAAACCGCTGGTGGACGGCCCGCAGAGCGCCGTCGTGACGGGACCGGCAGGCGAGGAAATCTTCTGTGATGAACATGGCCGCGTGCGGGTGAAATTTTACTGGGACCGCTATAACCCGTCCAACCAGGAAAGTTCATGCTGGATCCGTGTGGCACAGGCGTGGGCAGGCACCGGCTTTGGCAACCTGGCGATACCGCGAGTGGGTCAGGAGGTGATTGTGGACTTCCTCAACGGCGACACGGACCAGCCGATTATTATGGGGCGTACCTACCACCATGAAAACCGc                                                                                                                                            >  1:560700/1‑299 (MQ=25)
                                                                                                                                                                         ccGCTGGTGGACGGCCCGCAGAGCGCCGTCGGGACGGGACCGGCAGGCGAGGAAATCTTCTGTGATGAACATGGCCGCGTGCGGGTGAAATTTAACTGGGACCGCTATAACCCGTCAAACCAGGACAGTTCATGCTGGATCCGTGTGGCACAGGCGTGGGCAGGCACCGGCTTTGGCAACCTGGCGATACCGCGAGTGGGTCAGGAGGTGATTGTGGACTTCCTCAACGGCGACCCGGACCAGCCGATTATTATGGGGCGTACCTACCACCATGAAAACCGCACACCCGGCAGCCTGcc                                                                                                                           <  1:243335/299‑1 (MQ=37)
                                                                                                                                                                                tggACGGCCCGCAGAGCGCCGTCGTGACGGGACCGGCAGGCGAGGAAATCTTCTGTGATGAACATGGCCGCGTGCGGGTGAAATTTAACTGGGACCGCTATAACCCGTCAAACCAGGACAGTTCATGCTGGATCCGTGTGGCACAGGCGTGGGCAGGCACCGGCTTTGGCAACCTGGCGATACCGCGAGTGGGTCAGGAGGTGATTGTGGACTTCCTCAACGGCGACCCGGACCAGCCGATTATTATGGGGCGTACCTACCACCATGAAAACCGCACACCCGGCAGCCTGCCGGGGACg                                                                                                                    >  1:140541/1‑299 (MQ=37)
                                                                                                                                                                                                cgcCGTCGTGACGGGACCGGCAGGCGAGGAAATCTTCTGTGATGAACATGGCCGCGTGCGGGGGAAATGTAACTGGGACCGCTATAACCCGTCAAACCAGGACAGTTCATGCTGGATCCGTGTGGCCCCGGCGTGGGCCGGCACCGGCTTTGGCAACCTGGCGATTCCGCGCGTGGGTCAGGCGGGGATTGTGGACTTCCGCAACGGCCACCCGGCACAGCCGATTATTATGGGGCGTACCTACCACCTCCAAAAACGCACCCCCGGACGCCGCCCGGCGCCGACGCCGCAGATGACGa                                                                                                    >  2:215250/1‑299 (MQ=11)
                                                                                                                                                                                                  ccGTCGTGACGGGACCGGCAGGCGAGGAAATCTTCTGTGATGAACATGGCCGCGTGCGGGTGAAATTTAACTGGGACCGCTATAACCCGTCAAACCAGGACAGTTCATGCTGGATCCGTGTGGAACAGGCGTGGGCAGGCACCGGCTTTGGCAACCTGGCGATACCGCGAGTGGGTCAGGAGGTGATTGTGGACTTCCTCAACGGCGACCCGGACCAGCCGATTATTATGGGGCGTACCTACCACCATGAAAACCGCACACCCGGCAGCCTGCCGGGGACGAAGACGCAGATGACGATc                                                                                                  <  1:531637/299‑1 (MQ=37)
                                                                                                                                                                                                   cgtcgtGACGGGACCGGCAGGCGAGGAAATCTTCTGTGATGAACATGGCCGCGTGCGGGTGAAATTTAACTGGGACCGCTATAACCCGTCAAACCAGGACAGTTCATGCTGGATCCGTGTGGCACAGGCGTGGGCAGGCACCGGCTTTGGCAACCTGGCGATACCGCGAGTGGGTCAGGAGGTGATTGTGGACTTCCTCAACGGCGACCCGGACCAGCCGATTATTATGGGGCGTACCTACCACCATGAAAACCGCACACCCGGCAGCCTGCCGGGGACGAAGACGCAGATGACGATcc                                                                                                 <  1:260439/299‑1 (MQ=37)
                                                                                                                                                                                                                              aaTCTTCTGTGATGAACATGGCCGCGTGCGGGTGAAATTTAACTGGGACCGCTATAACCCGTCAAACCAGGACAGTTCATGCTGGATCCGTGTGGCACAGGCGTGGGCAGGCACCGGCTTTGCCAACCTGGCGATACCGCGAGTGGGTCAGGCGGTGATTGTGGACTTCCTCAACGGCGAACCGGACCAGCCGATTAGTATGGGGCGTACCTACCACCATGAAAACCGCACACCCGGCAGCCTGCCGGGGCCGAAGACGCAGATGACGATCCGCTCGAAAACCTATAAGGTCAGCggga                                                                      >  2:230086/1‑298 (MQ=255)
                                                                                                                                                                                                                               aTCTTCTGTGATGAACATGGCCGCGTGCGGGTGAAATTTAACTGGGACCGCTATAACCCGTCAAACCAGGACAGTTCATGCTGGATCCGTGTGGCACAGGCGTGGGCAGGCACCGGCTTTGGCAACCTGGCGATACCGCGAGTGGGTCAGGAGGTGATTGTGGACTTCCTCAACGGCGACCCGGACCAGCCGATTATTATGGGGCGTACCTACCACCATGAAAACCGCACACCCGGCAGCCTGCCGGGGACGAAGACGCAGATGACGATCCGCTCGAAAACCTATAAGGGCAGCGGGtt                                                                     >  1:547519/1‑299 (MQ=37)
                                                                                                                                                                                                                                    ctGTGATGAAAATGGCCGCGTGCGGGTGAAATTTAACTGGCACCCGTCTAACCCGTCAAGCCAGGACAGTTCATGCTGGATCCGTGTGGCCCAGGCGGGGGCAGGCACCGGCTTTGGCAACCTGGGGATACCGCGAGTGGGTCAGGAGGTGATTGTGGACTTCCTCAACGGCGACCCGGACCAGCCGATTATTATGGGGCGTACCTACCACCATGAAAACCGCACACCCGGCAGCCTGCCGGGGACGAAGACGCAGATGACGATCCGCTCGAAAACCTATAAGGGCAGCGGGTTTAATg                                                                <  2:241730/299‑1 (MQ=11)
                                                                                                                                                                                                                                       tgatgaACATGGCCGCGTGCGGGTGAAATTTAACTGGGACCGCTATAACCCGTCAAACCAGGACAGTTCATGCTGGATCCGTGTGGCACAGGCGTGGGCAGGCACCGGCTTTGGCAACCTGGCGATACCGCGAGTGGGTCAGGAGGTGATTGTGGACTTCCTCAACGGCGACCCGGACCAGCCGATTATTATGGGGCGTACCTACCACCATGAAAACCGCACACCCGGCAGCCTGCCGGGGACGAAGACGCAGATGACGATCCGCTCGAAAACCTATAAGGGCAGCGGGTTTAATGAAc                                                             >  1:72289/1‑299 (MQ=37)
                                                                                                                                                                                                                                               atgtcctggtccGGGTGGAATTTAAGTGGGACGGCGATGACCGGTCAGACCAGGACAGTTCATGCTGGGTCCGGGTGGAACAGGCGTGGGCAGGGAACGGCTTTGGCATCCTGGCGATACCGCTAGTGGGTCAGGAGGTGAGTGTGGACTTCCTCGACGGCGACCCGGACCAGCCGATTATTATGGGGCGGACCTGACACCATGAAAACCGCACACCCGGCAGCCTGCCGGGGACGAAGACGCAGATGACGATCCGCTCGAAAACCTATAAGGGCGGCGGGTTTAATGAACTGAAGttt                                                     <  2:108884/288‑1 (MQ=11)
                                                                                                                                                                                                                                                            ggTGAAATTTAACTGGGACCGCTATAACCCGTCAAACCAGGACAGTTCATGCTGGATCCGTGTGGCACAGGCGTGGGCAGGCACCGGCTTTGGCAACCTGGCGATACCGCGAGTGGGTCAGGAGGTGATTGTGGACTTCCTCAACGGCGACCCGGACCAGCCGATTATTATGGGGCGTACCTACCACCATGAAAACCGCACACCCGGCAGCCTGCCGGGGACGAAGACGCAGATGACGATCCGCTCGAAAACCTATAAGGGCAGCGGGTTTAATGAACTGAAGTTTGACGATGCGACAg                                       <  1:383157/299‑1 (MQ=37)
                                                                                                                                                                                                                                                               gAAATTTAACTGGGACCGCTATAACCCGTCAAACCAGGACAGTTCATGCTGGATCCGTGTGGCACAGGCGTGGGCAGGCACCGGCTTTGGCAACCTGGCGATACCGCGAGTGGGTCAGGAGGTGATTGTGGACTTCCTCAACGGCGACCCGGACCAGCCGATTATTATGGGGCGTACCTACCACCATGAAAACCGCACACCCGGCAGCCTGCCGGGGACGAAGACGCAGATGACGATCCGCTCGAAAACCTATAAGGGCAGCGGGTTTAATGAACTGAAGTTTGACGATGCGACAGGGa                                    <  1:585178/299‑1 (MQ=37)
                                                                                                                                                                                                                                                                        cTGGGACCGCTAGAACCCGTCAAGCCAGAACAGTTCATGCTGGATCCGTGTGGCACAGGCGTGTGCAGGCACCGGCTTTGGCAACCTGGCGATACCGCGAGTGGGTCAGGAGGTGATTGTGGACTTCCTCAACGGCGACCCGGACCAGCCGATTATTATGGGGCGTACCTACCACCATGAAAACCGCACACCCGGCAGCCTGCCGGGGACGAAGACGCAGATGACGATCCGCTCGAAAACCTATAAGGGCAGCGGGTTTAATGAACTGAAGTTTGACGATGCGACAGGGAAAGAACAgg                           <  1:621737/299‑1 (MQ=255)
                                                                                                                                                                                                                                                                          gggACCGCTATAACCCGTCAAACCAGGACAGTTCATGCTGGATCCGTGTGGCACAGGCGTGGGCAGGCACCGGCTTTGGCAACCTGGCGATACCGCGAGTGGGTCAGGAGGTGATTGTGGACTTCCTCAACGGCGACCCGGACCAGCCGATTATTATGGGGCGTACCTACCACCATGAAAACCGCACACCCGGCAGCCTGCCGGGGACGAAGACGCAGATGACGATCCGCTCGAAAACCTATAAGGGCAGCGGGTTTAATGAACTGAAGTTTGACGATGCGACAGGGAAAGAACAGGTc                         >  1:318583/1‑299 (MQ=37)
                                                                                                                                                                                                                                                                            gACCGCTATAACCCGTCAAACCAGGACAGGTCATGCTGGATCCGTGTGGCAAAGGCGTGGGCAGGCACCGGCTTTGGCAACCTGGCGATACCGCGAGTGGGTCAGGAGGTGATTGTGGACTTCCTCAACGGCGACCCGGACCAGCCGATTATTATGGGGCGTACCTACCACCATGAAAACCGCACACCCGGCAGCCTGCCGGGGACGAAGACGCAGATGACGATCCGCTCGAAAACCTATAAGGGCAGCGGGTTTAATGAACTGAAGTTTGACGATGCGACAGGGAAAGAACAGGTCTa                       <  1:95156/299‑1 (MQ=37)
                                                                                                                                                                                                                                                                               cGCTATAACCCGTCAAACCAGGACAGTTCATGCTGGATCCGTGTGGCACAGGCGTGGGCAGGCACCGGCTTTGGCAACCTGGCGATACCGCGAGTGGGTCAGGAGGTGATTGTGGACTTCCTCAACGGCGACCCGGACCAGCCGATTATTATGGGGCGTACCTACCACCATGAAAACCGCACACCCGGCAGCCTGCCGGGGACGAAGACGCCGCGGACGATCCGCCCGCAACCCCAAAAGGGCACCGGGCTTAATGAACTGAAGGTTGACGATGCCCAAGGGAAAGAAAAggacaacca                    >  2:267050/1‑292 (MQ=11)
                                                                                                                                                                                                                                                                                 cTATAACCCGTCAAACCAGGACAGTTCATGCTGGATCCGTGTGGCACAGGCGTGGGCAGGCACCGGCTTTGGCAACCTGGCGATACCGCGAGTGGGTCAGGAGGTGATTGTGGACTTCCTCAACGGCGACCCGGACCAGCCGATTATTATGGGGCGTACCTACCACCATGAAAACCGCACACCCGGCAGCCTGCCGGGGACGAAGACGCAGATGACGATCCGCTCGAAAACCTATAAGGGCAGCGGGTTTAATGAACTGAAGTTTGACGATGCGACAGGGAAAGAACAGGTCTACATcc                  >  1:510932/1‑299 (MQ=37)
                                                                                                                                                                                                                                                                                   catggctacaaaaCCAGACAGGTCAGGCTGGGTCCGGGTGGCACAGGCGTGGGCATGCACCGGCTTTGGGAACCTGGCGGTGCCGCGAGTGGGTCAGGAGGGGATGGTGGACTGCCTCAACGGCTACCCGGACCAGCCGATTATTATGGGGCGTACCTACCACCATGAAAACCGCACACCCGGCAGCCTGCCGGGGACGAAGACGCAGATGACGATCCGCTCGAAAACCTATAAGGGCAGCGGGTTTAATGAACTGAAGTTTGACGATGCGACAGGGAAAGAACAGGTCTACATCCAcg               <  2:546636/289‑1 (MQ=11)
                                                                                                                                                                                                                                                                                        ccGTCAAACCAGGACAGTTCATGCTGGATCCGTGTGGCACAGGCGTGGGCAGGCACCGGCTTTGGCAACCTGGCGATACCGCGAGTGGGTCAGGAGGTGATTGTGGACTTCCTCAACGGCGACCCGGACCAGCCGATTATTATGGGGCGTACCTACCACCATGAAAACCGCACACCCGGCAGCCTGCCGGGGACGAAGACGCAGATGACGATCCGCTCGAAAACCTATAAGGGCAGCGGGTTTAATGAACTGAAGTTTGACGATGCGACAGGGAAAGAACAGGTCTACATCCACGCGCa           >  1:226220/1‑299 (MQ=37)
                                                                                                                                                                                                                                                                                        ccGTCAAACCAGGACAGTTCATGCTGGATCCGTGTGGCACAGGCGTGGGCAGGCACCGGCTTTGGCAACCTGGCGATACCGCGAGTGGGTCAGGAGGTGATTGTGGACTTCCTCAACGGCGACCCGGACCAGCCGATTATTATGGGGCGTACCTACCACCATGAAAACCGCACACCCGGCAGCCTGCCGGGGACGAAGACGCAGATGACGATCCGCCCAAAACCCCACAAGGGCAGGGGGTTTAATGAACCGGAGTTTGACGATGCGACAGGGAGAGAACAGGCCGACATCCAcgaaca           >  2:267432/1‑295 (MQ=11)
                                                                                                                                                                                                                                                                                           tCAAACCAGGACAGTTCATGCTGGATCCGTGTGGCACAGGCGTGGGCAGGCACCGGCTTTGGCAACCTGGCGATACCGCGAGTGGCACAGCAGGTGATGGTGGACTCCCTCAACCGCGCCCCGGACGATCCGATTATAATGGGGCGTACCTACCCCCAGGAAAACAGCACCACAGGCACGCCACCGGTGACAAAGACGACGATGACGATCCCCCATAAAACCTATAAGGGCCCCGGGACTACTAAACTGAAGCGTAACGATCACACAGCGCAAACACGGATCTACACGCACGCCCagaa        >  2:461477/1‑299 (MQ=2)
                                                                                                                                                                                                                                                                                            cAAACCAGGACAGTTCATGCTGGATCCGTGTGGCACAGGCGTGGGCAGGCACCGGCTTTGGCAACCTGGCGATACCGCGAGTGGGTCAGGAGGTGATTGTGGACTTCCTCAACGGCGACCCGGACCAGCCGATTATTATGGGGCGTACCTACCACCATGAAAACCGCACACCCGGCAGCCTGCCGGGGACGAAGACGCAGATGACGATCCGCTCGAAAACCTATAAGGGCAGCGGGTTTAATGAACTGAAGTTTGACGATGCGACAGGGAAAGAACAGGTCTACATCCACGCGCagaag       <  1:137485/299‑1 (MQ=37)
                                                                                                                                                                                                                                                                                              aaCCAGGACAGTTCATGCTGGATCCGTGTGGCACAGGCGTGGGCAGGCACCGGCTTTGGCAACCTGGCGATACCGCGAGTGGGTCAGGAGGTGATTGTGGACTTCCTCAACGGCGACCCGGACCAGCCGATTATTATGGGGCGTACCTACCACCACGAAAACCGCACACCCGGCAGCCTGCCGGGGACGCAGACGAGGATGCCGATCCGCTCGAAACCCTAAAAGGGCACCGGGGTTAATGAACTGAAGTTTGACGATGCCAACGGGAAAAAAACAGACGACATCCACGGCCAgaagaa     >  2:437291/1‑299 (MQ=11)
                                                                                                                                                                                                                                                                                               aCCAGGACAGTTCATGCTTGATCCGTGTGGCACAGGCGTGGGCCGGCACCGTCTTTGGCAACCTGGCGATACCGCTAGTGGGTCAGGAGGTGATTGTGGACCTCCTCAACCGCGACCCGGACCAGCCGATTATTATGGGGCGTACCTACCACCATGAAAACCGCACACCCGGCAGCCTGCCGGGGACGAAGACGCAGATGACGATCCGCTCGAAAACCTATAAGGGCAGCGGGTTTAATGAACTGAAGTTTGACGATGCGACAGCGAACGAACAGGTCTACATCCACGCGCAGAAGAAc    <  1:392949/299‑1 (MQ=255)
                                                                                                                                                                                                                                                                                                 cAGGACAGTTCATGCTGGATCCGTGTGGCACAGGCGTGGGCAGGCACCGGCTTTGGCAACCTGGCGATACCGCGAGAGGGTCAGGAGGTGATTGTGGACTTCCTCAACGGCGACCCGGCCCAGCCGATTATTATGGGGCGTACCTACCACCACGAAAACCGCACACCCGGCTGCCTGCCGGGGACGAAGACGCAGATGACGCTCCGCTCGGCAACCTATCAGGGCCGGGGGTTGAAGTAACTGAAGTTGACGCATGCGAAAGGGAACGAACCGCTCTAAATACACGACCAGAAGAACaa  >  2:513167/1‑298 (MQ=11)
                                                                                                                                                                                                                                                                                                      |                                                                                                                                                                                                                                                                                                      
GGGGCATCCGCAGGCGAACCTGAACCGGGAATGGCAGGTGGTGGCGAGCGATCTTCACGGCGAACAGCCGCAGGCGGTACCGGGACGCAGGGGTTCAGGCACCACGCTGGATAACCATTTTGCGGTAATACCGGCAGACAGAACATGGCGACCACAGCCGTTGCTGAAACCGCTGGTGGACGGCCCGCAGAGCGCCGTCGTGACGGGACCGGCAGGCGAGGAAATCTTCTGTGATGAACATGGCCGCGTGCGGGTGAAATTTAACTGGGACCGCTATAACCCGTCAAACCAGGAAAGTTCATGCTGGATCCGTGTGGCACAGGCGTGGGCAGGCACCGGCTTTGGCAACCTGGCGATACCGCGAGTGGGTCAGGAGGTGATTGTGGACTTCCTCAACGGCGACCCGGACCAGCCGATTATTATGGGGCGTACCTACCACCATGAAAACCGCACACCCGGCAGCCTGCCGGGGACGAAGACGCAGATGACGATCCGCTCGAAAACCTATAAGGGCAGCGGGTTTAATGAACTGAAGTTTGACGATGCGACAGGGAAAGAACAGGTCTACATCCACGCGCAGAAGAACAT  >  CP006584/1587947‑1588534

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