Predicted mutation
evidence seq id position mutation freq annotation gene description
RA AE015451 2,051,065 T→A 9.3% intergenic (‑63/+23) PP_1825 ← / ← PP_1826 conserved exported protein of unknown function/Hydrolase, isochorismatase family

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*AE0154512,051,0650TA 9.3% 210.8 / 12.1 65intergenic (‑63/+23)PP_1825/PP_1826conserved exported protein of unknown function/Hydrolase, isochorismatase family
Reads supporting (aligned to +/- strand):  ref base T (26/32);  new base A (3/3);  total (29/36)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.37e-01

AGTTCCTGCAGGGCAGGGTAGTTGCCTAGCAGCCCGGCGACGCCGGCCATGGCCAGCAATACGCCCAGCGTCGCCAGCAGGGCCGACAGGCCGGCCGCCAGGCGGGCGCCCCAGTTACGCGGTTCGCGCGGGCGCAGGCGCTTGGCGTCGAAACTGCCTTTAAGCTTCATTCCGCTCTCCTCGGTGGATATCCGGAATTCGACCTGCTGCTGCCCATCGGGTTCCGCCCGGCTGCCGGGCGGTCCTGCAAGCCGCTCAGATCAGGCTGGCGGTGGGGGCCACGCAAGCGAAGTTGTCAGCCATCACGGCCATTTCGCACTGGTGAATCTGCGCGGCCGGAATGACGCCGTCCTTGAATGCCAGGTCACGGGTGGCCGACGCGTCTTCTACCAGGGTGCACCGATAACCATAGTCCTTGGCGCGGCGAACGGTGGTGCTGACGCTGGAATGGCTCATGAAACCAC  >  AE015451/2050833‑2051296
                                                                                                                                                                                                                                        |                                                                                                                                                                                                                                       
