Predicted mutation
evidence seq id position mutation freq annotation gene description
RA AE015451 4,475,626 A→G 7.4% G103G (GGT→GGC PP_3965 ← transposase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*AE0154514,475,6260AG 7.4% 191.3 / 4.6 54G103G (GGT→GGCPP_3965transposase
Reads supporting (aligned to +/- strand):  ref base A (22/28);  new base G (2/2);  total (24/30)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.38e-01

CTGTTCGATGAGGATTTCATCGTTCTGAATCTTCCTGCTCATCGCCTCGACCTGGGATTGCAATTGCAACGCCTGTGCCGCCAAGGCACGAAGCTGTTCCGGGGTCATCTGGTCGAGGTTGGGCGAGAAAGTCATGCCGCCGATTGTGCCAGAGCAGGCGCGGGGCGACGACAAGTAGACCGGCCAAATGGCCGGCGTTTACAGCATGCTAATCACACCGCCTGCACCGACGCGCTGCCACGGCAGGCCCAGCACCAACGCCTGGAGTTGTTCGCTGTCGAGTTCGACCTCGCAGCCGTGGCGAATGCCGGGCCAGTGGAACTTGCCTTGGTTCAGTCGACGCGCTGCCAGCCAGACACCGACGCCATCGTGCACCAGCACTTTCATCCGGTTAGCCCGGCGATTGGCGAACAGATAAGCGCAGTGCGGCTTCGCCGCACCGAACACCGCAATCACCCGG  >  AE015451/4475401‑4475860
                                                                                                                                                                                                                                 |                                                                                                                                                                                                                                          
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)
                            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)
                                                       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)
                                                                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)
                                                                         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)
                                                                                                             tGGTCGAGGTTGGGCGAGAAAGTCATGCCGCCGATTGTGCCAGAGCAGGCGCGGGGCGACGACAAGTAGACCGGCCAAATGGCCGGCGTTTACAGCATGCTAATCACACCGCCTGCACCGACGCGCTGCCACGGCAGGCCCAGCACCAACGCCTGGAGTTGTTCGCTGTCGAGTTCGACCTCGCAGCCGTGGCGAATGCCGGGc                                                                                                                                                     <  1:654322/204‑1 (MQ=255)
                                                                                                             tGGTCGAGGTTGGGCGAGAAAGTCATGCCGCCGATTGTGCCAGAGCAGGCGCGGGGCGACGACAAGTAGACCGGCCAAATGGCCGGCGTTTACAGCATGCTAATCACACCGCCTGCACCGACGCGCTGCCACGGCAGGCCCAGCACCAACGCCTGGAGTTGTTCGCTGTCGAGTTCGACCTCGCAGCCGTGGCGAATGCCGGGc                                                                                                                                                     >  2:654319/1‑204 (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                                                                                               >  2:925281/1‑182 (MQ=255)
                                                                                                                                                                                         aaaTGGCCGGCGTTTACAGCATGCTAATCACACCGCCTGCACCGACGCGCTGCCACGGCAGGCCCAGCACCAACGCCTGGAGTTGTTCGCTGTCGAGTTCGACCTCGCAGCCGTGGCGAATGCCGGGCCAGTGGAACTTGCCTTGGTTCAGTCGACGCGCTGCCAGCCAGACACCGACGCCa                                                                                               <  1:925284/182‑1 (MQ=255)
                                                                                                                                                                                                   cGTTTACAGCATGCTAATCACACCGCCTGCACCGACGCGCTGCCACGGCAGGCCCAGCACCAACGCCTGGAGTTGTTCGCTGTCGAGTTCGACCTCGCAGCCGTGGCGAATGCCGGGCCAGTGGAACTTGCCTTGGTTCAGTCGACGCGCTGCCAGCCAGACACCGa                                                                                                    >  1:131057/1‑167 (MQ=255)
                                                                                                                                                                                                                cTAATCACACCGCCTGCACCGACGCGCTGCCACGGCAGGCCCAGCACCAACGCCTGGAGTTGTTCGCTGTCGAGTTCGACCTCGCCGCCGTGGCGAATGCCGGGCCAGTGGAACTTGCCTTGGTTCAGTCgg                                                                                                                          >  2:216947/1‑131 (MQ=255)
                                                                                                                                                                                                                cTAATCACACCGCCTGCACCGACGCGCTGCCACGGCAGGCCCAGCACCAACGCCTGGAGTTGTTCGCTGTCGAGTTCGACCTCGCAGCCGTGGCGAATGCCGGGCCAGTGGAACTTGCCGTGGtt                                                                                                                                 >  2:299177/1‑125 (MQ=255)
                                                                                                                                                                                                                      ccacCGCCTGCACCGACGCGCTGCCACGGCAGGCCCAGCACCAACGCCTGGAGTTGTTCGCTGTCGAGTTCGACCTCGCAGCCGTGGCGAATGCCGGGCCAGTGGAACTTGCCTTGGTTCAGTCGACGCGCTGCCAGCCAGACACCGACGCCATCGTGCACCAGCACTTTCATCCGGTTAGCCCGGCGATTGGCGAACAGATAAGCGCAGTGCGGCt                               <  1:174267/216‑1 (MQ=255)
                                                                                                                                                                                                                        cccGCCTGCACCGCCGCGCTGCCACGGCAGGCCCAGCCCCAGCGCCTGGAGTTGTTCGCTGTCGAGTTCGACCTCGCAGCCGTGGCGAATGCCGGGCCAGTGGAACTTGCCTTGGTTCAGTCGACGCGCTGCCAGCCAGACACCGACGCCAt                                                                                              <  2:262822/151‑1 (MQ=255)
                                                                                                                                                                                                                           gCCGGCACCGACGCGCTGCCACGGCAGGCCCAGCCCCAACGCCCGGAgt                                                                                                                                                                                                  >  2:278092/1‑49 (MQ=255)
                                                                                                                                                                                                                             cTGCACCGACGCGCTGCCACGGCAGGCCCAGCACCAACGCCTGGAGTTGTTCGCTGTCGAGTTCGACCTCGCAGCCGTGGCGAATGCCGGGCCAGTGGACCTTGCCTTGGTTCAGTCGACGCGCTGCCAGCCCGACACCGACGCCATCGTGCACCACCACTCTCCTCCgtg                                                                      >  2:56776/1‑169 (MQ=255)
                                                                                                                                                                                                                              tGCACCGACGCGCTGCCACGGCAGGCCCAGCACCAACGCCTGGAGTTGTTCGCTGTCGAGTTCGACCTCGCAGCCGTGGCGAATGCCGGGCCAGTGGAACTTGCCTTGGTTCAGTCGACGCGCTGCCAGCCAGACACCGACGCCATCGTGCACCAGCACTTTCATCCGGTTAGCCCGGCGATTGGCGAACAGATAAGCGCAGTGCGGCTTCGCCGCACCGAACACCGCAATCAccc    >  1:314643/1‑236 (MQ=255)
                                                                                                                                                                                                                                cACCGACGCGCTGCCACGGCAGGCCCAGCACCAACGCCTGGAGTTGTTCGCTGTCTAGTTCGACCTCGCAGCCGTGGCGAATGCCGGGCCAGTGGAACTTGCCTTGGTTCAGTCGACGCGCTGCCAGCCAGACACCGACGCCATCGTGCAc                                                                                       >  2:309312/1‑151 (MQ=255)
                                                                                                                                                                                                                                cACCGACGCGCTGCCACGGCAGGCCCAGCACCAACGCCTGGAGTTGTTCGCTGTCGAGTTCGACCTCGCAGCCGTGGCGAATGCCGGGCCAGTGGAACTTGCCTTGGTTCAGTCGACGCGCTGCCAGCCAGACACCGACGCCATCGTGCACCAGCACTTTCATCCGGTTAGCCCGGCGATTGGCGAACAGATAAGCGCAGTGCGGCTTCGCCGCACCGAACACCGCAATCACCCgg  >  1:354046/1‑236 (MQ=255)
                                                                                                                                                                                                                                 |                                                                                                                                                                                                                                          
CTGTTCGATGAGGATTTCATCGTTCTGAATCTTCCTGCTCATCGCCTCGACCTGGGATTGCAATTGCAACGCCTGTGCCGCCAAGGCACGAAGCTGTTCCGGGGTCATCTGGTCGAGGTTGGGCGAGAAAGTCATGCCGCCGATTGTGCCAGAGCAGGCGCGGGGCGACGACAAGTAGACCGGCCAAATGGCCGGCGTTTACAGCATGCTAATCACACCGCCTGCACCGACGCGCTGCCACGGCAGGCCCAGCACCAACGCCTGGAGTTGTTCGCTGTCGAGTTCGACCTCGCAGCCGTGGCGAATGCCGGGCCAGTGGAACTTGCCTTGGTTCAGTCGACGCGCTGCCAGCCAGACACCGACGCCATCGTGCACCAGCACTTTCATCCGGTTAGCCCGGCGATTGGCGAACAGATAAGCGCAGTGCGGCTTCGCCGCACCGAACACCGCAATCACCCGG  >  AE015451/4475401‑4475860

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: