Predicted mutation
evidence seq id position mutation freq annotation gene description
RA AE015451 4,475,598 T→A 7.6% intergenic (‑63/+1) PP_3964 ← / ← PP_3965 transposase/transposase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*AE0154514,475,5980TA 7.6% 178.6 / 6.4 53intergenic (‑63/+1)PP_3964/PP_3965transposase/transposase
Reads supporting (aligned to +/- strand):  ref base T (16/32);  new base A (2/2);  total (19/34)
Fisher's exact test for biased strand distribution p-value = 6.02e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.05e-02

CGTTTGAGCAGGGCGATTTCGAACTTGAACTGTTCGATGAGGATTTCATCGTTCTGAATCTTCCTGCTCATCGCCTCGACCTGGGATTGCAATTGCAACGCCTGTGCCGCCAAGGCACGAAGCTGTTCCGGGGTCATCTGGTCGAGGTTGGGCGAGAAAGTCATGCCGCCGATTGTGCCAGAGCAGGCGCGGGGCGACGACAAGTAGACCGGCCAAATGGCCGGCGTTTACAGCATGCTAATCACACCGCCTGCACCGACGCGCTGCCACGGCAGGCCCAGCACCAACGCCTGGAGTTGTTCGCTGTCGAGTTCGACCTCGCAGCCGTGGCGAATGCCGGGCCAGTGGAACTTGCCTTGGTTCAGTCGACGCGCTGCCAGCCAGACACCGACGCCATCGTGCACCAGCACTTTCATCCGGTTAGCCCGGCGATTGGC  >  AE015451/4475372‑4475808
                                                                                                                                                                                                                                  |                                                                                                                                                                                                                  