aGTTCCTGCAGGGCAGGGTAGTTGCCTAGCAGCCCGGCGACGCCGGCCATGGCCAGCAATACGCCCAGCGTCGCCAGCAGGGCCGACAGGCCGGCCGCCAGGCGGGCGCCCCAGTTACGCGGTTCGCGCGGGCGCAGGCGCTTGGCGTCGAAACTGCCTTTAAGCTTCATTCCGCTCTCCTCGGTGGATATCCGGAATTCGACCTGCTGCTGCCCATCGGGTTCCGCCCGGCTGcc                                                                                                                                                                                                                                      <  1:1065743/236‑1 (MQ=255)
 gTTCCTGAAGGGCAGGGTAGTTGCCTAGCAGCCCGGCGACGCCGGCCATGGCCAGCAATACGCCCAGCGTCGCCAGCAGGGCCGACAGGCCGGCCGCCAGGCGGGCGCCCCAGTTACGCGGTTCGCGCGGGCGCAGGCGCTTGGCGTCGAAACTGCCTTTAAGCTTCATTCCGCTCTCCTCGGTGGATATCCGGAATTCGACCTGCTGCTGCCCATCGGGTTCCGCCCGGCTGcc                                                                                                                                                                                                                                      <  1:457160/235‑1 (MQ=255)
            gcagggTAGTTGCCTAGCAGCCCGGCGACGCCGGCCATGGCCAGCAATACGCCCAGCGTCGCCAGCAGGGCCGACAGGCCGGCCGCCAGGCGGGCGCCCCAGTTACGCGGTTCGCGCGGGCGCAGGCGCTTGGCGTCGAAACTGCCTTTAAGCTTCATTCCGCTCTCCTCGGTGGATATCCGGAATTCGACCTGCTGCTGCCCATCGGGTTCCGCCCGGCTGCCGGGCGGTCCTGc                                                                                                                                                                                                                          <  1:860958/236‑1 (MQ=255)
                ggTAGTTGCCTAGCAGCCCGGCGACGCCGGCCATGGCCAGCAATACGCCCAGCGTCGCCAGCAGGGCCGACAGGCCGGCCGCCAGGCGGGCGCCCCAGTTACGCGGTTCGCGCGGGCGCAGGCGCTTGGCGTCGAAACTGCCTTTAAGCTTCATTCCGCTCTCCTCGGTGGATATCCGGAATTCGACCTGCTGCTGCCCATCGGGTTCCGCCCGGCTGCCGGGCGGTCCTGCaa                                                                                                                                                                                                                        >  2:753316/1‑234 (MQ=255)
                        ccTAGCCGCCCGGCGACGCCGGCCATGGCCAGCAATACGCCCAGCGTCGCCAGCAGGGCCGACAGGCCGGCCGCCAGGCGGGCGCCCCAGTTACGCGGTTCGCGCGGGCGCAGGCGCTTGGCGTCGAAACTGCCTTTAAGCTTCATTCCGCTCTCCTCGGTGGATATCCGGAATTCGACCTGCTGCTGCCCATCGGGTTCCGCCCGGCt                                                                                                                                                                                                                                         <  1:511950/209‑1 (MQ=255)
                              agcCCGGCGACGCCGGCCATGGCCAGCAATACGCCCAGCGTCGCCAGCAGGGCCGACAGGCCGGCCGCCAGGCGGGCGCCCCAGTTACGCGGTTCGCGCGGGCGCAGGCGCTTGGCGTCGAAACTGCCTTTAAGCTTCATTCCGCTCTCCCCGGTGGATATCCGGAATTCGACCTGCTGCTGCCCATCGGGTTCCGCCCGGCTGCCGGGCGGTCCTGCAAGCCGCTCAGATCAgg                                                                                                                                                                                                         <  1:1038841/235‑1 (MQ=255)
                                  cGGCGACGCCGGCCATGGCCAGCAATACGCCCAGCGTCGCCAGCAGGGCCGACAGGCCGGCCGCAAGGCGGGCGCCCCAGTTACGCGGTTCGCGCGGGCGCAGGCACTTGGCGTCGAAACTGCCTTTAAGCTTCATTCCGCTCTCCTCGGTGGATATCCGGAATTCGACCTGCTGCTGCCCATCGGGTTCCGCCCGGCTGCCGGGCg                                                                                                                                                                                                                                 >  2:655901/1‑207 (MQ=255)
                                           cGGCCATGGCCAGCAATACGCCCAGCGTCGCCAGCAGGGCCGACAGGCCGGCCGCCAGGCGGGCGCCCCAGTTACGCGGTTCGCGCGGGCGCAGGCGCTTGGCGTCGAAACTGCCTTTAAGCTTCATTCCGCTCTCCTCGGTGGATATCCGGAATTCGACCTGCTGCTGCCCATCGGGTTCCGCCCGGCTGCCGGGCGGTCCTGCAAGCCGCTCAGATCAGGCTGGCGGTggggg                                                                                                                                                                                            <  1:1149026/235‑1 (MQ=255)
                                                   gCCAGCAATACGCCCAGCGTCGCCAGCAGGGCCGACAGGCCGGCCGCCAGGCGGGCGCCCCAGTTACGCGGTTCGCGCGGGCGCAGGCGCTTGGCGTCGAAACTGCCTTTAAGCTTCATTCCGCTCTCCTCGGTGGATATCCGGAATTCGACCTGCTGCTGCCCATCGGGTTCCGCCCGGCTGCCGGGCGGTCCTGCAAGCCGCTCAGATCAGGCTGGCGGTGGGGGCCACGCaa                                                                                                                                                                                    <  1:818468/235‑1 (MQ=255)
                                                              gCCCAGCGTCGCCAGCAGGGCCGACAGGCCGGCCGCCAGGCGGGCGCCCCAGTTACGCGGTTCGCGCGGGCGCAGGCGCTTGGCGTCGAAACTGCCTTTAAGCTTCATTCCGCTCTCCTCGGTGGATATCCGGAATTCGACCTGCTGCTGCCCATCGGGTTCCGCCCGGCTGCCGGGCGGt                                                                                                                                                                                                                               <  2:897305/181‑1 (MQ=255)
                                                               cccAGCGTCGCCAGCAGGGCCGACAGGCCGGCCGCCAGGCGGGCGCCCCAGTTACGCGGTTCGCGCGGGCGCAGGCGCTTGGCGTCGAAACTGCCTTTAAGCTTCATTCCGCTCTCCTCGGTGGATATCCGGAATTCGACCTGCTGCTGCCCATCGGGTTCCGCCCGGCTGCCGGGCgg                                                                                                                                                                                                                                <  1:555142/179‑1 (MQ=255)
                                                                ccAGCGTCGCCAGCAGGGCCGACAGGCCGGCCGCCAGGCGGGCGCCCCAGTTACGCGGTTCGCGCGGGCGCAGGCGCTTGGCGTCGAAACTGCCTTTAAGCTTCATTCCGCTCTCCTCGGTGGATATCCGGAATTCGACCTGCTGCTGCCCATCGGGTTCCGCCCGGCTGCCGGGCGGTCCTGCAAGCCGCTCAGATCAgg                                                                                                                                                                                                         >  1:160959/1‑201 (MQ=255)
                                                                ccAGCGTCGCCAGCAGGGCCGACAGGCCGGCCGCCAGGCGGGCGCCCCAGTTACGCGGTTCGCGCGGGCGCAGGCGCTTGGCGTCGAAACTGCCTTTAAGCTTCATTCCGCTCTCCTCGGTGGATATCCGGAATTCGACCTGCTGCTGCCCATCGGGTTCCGCCCGGCTGCCGGGCGGTCCTGCAAGCCGCTCAGATCAGGCTGGCGGTGGGGGCCACGCAAGCGAAGTTGTc                                                                                                                                                                         >  2:737619/1‑233 (MQ=255)
                                                                     gTCGCCAGCAGGGCCGACAGGCCGGCCGCCGGGCGGGCGCCCCAGTTACGCGGTTCGCGCGGGCGCAGGCGCTTGGCGTCGAAACTGCCTTTAAGCTTCATTCCGCTCTCCTCGGTGGATATCCGGAATTCGACCTGCTGCTGCCCATCGGGTTCCGCCCGGCTGCCGGGCGGTCCTGCAAGCCGCTCAGATCAGGCTGGCGGTGGGGGCCACGCAAGCGAAGTTGTCAGCCAt                                                                                                                                                                   >  1:1014054/1‑234 (MQ=255)
                                                                     gTCGCCAGCAGGGCCGACAGGCCGGCCGCCAGGCGGGCGCCCCAGTTACGCGGTTCGCGCGGGCGCAGGCGCTTGGCGTCGAAACTGCCTTTAAGCTTCATTCCGCTCTCCTCGGTGGATATCCGGAATTCGACCTGCTGCTGCCCATCGGGTTCCGCCCGGCTGCCGGGCGGTCCTGCAAGCCGCTCAGATCAGGCTGGCGGTGGGGGCCACGCAAGCGAAGTTGTCAGCCATCa                                                                                                                                                                 >  1:504126/1‑236 (MQ=255)
                                                                             cagGGCCGCCAGGCCGGCCGCCAGGCGGGCGCCCCAGTTACGCGGTTCGCGCGGGCGCAGGCGCTTGGCGTCGAAACTGCCTTTAAGCTTCATTCCGCTCTCCTCGGTGGATATCCGGAATTCGACCTGCTGCTGCCCATCGGGTTCCGCCCGGCTGCCGGGCGGTCCTGCAAGCCGCTCAGATCAGGCTGGc                                                                                                                                                                                                    <  2:1042337/193‑1 (MQ=255)
                                                                                       aGGCCGGCCGCCAGGCGGGCGCCCCAGTTACGCGGTTCGCGCGGGCGCAGGCGCTTGGCGTCGAAACTGCCTTTAAGCTTCATTCCGCTCTCCTCGGTGGATATCCGGAATTCGACCTGCTGCTGCCCATCGGGTTCCGCCCGGCTGCCGGGCGGTCCTGCAAGCCGCTCAGATCAGGCTGGCGGTGGGGGCCACGCAAGCGAAGTTGTCAGCCATCACGGcc                                                                                                                                                            <  1:955700/223‑1 (MQ=255)
                                                                                       aGGCCGGCCGCCAGGCGGGCGCCCCAGTTACGCGGTTCGCGCGGGCGCAGGCGCTTGGCGTCGAAACTGCCTTTAAGCTTCATTCCGCTCTCCTCGGTGGATATCCGGAATTCGACCTGCTGCTGCCCATCGGGTTCCGCCCGGCTGCCGGGCGGTCCTGCAAGCCGCTCAGATCAGGCTGGCGGTGGGGGCCACGCAAGCGAAGTTGTCAGCCATCACGGcc                                                                                                                                                            >  2:955697/1‑223 (MQ=255)
                                                                                         gccggccgCCAGGCGGGCGCCCCAGTTACGCGGTTCGCGCGGGCGCAGGCGCTTGGCGTCGAAACTGCCTTTAAGCTTCATTCCGCTCTCCTCGGTGGATATCCGGAATTCGACCTGCTGCTGCCCATCGGGTTCCGCCCGGCTGCCGGGCGGTCCTGCAAGCCGCTCAGATCAGGCTGGCGGTGGGGGCCACGCAAGCGAAGTTGTCAGCCATCACGGcc                                                                                                                                                            <  1:944309/221‑1 (MQ=255)
                                                                                           cggccgCCAGGCGGGCGCCCCAGTTACGCGGTTCGCGCGGGCGCAGGCGCTTGGCGTCGAAACTGCCTTTAAGCTTCATTCCGCTCTCCTCGGTGGATATCCGGAATTCGACCTGCTGCTGCCCATCGGGTTCCACCCGGCTGCCGGGCGGTCCTGCAAGCCGCTCAGATCAGGCTGGCGGt                                                                                                                                                                                                 >  2:441515/1‑182 (MQ=255)
                                                                                                 ccAGGCGGGCGCCCCAGTTACGCGGTTCGCGCGGGCGCAGGCGCTTGGCGTCGAAACTGCCTTTAAGCTTCATTCCGCTCTCCTCGGTGGATATCCGGAATTCGACCTGCTGCTGCCCATCGGGTTCCGCCCGGCTGCCGGGCGGTCCTGCAAGCCGCTCAGATCAGGCTGGCGGTGGGGGCCACGCAAGCGAAGTTGTCAGCCATCACGGCCATTTCGCACTGGTGAATCTgcgc                                                                                                                                     >  1:470789/1‑236 (MQ=255)
                                                                                                 ccAGGCGGGCGCCCCAGTTACGCGGTTCGCGCGGGCGCAGGCGCTTGGCGTCGAAACTGCCTTTAAGCTTCATTCCGCTCTCCTCGGTGGATATCCGGAATTCGACCTGCTGCTGCCCATCGGGTTCCGCCCGGCAGCCgg                                                                                                                                                                                                                                    >  2:401310/1‑141 (MQ=255)
                                                                                                 ccAGGCGGGCGCCCCAGTTACGCGGTTCGCGCGGGCGCAGGCGCTTGGCGTCGAAACTGCCTTTAAGCTTCATTCCGCTCTCCTCGGTGGATATCCGGAATTCGACCTGCTGCTGCCCATCGGGTTCCGCCCGGCAGCCgg                                                                                                                                                                                                                                    <  1:401313/141‑1 (MQ=255)
                                                                                                     gcgggcgCCCCAGTTACGCGGTTCGCGCGGGCGCAGGCGCTTGGCGTCGAAACTGCCTTTAAGCTTCATTCCGCTCTCCTCGGTGGATATCCGGAATTCGACCTGCTGCTGCCCATCGGGTTCCGCCCGGCTGcc                                                                                                                                                                                                                                      <  2:950763/135‑1 (MQ=255)
                                                                                                     gcgggcgCCCCAGTTACGCGGTTCGCGCGGGCGCAGGCGCTTGGCGTCGAAACTGCCTTTAAGCTTCATTCCGCTCTCCTCGGTGGATATCCGGAATTCGACCTGCTGCTGCCCATCGGGTTCCGCCCGGCAGcc                                                                                                                                                                                                                                      <  2:117392/135‑1 (MQ=255)
                                                                                                                  ttACGCGGTTCGCGCGGGCGCAGGCGCTTGGCGTCGAAACTGCCTTTAAGCTTCATTCCGCTCTCCTCGGTGGATATCCGGAATTCGACCTGCTGCTGCCCATCGGGTTCCGCCCGGCTGCCGGGCGGTCCTGCAAGCCGCTCAGATCAGGCTGGCGGTgg                                                                                                                                                                                               <  2:452277/161‑1 (MQ=255)
                                                                                                                             gcgcgGGCGCAGGCGCTTGGCGTCGAAACTGCCTTTAAGCTTCATTCCGCTCTCCTCGGTGGATATCCGGAATTCGACCTGCTGCTGCCCATCGGGTTCCGCCCGGCTGCCGGGCGGTCCTGCAAGCCGCTCAGATCAGGCTGGCGGTGGGGGcc                                                                                                                                                                                          <  1:285643/155‑1 (MQ=255)
                                                                                                                                     gcAGGCGCTTGGCGTCGAAACTGCCTTTAAGCTTCATTCCGCTCTCCCCGGTGGATATCCGGAATTCGACCTGCTGCTGCCCATCGGGTTCCGCCCGGCTGCCGGGCGGTCCTGCAAGCCGCTCa                                                                                                                                                                                                                <  2:292789/125‑1 (MQ=255)
                                                                                                                                       aGGCGCTTGGCGTCGAAACTGCCTTTAAGCTTCATTCCGCTCTCCTCGGTGGATATCCGGAATTCGACCTGCTGCTGCCCATCGGGTTCCGCCCGGCTGCCGGGCGGTCCTGCAAGCCGCTCAGATCAGGCTGGCGGTGGGGGCCACGCAAGCGAAGTTGTCAGCCATCACGGCCATTTCGCACTGGTGAATCTGCGCGGCCGGAATGACGCCGTc                                                                                                                   <  1:996169/216‑1 (MQ=255)
                                                                                                                                       aGGCGCTTGGCGTCGAAACTGCCTTTAAGCTTCATTCCGCTCTCCTCGGTGGATATCCGGAATTCGACCTGCTGCTGCCCATCGGGTTCCGCCCGGCTGCCGGGCGGTCCTGCAAGCCGCTCAGATCAGGCTGGCGGTGGGGGCCACGCAAGCGAAGTTGTCAGCCATCACGGCCATTTCGCACTGGTGAATCTGCGCGGCCGGAATGACGCCGTc                                                                                                                   >  2:996166/1‑216 (MQ=255)
                                                                                                                                       aGGCGCTTGGCGTCGAAACTGCCTTTAAGCTTCATTCCGCTCTCCTCGGTGGATATCCGGAATTCGACCTGCTGCTGCCCATCGGGTTCCGCCCGGCTGCCGGGCGGTCCTGCAAGCCGCTCAGATCAGGCTGGCGGTGGGGGCCACGCAAGCGAAGTTGTCAGCCATCACGGCCATTTCGCACTGGTGAATCTGCGCGGCCGGAATGACGCCGTc                                                                                                                   >  1:801098/1‑216 (MQ=255)
                                                                                                                                         gcgcTTGGCGTCGAAGCTTCCTTTCAGCTTCATTCCGCTCTCCTCGGTGGGTCTCCGGACTTCGCCCTGCTGCTGCCCATCGGGTTCCGCCCGGCTGCCGGGCGGTCCTGCAAGcc                                                                                                                                                                                                                     <  2:321994/116‑1 (MQ=255)
                                                                                                                                             ttGGCGTCGAAACTGCCTTTAAGCTTCATTCCGCTCTCCTCGGTGGATATCCGGAATTCGACCTGCTGCTGCCCATCGGGTTCCGCCCGGCTGCCGGGCGGTCCTGCAAGCCGCTCAGATCAGGCTGGCGGTGGGGGCCACGCAAGCGAAGTTGTCAGCCATCACGGCCATTTCGCACTGGTGAATCTGCGCGGCCGGAATGACGCCGTCCTTGAATGCCAGGTCACGGGTGGCCg                                                                                         <  1:1095147/236‑1 (MQ=255)
                                                                                                                                              tGGCGTCGAAACTGCCTTTAAGCTTCATTCCGCTCTCCTCGGTGGATATCCGGAATTCGACCTGCTGCTGCCCATCGGGTTCCGCCCGGCTGCCGGGCGGTCCTGCAAGCCGCTCAGATCAGGCTGGCGGTGGGGGCCACGCAAGCGAAGTTGTCAGCCATCACGGCCATTTCGCACTGGTGAATCTGCGCGGCCGGAATGACGCCGTCCTTGAATGCCAGGTCACGGGTGGCCg                                                                                         <  1:694586/235‑1 (MQ=255)
                                                                                                                                                      aaaCTGCCTTTAAGCTTCATTCCGCTCTCCTCGGTGGATATCCGGAATTCGACCTGCTGCTGCCCATCGGGTTCCGCCCGGCTGCCGGGCGGTCCTGCAAGCCGCTCAGATCAGGCTGGCGGTGGGGGCCACGCAAGCGAAGTTGTCAGCCATCACGGCCATTTCGCACTGGTGAATCTGCGCGGCCGGaa                                                                                                                             >  1:1039726/1‑191 (MQ=255)
                                                                                                                                                          tGCCTTTAAGCTTCATTCCGCTCTCCTCGGTGGATATCCGGAATTCGACCTGCTGCTGCCCATCGGGTTCCGCCCGGCTGCCGGGCGGTCCTGCAAGCCGCTCAGATCAGGCTGGCGGTGGGGGCCACGCAAGCGAAGTTGTCAGCCATCACGGCCAttt                                                                                                                                                        <  2:684629/160‑1 (MQ=255)
                                                                                                                                                            ccTTTAAGCTTCATTCCGCTCTCCTCGGTGGATATCCGGAATTCGACCTGCTGCTGCCCATCGGGTTCCGCCCGGCTGCCGGGCGGTCCTGCAAGCCGCTCAGATCAGGCTGGCGGTGGGGGCCACGCAAGCGAAGTTGTCAGCCATCACGGCCATTTCGCACTGGTGAATCTGCGCGGCCGGAATGACGCCGTCCTTGAATGCCAGGTCACGGGTGGCCGACGCGTCTTCTAcc                                                                           <  1:753618/235‑1 (MQ=255)
                                                                                                                                                                    cTTCATTCCGCTCTCCTCGGTGGATATCCGGAATTCGACTGGCTGCTGCCCATCGGGTTCCGCCCGGCTGCCGGGCGGTCCTGCAAGCCGCTCAGATCAGGCTGGCGGTGGGGGCCACGCAAGCGAAGTTGTCAGCCATCACGGCCATTTCGCACTGGTGAATCTGCGCGGCCGGAATGACGCCGTCCTTGAATGCCAGGTCACGGGTGGCCg                                                                                         <  1:101488/213‑1 (MQ=255)
                                                                                                                                                                     ttCATTCCGCTCTCCTCGGTGGATATCCGGAATTCGACCTGCTGCTGCCCATCGGGTTCCGCCCGGCTGCCGGGCGGTCCTGCAAGCCGCTCAGATCAGGCTGGCGGTGGGGGCCACGCAAGCGAAGTTGTCAGCCATCACGGCCATTTCGCACTGGTGAATCTGCGCGGCCGGAATGACGCCGTCCTTGAATGCCAGGTCACGGGTGGCCGACGCGTCTTCTACCAGGGTGCAcc                                                                 <  1:347067/236‑1 (MQ=255)
                                                                                                                                                                        aTTCCGCTCTCCTCGGTGGATATCCGGAATTCGACCTGCTGCTGCCCATCGGGTTCCGCCCGGCTGCCGGGCGGTCCTGCAAGCCGCTCAGATCAGGCTGGCGGTGGGGGCCACGCAAGCGAAGTTGTCAGCCATCACGGCCATTTCGCACTGGTGAATCTGCGCGGCCGGAATGACGCCGTCCTTGAATGCCAGGTCACGGGTGGCCGACGCGTCTTCTACCAGGGTGCAcc                                                                 <  2:871057/233‑1 (MQ=255)
                                                                                                                                                                         ttCCGCTCTCCTCGGTGGATATCCGGAATTCGACCTGCTGCTGCCCATCGGGTTCCGCCCGGCAGCCGGGCGGTCCTGCAAGCCGCTCAGATCAGGCTGGCg                                                                                                                                                                                                   >  2:285072/1‑102 (MQ=255)
                                                                                                                                                                         ttCCGCTCTCCTCGGTGGATATCCGGAAGTCGACCTGCTGCTGCCCATCGGGTTCCGCCCGGCTGCCGGGCGGTCCTGCAAGCCGCTCAGATCAGGCTGGCGGTGGGGGCCACGCAAGCGAAGTTGTCAGCCATCACGGCCATTTCGCACTGGTGAATCTGCGCGGCCGGAATGACGCCGTCCTTGAATGCCAGGTCACGGGTGGCCGAcgcgcc                                                                                  >  1:572346/1‑213 (MQ=255)
                                                                                                                                                                                    tcGGTGGATATCCGGAATTCGACCTGCTGCTGCCCATCGGGTTCCGCCCGGCTGCCGGGCGGTCCTGCACGCCGCTCAGATCAGGCTGGCGGTGGGGGCCACGCCCGCGACGtt                                                                                                                                                                            >  2:84381/1‑114 (MQ=255)
                                                                                                                                                                                    tcGGTGGATATCCGGAATTCGACCTGCTGCTGCCCATCGGGTTCCGCCCGGCTGCCGGGCGGTCCTGCAAGCCGCTCAGATCAGGCTGGCGGTGGGGGCCACGCAAGCg                                                                                                                                                                                 >  2:157614/1‑109 (MQ=255)
                                                                                                                                                                                    tcGGTGGATATCCGGAATTCGACCTGCTGCTGCCCATCGGGTTCCGCCCGGCTGCCGGGCGGTCCTGCAAGCCGCTCAGATCAGGCTGGCGGTGGGGGCCACGCAAGCg                                                                                                                                                                                 <  1:157614/109‑1 (MQ=255)
                                                                                                                                                                                    tcGGTGGATATCCGGAATTCGACCTGCTGCTGCCCATCGGGTTCCGCCCGGCAGCCGGGCGGTCCTGCAAGCCGCTCAGATCAGGCTGGCGGTGGGGGCCACGCAAGCGAAGTTGTCAGCCATCACGGCCATTTCGCACTGGTGAATCTGCGCGGCCGGAATGACGCCGTCCTTGAATGCCAGGTCAc                                                                                                  >  1:609360/1‑188 (MQ=255)
                                                                                                                                                                                      ggTGGATATCCGGAATTCGACCTGCTGCTGCCCATCGGGTTCCGCCCGGCTGCCGGGCGGTCCTGCAAGCCGCTCAGATCAGGCTGGCGGTGGGGGCCACGCAAGCGAAGTTGTCAGCCATCACGGCCATTTCGCACTGGTGAATCTGCGCGGCCGGAATGACGCCGTCCTTGAATGCCAGGTCACGGGTGGCCGACGCGTCTTCTACCAgg                                                                        >  2:793919/1‑212 (MQ=255)
                                                                                                                                                                                      ggTGGATATCCGGAATTCGACCTGCTGCTGCCCATCGGGTTCCGCCCGGCTGCCGGGCGGTCCTGCAAGCCGCTCAGATCAGGCTGGCGGTGGGGGCCACGCAAGCGAAGTTGTCAGCCATCACGGCCATTTCGCACTGGTGAATCTGCGCGGCCGGAATGACGCCGTCCTTGAATGCCAGGTCACGGGTGGCCGACGCGTCTTCTACCAgg                                                                        <  1:793922/212‑1 (MQ=255)
                                                                                                                                                                                      ggTGGATATCCGGAATTCGACCTGCTGCTGCCCATCGGGTTCCGCCCGGCTGCCGGGCGGTCCTGCAAGCCGCTCAGATCAGGCTGGCGGTGGGGGCCACGCAAGCGAAGTTGTCAGCCATCACGGCCATTTCGCACTGGGGAATCTGCGCGGCCGGAATGACGCCGTCCTTGAATGCCAGGTCACGGGGGGCCGACGCGTCTTCTACCAGGGTGCACCGATAACCAtc                                                       >  1:340309/1‑228 (MQ=255)
                                                                                                                                                                                         ggATATCCGGAATTCGACCTGCTGCTGCCCATCGGGTTCCGCCCGGCTGCCGGGCGGTCCTGCAAGCCGCTCAGATCAGGCTGGCGGTGGGGGCCACGCAAGCGAAGTTGTCAGCCATCACGGCCATTTCGCACTGGTGAATCTGCGCGGCCGGAATGACGCCGTCCTTGCATGCCAGGTCACGGGGGGCCGAcg                                                                                      >  1:142894/1‑195 (MQ=255)
                                                                                                                                                                                              tCCGGAATTCGACCTGCTGCTGCCCATCGGGTTCCGCCCGGCTGCCGGGCGGTCCTGCAAGCCGCTCAGATCAGGCTgg                                                                                                                                                                                                     >  2:1085773/1‑79 (MQ=255)
                                                                                                                                                                                                cGGAATTCGACCTGCTGCTGCCCATCGGGTTCCGCCCGGCTGCCGGGCGGTCCTGCAAGCCGCTCAGATCAGGCTGGCGGTGGGGGCCACGCAAGCGAAGTTGTCAGCCATCACGGCCATTTCGCACTGGTGAATCTGCTCGGCCGGAATGACGCCGTcc                                                                                                                  >  2:553379/1‑160 (MQ=255)
                                                                                                                                                                                                cGGAATTCGACCTGCTGCTGCCCAGCGGGTTCCGCCCGGCTGCCGGGCGGTCCTGCAAGCCGCTCAGATCAGGCTGGCGGTGGGGGCCACGCAAGCGAAGTTGTCAGCCATCACGGCCATTTCGCACTGGTGAATCTGCGCGGCCGGAATGACGCCGTCCTTGAATGCCa                                                                                                        <  2:626342/170‑1 (MQ=255)
                                                                                                                                                                                                        gcccTGCTGCTGCCGATCGGGTTCCGCCCGGCTGCCGGGCGGTCCTGCAAGCCGATCAGATCAGGCTGGCGGTGGGGGCCACGCAAGCGAAGTTGTCAGCCATCACGGCCATTTCGCACTGGTGAATCTGCGCGGCCGGAATGACGCCGTCCTTGAATGCCAGGTCACGGGTGGCCGAcgc                                                                                     <  1:644729/179‑1 (MQ=255)
                                                                                                                                                                                                        gACCTGCTGCTGCCCATCGGGTTCCGCCCGGCTGCCGGGCGGTCCTGCAAGCCGCTCAGATCAGGCTGGCGGTGGGGGCCACGCAAGCGAAGTTGTCAGCCATCACGGCCATTTCGCACTGGTGAATCTGCGCGGCCGGGATGACGCCGTCCTTGAATGCCAGGTCACGGGTGGCCGACGCGTCTTCTACCAGGGt                                                                      <  1:872290/196‑1 (MQ=255)
                                                                                                                                                                                                                      cATCGGGTTCCGCCCGGCTGCCGGGCGGTCCTGCAAGCCGCTCAGATCAGGCTGGCGGTGGGGGCCACGCAAGCGAAGTTGTCAGCCATCACGGCCATTTCGCACTGGTGAATCTGCGCGGCCGGAATGACGCCGTCCTTGAATGCCAGGTCACGGGTGGCCGACGCGTCTTCTACCAGGGTGCACCGATAACCATAGTCCTTGGCGCGGCGAACGGTGGTGCTGACGCTGGAATg                >  1:483586/1‑236 (MQ=255)
                                                                                                                                                                                                                       aTCGGGTTCCGCCCGGCGGCCGGGCGGTCCTGCAAGCCGCGCAGATCAGGCTGGCGGTGGGGGTCACGCAAGCGAAGTTGTCAGCCATCACGGCCATTTCGCACTGGTGAATCTGCGCGGCCGGAATGACGc                                                                                                                       <  2:9720/132‑1 (MQ=255)
                                                                                                                                                                                                                         cGGGTTCCGCCCGGCTGCCGGGCGGTCCTGCAAGCCGCTCAGATCAGGCTGGCGGTGGGGGCCACGCAAGCGAAGTTGTCAGCCATCACGGCCATTTCGCACTGGTGAATCTGCGCGGCCGGAATGACGCCGTCCTTGAATGCCAGGTCACGGGTGGCCGCCGCGt                                                                                   >  2:16208/1‑166 (MQ=255)
                                                                                                                                                                                                                         cGGGTTCCGCCCGGCTGCCGGGCGGTCCTGCAAGCCGCTCAGATCAGGCTGGCGGTGGGGGCCACGCAAGCGAAGTTGTCAGCCATCACGGCCATTTCGCACTGGTGAATCTGCGCGGCCGGAATGACGCCGTCCTTGAATGCCAGGTCACGGGTGGCCGACGCCTCTTCTACCAGGGTGCACCGATAACCATAGTCCTTGGCGCGGCGAACGgtggtg                              >  2:772920/1‑219 (MQ=255)
                                                                                                                                                                                                                           ggTTCCGCCCGGCTGCCGGGCGGTCCTGCAAGCCGCTCAGATCAGGCTGGCGGTGGGGGCCACGCAAGCGAAGTTGTCAGCCATCACGGCCATTTCGCACTGGTGAATCTGCGCGGCCGGAATGACGCCGTCCTTGAATGCCAGGTCACGGGTGGCCGACGCGTCTTCTACCAGGGTGCACCGATAACCATAGTCCTTGGCGCGGCGAACGGTGGTGCTGACGCTGGAATGGCTCa           <  1:299794/236‑1 (MQ=255)
                                                                                                                                                                                                                                cGCCCGGCTGCCGGGCGGTCCTGCAAGCCGCTCAGATCAGGCTGGCGGTGGGGGCCACGCAAGCGAAGTTGTCAGCCATCACGGCCATTTCGCACTGGTGAAGCTGCGCGGCCGGAATGACGCCGTCCTTGAATGCCAGGTCACGGGTGGCCg                                                                                         <  2:364537/153‑1 (MQ=255)
                                                                                                                                                                                                                                  cccGTCTGCCGGGCGGTCCTGCAAGCCGCTCAGATCAGGCTGGCGGTGGGGGCCGCGCAAGCGAAGTTGCCAGCCATCACGGCCATTTCGCACTGGTGAATCTGCGCGGCCGGAATGACGCCGTCCTTGAATGCCAGGTCACGGGTGGCCGACGCGTCTTCTACCAGGGTGCACCGATAACCATAGTCCTTGGCGCGGCGa                                       <  1:309451/201‑1 (MQ=255)
                                                                                                                                                                                                                                   ccGGCTGCCGGGCGGTCCTGCAAGCCGCTCAGATCAGGCTGGCGGTGGGGGCCACGCAAGCGAAGTTGTCAGCCATCACGGCCATTTCGCACTGGTGAATCTGCGCGGCCGGAATGACGCCGTCCTTGAAt                                                                                                            >  2:1260294/1‑131 (MQ=255)
                                                                                                                                                                                                                                   ccGGCAGCCGGGCGGTCCTGCAAGCCGCTCAGATCAGGCTGGCGGTGGGGGCCACGCAAGCGAAGTTGTCAGCCATCACGGCCATTTCGCACTGGTGAATCTGCGCGGCCGGAATGACGCCGTCCTTGAATGCCAGGTCACGGGTGGCCGACGCGTCTTCTACCAGGGTGCACCGATAACCATAGTCCTTGGCGCGGCGAACGGTGGTGCTGACGCTGGAATGGCTCATGAAAcc    <  1:1229520/235‑1 (MQ=255)
                                                                                                                                                                                                                                     tgCTGCCGGGCGGTCCTGCAAGCCGCTCAGATCAGGCTGGCGGCGGGGGCCACGCAAGCGAAGTTGTCAGCCATCACGGCCATTTCTCACTGGGGAATCTGCGCGGCCg                                                                                                                                >  1:341110/2‑109 (MQ=255)
                                                                                                                                                                                                                                     ggCTGCCGGGCGGTCCTGCAAGCCGCTCAGATCAGGCTGGCGGTGGGGGCCACGCAAGCGAAGTTGTCAGCCATCACGGCCATTTCGCACTGGTGAATCTGCGCGGCCGGAATGACGCCGTCCTTGAATGCCAGGTCACGGGTGGCCGACGCGTCTTCTACCAGGGTGCACCGATAACCATAGTCCTTGGCGCGGCGAACGGTGGTGCTGACGCTGGAATGGCTCATGAAACcac  >  1:786694/1‑235 (MQ=255)
                                                                                                                                                                                                                                        |                                                                                                                                                                                                                                       
AGTTCCTGCAGGGCAGGGTAGTTGCCTAGCAGCCCGGCGACGCCGGCCATGGCCAGCAATACGCCCAGCGTCGCCAGCAGGGCCGACAGGCCGGCCGCCAGGCGGGCGCCCCAGTTACGCGGTTCGCGCGGGCGCAGGCGCTTGGCGTCGAAACTGCCTTTAAGCTTCATTCCGCTCTCCTCGGTGGATATCCGGAATTCGACCTGCTGCTGCCCATCGGGTTCCGCCCGGCTGCCGGGCGGTCCTGCAAGCCGCTCAGATCAGGCTGGCGGTGGGGGCCACGCAAGCGAAGTTGTCAGCCATCACGGCCATTTCGCACTGGTGAATCTGCGCGGCCGGAATGACGCCGTCCTTGAATGCCAGGTCACGGGTGGCCGACGCGTCTTCTACCAGGGTGCACCGATAACCATAGTCCTTGGCGCGGCGAACGGTGGTGCTGACGCTGGAATGGCTCATGAAACCAC  >  AE015451/2050833‑2051296

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