Predicted mutation
evidence seq id position mutation freq annotation gene description
RA AE015451 4,083,052 T→A 5.2% G17G (GGT→GGA) ‡ PP_3593 → Amino acid ABC transporter, periplasmic binding protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*AE0154514,083,0520TA 5.2% 263.2 / 5.8 77G17G (GGT→GGA) ‡PP_3593Amino acid ABC transporter, periplasmic binding protein
Reads supporting (aligned to +/- strand):  ref base T (40/31);  new base A (2/2);  total (44/33)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.25e-02

TGGGAATCCACCTCGACGATGCTGTCGTGGGTGCATCTGATGGTCAACGAGGTCTGCAAGAAGCTCGGACCGGATGTTGAAAAACGCTTGGAAGCAACTGCCGCTCTACATAACGAATTCGTCGGCTACACCTCGTGGTCGACGGGCAAACAACAATAGCAAGAGTGCAATGAGGTGCGAGATGAACAAGTCCATGGCCTTGCTGGGTGCGTGCGCCCTGCTGCTGGCGGGTGCCGCCAGCGCTGAAACCCTGCGGTTCGCCACCGAGGGTGCCTACCCGCCCTTCAACTATGTAGACGCCGATAACAAGCTGCACGGCTTTGATGTCGACATCACCCATGCCCTGTGCGAACAGATGAAAGTCGAATGCACGCTGGTGGCGCAGGACTGGGAAGGCATCATCCCCGCCCTGATGGCGCGCAAGTATGACGCGATAGTGGCATCGATGATCGACAC  >  AE015451/4082821‑4083276
                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                
