Predicted mutation
evidence seq id position mutation freq annotation gene description
RA AE015451 4,475,554 G→C 6.1% intergenic (‑19/+45) PP_3964 ← / ← PP_3965 transposase/transposase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*AE0154514,475,5540GC 6.1% 225.6 / 6.4 66intergenic (‑19/+45)PP_3964/PP_3965transposase/transposase
Reads supporting (aligned to +/- strand):  ref base G (18/44);  new base C (2/2);  total (20/46)
Fisher's exact test for biased strand distribution p-value = 5.79e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.37e-01

CAAGTTGCTTTGCGCCGGGCTGATTTGCTCGCTGCGCTTGGCAAACTTGTGCCGTTTGAGCAGGGCGATTTCGAACTTGAACTGTTCGATGAGGATTTCATCGTTCTGAATCTTCCTGCTCATCGCCTCGACCTGGGATTGCAATTGCAACGCCTGTGCCGCCAAGGCACGAAGCTGTTCCGGGGTCATCTGGTCGAGGTTGGGCGAGAAAGTCATGCCGCCGATTGTGCCAGAGCAGGCGCGGGGCGACGACAAGTAGACCGGCCAAATGGCCGGCGTTTACAGCATGCTAATCACACCGCCTGCACCGACGCGCTGCCACGGCAGGCCCAGCACCAACGCCTGGAGTTGTTCGCTGTCGAGTTCGACCTCGCAGCCGTGGCGAATGCCGGGCCAGTGGAACTTGCCTTGGTTCAGTCGACGCGCTGCCAGCCAGACACCGACGCCATCGTGCACCAGC  >  AE015451/4475320‑4475779
                                                                                                                                                                                                                                          |                                                                                                                                                                                                                                 