cGTTTGAGCAGGGCGATTTCGAACTTGAACTGTTCGATGAGGATTTCATCGTTCTGAATCTTCCTGCTCATCGCCTCGACCTGGGATTGCAATTGCAACGCCTGTGCCGCCAAGGCACGAAGCTGTTCCGGGGTCATCTGGTCGAGGTTGGGCGAGAAAGTCATGCCGCCGATTGTGCCAGAGCAGGCGCGGGGCGACGACAAGTAGACCGGCCAAATGGCCGGCGtt                                                                                                                                                                                                                   >  1:262825/1‑228 (MQ=255)
cGTTTGAGCAGGGCGATTTCGAACTTGAACTGTTCGATGAGGATTTCATCGTTCTGAATCTTCCTGCTCATCGCCTCGACCTGGGATTGCAATTGCAACGCCTGTGCCGCCAAGGCACGAAGCTGTTCCGGGGTCATCTGGTCGAGGTTGGGCGAGAAAGTCATGCCGCCGATTGTGCCAGAGCAGGCGCGGGGCGACGACAAGTAGACCGGCCAAATGGCCGGCGTTTACAGCAt                                                                                                                                                                                                           >  1:21414/1‑236 (MQ=255)
                             cTGTTCGATGAGGATTTCATCGTTCTGAATCTTCCTGCTCATCGCCTCGACCTGGGATTGCAATTGCAACGCCTGTGCCGCCAAGGCACGAAGCTGTTCCGGGGTCATCTGGTCGAGGTTGGGCGAGAAAGTCATGCCGCCGATTGTGCCAGAGCAGGCGCGGGGCGACGACAAGTAGACCGGCCAAATGGCCGGCGTTTACAGCATGCTAATCACACCGCCTGCACCGa                                                                                                                                                                                    <  2:587541/230‑1 (MQ=255)
                              tGTTCGATGAGGATTTCATCGTTCTGAATCTTCCTGCTCATCGCCTCGACCTGGGATTGCAATTGCAACGCCTGTGCCGCCAAGGCACGAAGCTGTTCCGGGGTCATCTGGTCGAGGTTGGGCGAGAAAGTCATGCCGCCGATTGTGCCAGAGCAGGCGCGGGGCGACGACAAGTAGACCGGCCAAATGGCCGGCGTTTACAGCATGCTAATCACACCGCCTGCACCGACGCGCt                                                                                                                                                                              <  1:925664/235‑1 (MQ=255)
                                       aGGATTTCATCGTTCTGAATCTTCCTGCTCATCGCCTCGACCTGGGATTGCAATTGCAACGCCTGTGCCGCCAAGGCACGAAGCTGTTCCGGGGTCATCTGGTCGAGGTTGGGCGAGAAAGTCATGCCGCCGATTGTGCCAGAGCAGGCGCGGGGCGACGACAAGTAGACCGGCCAAATGGCCGGCGTTTACAGCATGCTAATCACACCGCCTGCACCGACGCGCTGCCACGGCAg                                                                                                                                                                    <  1:45216/236‑1 (MQ=255)
                                              cATCGTTCTGAATCTTCCTGCTCATCGCCTCGACCTGGGATTGCAATTGCAACGCCTGTGCCGCCAAGGCACGAAGCTGTTCCGGGGTCATCTGGTCGAGGTTGGGCGAGAAAGTCATGCCGCCGATTGTGCCAGAGCAGGCGCGGGGCGACGACAAGTAGACCGGCCAAATGGCCGGCGTTTACAGCATGCTAATCACACCGCCTGCACCGACGCGCTGCCACGGCAGGCCCAGc                                                                                                                                                             <  1:904666/236‑1 (MQ=255)
                                                       gAATCGTCCTGCTCATCGCCTCGACCTGGGATGGCAATTGCAACGCCTGTGCCGCCAAGGCACGAAGCTGTTCCGGGGTCATCTGGTCCAGGTTGGGCGAGAAAGTCATGCCGCCGATTGTGCCAGAGCAGGCGCGGGGCGACGACAAGTAGCCCGGCCAAATGGCCGGCGTTTACAg                                                                                                                                                                                                              <  2:580586/178‑1 (MQ=255)
                                                         aTCTTCCTGCTCATCGCCTCGACCTGGGATTGCAATTGCAACGCCTGTGCCGCCAAGGCACGAAGCTGTTCCGGGGTCATCTGGTCGAGGTTGGGCGAGAAAGTCATGCCGCCGATTGTGCCAGAGCAGGCGCGGGGCGACGACAAGTAGACCGGCCAAATGGCCGGCGTTTACAGCATGCTAATCACACCGCCTGCACCGACgcgc                                                                                                                                                                               <  1:46471/207‑1 (MQ=255)
                                                               cTGCTCATCGCCTCGACCTGGGATTGCAATTGCAACGCCTGTGCCGCCAAGGCACGAAGCTGTTCCGGGGTCATCTGGTCGAGGTTGGGCGAGAAAGTCATGCCGCCGATTGTGCCAGAGCAGGCGCGGGGCGACGACAAGTAGACCGGCCAAATGGCCGGCGTTTACAGCATGCTAATCACACCGCCTGCAc                                                                                                                                                                                       <  2:21414/193‑1 (MQ=255)
                                                                tGCTCATCGCCTCGACCTGGGATTGCAATTGCAACGCCTGTGCCGCCAAGGCACGAAGCTGTTCCGGGGTCATCTGGTCGAGGTTGGGCGAGAAAGTCATGCCGCCGATTGTGCCAGAGCAGGCGCGGGGCGACGACAAGTAGACCGGCCAAATGGCCGGCGTTTACAGCATGCTAATCACACCGCCTGCACCGACGCGCTGCCACGGCAGGCCCAGCACCAACGCCTGGAgttgt                                                                                                                                           <  1:504678/236‑1 (MQ=255)
                                                                      tCGCCTCGACCTGGGATTGCAATTGCAACGCCTGTGCCGCCAAGGCACGAAGCTGTTCCGGGGTCATCTGGTCGAGGTTGGGCGAGAAAGTCATGCCGCCGATTGTGCCAGAGCAGGCGCGGGGCGACGACAAGTAGACCGGCCAAATGGCCGGCGTTTACAGCATGCTAATCACACCGCCTGCACCGACGCGCTGCCACGGCAGGCCCAGCACCAACGCCTGGAGTTGTTCGc                                                                                                                                       <  1:758970/234‑1 (MQ=255)
                                                                        gCCTCGACCTGGGATTGCAATTGCAACGCCTGTGCCGCCAAGGCACGAAGCTGTTCCGGGGTCATCTGGTCGAGGTTGGGCGAGAAAGTCATGCCGCCGATTGTGCCAGAGCAGGCGCGGGGCGACGACAAGTAGACCGGCCAAATGGCCGGCGTTTACAGCATGCTAATCACACCGCCTGCACCGACGCGCTGCCAc                                                                                                                                                                         <  2:556573/198‑1 (MQ=255)
                                                                        gCCTCGACCTGGGATTGCAATTGCAACGCCTGTGCCGCCAAGGCACGAAGCTGTTCCGGGGTCATCTGGTCGAGGTTGGGCGAGAAAGTCATGCCGCCGATTGTGCCAGAGCAGGCGCGGGGCGACGACAAGTAGACCGGCCAAATGGCCGGCGATTACAGCATGCTGATCACGCCACCCGCGCCGACTCGCt                                                                                                                                                                              <  2:745361/193‑1 (MQ=11)
                                                                        gCCTCGACCTGGGATTGCAATTGCAACGCCTGTGCCGCCAAGGCACGAAGCTGTTCCGGGGTCATCTGGTCGAGGTTGGGCGAGAAAGTCATGCCGCCGATTGTGCCAGAGCAGGCGCGGGGCGACGACAAGTAGACCGGCCAAATGGCCGGCGATTACAGCATGCTGATCACGCCACCCGCGCCGACTCGCt                                                                                                                                                                              >  1:745364/1‑193 (MQ=11)
                                                                            cGACCTGGGATTGCAATTGCAACGCCTGTGTCGCCAAGGCACGAAGCTGTTCCGGGGTCATCTGGTCGAGGTTGGGCGAGAAAGTCATGCCGCCGATTGTGCCAGAGCAGGCGCGGGGCGACGACAAGTAGACCGGCCAAATGGCCGGCGt                                                                                                                                                                                                                    <  2:345750/151‑1 (MQ=255)
                                                                               ccggggATTGCAATTGCAACGCCTGTGCCGCCAAGGCACGAAGCTGTTCCGGGGTCATCTGGTCGAGGTTGGGCGAGAAAGTCATGCCGCCGATTGTGCCAGAGCAGGCGCGGGGCGACGACAAGTAGACCGGCCAAATGGCCGGCGTTTACAg                                                                                                                                                                                                              <  2:467110/151‑1 (MQ=255)
                                                                                 tGGGATTGCAATTGCAACGCCTGTGCCGCCAAGGCACGAAGCTGTTCCGGGGTCATCTGGTCGAGGTTGGGCGAGAAAGTCATGCCGCCGATTGTGCCAGAGCAGGCGCGGGGCGACGACAAGTAGACCGGCCAAATGGCCGGCGttt                                                                                                                                                                                                                  <  2:647984/148‑1 (MQ=255)
                                                                                    gATTGCAATTGCAACGCCTGTGCCGCCAAGGCACGAAGCTGTTCCGGGGTCATCTGGTCGAGGTTGGGCGAGAAAGTCATGCCGCCGATTGTGCCAGAGCAGGCGCGGGGCGACGACAAGTAGACCGGCCAAATGGCCGGCGTTTACAGCATGCTAATCACACCGCCTGCACCGACGCGCTGCCAc                                                                                                                                                                         <  2:858882/186‑1 (MQ=255)
                                                                                           attgcaaCGCCGGTGCCGCCAAGGCACGAAGCTGTTCCGGGGTCATCTGGTCGAGGTTGGGCGAGAAAGTCATGCCGCCGATTGTGCCAGAGCAGGCGCGGGGCGACGACAAGTAGACCGGCCAAATGGCCGGCGTTTACAGCATGCTAATCACACCGCCTGCACCGACGCGCTGCCACGGCAGGCCCAGCACCAACGc                                                                                                                                                     <  2:22607/199‑1 (MQ=255)
                                                                                            ttgcaaCGCCTGTGCCGCCAAGGCACGAAGCTGTTCCGGGGTCATCTGGTCGAGGTTGGGCGAGAAAGTCATGCCGCCGATTGTGCCAGAGCAGGCGCGGGGCGACGACAAGTAGACCGGCCAAATGGCCGGCGTTTACAGCATGCTGATCACGCCACCCGCGCCGACTCGCt                                                                                                                                                                              <  2:812989/173‑1 (MQ=11)
                                                                                            ttgcaaCGCCTGTGCCGCCAAGGCACGAAGCTGTTCCGGGGTCATCTGGTCGAGGTTGGGCGAGAAAGTCATGCCGCCGATTGTGCCAGAGCAGGCGCGGGGCGACGACAAGTAGACCGGCCAAATGGCCGGCGTTTACAGCATGCTGATCACGCCACCCGCGCCGACTCGCt                                                                                                                                                                              >  1:812992/1‑173 (MQ=11)
                                                                                            ttgcaaCGCCTGTGCCGCCAAGGCACGAAGCTGTTCCGGGGTCATCTGGTCGAGGTTGGGCGAGAAAGTCATGCCGCCGATTGTGCCAGAGCAGGCGCGGGGCGACGACAAGTAGACCGGCCAAATGGCCGGCGATTAc                                                                                                                                                                                                                >  2:818430/1‑139 (MQ=255)
                                                                                            ttgcaaCGCCTGTGCCGCCAAGGCACGAAGCTGTTCCGGGGTCATCTGGTCGAGGTTGGGCGAGAAAGTCATGCCGCCGATTGTGCCAGAGCAGGCGCGGGGCGACGACAAGTAGACCGGCCAAATGGCCGGCGATTAc                                                                                                                                                                                                                <  1:818433/139‑1 (MQ=255)
                                                                                             tgcaaCGCCTGTGCCGCCAAGGCACGAAGCTGTTCCGGGGTCATCTGGTCGAGATTGGGCGAGAAAGTCATGCCGCCGATTGTGCCAGACCAGGCGCGATGCGACGACAAGTAGACCGGCCAAATGGCCGGCGTTTACAGCATGCTAATCACACCGCCTGCACCGACGCGCTGCCACGGCAGGCCCAGCACCAACGCCTGGAGTTGTTCGCTGTCGAGTTCGACCTCGCAGCCg                                                                                                                >  2:382549/1‑234 (MQ=255)
                                                                                                aaCGCCTGTGCCGCCAAGGCACGAAGCTGTTCCGGGGTCATCTGGTCGAGGTTGGGCGAGAAAGTCATGCCGCCGATTGTGCCAGAGCAGGCGCGGGGCGACGACAAGTAGACCGGCCAAATGGCCGGCGTTTACAGCATGCTAATCACACCGCCTGCACCGACGCGCTGCCACGGCAGGCCCAGCACCAACGCCTGGAGTTGTTCGCTGTCGAGt                                                                                                                               <  1:430110/216‑1 (MQ=255)
                                                                                                aaCGCCTGTGCCGCCAAGGCACGAAGCTGTTCCGGGGTCATCTGGTCGAGGTTGGGCGAGAAAGTCATGCCGCCGATTGTGCCAGAGCAGGCGCGGGGCGACGACAAGTAGACCGGCCAAATGGCCGGCGTTTACAGCATGCTAATCACACCGCCTGCACCGACGCGCTGCCACGGCAGGCCCAGCACCAACGCCTGGAGTTGTTCGCTGTCGAGt                                                                                                                               >  2:430107/1‑216 (MQ=255)
                                                                                                   gCCTGTGCCGCCAAGGCACGAAGCTGTTCCGGGGTCATCTGGTCGAGGTTGGGCGAGAAAGTCATGCCGCCGATTGTGCCAGAGCAGGCGCGGGGCGACGACAAGTAGACCGGCCAAATGGCCGGCGTTTACAGCATGCTAATCACACCGCCTGCACCGACGCGCTGCCACg                                                                                                                                                                        <  2:741256/172‑1 (MQ=255)
                                                                                                    ccTGTGCCGCCAAGGCACGAAGCTGTTCCGGGGTCATCTGGGCGAGGTTGGGCGAGAAAGTCATGCCGCCGATTGTGCCAGAGCAGGCGCGGGGCGACGACAAGTAGACCGGCCAAATGGCCGGCGTTTACAg                                                                                                                                                                                                              <  2:125412/133‑1 (MQ=255)
                                                                                                     cggtgCCGCCAAGGCACGAAGCTGTTCCGGGGTCATCTGGCAGAGGTTGGGCGAGAAAGTCATGCCGCCGATGGTGCCAGAGCAGGCGCGGGGCGACGACAAGTAGACCGGCCAAATGGCCGGCGTTTACAGCATGCTAATaac                                                                                                                                                                                                  <  1:277442/142‑1 (MQ=255)
                                                                                                      tgtgCCGCCAAGGCACGAAGCTGTTCCGGGGTCATCTGGTCGAGGTTGGGCGAGAAAGTCATGCCGCCGATTGTGCCAGAGCAGGCGCGGGGCGACGACAAGTAGACCGGCCAAATGGCCGGCGGTTACAGCATGCTAATCACACCGCCTGCACCGACGCGCTGCCACGCCAGGccccg                                                                                                                                                              >  2:179143/1‑177 (MQ=255)
                                                                                                       gtgCCGCCAAGGCACGAAGCTGTTCCGGGGTCATCTGGTCGAGATTGGGCGAGAAAGTCATGCCGCCGATTGTGCCAGACCAGGCGCGATGCGACGACAAGTAGACCGGCCAAATGGCCGGCGTTTACAGCATGCTAATCACACCGCCTGCACCGACGCGCTGCCACGGCAGGCCCAGCACCAACGCCTGGAGTTGTTCGCTGTCGAGTTCGACCTCGCAGCCGTGGCGAATGcc                                                                                                     <  1:382552/235‑1 (MQ=255)
                                                                                                          ccgccAAGGCACGAAGCTGTTCCGGGGTCATCTGGTCGAGGTTGGGCGAGAAAGTCATGCCGCCGATTGTGCCAGAGCAGGCGCGGGGCGACGACAAGTAGACCGGCCAAATGGCCGGCGTTTACAGCATGCTAATCACACCGCCTGCACCGACGCGCTGCCACGGCAGGCCCAGCACCAACGCCTGGAGTTGTTCGCTGTCGAGTTCGACCTCGCAGCCGt                                                                                                               >  1:459653/1‑222 (MQ=255)
                                                                                                                 ggCACGAAGCTGTTCCGGGGTCATCTGGTCGAGGTTGGGCGAGAAAGTCATGCCGCCGATTGTGCCAGAGCAGGCGCGGGGCGACGACAAGTAGACCGGCCAAATGGCCGGCGTTTACAGCATGCTAATCACACCGCCTGCACCg                                                                                                                                                                                     <  2:389670/145‑1 (MQ=255)
                                                                                                                                cGGGGTCATCTGGTCGAGGTTGGGCGAGAAAGTCATGCCGCCGATTGTGCCAGAGCAGGCGCGGGGCGACGACAAGTAGACCGGCCAAATGGCCGGCGTTTACAGCATGCTAATCACACCGCCTGCACCGACGCGCTGCCACGGCAGGCCCAGCACCAACGCCTGGAGTTGTTCGc                                                                                                                                       <  1:552577/176‑1 (MQ=255)
                                                                                                                                cGGGGTCATCTGGTCGAGGTTGGGCGAGAAAGTCATGCCGCCGATTGTGCCAGAGCAGGCGCGGGGCGACGACAAGTAGACCGGCCAAATGGCCGGCGTTTACAGCATGCTAATCACACCGCCTGCACCGACGCGCTGCCACGGCAGGCCCAGCACCAACGCCTGGAGTTGTTCGCTGTCGAGTTCGACCTCGCAGCCGTGGCGAATGCCGGGCCAGTGGAACTTGCCTTGGTTCa                                                                           >  1:740135/1‑236 (MQ=255)
                                                                                                                                  gggTCATCTGGTCGGGGTTGGGCGAGAAAGTCATGCCGCCGATTGTGCCAGAGCAGGCGCGGGGCGACGACAAGTAGACCGGCCAAATGGCCGGCGTTTACAGCATGCt                                                                                                                                                                                                        <  2:30382/109‑1 (MQ=255)
                                                                                                                                          tGGTCGAGGTTGGGCGAGAAAGTCATGCCGCCGATTGTGCCAGAGCAGGCGCGGGGCGACGACAAGTAGACCGGCCAAATGGCCGGCGTTTACAGCATGCTAATCACACCGCCTGCACCGACGCGCTGCCACGGCAGGCCCAGCACCAACGCCTGGAGTTGTTCGCTGTCGAGTTCGACCTCGCAGCCGTGGCGAATGCCGGGc                                                                                                 >  2:654319/1‑204 (MQ=255)
                                                                                                                                          tGGTCGAGGTTGGGCGAGAAAGTCATGCCGCCGATTGTGCCAGAGCAGGCGCGGGGCGACGACAAGTAGACCGGCCAAATGGCCGGCGTTTACAGCATGCTAATCACACCGCCTGCACCGACGCGCTGCCACGGCAGGCCCAGCACCAACGCCTGGAGTTGTTCGCTGTCGAGTTCGACCTCGCAGCCGTGGCGAATGCCGGGc                                                                                                 <  1:654322/204‑1 (MQ=255)
                                                                                                                                          tGGTCGAGGTTGGGCGAGAAAGTCATGCCGCCGATTGTGCCAGAGCAGGCGCGGGGCGACGACAAGTAGACCGGCCAAATGGCCGGCGTTTACAGCATGCTAATCACACCGCCTGCACCGACGCGCTGCCACGGCAGGCCCAGCACCAACGCCTGGAGTTGTTCGCTGTCGAGTTCGACCTCGCAGCCGTGGCGAATGCCGGGCCAGTGGAACTTGCCTTg                                                                                >  2:756772/1‑221 (MQ=255)
                                                                                                                                          tGGTCGAGGTTGGGCGAGAAAGTCATGCCGCCGATTGTGCCAGAGCAGGCGCGGGGCGACGACAAGTAGACCGGCCAAATGGCCGGCGTTTACAGCATGCTAATCACACCGCCTGCACCGACGCGCTGCCACGGCAGGCCCAGCACCAACGCCTGGAGTTGTTCGCTGTCGAGTTCGACCTCGCAGCCGTGGCGAATGCCGGGCCAGTGGAACTTGCCTTGGTTCAGTCGa                                                                      >  2:327344/1‑231 (MQ=255)
                                                                                                                                          tGGTCGAGGTTGGGCGAGAAAGTCATGCCGCCGATTGTGCCAGAGCAGGCGCGGGGCGACGACAAGTAGACCGGCCAAATGGCCGGCGTTTACAGCATGCTAATCACACCGCCTGCACCGACGCGCTGCCACGGCAGGCCCAGCACCAACGCCTGGAGTTGTTCGCTGTCGAGTTCGACCTCGCAGCCGTGGCGAATGCCGGGCCAGTGGAACTTGCCTTGGTTCAGTCGAcg                                                                    >  2:792757/1‑233 (MQ=255)
                                                                                                                                                     gggCGAGAAAGTCATGCCGCCGATTGTGCCAGAGCAGGCGCGGGGCGACGACAAGTAGACCGGCCAAATGGCCGGCGTTTACAGCATGCTAATCACACCGCCTGCACCGACGCGCTGCCACGGCAGGCCCAGCACCa                                                                                                                                                         >  2:108530/1‑137 (MQ=255)
                                                                                                                                                          agaAAGTCATGCCGCCGATTGTGCCAGACCAGGCGCGATGCGACGACAAGTAGACCGGCCAAATGGCCGGCGTTTACAGCATGCTAATCACACCGCCTGCACCGACGCGCTGCCACGGCAGGCCCAGCGCCAACGCCTGGAGTTGTTCGCTGTCGAGTTCGACCTc                                                                                                                       <  1:559341/166‑1 (MQ=255)
                                                                                                                                                          agaAAGTCATGCCGCCGATTGTGCCAGACCAGGCGCGATGCGACGACAAGTAGACCGGCCAAATGGCCGGCGTTTACAGCATGCTAATCACACCGCCTGCACCGACGCGCTGCCACGGCAGGCCCAGCACCAACGCCTGGAGTTGTTCGCTGTCGAGTTCGACCTc                                                                                                                       >  2:559338/1‑166 (MQ=255)
                                                                                                                                                             aaGTCATGCCGCCGATTGTGCCAGAGCAGGCGCGGGGCGACGACAAGTAGACCGGCCAAATGGCCGGCGTTTACAGCATGCTAATCACACCGCCTGCACCGACGCGCTGCCACGGCAGGCCCAGCACCAACGCCTGGAGTTGTTCGCTGTCGAGTTCGACCTCGCAGCCGTGGCGAATGCCGGGCCAGTGGAACTTGCCTTGGTTCAGTCGACGCGCTGCCAGCCAGACACCGACg                                              <  1:786641/236‑1 (MQ=255)
                                                                                                                                                                tCATGCCGCCGATTGTGCCAGAGCAGGCGCGGGGCGACGACAAGTAGACCGGCCAAATGGCCGGCGTTTACAGCATGCTAATCACACCGCCTGCACCGACGCGCTGCCACGGCAGGCCCAGCACCAACGCCTGGAGTTGTTCGCTGTCGAGTTCGACCTCGCAGCCGTGGCGAATGCCGGGCCAGTGGAACTTGCCTTGGTTCAGTCGACGCGCTGCCAGCCAGACACCGACg                                              <  2:365808/233‑1 (MQ=255)
                                                                                                                                                                        ccgATTGTGCCAGAGCAGGCGCGGGGCGACGACAAGTAGACCGGCCAAATGGCCGGCGTTTACAGCATGCTAATCACACCGCCTGCACCGACGCGCTGCCACGGCAGGCCCAGCACCAACGCCTGGAGTTGTTCGCTGTCGAGTTCGACCTCGCAGCCGTGGCGAATGCCGGGCCAGTGGAACTTGCCTTGGTTCAGTCGACGCGCTGCCAGCCAGACACCGACGCCATCGTGCAc                                   <  1:792760/236‑1 (MQ=255)
                                                                                                                                                                           aTTGTGCCAGAGCAGGCGCGGGGCGACGACAAGTAGACCGGCCAAATGGCCGGCGTTTACAGCATGCTAATCACACCGCCTGCACCGACGCGCTGCCACGGCAGGCCCAGCa                                                                                                                                                            >  2:381065/1‑112 (MQ=255)
                                                                                                                                                                           aTTGTGCCAGAGCAGGCGCGGGGCGACGACAAGTAGACCGGCCAAATGGCCGGCGTTTACAGCATGCTAATCACACCGCCTGCACCGACGCGCTGCCACGGCAGGCCCAGCa                                                                                                                                                            <  1:381068/112‑1 (MQ=255)
                                                                                                                                                                                 ccAGAGCAGGCGCGGGGCGACGACAAGTAGACCGGCCAAATGGCCGGCGTTTACAGCATGCTAATCACACCGCCTGCACCGACGCGCTGCCACGGCAGGCCCAGCACCAACGCCTGGAGTTGTTCGCTGTCGAGTTCGACCTCGCAGCCGTGGCGAATGCCGGGCCAGTGGAACTTGCCTTGGTTCAGTCGACGCGCTGCCAGCCAGACACCGACGCCATCGTGCACCAGc                               <  2:548330/231‑1 (MQ=255)
                                                                                                                                                                                                          aGTAGACCGGCCAAATGGCCGGCGTTTACAGCATGCTAATCACACCGCCTGCACCGACGCGCTGCCACGGCAGGCCCAGCACCAACGCCTGGAGTTGTTCGCTGTCGAGTTCGACCTCGCAGCCGTGGCGAATGCCGGGCCAGTGGAACTTGCCTTGGTTCAGTCGACGCGCTGCCAGCCAGACACCGACGCCATCGTGCACCAGCACTTTCATCCGGTTAGCCCGGCGATTGGc  <  1:756775/235‑1 (MQ=255)
                                                                                                                                                                                                           gTAGACCGGCCAAATGGCCGGCGTTTACAGCATGCTAATCACACCGCCTGCACCGACGCGCTGCCACGGCAGGCCCAGCACCAACGCCTGGAGTTGTTCGCTGTCGAGTTCGACCTCGCAGCCGTGGCGAATGCCGGGCCAGTGGAACTTGCCTTGGTTCAGTCGACGCGCTGCCAGCCAGACACCg                                                 <  2:924358/187‑1 (MQ=255)
                                                                                                                                                                                                                    ccAAATGGCCGGCGTTTACAGCATGCTAATCACACCGCCTGCACCGACGCGCGGCCACGGCAGGCCCAGCACCAACGCACGGAgttgtt                                                                                                                                          >  2:176565/1‑89 (MQ=255)
                                                                                                                                                                                                                      aaaTGGCCGGCGTTTACAGCATGCTAATCACACCGCCTGCACCGACGCGCTGCCACGGCAGGCCCAGCACCAACGCCTGGAGTTGTTCGCTGTCGAGTTCGACCTCGCAGCCGTGGCGAATGCCGGGCCAGTGGAACTTGCCTTGGTTCAGTCGACGCGCTGCCAGCCAGACACCGACGCCa                                           <  1:925284/182‑1 (MQ=255)
                                                                                                                                                                                                                      aaaTGGCCGGCGTTTACAGCATGCTAATCACACCGCCTGCACCGACGCGCTGCCACGGCAGGCCCAGCACCAACGCCTGGAGTTGTTCGCTGTCGAGTTCGACCTCGCAGCCGTGGCGAATGCCGGGCCAGTGGAACTTGCCTTGGTTCAGTCGACGCGCTGCCAGCCAGACACCGACGCCa                                           >  2:925281/1‑182 (MQ=255)
                                                                                                                                                                                                                                cGTTTACAGCATGCTAATCACACCGCCTGCACCGACGCGCTGCCACGGCAGGCCCAGCACCAACGCCTGGAGTTGTTCGCTGTCGAGTTCGACCTCGCAGCCGTGGCGAATGCCGGGCCAGTGGAACTTGCCTTGGTTCAGTCGACGCGCTGCCAGCCAGACACCGa                                                >  1:131057/1‑167 (MQ=255)
                                                                                                                                                                                                                                  |                                                                                                                                                                                                                  
CGTTTGAGCAGGGCGATTTCGAACTTGAACTGTTCGATGAGGATTTCATCGTTCTGAATCTTCCTGCTCATCGCCTCGACCTGGGATTGCAATTGCAACGCCTGTGCCGCCAAGGCACGAAGCTGTTCCGGGGTCATCTGGTCGAGGTTGGGCGAGAAAGTCATGCCGCCGATTGTGCCAGAGCAGGCGCGGGGCGACGACAAGTAGACCGGCCAAATGGCCGGCGTTTACAGCATGCTAATCACACCGCCTGCACCGACGCGCTGCCACGGCAGGCCCAGCACCAACGCCTGGAGTTGTTCGCTGTCGAGTTCGACCTCGCAGCCGTGGCGAATGCCGGGCCAGTGGAACTTGCCTTGGTTCAGTCGACGCGCTGCCAGCCAGACACCGACGCCATCGTGCACCAGCACTTTCATCCGGTTAGCCCGGCGATTGGC  >  AE015451/4475372‑4475808

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