tGGGAATCCACCTCGACGATGCTGTCGTGGGTGCATCTGATGGTCAACGAGGTCTGCAAGAAGCTCGGACCGGATGTTGAAAAACGCTTGGAAGCAACTGCCGCTCTACATAACGAATTCGTCGGCTACACCTCGTGGTCGACGGGCAAACAACAATAGCAAGAGTGCAATGAGGTGCGAGATGAACAAGTCCATGGCCTTGCTGGGTGCGTGCGCCCTGCTGCTGGCGGGTg                                                                                                                                                                                                                                 >  2:632705/1‑233 (MQ=255)
 gggAATCCACCTCGACGATGCTGTCGTGGGTGCATCTGATGGTCAACGAGGTCTGCAAGAAGCTCGGACCGGATGTTGAAAAACGCTTGGAAGCAACTGCCGCTCTACATAACGAATTCGTCGGCTACACCTCGTGGTCGACGGGCAAACAACAATAGCAAGAGTGCAATGAGGTGCGAGATGAACAAGTCCATGGCCTTGCTGGGTGCGTGCGCCCTGCTGCTGGCGGGTgccgc                                                                                                                                                                                                                             >  1:1072017/1‑236 (MQ=255)
   gAATCCACCTCGACGATGCTGTCGTGGGTGCATCTGATGGTCAACGAGGTCTGCAAGAAGCTCGGACCGGATGTTGAAAAACGCTTGGAAGCAACTGCCGCTCTACATAACGAATTCGTCGGCTACACCTCGTGGTCGACGGGCAAACAACAATAGCAAGAGTGCAATGAGGTGCGAGATGAACAAGTCCATGGCCTTGCTGGGTGCGTGCGCCCTGCTGCTGGCGGGTGCTGCCa                                                                                                                                                                                                                           <  1:632708/236‑1 (MQ=255)
          ccTCGACGATGCTGTCGTGGGTGCATCTGATGGTCAACGAGGTCTGCAAGAAGCTCGGACCGGATGTTGAAAAACGCTTGGAAGCAACTGCCGCTCTACATAACGAATTCGTCGGCTACACCTCGTGGTCGACGGGCAAACAACAATAGCAAGAGTGCAATGAGGTGCGAGATGAACAAGTCCATGGCCTTGCTGGGTGCGTGCGCCCTGCTGCTGGCGGGTg                                                                                                                                                                                                                                 >  1:372402/1‑223 (MQ=255)
          ccTCGACGATGCTGTCGTGGGTGCATCTGATGGTCAACGAGGTCTGCAAGAAGCTCGGACCGGATGTTGAAAAACGCTTGGAAGCAACTGCCGCTCTACATAACGAATTCGTCGGCTACACCTCGTGGTCGACGGGCAAACAACAATAGCAAGAGTGCAATGAGGTGCGAGATGAACAAGTCCATGGCCTTGCTGGGTGCGTGCGCCCTGCTGCTGGCGGGTGCCGCCAGCGCTGa                                                                                                                                                                                                                    >  1:1264705/1‑236 (MQ=255)
             cgacgaTGCTGTCGTGGGTGCATCTGATGGTCAACGAGGTCTGCAAGAAGCTCGGACCGGATGTTGAAAAACGCTTGGAAGCAACTGCCGCTCTACATAACGAATTCGTCGGCTACACCTCGTGGTCGACGGGCAAACAACAATAGCAAGAGTGCAATGAGGTGCGAGATGAACAAGTCCATGGCCTTGCTGGGTGCGTGCGCCCTGCTGCTGGCGGGTGCCGCCAGCGCTg                                                                                                                                                                                                                     >  2:1091209/1‑232 (MQ=255)
                cgaTGCTGTCGTGGGTGCATCTGATGGTCAACGAGGTCTGCAAGAAGCTCGGACCGGATGTTGAAAAACGCTTGGAAGCAACTGCCGCTCTACATAACGAATTCGTCGGCTACACCTCGTGGTCGACGGGCAAACAACAATAGCAAGAGTGCAATGAGGTGCGAGATGAACAAGTCCATGGCCTTGCTGGGTGCGTGCGCCCTGCTGCTGGCGGGTGCCGCCAGCGCTGAAAcc                                                                                                                                                                                                                <  2:694099/234‑1 (MQ=255)
                              gTGCATCTGATGGTCAACGAGGTCTGCAAGAAGCTCGGACCGGATGTTGAAAAACGCTTGGAAGCAACTGCCGCTCTACATAACGAATTCGTCGGCTACACCTCGTGGTCGACGGGCAAACAACAATAGCAAGAGTGCAATGAGGTGCGAGATGAACAAGTCCATGGCCTTGCTGGGTGCGTGCGCCCTGCTGCTGGCGGGTGCCGCCAGCGCTGAAACCCTGCGGTTCGCCAc                                                                                                                                                                                                  <  2:1072014/234‑1 (MQ=255)
                                gCATCTGATGGTCAACGAGGTCTGCAAGAAGCTCGGACCGGATGTTGAAAAACGCTTGGAAGCAACTGCCGCTCTACATAACGAATTCGTCGGCTACACCTCGTGGTCGACGGGCAAACAACAATAGCAAGAGTGCAATGAGGTGCGAGATGAACAAGTCCATGGCCTTGCTGGGTGCGTGCGCCCTGCTGCTGGCGGGTGCCGCCAGCGCTGAAACCCTGCGGTTCGCCACCGAg                                                                                                                                                                                              >  1:1078427/1‑236 (MQ=255)
                                          gTCAACGAGGTCTGCAAGAAGCTCGGACCGGATGTTGAAAAACGCTTGGAAGCAACTGCCGCTCTACATAACGAATTCGTCGGCTACACCTCGTGGTCGACGGGCAAACAACAATAGCAAGAGTGCAATGAGGTGCGAGATGAACAAGTCCATGGCCTTGCTGGGTGCGTGCGCCCTGCTGCTGGCGGGTGCCGCCAGCGCTGAAACCCTGCGGTTCGCCACCGAGGGTGCCt                                                                                                                                                                                       <  2:1264702/233‑1 (MQ=255)
                                          gTCAACGAGGTCTGCAAGAAGCTCGGACCGGATGTTGAAAAACGCTTGGAAGCAACTGCCGCTCTACATAACGAATTCGTCGGCTACACCTCGTGGTCGACGGGCAAACAACAATAGCAAGAGTGCAATGAGGTGCGAGATGAACAAGTCCATGGCCTTGCTGGGTGCGTGCGCCCTGCTGCTGGCGGGTGCCGCCAGCGCTGAAACCCTGCGGTTCGCCACCGAGGGTGCCTAc                                                                                                                                                                                     <  1:1044844/235‑1 (MQ=255)
                                           tCAACGAGGTCTGCAAGAAGCTCGGACCGGATGTTGAAAAACGCTTGGAAGCAACTGCCGCTCTACATAACGAATTCGTCGGCTACACCTCGTGGTCGACGGGCAAACAACAATAGCAAGAGTGCAATGAGGTGCGAGATGAACAAGTCCATGGCCTTGCTGGGTGCGTGCGCCCTGCTGCTGGCGGGTGCCGCCAGCGCTGAAACCCTGCGGTTCGCCAc                                                                                                                                                                                                  <  1:444453/221‑1 (MQ=255)
                                                               cTCGGACCGGATGTTGAAAAACGCTTGGAAGCAACTGCCGCTCTACATAACGAATTCGTCGGCTACACCTCGTGGTCGACGGGCAAACAACAATAGCAAGAGTGCAATGAGGTGCGAGATGAACAAGTCCATGGCCTTGCTGGGTGCGTGCGCCCTGCTGCTGGCGGGTGCCGCCAGCGCTGAAACCCTGCGGTTCGCCACCGAGGGTGCCTACCCGCCCTTCAACTATGTAGa                                                                                                                                                                 <  1:337936/234‑1 (MQ=255)
                                                                 cGGACCGGATGTTGAAAAACGCTTGGAAGCAACTGCCGCTCTACATAACGAATTCGTCGGCTACACCTCGTGGTCGACGGGCAAACAACAATAGCAAGAGTGCAATGAGGTGCGAGATGAACAAGTCCATGGCCTTGCTGGGTGCGTGCGCCCTGCTGCTGGCGGGTGCCGCCAgcg                                                                                                                                                                                                                        >  2:813503/1‑177 (MQ=255)
                                                                 cGGACCGGATGTTGAAAAACGCTTGGAAGCAACTGCCGCTCTACATAACGAATTCGTCGGCTACACCTCGTGGTCGACGGGCAAACAACAATAGCAAGAGTGCAATGAGGTGCGAGATGAACAAGTCCATGGCCTTGCTGGGTGCGTGCGCCCTGCTGCTGGCGGGTGCCGCCAGCGCTGAAACCCTGCGGTTCGCCACCGAGGGTGCCTACCCGCCCTTc                                                                                                                                                                            <  1:16248/221‑1 (MQ=255)
                                                                 cGGACCGGATGTTGAAAAACGCTTGGAAGCAACTGCCGCTCTACATAACGAATTCGTCGGCTACACCTCGTGGTCGACGGGCAAACAACAATAGCAAGAGTGCAATGAGGTGCGAGATGAACAAGTCCATGGCCTTGCTGGGTGCGTGCGCCCTGCTGCTGGCGGGTGCCGCCAGCGCTGAAACCCTGCGGTTCGCCACCGAGGGTGCCTACCCGCCCTTCAACTATGTAGACGcc                                                                                                                                                             >  1:852407/1‑236 (MQ=255)
                                                                        gATGTTGAAAAACGCTTGGAAGCAACTGCCGCTCTACATAACGAATTCGTCGGCTACACCTCGTGGTCGACGGGCAAACAACAATAGCAAGAGTGCAATGAGGTGCGAGATGAACAAGTCCATGGCCTTGCTGGGTGCGTGCGCCCTGCTGCTGGCGGGTGCCGCCAGCGCTGAAACCCTGCGGTTCGCCACCGAGGGTGCCTACCCGCCCTTCAACTATGTAGACGCCGATAAc                                                                                                                                                       <  1:1277276/235‑1 (MQ=255)
                                                                            ttGAAAAACGCTTGGAAGCAACTGCCGCTCTACATAACGAATTCGTCGGCTACACCTCGTGGTCGACGGGCAAACAACAATAGCAAGAGTGCAATGAGGTGCGAGATGAACAAGTCCATGGCCTTGCTGGGTGCGTGCGCCCTGCTGCTGGCGGGTGCCGCCAGCGCTGAAACCCTGCGGTTCGCCACCGAGGGTGCCTACCCGACCTTCAACTATGTAGACGCCGAt                                                                                                                                                          >  2:909747/1‑228 (MQ=255)
                                                                             tGAAAAACGCTTGGAAGCAACTGCCGCTCTACATAACGAATTCGTCGGCTACACCTCGTGGTCGACGGGCAAACAACAATAGCAAGAGTGCAATGAGGTGCGAGATGAACAAGTCCATGGCCTTGCTGGGTGCGTGCGCCCTGCTGCTGGCGGGTGCCGCCAGCGCTGAAACCCTGCGGTTCGCCACCGAGGGTGCCTACCCGCCCTTCAACTATGTAGACGCCGATAACAAGCt                                                                                                                                                  <  2:852404/235‑1 (MQ=255)
                                                                                    cGCTTGGAAGCAACTGCCGCTCTACATAACGAATTCGTCGGCTACACCTCGTGGTCGACGGGCAAACAACAATAGCAAGAGTGCAATGAGGTGCGAGATGAACAAGTCCATGGCCTTGCTGGGTGCGTGCGCCCTGCTGCTGGCGGGt                                                                                                                                                                                                                                  <  1:1137747/148‑1 (MQ=255)
                                                                                     gCTTGGAAGCAACTGCCGCTCTACATAACGAATTCGTCGGCTACACCTCGTGGTCGACGGGCAAACAACAATAGCAAGAGTGCAATGAGGTGCGAGATGAACAGGTCCATGGCCTTGCTGGGTGCGTGCGCCCTGCTGCTGGCGGGTGCCGCCAGCGCTGAAACCCTGCGGTTCGCCACCGAg                                                                                                                                                                                              <  2:448915/183‑1 (MQ=255)
                                                                                         ggAAGCAACTGCCGCTCTACATAACGAATTCGTCGGCTACACCTCGTGGTCGACGGGCAAACAACAATAGCAAGAGTGCAATGAGGTGCGAGATGAACAAGTCCATGGCCTTGCTGGGTGCGTGCGCCCTGCTGCTGGCGGGTGCCGCCac                                                                                                                                                                                                                          >  2:287072/1‑150 (MQ=255)
                                                                                         ggAAGCAACTGCCGCTCTACATAACGAATTCGTCGGCTACACCTCGTGGTCGACGGGCAAACAACAATAGCAAGAGTGCAATGAGGTGCGAGATGAACAAGTCCATGGCCTTGCTGGGTGCGTGCGCCCTGCTGCTGGCGGGTGCCGCCAGCGCTGAAAccc                                                                                                                                                                                                               >  1:236274/1‑162 (MQ=255)
                                                                                         ggAAGCAACTGCCGCTCTACATAACGAATTCGTCGGCTACACCTCGTGGTCGACGGGCAAACAACAATAGCAAGAGTGCAATGAGGTGCGAGATGAACAAGTCCATGGCCTTGCTGGGTGCGTGCGCCCTGCTGCTGGCGGGTGCCGCCAGCGCTGAAACCCTGCGGTTCGCCACCGAGGGTGCCTACCCGCCCTTCAACTATGTAGACGCCGATAACaa                                                                                                                                                     >  2:872429/1‑220 (MQ=255)
                                                                                                     cGCTCTACATAACGAATTCGTCGGCTACACCTCGTGGTCGACGGGCAAACAACAATAGCAAGAGTGCAATGAGGTGCGAGATGAACAAGTCCATGGCCTTGCTGGGTGCGTGCGCCCTGCTGCTGGCGGGTGCCGCCAGCGCTGAAACCCTGCGGTTCGCCACCGAGGGTGCCTACCCGccc                                                                                                                                                                               <  2:979317/182‑1 (MQ=255)
                                                                                                           aCATAACGAATTCGTCGGCTACACCTCGTGGTCGACGGGCAAACAACAATAGCAAGAGTGCAATGAGGTGCGAGATGAACAAGTCCATGGCCTTGCTGGGTGCGTGCGCCCTGCTGCTGGCGGGTGCCGCCAGCGCTGAAACCCTGCGGTTCGCCACCGAGGGTGCCTACCCGCCCTTCAACTATGTAGACGCCGATAACAAGCTGCACGGCTTTGATGTCGACATCACCCATg                                                                                                                     <  1:781433/234‑1 (MQ=255)
                                                                                                                aCGAATTCGTCGGCTACACCTCGTGGTCGACGGGCAAACAACAATAGCAAGAGTGCAATGAGGTGCGAGATGAACAAGTCCATGGCCTTGCTGGGTGCGTGCGCCCTGCTGCTGGCGGGTGCCGCCAGCGCTGAAACCCTGCGGTTCGCCACCGAGGGTGCCTAc                                                                                                                                                                                     <  2:836567/165‑1 (MQ=255)
                                                                                                                    aTTCGTCGGCTACACCTCGTGGTCGACGGGCAAACAACAATAGCAAGAGTGCAATGAGGTGCGAGATGAACAAGTCCATGGCCTTGCTGGGTGCGTGCGCCCTGCTGCTGGCGGGTGCCGCCAGCGCTGAAACCCTGCGGTTCGCCACCGAGGGTGCCTACCCGCCCTTCAACTATGTAGACGCCGATAACAAGCTGCACGGCTTTGATGTCGACATCACCCATGCCCTGTGc                                                                                                             >  1:1026190/1‑233 (MQ=255)
                                                                                                                       cgtcgGCTACACCTCGTGGTCGACGGGCAAACAACAATAGCAAGAGTGCAATGAGGTGCGAGATGAACAAGTCCATGGCCTTGCTGGGTGCGTGCGCCCTGCTGCTGGCGGCACCCGCCAGCGCTGAAACCCTGCGGTTCGCCa                                                                                                                                                                                                   <  2:1100077/144‑1 (MQ=255)
                                                                                                                       cgtcgGCTACACCTCGTGGTCGACGGGCAAACAACAATAGCAAGAGTGCAATGAGGTGCGAGATGAACAAGTCCATGGCCTTGCTGGGTGCGTGCGCCCTGCTGCTGGCGGCACCCGCCAGCGCTGAAACCCTGCGGTTCGCCa                                                                                                                                                                                                   >  1:1100080/1‑144 (MQ=255)
                                                                                                                         tcgGCTACACCTCGTGGTCGACGGGCAAACAACAATAGCAAGAGTGCAATGAGGTGCGAGATGAACAAGTCCATGGCCTTGCTGGGTGCGTGCGCCCTGCTGCTGGCGGGTGCCGCCCGCGCTGAAACCCTGCGGTTCGCCACCGAgg                                                                                                                                                                                             >  2:781430/1‑148 (MQ=255)
                                                                                                                         tcgGCTACACCTCGTGGTCGACGGGCAAACAACAATAGCAAGAGTGCAATGAGGTGCGAGATGAACAAGTCCATGGCCTTGCTGGGTGCGTGCGCCCTGCTGCTGGCGGGTGCCGCCAGCGCTGAAACCCTGCGGt                                                                                                                                                                                                         >  1:288015/1‑136 (MQ=255)
                                                                                                                                   cTCGTGGTCGACGGGCAAACAACAATAGCAAGAGTGCAATGAGGTGCGAGATGAACAAGTCCATGGCCTTGCTGGGTGCGTGCGCCCTGCTGCTGGCGGGTGCCGCCAGCGCTGAAACCCTGCGGTTCg                                                                                                                                                                                                      <  1:1103735/129‑1 (MQ=255)
                                                                                                                                   cTCGTGGTCGACGGGCAAACAACAATAGCAAGAGTGCAATGAGGTGCGAGATGAACAAGTCCATGGCCTTGCTGGGTGCGTGCGCCCTGCTGCTGGCGGGTGCCGCCAGCGCTGAAACCCTGCGGTTCg                                                                                                                                                                                                      >  2:1103732/1‑129 (MQ=255)
                                                                                                                                   cTCGTGGTCGACGGGCAAACAACAATAGCAAGAGTGCAATGAGGTGCGAGATGAACAAGTCCATGGCCTTGCTGGGTGCGTGCGCCCTGCTGCTGGCGGGTGCCGCCAGCGCTGAAACCCTGCGGTTCGCCACCGAGGGTGCCTACCCGCCCTTCAACTATGTAGACGCCGATAACAAGCTGCACGGCTTTGATGTCGACATCACCCATGccc                                                                                                                  >  1:439324/1‑213 (MQ=255)
                                                                                                                                   cTCGTGGTCGACGGGCAAACAACAATAGCAAGAGTGCAATGAGGTGCGAGATGAACAAGTCCATGGCCTTGCTGGGTGCGTGCGCCCTGCTGCTGGCGGGTGCCGCCAGCGCTGAAACCCTGCGGTTCGCCACCGAGGGTGCCTACCCGCCCTTCAACTATGTAGACGCCGATAACAAGCTGCACGGCTTTGATGTCGACATCACCCATGCCCTGTGCGAACAGATGAAAGTCGa                                                                                            <  1:13954/235‑1 (MQ=255)
                                                                                                                                    tCGTGGTCGACGGGCAAACAACAATAGCAAGAGTGCAATGAGGTGCGAGATGAACAAGTCCATGGCCTTGCTGGGTGCGTGCGCCCTGCTGCTGGCGGGTGCCGCCAGCGCTGAAACCCTGCGGTTCGCCACCGAGGGTGCCTACCCGCCCTTCAACTATGTAGACGCCGATAACAAGCTGCACGGCTTTGATGTCGACATCACCCATGCCCTGTGCGAACAGATGAAAGTCGAAt                                                                                          <  1:452051/236‑1 (MQ=255)
                                                                                                                                      gTGGTCGACGGGCAAACAACAATAGCAAGAGTGCAATGAGGTGCGAGATGAACAAGTCCATGGCCTTGCTGGGTGCGTGCGCCCTGCTGCTGGCGGGTgcc                                                                                                                                                                                                                               <  2:951155/101‑1 (MQ=255)
                                                                                                                                      gTGGTCGACGGGCAAACAACAATAGCAAGAGTGCAATGAGGTGCGAGATGAACAAGTCCATGGCCTTGCTGGGTGCGTGCGCCCTGCTGCTGGCGGGTgcc                                                                                                                                                                                                                               >  1:951158/1‑101 (MQ=255)
                                                                                                                                      gTGGTCGACGGGCAAACAACAATAGCAAGAGTGCAATGAGGTGCGAGATGAACAAGTCCATGGCCTTGCTGGGTGCGTGCGCCCTGCTGCTGGCGGGTGCCGCCAGCGCTGAAACCCTGCGGTTCGCCACCGAGGGTGCCTACCCGCCCTTCAACTATGTAGACg                                                                                                                                                               <  2:559819/165‑1 (MQ=255)
                                                                                                                                       tGGTCGACGGGCAAACAACAATAGCAAGAGTGCAATGAGGTGCGAGATGAACAAGTCCATGGCCTTGCTGGGTGCGTGCGCCCTGCTGCTGGCGGGTGCCGCCAGCGCTGAAACCCTGCGGTTCGCCACCGAGGGTGCCt                                                                                                                                                                                       <  2:589356/140‑1 (MQ=255)
                                                                                                                                        ggTCGACGGGCAAACAACAATAGCAAGAGTGCAATGAGGTGCGAGATGAACAAGTCCATGGCCTTGCTGGGTGCGTGCGCCCTGCTGCTGGCGGGTGCCGCCAGCGCTGAAACCCTGCGGTTCGCCACCGAGGGTGCCTACCCGCCCTTCAACTATGTAGACGCCGATAACAAGCTGc                                                                                                                                                >  1:1237473/1‑178 (MQ=255)
                                                                                                                                         gTCGACGGGCAAACAACAATAGCAAGAGTGCAATGAGGTGCGAGATGAACAAGTCCATGGCCTTGCTGGGTGCGTGCGCCCTGCTGCTGGCGGGGGCCGCCac                                                                                                                                                                                                                          >  2:287425/1‑102 (MQ=255)
                                                                                                                                          tCGACGGGCAAACAACAATAGCAAGAGTGCAATGAGGTGCGAGATGAACAAGTCCATGGCCTTGCTGGGTGCGTGCGCCCTGCTGCTGGCGGGTGCCGCCAGCGCTGAAACCCTGCGGTTCGCCACCGAGGGTGCCTACCCGCCCTTCAACTATGTAGACGCCGATAACAAGCTGCACGGCTTTGATGTCGACATCACCCATGCCCTGTGCGAACAGATGAAAGTCGAATGCACGc                                                                                    <  1:1046480/236‑1 (MQ=255)
                                                                                                                                                ggCAAACAACAATAGCAAGAGTGCAATGAGGTGCGAGATGAACAAGTCCATGGCCTTGCTGGGTGCGTGCGCCCTGCTGCTGGCGGGTGCCGCCAGCGCTGAAACCCTGCGGTTCGCCACCGAGGGTGCCTACCCGCCCTTCAACTATGTAGACGCCGATAACAAGCTGCACGGCTTTGATGTCGACATCACCCATGc                                                                                                                    >  2:1085201/1‑198 (MQ=255)
                                                                                                                                                    aacaacAATAGCAAGAGTGCAATGAGGTGCGAGATGAACAAGTCCATGGCCTTGCTGGGTGCGTGCGCCCTGCTGCTGGCGGGTGCCGCCAGCGCTGAAACCCTg                                                                                                                                                                                                             >  1:134976/1‑105 (MQ=255)
                                                                                                                                                    aacaacAATAGCAAGAGTGCAATGAGGTGCGAGATGAACAAGTCCATGGCCTTGCTGGGTGCGTGCGCCCTGCTGCTGGCGGGTGCCGCCAGCGCTGAAACCCTGCGGTTCGCCACCGAGGGTGCCTACCCGCCCTTCAACTATGTAGACGCCGATAACAAGCTGCACGGCTTTGATGTCGACATCACCCATGCCCTGTGCGAACAGATGAAAGTCGAATGCACGCTGGTGGCGCa                                                                          <  1:19156/236‑1 (MQ=255)
                                                                                                                                                            tAGCAAGAGTGCAATGAGGTGCGAGATGAACAAGTCCATGGCCTTGCTGGGTGCGTGCGCCCTGCTGCTGGCGGGTGCCGCCAGCGCTGAAACCCTGCGGTTCGCCACCGAGGGTGCCTACCCGCCCTTCAACTATGTAGACGCCGATAACAAGCt                                                                                                                                                  <  2:894642/156‑1 (MQ=255)
                                                                                                                                                              gCAAGAGTGCAATGAGGTGCGAGATGAACAAGTCCATGGCCTTGCTGGGTGCGTGCGCCCTGCTGCTGGCGGGTgccgc                                                                                                                                                                                                                             >  2:346263/1‑79 (MQ=255)
                                                                                                                                                              gCAAGAGTGCAATGAGGTGCGAGATGAACAAGTCCATGGCCTTGCTGGGTGCGTGCGCCCTGCTGCTGGCGGGTGCCGCCAGCGCTGAAACCCTGCGGTTCGCCACCGAGGGTGCCTACCCGCCCTTCAACTATGTAGACGCCGATAACAAGCTGCACGGCTTTGATGTCGACATCACCCATGCCCTGTGCGAACAGATGAAAGTCGAATGCACGCTGGTGGCGCAGGACTGGGa                                                                 >  1:854485/1‑235 (MQ=255)
                                                                                                                                                                    gTGCAATGAGGTGCGAGATGAACAAGTCCATGGCCTTGCTGGGTGCGTGCGCCCTGCTGCTGGCGGGTGCCGCCAGCGCTGAAACCCTggggg                                                                                                                                                                                                         >  1:359144/1‑92 (MQ=255)
                                                                                                                                                                    gTGCAATGAGGTGCGAGATGAACAAGTCCATGGCCTTGCTGGGTGCGTGCGCCCTGCTGCTGGCGGGTGCCGCCAGCGCTGAAACCCTGCGGTTCGCCACCGAGGGTGCCTACCCGCCCTTCAACTATGTAGACGCCGATAACAAGCTGCACGGCTTTGAt                                                                                                                                     >  2:414473/1‑161 (MQ=255)
                                                                                                                                                                      gCAATGAGGTGCGAGATGAACAAGTCCATGGCCTTGCTGGGTGCGTGCGCCCTGCTGCTGGCGGGTGGCGCCAGCGCTGAAACCCCGCGGGTCGCCACCCGGGGGGCCCACCCg                                                                                                                                                                                  >  1:227505/1‑114 (MQ=255)
                                                                                                                                                                      gCAATGAGGTGCGAGATGAACAAGTCCATGGCCTTGCTGGGTGCGTGCGCCCTGCTGCTGGCGGGTGCCGCCAGCGCTGAAACCCTGCGGTTCGCCACCGAGGGTGCCTACCCGCCCTTCAACTATGTAGACGCCGATAACAAGCTGCACGGCTTTGATGTCGACATCACCCATGCCCTGTGCGAACAGATGAAAGTCGAATGCACGCTGGTGGCGCAGGACTGGGAAGGcatca                                                         >  1:929785/1‑235 (MQ=255)
                                                                                                                                                                       cAATGAGGTGCGAGATGAACAAGTCCATGGCCTTGCTGGGTGCGTGCGCCCTGCTGCTGGCGGGTGCCGCCAGCGCTGAAACCCTGCGGTTCGCCACCGAGGGTGCCTACCCGCCCTTCAACTATGTAGACGCCGATAACAAGCTGCACGGCTTTGATGTCGa                                                                                                                                >  2:1218129/1‑163 (MQ=255)
                                                                                                                                                                         aTGAGGTGCGAGATGAACAAGTCCATGGCCTTGCTGGGTGCGTGCGCCCTGCTGCTGGCGGGTGCCGCCAGCGCTGAAACCCTGCGGTTCGCCACCGAGGGTGCCTACCCGCCCTTCAACTATGTAGACGCCGATAACAAGCTGCACGGCTTTGATGTCGa                                                                                                                                >  1:1003891/1‑161 (MQ=255)
                                                                                                                                                                         aTGAGGTGCGAGATGAACAAGTCCATGGCCTTGCTGGGTGCGTGCGCCCTGCTGCTGGCGGGTGCCGCCAGCGCTGAAACCCTGCGGTTCGCCACCGAGGGTGCCTACCCGCCCTTCAACTATGTAGACGCCGATAACAAGCTGCACGGCTTTGATGTCGa                                                                                                                                <  2:1003888/161‑1 (MQ=255)
                                                                                                                                                                                  gagaTGAACAAGTCCATGGCCTTGCTGGGTGCGTGCGCCCTGCTGCTGGCGGGTGCCGCCAGCGCTGAAAcc                                                                                                                                                                                                                <  2:487359/72‑1 (MQ=255)
                                                                                                                                                                                  gagaTGAACAAGTCCATGGCCTTGCTGGGTGCGTGCGCCCTGCTGCTGGCGGGTGCCGCCAGCGCTGAAAcc                                                                                                                                                                                                                >  1:487362/1‑72 (MQ=255)
                                                                                                                                                                                        aaCAATTCCATGGCCTTGCTGGGTGCGTGCTCCCGGCTGCTGGCGGGTGCCGCCAGCGCTGAAACCCTGCGGTTCGCCACCGAGGGTGTCTACCCGCCCTTCAACTATGTAGACGCCGATAACAAGCTGCACGGCTTTg                                                                                                                                       <  2:698380/139‑1 (MQ=255)
                                                                                                                                                                                            ggTCCGTGGCCTTGCTGTGTGCGTGCGCCCTGCTGCTGGCGGGTGCCGCCAGCGCTGTAACCCTGCGGTTCGCCACCGAGGGTGCCTACCCGCCCTTCAACTAt                                                                                                                                                                      <  2:67594/103‑1 (MQ=255)
                                                                                                                                                                                                cATGGCCTTGCTGGGTGCGTGCGCCCTGCTGCTGGCGGGTGCCGCCCGCGCTGAAAACCTGCGGTTCGCCACCGAGGGTGCCTACCCGCCCCTCAAcca                                                                                                                                                                       >  1:231943/1‑97 (MQ=255)
                                                                                                                                                                                                     ccTTGCTGGGTGCGTGCGCCCTGCTGCTGGCGGGGgccgc                                                                                                                                                                                                                             >  2:161885/1‑40 (MQ=255)
                                                                                                                                                                                                          cTGGGTGCGTGCGCCCTGCTGCTGGCGGGTGCCGCCAGCGCTGAAACCCTGCGGTTCGCCACCGAGGGTGCCTACCCGCCCTTCAACTATGTAGACGCCGATAACAAGCTGCACGGCTTTGATGTCGACATCACCCATGCCCTGTGCGAACAGATGAAAGTCGAATGCACGCtggtg                                                                               <  1:590515/177‑1 (MQ=255)
                                                                                                                                                                                                           tGGGTGCGTGCGCCCTGCTGCTGGCGGCACCCGCCATCGCTGAAACCCTGCGGTTCGCCACCGAGGGTGCCTACCCGCCCTTCAACTATGTAGACGcc                                                                                                                                                             <  2:675640/98‑1 (MQ=255)
                                                                                                                                                                                                           tGGGTGCGTGCGCCCTGCTGCTGGCGGCACCCGCCAGCGCTGAAACCCTGCGGTTCGCCACCGAGGGTGCCTACCCGCCCTTCAACTATGTAGACGcc                                                                                                                                                             >  1:675643/1‑98 (MQ=255)
                                                                                                                                                                                                            gggTGCGTGCGCCCTGCTGCTGGCGGGTGCCGCCAGCGCTGAAACCCTGCGGTTCGCCACCGAGGGTGCCTACCCGCCCTTCAACTATGTAGACGCCGATAACAAGCTGCACGGCTTTGATGTCGACATCACCCATGCCCTGTGCGa                                                                                                           >  1:1245942/1‑147 (MQ=255)
                                                                                                                                                                                                            gggTGCGTGCGCCCTGCTGCTGGCGGGTGCCGCCAGCGCTGAAACCCTGCGGTTCGCCACCGAGGGTGCCTACCCGCCCTTCAACTATGTAGACGCCGATAACAAGCTGCACGGCTTTGATGTCGACATCACCCATGCCCTGTGCGa                                                                                                           <  2:1245939/147‑1 (MQ=255)
                                                                                                                                                                                                                  gtgcgCCCTGCTGCTGGCGGGTGCCGCCAGCGCTGAAACCCTGCGGTTCGCCACCGAGGGTGCCTACCCGCCCTTCAACTATGTAGCCGCCGATAACAAGCTGCACGGCTTTGATGTCGACATCACCCATGCCCTGTGCGAACCGATg                                                                                                    >  2:54612/1‑148 (MQ=255)
                                                                                                                                                                                                                  gtgcgCCCTGCTGCTGGCGGGTGCCGCCAGCGCTGAAACCCTGCGGTTCGCCACCGAGGGTGCCTACCCGCCCTTCAACTATGTAGACGCCGATAACAAGCTGCACGGCTTTGATGTCGa                                                                                                                                >  2:526234/1‑120 (MQ=255)
                                                                                                                                                                                                                  gtgcgCCCTGCTGCTGGCGGGTGCCGCCAGCGCTGAAACCCTGCGGTTCGCCACCGAGGGTGCCTACCCGCCCTTCAACTATGTAGACGCCGATAACAAGCTGCACGGCTTTGATGTCGa                                                                                                                                <  1:526237/120‑1 (MQ=255)
                                                                                                                                                                                                                  gtgcgCCCTGCTGCTGGCGGGTGCCGCCAGCGCTGAAACCCTGCGGTTCGCCACCGAGGGTGCCTACCCGCCCTTCAACTATGTAGACGCCGATAACAAGCTGCACGGCTTTGATGTCGACATCACCCATGCCCTGTGCGAACAGATGAAAGTCGAATGCACGCTGGTGGCGCAGGACTGGGAAGGCATCATCCCCGCCCTGATGGCGCGCAAGTATGACGCGATAGTGGCATCga            >  1:495633/1‑236 (MQ=255)
                                                                                                                                                                                                                      gcCCTGCTGCTGGCGGGTGCCGCCAGCGCTGAAACCCTGCGGTTCGCCACCGAGGGTGCCTACCCGCCCTTCAACTATGTAGACGCCGATAACAAGCTGCACGGCTTTGATGTCGACATCACCCATGCCCTGTGCGAACAGATGAAAGTCGAATGCACGCTGGTGGCGCAGGACTGGGAAGGCATCATCCCCGCCCTGATGGCGCGCAAGTATGACGCGATAGTGGCATCgatgat        >  1:711112/1‑236 (MQ=255)
                                                                                                                                                                                                                       cccTGCTGCTGGCGGGTGCCGCCAGCGCTGAAACCCTGCGGTTCGCCACCGAGGGTGCCTACCCGCCCTTCAACTATGTAGACGCCGATAACAAGCTGCACGGCt                                                                                                                                          >  2:590512/1‑105 (MQ=255)
                                                                                                                                                                                                                        ccTGCTGCTGGCGGGTGCCGCCAGCGCTGAAACCCTGCGGTTCGCCACCGAGGGTGCCTACCCGCCCTTCAACTATGTAGACGCCGATAACAAGCTGCACGGCTTTGATGTCGACATCACCCATGCCCTGTGCGAACAGATGAAAGTCGAATGCACGCTGGTGGCGCAGGACTGGGAAgg                                                              <  1:117757/180‑1 (MQ=255)
                                                                                                                                                                                                                            ctgctgGCGGGTGCCGCCAGCGCTGAAACCCTGCGGTTCGCCACCGAGGGTGCCTACCCGCCCTTCAACTATGTAGACGCCGATAACAAGCTGCACGGCTTTGATGTCGACATCACCCATGccc                                                                                                                  >  2:602366/1‑124 (MQ=255)
                                                                                                                                                                                                                            ctgctgGCGGGTGCCGCCAGCGCTGAAACCCTGCGGTTCGCCACCGAGGGTGCCTACCCGCCCTTCAACTATGTAGACGCCGATAACAAGCTGCACGGCTTTGATGTCGACATCACCCATGCCCTGTGCGAACAGATGAAAGTCGAATGCACGCTGGTGGCGCAGGACTGGGAAGGCATCATCCCCGCCCTGATGGCGCGCAAGTATGACGCGATAGTGGCATCGATGATCGacac  >  1:633596/1‑236 (MQ=255)
                                                                                                                                                                                                                                  gCGGGTGCCGCCAGCGCTGAAACCCTGCGGTTCGCCACCGAGGGTGCCTACCCGCCCTTCAACTATGTAGACGCCGATAACAAGCTGCACGGCTTTGATGTCGACATCACCCATGCCCTGTGCGAACAGATGAAAGTCGAATGCACGCTGGTGGCGCAGGACTGGGAAGGCATCATCCCCGccc                                                >  1:921777/1‑184 (MQ=255)
                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                
TGGGAATCCACCTCGACGATGCTGTCGTGGGTGCATCTGATGGTCAACGAGGTCTGCAAGAAGCTCGGACCGGATGTTGAAAAACGCTTGGAAGCAACTGCCGCTCTACATAACGAATTCGTCGGCTACACCTCGTGGTCGACGGGCAAACAACAATAGCAAGAGTGCAATGAGGTGCGAGATGAACAAGTCCATGGCCTTGCTGGGTGCGTGCGCCCTGCTGCTGGCGGGTGCCGCCAGCGCTGAAACCCTGCGGTTCGCCACCGAGGGTGCCTACCCGCCCTTCAACTATGTAGACGCCGATAACAAGCTGCACGGCTTTGATGTCGACATCACCCATGCCCTGTGCGAACAGATGAAAGTCGAATGCACGCTGGTGGCGCAGGACTGGGAAGGCATCATCCCCGCCCTGATGGCGCGCAAGTATGACGCGATAGTGGCATCGATGATCGACAC  >  AE015451/4082821‑4083276

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