cAAGTTGCTTTGCGCCGGGCTGATTTGCTCGCTGCGCTTGGCAAACTTGTGCCGTTTGAGCAGGGCGATTTCGAACTTGAACTGTTCGATGAGGATTTCATCGTTCTGAATCTTCCTGCTCATCGCCTCGACCTGGGATTGCAATTGCAACGCCTGTGCCGCCAAGGCACGAAGCTGTTCCGGGGTCATCTGGTCGAGGTTGGGCGAGAAAGTCATGCCGCCGATTGTGCCagag                                                                                                                                                                                                                                   <  1:880971/235‑1 (MQ=255)
  aGTTGCTTTGCGCCGGGCTGATTTGCTCGCTGCGCTTGGCAAACTTGTGCCGTTTGAGCAGGGCGATTTCGAACTTGAACTGTTCGATGAGGATTTCATCGTTCTGAATCTTCCTGCTCATCGCCTCGACCTGGGATTGCAATTGCAACGCCTGTGCCGCCAAGGCACGAAGCTGTTCCGGGGTCATCTGGTCGAGGTTGGGCGAGAAAGTCATGCCGCCGATTGTGCCAGAGCAg                                                                                                                                                                                                                                <  1:800689/236‑1 (MQ=255)
           gcgcCGGGCTGATTTGCTCGCTGCGCTTGGCAAACTTGTGCCGTTTGAGCAGGGCGATTTCGAACTTGAACTGTTCGATGAGGATTTCATCGTTCTGAATCTTCCTGCTCATCGCCTCGACCTGGGATTGCAATTGCAACGCCTGTGCCGCCAAGGCACGAAGCTGTTCCGGGGTCATCTGGTCGAGGTTGGGCGAGAAAGTCATGCCGCCGATTGTGCCAGAGCAGGCGCGGGGc                                                                                                                                                                                                                       <  1:914093/236‑1 (MQ=255)
               cGGGCTGATTTGCTCGCTGCGCTTGGCAAACTTGTGCCGTTTGAGCAGGGCGATTTCGAACTTGAACTGTTCGATGAGGATTTCATCGTTCTGAATCTTCCTGCTCATCGCCTCGACCTGGGATTGCAATTGCAACGCCTGTGCCGCCAAGGCACGAAGCTGTTCCGGGGTCATCTGGTCGAGGTTGGGCGAGAAAGTCATGCCGCCGATTGTGCCAGAGCAGGCGCGGGGcgacg                                                                                                                                                                                                                   <  1:887985/236‑1 (MQ=255)
                              gCTGCGCTTGGCAAACTTGTGCCGTTTGAGCAGGGCGATTTCGAACTTGAACTGTTCGATGAGGATTTCATCGTTCTGAATCTTCCTGCTCATCGCCTCGACCTGGGATTGCAATTGCAACGCCTGTGCCGCCAAGGCACGAAGCTGTTCCGGGGTCATCTGGTCGAGGTTGGGCGAGAAAGTCATGCCGCCGATTGTGCCAGAGCAGGCGCGGGGCGACGACAAGTAGACCGGc                                                                                                                                                                                                     <  1:916951/235‑1 (MQ=255)
                               cTGCGCTTGGCAAACTTGTGCCGTTTGGGCAGGGCGATTTCGAACTTGAACTGTTCGATGAGGATTTCATCGTTCTGAATCTTCCTGCTCATCGCCTCGACCTGGGATTGCAATTGCAACGCCTGTGCCGCCAAGGCACGAAGCTGTTCCGGGGTCATCTGGTCGAGGTTGGGCGAGAAAGTCATGCCGCCGATTGTGCCAGAGCAGGCGCGGGGc                                                                                                                                                                                                                       <  2:797416/216‑1 (MQ=255)
                                  cgcTTGGCAAACTTGTGCCGTTTGAGCAGGGCGATTTCGAACTTGAACTGTTCGATGAGGATTTCATCGTTCTGGATCTTCCTGCTCATCGCCTCGACCTGGGATTGCAATTGCAACGCCTGTGCCGCCAAGGCACGAAGCTGTTCCGGGGTCATCTGGTCGAGGTTGGGCGAGAAAGTCATGCCGCCGATTGTGCCAGAGCAGGCGCGGGGCGACGACAAGTAGACCGGCCaaa                                                                                                                                                                                                 <  1:773704/235‑1 (MQ=255)
                                    cTTGGCAAACTTGTGCCGTTTGAGCAGGGCGATTTCGAACTTGAACTGTTAGATGAGGATTTCATCGTTCTGAATCTTCCTGCTCATCGCCTCGACCTGGGATTGCAATTGCAACGCCTGTGCCGCCAAGGCACGAAGCTGTTCCGGGGTCATCTGGTCGAGGTTGGGCGAGAAAGTCATGCCGCCGATTGTGCCAGAGCAGGCGCGGGGCGACGACAAGTAGACCGGCCAAATg                                                                                                                                                                                               <  1:739432/235‑1 (MQ=255)
                                      tagCAAACTTGTGCCGTTTGAGCAGGGCGATTTCGAACTTGAACTGTTCGATGAGGATTTCATCGTTCTGAATCTTCCTGCTCATCGCCTCGACCTGGGATTGCAATTGCAACGCCTGTGCCGCCAAGGCACGAAGCTGTTCCGGGGTCATCTGGTCGAGGTTGGGCGAGAAAGCCATGCCGCCGATTGTGCCAGAGCAGGCGCGGGGCGACGACAAGTAGACCGGCCaaa                                                                                                                                                                                                 <  1:912137/229‑1 (MQ=255)
                                                    cGTTTGAGCAGGGCGATTTCGAACTTGAACTGTTCGATGAGGATTTCATCGTTCTGAATCTTCCTGCTCATCGCCTCGACCTGGGATTGCAATTGCAACGCCTGTGCCGCCAAGGCACGAAGCTGTTCCGGGGTCATCTGGTCGAGGTTGGGCGAGAAAGTCATGCCGCCGATTGTGCCAGAGCAGGCGCGGGGCGACGACAAGTAGACCGGCCAAATGGCCGGCGtt                                                                                                                                                                                      >  1:262825/1‑228 (MQ=255)
                                                    cGTTTGAGCAGGGCGATTTCGAACTTGAACTGTTCGATGAGGATTTCATCGTTCTGAATCTTCCTGCTCATCGCCTCGACCTGGGATTGCAATTGCAACGCCTGTGCCGCCAAGGCACGAAGCTGTTCCGGGGTCATCTGGTCGAGGTTGGGCGAGAAAGTCATGCCGCCGATTGTGCCAGAGCAGGCGCGGGGCGACGACAAGTAGACCGGCCAAATGGCCGGCGTTTACAGCAt                                                                                                                                                                              >  1:21414/1‑236 (MQ=255)
                                                       ttGAGCAGGGCGATTTCGAACTTGAACTGTTCGATGAGGATTTCATCGTTCTGAATCTTCCTGCTCATCGCCTCGACCTGGGATTGCAATTGCAACGCCTGTGCCGCCAAGGCACGAAGCTGTTCCGGGGTCATCTGGTCGAGGTTGGGCGAGAAAGTCATGCCGCCGATTGTGCCAGAGCAGGCGCGGGGCGACGACAAGTAGACCGGCCAAAt                                                                                                                                                                                                <  2:731512/215‑1 (MQ=255)
                                                        tGAGCAGGGCGATTTCGAACTTGAACTGTTCGATGAGGATTTCATCGTTCTGAATCTTCCTGCGCATCGCCTCGACCTGGGATTGCAATTGCAACGCCTGTGCCGCCAAGGCACGAAGCTGTTCCGGGGTCATCTGGTCGAGGTTGGGCGAGAAAGTCATGCCGCCGATTGTGCCAGAGCAGGCGCGGGGc                                                                                                                                                                                                                       <  2:215692/191‑1 (MQ=255)
                                                         gAGCAGGGCGATTTCGAACTTGAACTGTTCGATGAGGATTTCATCGTTCTGAATCTTCCTGCTCATCGCCTCGACCTGGGATTGCAATTGCAACGCCTGTGCCGCCAAGGCACGAAGCTGTTCCGGGGTCATCTGGTCGCGGTTGGGCGAGAAAGTCATGCCGCCGATTGTGCCAGAGCAGGCGCGGGGCGACGACacg                                                                                                                                                                                                              >  2:206888/1‑197 (MQ=255)
                                                                gCGATTTCGAACTTGAACTGTTCGATGAGGATTTCATCGTTCTGAATCTTCCTGCTCATCGCCTCGACCTGGGATTGCAATTGCAACGCCTGTGCCGCCAAGGCACGAAGCTGTTCCGGGGTCATCTGGTCGAGGTTGGGCGAGAAAGTCATGCCGCCGATTGTGCCAGAGCAGGCGCGGGGCGACGACAAGTAGACCGGCCAAATg                                                                                                                                                                                               <  1:307783/207‑1 (MQ=255)
                                                                                 cTGTTCGATGAGGATTTCATCGTTCTGAATCTTCCTGCTCATCGCCTCGACCTGGGATTGCAATTGCAACGCCTGTGCCGCCAAGGCACGAAGCTGTTCCGGGGTCATCTGGTCGAGGTTGGGCGAGAAAGTCATGCCGCCGATTGTGCCAGAGCAGGCGCGGGGCGACGACAAGTAGACCGGCCAAATGGCCGGCGTTTACAGCATGCTAATCACACCGCCTGCACCGa                                                                                                                                                       <  2:587541/230‑1 (MQ=255)
                                                                                  tGTTCGATGAGGATTTCATCGTTCTGAATCTTCCTGCTCATCGCCTCGACCTGGGATTGCAATTGCAACGCCTGTGCCGCCAAGGCACGAAGCTGTTCCGGGGTCATCTGGTCGAGGTTGGGCGAGAAAGTCATGCCGCCGATTGTGCCAGAGCAGGCGCGGGGCGACGACAAGTAGACCGGCCAAATGGCCGGCGTTTACAGCATGCTAATCACACCGCCTGCACCGACGCGCt                                                                                                                                                 <  1:925664/235‑1 (MQ=255)
                                                                                       gATGAGGATTTCATCGTTCTGAATCTTCCTGCTCATCGCCTCGACCTGGGATTGCAATTGCAACGCCTGTGCCGCCAAGGCACGAAGCTGTTCCGGGGTCATCTGGTCGAGGTTGGGCGAGAAAGTCATGCCGCCGATTGTGCCAGAGCAGGCGCGGGGCGACGACAAGTAGACCGGCCAAATGGcc                                                                                                                                                                                            >  1:14985/1‑187 (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)
                                                                                                                    tGCTCATCGCCTCGACCTGGGATTGCAATTGCAACGCCTGTGCCGCCAAGGCACGAAGCTGTTCCGGGGTCATCTGGTCGAGGTTGGGCGAGAAAGTCATGCCGCCGATTGTGCCAGAGCAGGCGCGGGGcgacg                                                                                                                                                                                                                   >  2:501660/1‑135 (MQ=255)
                                                                                                                    tGCTCATCGCCTCGACCTGGGATTGCAATTGCAACGCCTGTGCCGCCAAGGCACGAAGCTGTTCCGGGGTCATCTGGTCGAGGTTGGGCGAGAAAGTCATGCCGCCGATTGTGCCAGAGCAGGCGCGGGGCGACGACAAGTAGACCGGCCAAATGGCCGGCGTTTACAGCATGCTAATCACACCGCCTGCACCGACGCGCTGCCACGGCAGGCCCAGCACCAACGCCTGGAgttgt                                                                                                              <  1:504678/236‑1 (MQ=255)
                                                                                                                          tCGCCTCGACCTGGGATTGCAATTGCAACGCCTGTGCCGCCAAGGCACGAAGCTGTTCCGGGGTCATCTGGTCGAGGTTGGGCGAGAAAGTCATGCCGCCGATTGTGCCAGAGCAGGCGCGGGGCGACGACAAGTAGACCGGCCAAATGGCCGGCGTTTACAGCATGCTAATCACACCGCCTGCACCGACGCGCTGCCACGGCAGGCCCAGCACCAACGCCTGGAGTTGTTCGc                                                                                                          <  1:758970/234‑1 (MQ=255)
                                                                                                                          tCGCCTCGACCTGGGATTGCAATTGCAACGCCTCTGCCGCCAAGGCACGAAGCTGTGCCGGGGTCATCTGGTCGAGGTTGGGCGAGAAAGTCATGCCGCCGATTGTGCCCGCGCAGGCGCGGGGc                                                                                                                                                                                                                       <  2:626313/125‑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)
                                                                                                                                                      cGCCTGTGCCGCCAAGGCACGAAGCTGTTCCGGGGTCATCTGGTCGAGGTTGGGCGAGAAAGTCATGCCGCCGATTGTGCCAGAGCAGGCGCGGGGc                                                                                                                                                                                                                       <  2:492211/97‑1 (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                                                                    <  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                                                                                                                               <  1:381068/112‑1 (MQ=255)
                                                                                                                                                                                                                               aTTGTGCCAGAGCAGGCGCGGGGCGACGACAAGTAGACCGGCCAAATGGCCGGCGTTTACAGCATGCTAATCACACCGCCTGCACCGACGCGCTGCCACGGCAGGCCCAGCa                                                                                                                               >  2:381065/1‑112 (MQ=255)
                                                                                                                                                                                                                                     ccAGAGCAGGCGCGGGGCGACGACAAGTAGACCGGCCAAATGGCCGGCGTTTACAGCATGCTAATCACACCGCCTGCACCGACGCGCTGCCACGGCAGGCCCAGCACCAACGCCTGGAGTTGTTCGCTGTCGAGTTCGACCTCGCAGCCGTGGCGAATGCCGGGCCAGTGGAACTTGCCTTGGTTCAGTCGACGCGCTGCCAGCCAGACACCGACGCCATCGTGCACCAGc  <  2:548330/231‑1 (MQ=255)
                                                                                                                                                                                                                                          |                                                                                                                                                                                                                                 
CAAGTTGCTTTGCGCCGGGCTGATTTGCTCGCTGCGCTTGGCAAACTTGTGCCGTTTGAGCAGGGCGATTTCGAACTTGAACTGTTCGATGAGGATTTCATCGTTCTGAATCTTCCTGCTCATCGCCTCGACCTGGGATTGCAATTGCAACGCCTGTGCCGCCAAGGCACGAAGCTGTTCCGGGGTCATCTGGTCGAGGTTGGGCGAGAAAGTCATGCCGCCGATTGTGCCAGAGCAGGCGCGGGGCGACGACAAGTAGACCGGCCAAATGGCCGGCGTTTACAGCATGCTAATCACACCGCCTGCACCGACGCGCTGCCACGGCAGGCCCAGCACCAACGCCTGGAGTTGTTCGCTGTCGAGTTCGACCTCGCAGCCGTGGCGAATGCCGGGCCAGTGGAACTTGCCTTGGTTCAGTCGACGCGCTGCCAGCCAGACACCGACGCCATCGTGCACCAGC  >  AE015451/4475320‑4475779

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