Predicted mutation
evidence seq id position mutation freq annotation gene description
RA AE015451 4,475,632 G→T 7.8% R101R (CGC→CGA)  PP_3965 ← transposase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*AE0154514,475,6320GT 7.8% 182.2 / 7.0 51R101R (CGC→CGA) PP_3965transposase
Reads supporting (aligned to +/- strand):  ref base G (22/25);  new base T (2/2);  total (24/27)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.08e-01

TGTTCGATGAGGATTTCATCGTTCTGAATCTTCCTGCTCATCGCCTCGACCTGGGATTGCAATTGCAACGCCTGTGCCGCCAAGGCACGAAGCTGTTCCGGGGTCATCTGGTCGAGGTTGGGCGAGAAAGTCATGCCGCCGATTGTGCCAGAGCAGGCGCGGGGCGACGACAAGTAGACCGGCCAAATGGCCGGCGTTTACAGCATGCTAATCACACCGCCTGCACCGACGCGCTGCCACGGCAGGCCCAGCACCAACGCCTGGAGTTGTTCGCTGTCGAGTTCGACCTCGCAGCCGTGGCGAATGCCGGGCCAGTGGAACTTGCCTTGGTTCAGTCGACGCGCTGCCAGCCAGACACCGACGCCATCGTGCACCAGCACTTTCATCCGGTTAGCCCGGCGATTGGCGAACAGATAAGCGCAGTGCGGCTTCGCCGCACCGAACACCGCAATCACCCGG  >  AE015451/4475402‑4475860
                                                                                                                                                                                                                                      |                                                                                                                                                                                                                                    
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)
                                  tGCTCATCGCCTCGACCTGGGATTGCAATTGCAACGCCTGTGCCGCCAAGGCACGAAGCTGTTCCGGGGTCATCTGGTCGAGGTTGGGCGAGAAAGTCATGCCGCCGATTGTGCCAGAGCAGGCGCGGGGCGACGACAAGTAGACCGGCCAAATGGCCGGCGTTTACAGCATGCTAATCACACCGCCTGCACCGACGCGCTGCCACGGCAGGCCCAGCACCAACGCCTGGAgttgt                                                                                                                                                                                               <  1:504678/236‑1 (MQ=255)
                                        tCGCCTCGACCTGGGATTGCAATTGCAACGCCTGTGCCGCCAAGGCACGAAGCTGTTCCGGGGTCATCTGGTCGAGGTTGGGCGAGAAAGTCATGCCGCCGATTGTGCCAGAGCAGGCGCGGGGCGACGACAAGTAGACCGGCCAAATGGCCGGCGTTTACAGCATGCTAATCACACCGCCTGCACCGACGCGCTGCCACGGCAGGCCCAGCACCAACGCCTGGAGTTGTTCGc                                                                                                                                                                                           <  1:758970/234‑1 (MQ=255)
                                          gCCTCGACCTGGGATTGCAATTGCAACGCCTGTGCCGCCAAGGCACGAAGCTGTTCCGGGGTCATCTGGTCGAGGTTGGGCGAGAAAGTCATGCCGCCGATTGTGCCAGAGCAGGCGCGGGGCGACGACAAGTAGACCGGCCAAATGGCCGGCGTTTACAGCATGCTAATCACACCGCCTGCACCGACGCGCTGCCAc                                                                                                                                                                                                                             <  2:556573/198‑1 (MQ=255)
                                          gCCTCGACCTGGGATTGCAATTGCAACGCCTGTGCCGCCAAGGCACGAAGCTGTTCCGGGGTCATCTGGTCGAGGTTGGGCGAGAAAGTCATGCCGCCGATTGTGCCAGAGCAGGCGCGGGGCGACGACAAGTAGACCGGCCAAATGGCCGGCGATTACAGCATGCTGATCACGCCACCCGCGCCGACTCGCt                                                                                                                                                                                                                                  >  1:745364/1‑193 (MQ=11)
                                          gCCTCGACCTGGGATTGCAATTGCAACGCCTGTGCCGCCAAGGCACGAAGCTGTTCCGGGGTCATCTGGTCGAGGTTGGGCGAGAAAGTCATGCCGCCGATTGTGCCAGAGCAGGCGCGGGGCGACGACAAGTAGACCGGCCAAATGGCCGGCGATTACAGCATGCTGATCACGCCACCCGCGCCGACTCGCt                                                                                                                                                                                                                                  <  2:745361/193‑1 (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)
                                                                                                  cGGGGTCATCTGGTCGAGGTTGGGCGAGAAAGTCATGCCGCCGATTGTGCCAGAGCAGGCGCGGGGCGACGACAAGTAGACCGGCCAAATGGCCGGCGTTTACAGCATGCTAATCACACCGCCTGCACCGACGCGCTGCCACGGCAGGCCCAGCACCAACGCCTGGAGTTGTTCGc                                                                                                                                                                                           <  1:552577/176‑1 (MQ=255)
                                                                                                  cGGGGTCATCTGGTCGAGGTTGGGCGAGAAAGTCATGCCGCCGATTGTGCCAGAGCAGGCGCGGGGCGACGACAAGTAGACCGGCCAAATGGCCGGCGTTTACAGCATGCTAATCACACCGCCTGCACCGACGCGCTGCCACGGCAGGCCCAGCACCAACGCCTGGAGTTGTTCGCTGTCGAGTTCGACCTCGCAGCCGTGGCGAATGCCGGGCCAGTGGAACTTGCCTTGGTTCa                                                                                                                               >  1:740135/1‑236 (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                                                                                                                                                                                                                <  1:381068/112‑1 (MQ=255)
                                                                                                                                             aTTGTGCCAGAGCAGGCGCGGGGCGACGACAAGTAGACCGGCCAAATGGCCGGCGTTTACAGCATGCTAATCACACCGCCTGCACCGACGCGCTGCCACGGCAGGCCCAGCa                                                                                                                                                                                                                >  2:381065/1‑112 (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)
                                                                                                                                                                                                                                 ccGACGCGCTGCCACGGCAGGCCCAGCACCAACGCCTGGAGTTGTTCGCTGTCGAGTTCGACCTCGCAGCCGTGGCGAATGCCGGGCCAGTGGAACTTGCCTTGGTTCAGTCGACGCGCTGCCAGCCAGACACCGACGCCATCGTGCACCAGCACTTTCATCCGGTTAGCCCGGCGATTGGCGAACAGATAAGCGCAGTGCGGCTTCGCCGCACCGAACAcc              <  2:740132/222‑1 (MQ=255)
                                                                                                                                                                                                                                      gcgcTGCCACGGCAGGCCCAGCACCAACGCCTGGAGTTGTTCGCTGTCGAGTTCGACCTCGCAGCCGTGGCGAATGCCGGGCCAGTGGAACTTGCCTTGGTTCAGTCGACGCGCTGCCAGCCAGACACCGACGCCATCGTGCACCAGCACTTTCATCCGGTTAGCCCGGCGATTGGCGAACAGATAAGCGCAGTGCGGCTTCGCCGCACCGAACACCGCAATc        <  2:498661/223‑1 (MQ=255)
                                                                                                                                                                                                                                      |                                                                                                                                                                                                                                    
TGTTCGATGAGGATTTCATCGTTCTGAATCTTCCTGCTCATCGCCTCGACCTGGGATTGCAATTGCAACGCCTGTGCCGCCAAGGCACGAAGCTGTTCCGGGGTCATCTGGTCGAGGTTGGGCGAGAAAGTCATGCCGCCGATTGTGCCAGAGCAGGCGCGGGGCGACGACAAGTAGACCGGCCAAATGGCCGGCGTTTACAGCATGCTAATCACACCGCCTGCACCGACGCGCTGCCACGGCAGGCCCAGCACCAACGCCTGGAGTTGTTCGCTGTCGAGTTCGACCTCGCAGCCGTGGCGAATGCCGGGCCAGTGGAACTTGCCTTGGTTCAGTCGACGCGCTGCCAGCCAGACACCGACGCCATCGTGCACCAGCACTTTCATCCGGTTAGCCCGGCGATTGGCGAACAGATAAGCGCAGTGCGGCTTCGCCGCACCGAACACCGCAATCACCCGG  >  AE015451/4475402‑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: ‑