Predicted mutation
evidence seq id position mutation freq annotation gene description
RA AE015451 135,216 A→G 5.1% Q371R (CAG→CGG)  PP_0129 → GGDEF domain protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*AE015451135,2160AG 5.1% 273.0 / 4.1 79Q371R (CAG→CGG) PP_0129GGDEF domain protein
Reads supporting (aligned to +/- strand):  ref base A (44/31);  new base G (2/2);  total (46/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 = 2.98e-02

AGCCGCCTTACAGCCAGGTTGCCGCGCACATCGAGCAGACCCTGATCGGTCTGCTCGATGACCTCAGCCTGCCCGAGCGACACAAGGCTCAGGCGCAGGAAATGCGCGAGCGGGTTGCTCGCGGCTTGAACTGGTATGAACTGATCCCGGTGCTGGATGACTTGGCAGTACTCATGCTGGCGATCACCGACAGCGGGCAGCACGAATTCGAAACCTACCTGCAGCAGCTCAATGAACGCCTGGAAGGTTTCCAGAGCCACCTGCACGAGGCCAGCGCCGGCCACGCGGACAATAGCACCGCTGCCCGCGAGCTGGATAGCCAACTGCGTGAGCACGTCGACGGCCTGCAGAGCAGCGTGCAAGGTGCCGCCGATGTAGACAGCCTCAAGCACATTCTGGAAAACCGCCTGGAAGGCTTGCTGGTGACGATGGAC  >  AE015451/134991‑135424
                                                                                                                                                                                                                                 |                                                                                                                                                                                                                
aGCCGCCTTACAGCCAGGTTGCCGCGCACATCGAGCAGACCCTGATCGGTCTGCTCGATGACCTCAGCCTGCCCGAGCGACACAAGGCTCAGGCGCAGGAAATGCGCGAGCGGGTTGCTCGCGGCTTGAACTGGTATGAACTGATCCCGGTGCTGGATGACTTGGCAGTACTCATGCTGGCGATCACCGACAGCGGGCAGCACGAATTCGAAACCTACCTGCAGCAGCTCAATGaa                                                                                                                                                                                                        >  1:126464/1‑236 (MQ=255)
     ccTTACAGCCAGGTTGCCGCGCACATCGAGCAGACCCTGATCGGTCTGCTCGATGACCTCAGCCTGCCCGAGCGACACAAGGCTCAGGCGCAGGAAATGCGCGAGCGGGTTGCTCGCGGCTTGAACTGGTATGAACTGATCCCGGTGCTGGATGACTTGGCAGTACTCATGCTGGCGATCACCGACAGCGGGCAGCACGAATTCGAAACCTACCTGCAGCAGCTCAATGAACGc                                                                                                                                                                                                     >  2:400379/1‑234 (MQ=255)
                 gTTGCCGCGCACATCGAGCAGACCCTGATCGGTCTGCTCGATGACCTCAGCCTGCCCGAGCGACACAAGGCTCAGGCGCAGGAAATGCGCGAGCGGGTTGCTCGCGGCTTGAACTGGTATGAACTGATCCCGGTGCTGGATGACTTGGCAGTACTCATGCTGGCGATCACCGACAGCGGGCAGCACGAATTCGAAACCTACCTGCAGCAGCTCAATGAACGCCTGGAAGGTTTCCa                                                                                                                                                                                       <  1:734945/236‑1 (MQ=255)
                                 aGCAGACCCTGATCGGTCTGCTCGATGACCTCAGCCTGCCCGAGCGACACAAGGCTCAGGCGCAGGAAATGCGCGAGCGGGTTGCTCGCGGCTTGAACTGGTATGAACTGATCCCGGTGCTGGATGACTTGGCAGTACTCATGCTGGCGATCACCGACAGCGGGCAGCACGAATTCGAAACCTACCTGcagcag                                                                                                                                                                                                                 >  2:1136198/1‑194 (MQ=255)
                                          tGATCGGTCTGCTCGATGACCTCAGCCTGCCCGAGCGACACAAGGCTCAGGCGCAGGAAATGCGCGAGCGGGTTGCTCGCGGCTTGAACTGGTATGAACTGATCCCGGTGCTGGATGACTTGGCAGTACTCATGCTGGCGATCACCGACAGCGGGCAGCACGAATTCGAAACCTACCTGCAGCAGCTCAATGAACGCCTGGAAGGTTTCCAGAGCCACCTGCACGAGGcc                                                                                                                                                                    >  1:512081/1‑230 (MQ=255)
                                            aTCGGTCTGCTCGATGACCTCAGCCTGCCCGAGCGACACAAGGCTCAGGCGCAGGAAATGCGCGAGCGGGTTGCTCGCGGCTTGAACTGGTATGAACTGATCCCGGTGCTGGATGACTTGGCAGTACTCATGCTGGCGATCACCGACAGCGGGCAGCACGAATTCGAAACCTACCTGCAGCAGCTCAATGAACGCCTGGAAGGTTTCCAGAGCCACCTGCACGAGGCCAGCGCCg                                                                                                                                                             <  1:662097/235‑1 (MQ=255)
                                                   tGCTCGATGACCTCAGCCTGCCCGAGCGACACAAGGCTCAGGCGCAGGAAATGCGCGAGCGGGTTGCTCGCGGCTTGAACTGGTATGAACTGATCCCGGTGCTGGATGACTTGGCAGTACTCATGCTGGCGATCACCGACAGCGGGCAGCACGAATTCGAAACCTACCTGCAGCAGCTCAATGAACGCCTGGAAGGTTTCCAGAGCCACCTGCACGAGGCCAGCGCCGGCCAcgcg                                                                                                                                                     >  1:1018836/1‑236 (MQ=255)
                                                                   ccTGCCCGAGCGACACTAGGCTCAGGCGCAGGAAATGCGCGTGCGGGTTGCTTGCGGCTTGAACTGGTATGAACTGATCCCGGTGCTGGATGACTTGGCAGTACTCATGCTGGCGATCACCGACAGCGGGCAGCACGAATTCGAAACCTACCTGCAGCAGCTCAATGAACGCCTGGAAGGTTTCCAGAGCCATCTGCACGAGGCCAGCGCCGGCCAc                                                                                                                                                        <  1:675213/217‑1 (MQ=255)
                                                                     tGCCCGAGCGACACAAGGCTCAGGCGCAGGAAATGCGCGAGCGGGTTGCTCGCGGCTTGAACTGGTATGAACTGATCCCGGTGCTGGATGACTTGGCAGTACTCATGCTGGCGATCACCGACAGCGGGCAGCACGAATTCGAAACCTACCTGCAGCAGCTCAATGAACGCCTGGAAGGTTTCCAGAGCCACCTGCACGAGGCCAGCGCCGGCCAc                                                                                                                                                        <  2:448313/215‑1 (MQ=255)
                                                                          gAGCGACACAAGGCTCAGGCGCAGGAAATGCGCGAGCGGGTTGCTCGCGGATTGAACTGGTATGAACTGATCCCGGTGCTGGATGACTTGGCAGTACTCATGCTGGCGATCACCGACAGCGGGCAGCACGAATTCGAAACCTACCTGCAGCAGCTCAATGAACGCCTGGAAGGTTTCCAGAGCCACCTGCACGAGGCCAGCGcc                                                                                                                                                              <  1:375266/204‑1 (MQ=255)
                                                                           aGCGACACAAGGCTCAGGCGCAGGAAATGCGCGAGCGGGTTGCTCGCGGCTTGAACTGGTATGAACTGATCCCGGTGCTGGATGACTTGGCAGTACTCATGCTGGCGATCACCGACAGCGGGCAGCACGAATTCGAAACCTACCTGCAGCAGCTCa                                                                                                                                                                                                             <  2:213845/156‑1 (MQ=255)
                                                                           aGCGACACAAGGCTCAGGCGCAGGAAATGCGCGAGCGGGTTGCTCGCGGCTTGAACTGGTATGAACTGATCCCGGTGCTGGATGACTTGGCAGTACTCATGCTGGCGATCACCGACAGCGGGCAGCACGAATTCGAAACCTACCTGCAGCAGCTCa                                                                                                                                                                                                             >  1:213845/1‑156 (MQ=255)
                                                                              gACACAAGGCTCAGGCGCAGGAAATGCGCGAGCGGGTTGCTCGCGGCTTGAACTGGTATGAACTGATCCCGGTGCTGGATGACTTGGCAGTACTCATGCTGGCGATCACCGACAGCGGGCAGCACGAATTCGAAACCTACCTGCAGCAGCTCAATg                                                                                                                                                                                                          >  2:77644/1‑156 (MQ=255)
                                                                                      gCTCAGGCGCAGGAAATGCGCGAGCGGGTTGCTCGCGGCTTGACCTGGTATGAACTGATCCCGGTGCTGGATGACTTGGCAGTACTCATGCTGGCGATCACCGACAGCGGGCAGCACGAATTCGAAACCTACCTGCAGCAGCTCAATGAACGCCTGGAAGGTTTCCAGAGCCACCTGCACGAGGCCAGCGCCGGCCACGCGGACAATAGCACCGCTGCCCGCGAGCTGGATAGc                                                                                                                    <  1:588777/234‑1 (MQ=255)
                                                                                        tCAGGCGCAGGAAATGCGCGAGCGGGTTGCTCGCGGCTTGAACTGGTATGAACTGATCCCGGTGCTGGATGACTTGGCAGTACTCATGCTGGCGATCACCGACAGCGGGCAGCACGAATTCGAAACCTACCTGCAGCAGCTCAATGAACGCCTGGAAGGTTTCCAGAGCCACCTGCACGAGGCCAGCGCCGGCCACGCGGACAACAGCACCGCTGCCCGCGAGCTGGAGAGCCa                                                                                                                  >  1:953071/1‑234 (MQ=255)
                                                                                           ggCGCAGGAAATGCGCGAGCGGGTTGCTCGCGGCTTGAACTGGTATGAACTGATCCCGGTGCTGGATGACTTGTCAGTACTCATGCTGGCGATCACCGACAGCGGGCAGCACGAATTCGAAACCTACCTGCAGCAGCTCAATGAACGCCTGGAAGGTTTCCAGAGCCACCTGCACGAGGCCAGCGCCGGCCACGCGGACa                                                                                                                                                 <  2:387260/200‑1 (MQ=255)
                                                                                             cgcAGGAAATGCGCGAGCGGGTTGCTCGCGGCTTGAACTGGTATGAACTGATCCCGGTGCTGGATGACTTGGCAGTACTCATGCTGGCGATCACCGACAGCGGGCAGCACGAATTCGAAACCTACCTGCAGCAGCTCAATGAACGCCTGGAAGGTTTCCAGAGCCACCTGCACGAGGCCAGCGCCGGCCACGCGGACAATAGCACCGCTGCCCGCGAGCTGGATAGCCAACTGCGt                                                                                                           >  1:1075202/1‑236 (MQ=255)
                                                                                             cgcAGGAAATGCGCGAGCGGGTTGCTCGCGGCTTGAACTGGTATGAACTGATCCCGGTGCTGGATGACTTGGCAGTACTCATGCTGGCGATCACCGACAGCGGGCAGCACGAATTCGAAAACTACCTGCAGCAGATCAATGAACGCCTGGAAGGGTTCCAGAGCCACCTGCACGAgg                                                                                                                                                                      >  2:129140/1‑177 (MQ=255)
                                                                                                  gAAATGCGCGAGCGGGTTGCTCGCGGCTTGAACTGGTATGAACTGATCCCGGTGCTGGATGACTTGGCAGTACTCATGCTGGCGATCACCGACAGCGGGCAGCACGAATTCGAAACCTACCTGCAGCAGCTCAATGAACGCCTGGAAGGTTTCCAGAGCCACCTGCACGAGGCCAGCGCCGGCCACGCGGACAATAGCACCGCTGcc                                                                                                                                   >  2:670418/1‑207 (MQ=255)
                                                                                                  gAAATGCGCGAGCGGGTTGCTCGCGGCTTGAACTGGTATGAACTGATCCCGGTGCTGGATGACTTGGCAGTACTCATGCTGGCGATCACCGACAGCGGGCAGCACGAATTCGAAACCTACCTGCAGCAGCTCAATGAACGCCTGGAAGGTTTCCAGAGCCACCTGCACGAGGCCAGCGCCGGCCACGCGGACAATAGCACCGCTGCCCGCGAGCTGGATAGCCAACTGCGTGAg                                                                                                        >  1:97812/1‑234 (MQ=255)
                                                                                                        cgcgAGCGGGTTGCTCGCGGCTTGAACTGGTATGAACTGATCCCGGTGCTGGATGACTTGGCAGTACTCATGCTGGCGATCACCGACAGCGGGCAGCACGAATTCGAAACCTACCTGCAGCAGCTCAATg                                                                                                                                                                                                          >  2:704664/1‑130 (MQ=255)
                                                                                                         gcgAGCGGGTTGCTCGCGGCTTGCACTGGTATGAACTGATCCCGGTGCTGGATGACTTGGCAGTACCCATGCTGGCGATCACCGACAGCGGGCAGCACGAATTCGAAACCTACCTGCAGCAGCTCAATGAACGCCTGGAAGGTTTCCAGAGCCACCTGCACGAGGCCAGCGCCGGCCAc                                                                                                                                                        <  2:368328/179‑1 (MQ=255)
                                                                                                           gAGCGGGTTGCTCGCGGCTTGAACTGGTATGAACTGATCCCGGTGCTGGATGACTTGGCAGTACTCATGCTGGCGATCACCGACAGCGGGCAGCACGAATTCGAAACCTACCTGCAGCAGCTCAATGAACGCCTGGAAGGTTTCCAGAGCCACCTGCACGAGGCCAGCGCCGGCCACGCGGACAATAGCACCGCTGCCCGCGAGCTGGATAGCCAACTGCGTGAGCACGTCGACg                                                                                              >  1:545037/1‑235 (MQ=255)
                                                                                                                  ttGCTCGCGGCTTGAACTGGTATGAACTGATCCCGGTGCTGGATGACTTGGCAGTACTCATGCTGGCGATCACCGACAGCGGGCAGCACGAATTCGAAACCTACCTGCAGCAGCTCAAt                                                                                                                                                                                                           <  1:245713/119‑1 (MQ=255)
                                                                                                                          ggCTTGAACTGGTATGAACTGATCCCGGTGCTGGATGACTTGGCAGTACTCATGCTGGCGATCACCGACAGCGGGCAGCACGAATTCGAAACCTACCTGCAGCAGCTCAATGAACGCCTGGAAGGTTTCCAGAGc                                                                                                                                                                                   >  2:1100371/1‑135 (MQ=255)
                                                                                                                          ggCTTGAACTGGTATGAACTGATCCCGGTGCTGGATGACTTGGCAGTACTCATGCTGGCGATCACCGACAGCGGGCAGCACGAATTCGAAACCTACCTGCAGCAGCTCAATGAACGCCTGGAAGGTTTCCAGAGc                                                                                                                                                                                   <  1:1100374/135‑1 (MQ=255)
                                                                                                                                        ttaaCTGATCCCGGTGCTGGATGACTTGGCAGTACTCATGCTGGCGATCACCGACAGCGGGCAGCACGAATTCGAAACCTACCTGCAGCAGCTCAATGAACGCCTGGAAGGTTTCCAGAGCCACCTGCACGAGGCCAGCGCCGGCCACGCGGACAATAGCACCGCTGCCCGCGAGCTGGATAGCCAACTGCGTGAGCACGTCGAc                                                                                               <  1:848682/203‑1 (MQ=255)
                                                                                                                                         gAACTGATCCCGGTGCTGGATGACTTGGCAGTACTCATGCTGGCGATCACCGACAGCGGGCAGCACGAATTCGAAACCTACCTGCAGCAGCTCAATGAACGCCTGGAAGGTTTCCAGAGCCACCTGCACGAGGCCAGCGCCGGCCACGCGGACAATAGCACCGCTGCCCGCGAGCTGGa                                                                                                                        >  2:614781/1‑179 (MQ=255)
                                                                                                                                         gAACTGATCCCGGTGCTGGATGACTTGGCAGTACTCATGCTGGCGATCACCGACAGCGGGCAGCACGAATTCGAAACCTACCTGCAGCAGCTCAATGAACGCCTGGAAGGTTTCCAGAGCCACCTGCACGAGGCCAGCGCCGGCCACGCGGACAATAGCACCGCTGCCCGCGAGCTGGATAGCCAACTGCGTGAGCACGTc                                                                                                  >  2:778763/1‑201 (MQ=255)
                                                                                                                                             tGATCCCGGTGCTGGATGACTTGGCAGTACTCATGCTGGCGATCACCGACAGCGGGCAGCACGAATTCGAAACCTACCTGCAGCAGCTCAATGAACGCCTGGAAGGTTTCCAGAGCCACCTGCACGAGGCCAGCGCCGGCCACGCGGACAATAGCACCGCTGCCCGCGAGCTGGATAGCCAACTGCGTGAGCACGTCGACGGCCTGCAGAGCAGCGTGCAAGGTGCCGCCGATGTa                                                           >  1:458168/1‑236 (MQ=255)
                                                                                                                                                 cccGGTGCTGGATGACTTGGCAGTACTCATGCTGGCGATCACCGACAGCGGGCAGCCCGAATTCGAAACCTACCTGCAGCAGCTCAATGAACGCCTGGAAGGTTTCCAGAGCCACCTGCACGAGGCCAGCGCCGGCCACGCGGACAATAGCACCGCTGCCCGCGAGCTGGATAGCCAACTGCGTGAGCACGTCGACGGCCTGCAGAGCAGCGTGCAAGGTGCCGCCGATGt                                                            <  1:221134/231‑1 (MQ=255)
                                                                                                                                                 cccGGTGCTGGATGACTTGGCAGTACTCATGCTGGCGATCACCGACAGCGGGCAGCACGAATTCGAAACCTACCTGCAGCAGCTCAATGAACGCCTGGAAGGTTTCCAGAGCCACCTGCACGAGGCCAGCGCCGGCCACGCGGACAATAGCACCGCTGCCCGCGAGCTGGATAGCCAACTGCGTGAGCACGTCGACGGCCTGCAGAGCAGCGTGCAAGGTGCCGCCGATGTAGAc                                                        >  1:995715/1‑235 (MQ=255)
                                                                                                                                                  ccGGTGCTGGATGACTTGGCAGTACTCATGCTGGCGATCACCGACAGCGGGCAGCACGAATTCGAGACCTACCTGCAGCAGCTCAATGAACGCCTGGAAGGTTTCCAGAGCCACCTGCACGAGGCCAGCGCCGGCCACGCGGACAATAGCACCGCTGCCCGCGAGCTGGATAGCCAACTGCGTGAGCACGTCg                                                                                                 <  1:1284507/193‑1 (MQ=255)
                                                                                                                                                  ccGGTGCTGGATGACTTGGCAGTACTCATGCTGGCGATCACCGACAGCGGGCAGCACGAATTCGAAACCTACCTGCAGCAGCTCAATGAACGCCTGGAAGGTTTCCAGAGCCACCTGCACGAGGCCAg                                                                                                                                                                  >  2:633984/1‑128 (MQ=255)
                                                                                                                                                  ccGGTGCTGGATGACTTGGCAGTACTCATGCTGGCGATCACCGACAGCGGGCAGCACGAATTCGAAACCTACCTGCAGCAGCTCAATGAACGCCTGGAAGGTTTCCAGAGCCACCTGCACGAGGCCAGCGCCg                                                                                                                                                             >  2:1125328/1‑133 (MQ=255)
                                                                                                                                                  ccGGTGCTGGATGACTTGGCAGTACTCATGCTGGCGATCACCGACAGCGGGCAGCACGAATTCGAAACCTACCTGCAGCAGCTCAATGAACGCCTGGAAGGTTTCCAGAGCCACCTGCACGAGGCCAGCGCCg                                                                                                                                                             <  1:1125331/133‑1 (MQ=255)
                                                                                                                                                  ccGGTGCTGGATGACTTGGCAGTACTCATGCTGGCGATCACCGACAGCGGGCAGCACGAATTCGAAACCTACCTGCAGCAGCTCAATGAACGCCTGGAAGGTTTCCAGAGCCACCTGCACGAGGCCAGCGCCGGcc                                                                                                                                                          >  2:1269543/1‑136 (MQ=255)
                                                                                                                                                  ccGGTGCTGGATGACTTGGCAGTACTCATGCTGGCGATCACCGACAGCGGGCAGCACGAATTCGAAACCTACCTGCAGCAGCTCAATGAACGCCTGGAAGGTTTCCAGAGCCACCTGCACGAGGCCAGCGCCGGcc                                                                                                                                                          <  1:1269546/136‑1 (MQ=255)
                                                                                                                                                  ccGGTGCTGGATGACTTGGCAGTACTCATGCTGGCGATCACCGACAGCGGGCAGCACGAATTCGAAACCTACCTGCAGCAGCTCAATGAACGCCTGGAAGGTTTCCAGAGCCACCTGCACGAGGCCAGCGCCGGCCACGCGGACAATAGCAcc                                                                                                                                         >  2:411305/1‑153 (MQ=255)
                                                                                                                                                  ccGGTGCTGGATGACTTGGCAGTACTCATGCTGGCGATCACCGACAGCGGGCAGCACGAATTCGAAACCTACCTGCAGCAGCTCAATGAACGCCTGGAAGGTTTCCAGAGCCACCTGCACGAGGCCAGCGCCGGCCACGCGGACAATAGCACCGCTGCCCGCGAGcgg                                                                                                                          >  2:5878/1‑166 (MQ=255)
                                                                                                                                                  ccGGTGCTGGATGACTTGGCAGTACTCATGCTGGCGATCACCGACAGCGGGCAGCACGAATTCGAAACCTACCTGCAGCAGCTCAATGAACGCCTGGAAGGTTTCCAGAGCCACCTGCACGAGGCCAGCGCCGGCCACGCGGACAATAGCACCGCTGCCCGCGAGCTGGATAg                                                                                                                     >  2:627848/1‑173 (MQ=255)
                                                                                                                                                  ccGGTGCTGGATGACTTGGCAGTACTCATGCTGGCGATCACCGACAGCGGGCAGCACGAATTCGAAACCTACCTGCAGCAGCTCAATGAACGCCTGGAAGGTTTCCAGAGCCACCTGCACGAGGCCAGCGCCGGCCACGCGGACAATAGCACCGCTGCCCGCGAGCTGGATAGCCAACTGCGTGAGCACGTCGACGGCCTGCAGAGCAGCGTGCAAGGTGccgccg                                                                >  1:424228/1‑226 (MQ=255)
                                                                                                                                                  ccGGTGCTGGATGACTTGGCAGTACTCATGCTGGCGATCACCGACAGCGGGCAGCACGAATTCGAAACCTACCTGCAGCAGCTCAATGAACGCCTGGAAGGTTTCCAGAGCCACCTGCACGAGGCCAGCGCCGGCCACGCGGACAATAGCACCGCTGCCCGCGAGCTGGATAGCCAACTGCGTGAGCACGTCGACGGCCTGCAGAGCAGCGTGCAAGGTGCCGCCGATGTAGAc                                                        >  1:753823/1‑234 (MQ=255)
                                                                                                                                                  ccGGTGCTGGATGACTTGGCAGTACTCATGCTGGCGATCACCGACAGCGGGCAGCACGAATTCGAAACCTACCTGCAGCAGCTCAATGAACGCCTGGAAGGTTTCCAGAGCCACCTGCACGAGGCCAGCGCCGGCCACGCGGACAATAGCACCGCTGCCCGCGAGCTGGATAGCCAACTGCGTGAGCACGTCGACGGCCTGCAGAGCAGCGTGCAAGGTGCCGCCGATGTAGACa                                                       >  1:1028778/1‑235 (MQ=255)
                                                                                                                                                  ccGGTGCTGGATGACTTGGCAGTACTCATGCTGGCGATCACCGACAGCGGGCAGCACGAATTCGAAACCTACCTGCAGCAGCTCAATGAACGCCTGGAAGGTTTCCAGAGCCACCTGCACGAGGCCAGCGCCGGCCACGCGGACAATAGCACCGCTGCCCGCGAGCTGGATAGCCAACTGCGTGAGCACGTCGACGGCCTGCAGAGCAGCGTGCAAGGTGCCGCCGATGTAGACAg                                                      >  1:627540/1‑236 (MQ=255)
                                                                                                                                                  ccGGTGCTGGATGACTTGGCAGTACTCATGCTGGCGATCACCGACAGCGGGCAGCACGAATTCGAAACCTACCTGCAGCAGCTCAATGAACGCCTGGAAGGTTTCCAGAGCCACCTGCACGAGGCCAGCGCCGGCCACGCGGACAATAGCACCGCTGCCCGCGAGCTGGAGAGCCAACTGCGTGAGCACGTCGACGGCCTGCAGAGCAGCGTGCAAGGTGCCGCCGAt                                                              >  1:372994/1‑228 (MQ=255)
                                                                                                                                                       gCTGGATGACTTGGCAGTACTCATGCTGGCGATCACCGACAGCGGGCAGCACGAATTCGAAACCTACCTGCAGCAGCTCAATGAACGCCTGGAAGGTTTCCAGAGCCACCTGCACGAGGCCAGCGCCGGCCACGCGGACAATAGCACCGCTGCCCGCGAGCTGGATAGCCAACTGCGTGAGCACGTCGAc                                                                                               >  2:414093/1‑190 (MQ=255)
                                                                                                                                                        cTGGATGACTTGGCAGTACTCATGCTGGCGATCACCGACAGCGGGCAGCACGAATTCGAAACCTACCTGCAGCAGCTCAATGAACGCCTGGAAGGTTTCCAGAGCCACCTGCACGAGGCCAGCGCCGGCCACGCGGACAATAGCACCGCTGCCCGCGAGCTGGATAGCCAACTGCGTGAGCACGTCGACGGCCTGCa                                                                                       >  2:1151536/1‑197 (MQ=255)
                                                                                                                                                              gACTTGGCAGTACTCATGCTGGCGATCACCGACAGCGGGCAGCACGAATTCGAAACCTACCTGCAGCAGCTCAATGATCGCCTGGAAGGTTTCCAGAGCCACCTGCACGAGGCCAGCGCCGGCCACGCGGACAATAGCACCGCTGCCCGCGAGCTGGATAGCCAACTGCGt                                                                                                           <  1:848344/171‑1 (MQ=255)
                                                                                                                                                                   ggCAGTACTCATGCTGGCGATCACCGACAGCGGGCAGCACGAATTCGAAACCTACCTGCAGCGGCTCAATGAACGCCTGGAAGGTTTCCAGAGCCACCTGCACGAGGCCAGCGCCGGcc                                                                                                                                                          <  2:529914/119‑1 (MQ=255)
                                                                                                                                                                   ggCAGTACTCATGCTGGCGATCACCGACAGCGGGCAGCACGAATTCGAAACCTACCTGCAGCGGCTCAATGAACGCCTGGAAGGTTTCCAGAGCCACCTGCACGAGGCCAGCGCCGGcc                                                                                                                                                          >  1:529917/1‑119 (MQ=255)
                                                                                                                                                                   ggCAGTACTCATGCTGGCGATCACCGACAGCGGGCAGCACGAATTCGAAACCTACCTGCAGCAGCTCAATGAACGCCTGGAAGGTTTCCAGAGCCACCTg                                                                                                                                                                             <  2:906721/100‑1 (MQ=255)
                                                                                                                                                                   ggCAGTACTCATGCTGGCGATCACCGACAGCGGGCAGCACGAATTCGAAACCTACCTGCAGCAGCTCAATGAACGCCTGGAAGGTTTCCAGAGCCACCTg                                                                                                                                                                             >  1:906724/1‑100 (MQ=255)
                                                                                                                                                                   ggCAGTACTCATGCTGGCGATCACCGACAGCGGGCAGCACGAATTCGAAACCTACCTGCAGCAGCTCAATGAACGCCTGGAAGGTTTCCAGAGCCACCTGCACGAGGCCAGCGCCGGCCACGCGGACAATTGCACCGCTGCCCGCGAGc                                                                                                                            >  2:854034/1‑149 (MQ=255)
                                                                                                                                                                      aGTACTCATGCTGGCGATCACCGACAGCGGGCAGCACGAATTCGAAACCTACCTGCAGCAGCTCAATGAACGCCTGGAAGGTTTCCAGAGCCACCTGCACg                                                                                                                                                                         >  2:1150761/1‑101 (MQ=255)
                                                                                                                                                                      aGTACTCATGCTGGCGATCACCGACAGCGGGCAGCACGAATTCGAAACCTACCTGCAGCAGCTCAATGAACGCCTGGAAGGTTTCCAGAGCCACCTGCACGAGGCCAGCGCCGGCCACGCGGACAATAGCACCGCTGCCCGCGAGCTGGATAGCCAACTGCGTGAGCACGTCGACGGCCTGCAGAGCAGCGTGCAAGGTGCCGCCGATGTAGACAGCCTCAAGCACATTCTGGaa                                   <  1:326197/235‑1 (MQ=255)
                                                                                                                                                                             aagCTGGCGATCACCGACAGCGGGCAGCACGAATTCGAAACCTACCTGCAGCAGCTCAATGAACGCCTGGAAGGTTTCCAGAGCCACCTGCACGAGGCCAGCGCCGGCCAAGCGGACAATAGCACCGCTGCCCGCGAg                                                                                                                             >  2:775479/3‑138 (MQ=255)
                                                                                                                                                                             aTGCTGGCGATCACCGACAGCGGGCAGCACGAATTCGAAACCTACCTGCAGCAGCTCAATGAACGCCTGGAAGGTTTCCAGAGCCACCTGCACGAGGCCAGCGcc                                                                                                                                                              >  2:743263/1‑105 (MQ=255)
                                                                                                                                                                             aTGCTGGCGATCACCGACAGCGGGCAGCACGAATTCGAAACCTACCTGCAGCAGCTCAATGAACGCCTGGAAGGTTTCCAGAGCCACCTGCACGAGGCCAGCGCCGGCCACGCGGACAATAGCACCGCTGCCCGCGAGCTGGATAGCCAACTGCGTGAGCACGTCGACGGCCTGCAGAGCAGCGTGCAAGGTGCCGCCGATGTAGACAGCCTCAAGCACATTCTGGAAAACCGc                             >  2:649094/1‑234 (MQ=255)
                                                                                                                                                                                cTGGCGATCACCGACAGCGGGCAGCACGAATTCGAAACCTACCTGCAGCAGCTCAATGAACGCCTGGAAGGTTTCCAGAGCCACCTGCACGAGGCCAGCGCCGGCCACGCGGACAATAGCACCGCTGCCCGCGAGCTGGATAGCCAACTGCGTGAGCACGTCGACGGCCTGCAGAGCAGCGTGCAAGGTGCCGCCGATGt                                                            <  2:1018833/200‑1 (MQ=255)
                                                                                                                                                                                 tGGCGATCACCGACAGCGGGAAGCACGGATTCGAAACCTCCCTGCAGCAGCTCAATGAACGCCTGGAAGGTTTCCAGAGCCACCTg                                                                                                                                                                             <  2:247154/86‑1 (MQ=255)
                                                                                                                                                                                 tGGCGACCACCGACAGCGGGCAGCACGAATTCGAAACCTACCTGCAGCAGCTCAATGAACGCCTGGAAGGTTTCCAGAGCCACCTGCACGAGGCCAGCGCCGGCCACGCGGACAATAGCACCGCTGCCCGCGAGCTGGATAGCCAACTGCGTGAGCACg                                                                                                    <  2:522071/159‑1 (MQ=255)
                                                                                                                                                                                   gCGACCACCGACAGCGGGCAGCACGAATTCGAAACCTACCTGCAGCAGCTCAATGAACGCCTGGAAGGTTTCCAGAGCCACCTGCACGAGGCCAGCGCCGGCCAc                                                                                                                                                        <  2:464372/105‑1 (MQ=255)
                                                                                                                                                                                   gCGACCACCGACAGCGGGCAGCACGAATTCGAAAACTACCTGCAGCAGCTCAATGAACGCCTGGAAGGTTTCCAGAGCCACCTGCACGAGGCCAGCGCCGGCCACGCGGACAATAGCACCGCTGCCCGCGAGCTGGATAGCCAACTGCGTGAGCACGTCg                                                                                                 <  2:553768/160‑1 (MQ=255)
                                                                                                                                                                                     gATCACCGACAGCGGGCAGCACGAATTCGAAACCTACCTGCAGCGGCTCAATGAACGCCTGGAAGGTTTCCAGAGCCACCTGCACGAGGCCAGCGCCGGCCACGCGGACAATAGCACCGCTGCCCGCGAGCTGGATAGCCAACTGCGTGAGCACGTCGACGGCCTGCAGAGCAGCGTGCAAGGTGCCGCCGATGTAGACAGCCTCa                                                 <  2:1061309/206‑1 (MQ=255)
                                                                                                                                                                                     gATCACCGACAGCGGGCAGCACGAATTCGAAACCTACCTGCAGCGGCTCAATGAACGCCTGGAAGGTTTCCAGAGCCACCTGCACGAGGCCAGCGCCGGCCACGCGGACAATAGCACCGCTGCCCGCGAGCTGGATAGCCAACTGCGTGAGCACGTCGACGGCCTGCAGAGCAGCGTGCAAGGTGCCGCCGATGTAGACAGCCTCa                                                 >  1:1061312/1‑206 (MQ=255)
                                                                                                                                                                                     gATCACCGACAGCGGGCAGCACGAATTCGAAACCTACCTGCAGCAGCTCAATGAACGCCTGGAAGGTTTCCAGAGCCACCTGCACGAGGCCAGCGCCGGCCACGCGGACAATAGCACCGCTGCCCGCGAGCTGGATAGCCAACTGCGTGAGCACGTCGACGGCCTGCAGAGCAGCGTGCAAGGTGCCGCCGATGTAGACAGCCTCAAGCACATTCTGGAAAACCGCCTGGAAg                      >  2:737823/1‑233 (MQ=255)
                                                                                                                                                                                      aTCACCGACAGCGGGCAGCACGAATTCGAAACCTACCTGCAGCAGCTCAATGAACGCCTGGAAGGTTTCCAGAGCCACCTGCACGAGGCCAGCGCCGGCCACGCGGACAATAGCACCGCTGCCCGCGAGCTGGATAGCCAACTGc                                                                                                             >  2:916590/1‑145 (MQ=255)
                                                                                                                                                                                         aCCGACAGCGGGCAGCACGAATTCGAAACCTACCTGCAGCAGCTCAATGAACGCCTGGAAGGTTTCCAGAGCCACCTGCACGAGGCCAGCGCCGGCCACGCGGACAATAGCACCGCTGCCCGCGAGCTGGATAGCCAACTGCGt                                                                                                           <  1:496862/144‑1 (MQ=255)
                                                                                                                                                                                          ccGACAGCGGGCAGCACGAATTCGAAACCTACCTGCAGCAGCTCAATGAACGCCTGGAAGGTTTCCAGAGCCACCTGCACGAGGCCAGCGCCGGCCACGCGGACAATAGCACCGCTGCCCGCGAGCTGGATAGCCAACTGCGTGAGCACGTCGACGGCCTGCAGAGCAGCGTGCAAGGTGCCGCCGATGTAGACAGCCTCAAGCACATTCTGGAAAACCGCCTGGAAGGCTTGCt               <  2:1028775/235‑1 (MQ=255)
                                                                                                                                                                                              cAGCGGGCAGCACGAATTCGAAACCTACCTGCAGCAGCTCAATGAACGCCTGGAAGGTTTCCAGAGCCACCTGCACGAGGCCAGCGCCGGCCACGCGGACAATAGCACCGCTGCCCGCGAGCTGGATAGCCAACTGCGTGAGCACGTCGACGGCCTGCAGAGCAGCGTGCAAGGTGCCGCCGATGTAGACAGCCTCAAGCACATTCTGGAAAACCGCCTGGAAGGCTTGCTGGTGa          <  1:418280/236‑1 (MQ=255)
                                                                                                                                                                                                  gggCAGCACGAATTCGAAACCTACCTGCAGCAGCTCAATGAACGCCTGGAAGGTTTCCAGAGCCACCTGCACGAGGCCAGCTCCGGCCACGCGGACAATAGCACCGCTGCCCGCGAGCTGGATAGCCAACTGCGTGAGCACGTCGACGGcc                                                                                           >  1:1167686/1‑151 (MQ=255)
                                                                                                                                                                                                   ggCAGCACGAATTCGAAACCTACCTGCAGCAGCTCAATGAACGCCTGGAAGGTTTCCAGAGCCACCTGCACGAGGCCAGCGCCGGCCACGCGGACAATAGCACCg                                                                                                                                        >  1:221518/1‑105 (MQ=255)
                                                                                                                                                                                                    gcagcaCGAATTCGAAACCTACCTGCAGCAGCTCAATGAACGCCTGGAAGGTTTCCAGAGCCACCTGCACGAGGCCAGCGCCGGCCACGCGGACAATAGCACCGCTGCCCGCGAGc                                                                                                                            <  2:514442/116‑1 (MQ=255)
                                                                                                                                                                                                    gcagcaCGAATTCGAAACCTACCTGCAGCAGCTCAATGAACGCCTGGAAGGTTTCCAGAGCCACCTGCACGAGGCCAGCGCCGGCCACGCGGACAATAGCACCGCTGCCCGCGAGCTGGATAGCCAACTGCGTGAGCACGTCGACGGCCTGCAGAGCAGCGTGCAAGGTGCCGCCGATGTAGACAGCCTCAAGCACATTCTGGAAAACCGCCTGGAAGGCTTGCTGGTGACGATg     >  1:944315/1‑235 (MQ=255)
                                                                                                                                                                                                     cagcaCGAATTCGAAACCTACCTGCAGCAGCTCAATGAACGCCTGGAAGGTTTCCAGAGCCACCTGCACGAGGCCAGCGCCGGCCACGCGGACAATAGCACCGCTGCCCGCGAGCTGGATAGCCAACTGCGTGAGCACGTCGACGGCCTGCAGAGCAGCGTGCAAGGTGCCGCCGATGTAGACAGCCTCAAGCACATTCTGGAAAACCGCCTGGAAGGCTTGCTGGTGACGATg     <  2:974998/234‑1 (MQ=255)
                                                                                                                                                                                                      agcaCGAATTCGAAACCTACCTGCAGCAGCTCAATGAACGCCTGGAAGGTTTCCAGAGCCACCTGCACGAGGCCAGCGCCGGCCACGCGGACAATAGCACCGCTGCCCGCGAGCTGGATAGCCAACTGCGTGAGCACGTCGACGGCCTGCAGAGCAGCGTGCAAGGTGCCGCCGATGTAGACAGCCTCAAGCACATTCTGGAAAACCGCCTGGAAGGCTTGCTGGTGACGATGGAc  >  1:443147/1‑236 (MQ=255)
                                                                                                                                                                                                                cGAAACCTACCTGCAGCAGCTCAATGAACGCCTGGAAGGTTTCCAGAGCCACCTGCACGAGGCCAGCGCCGGCCACGCGGACAATAGCACCGCTGCCCGCGAGCTGGATAg                                                                                                                     >  2:398948/1‑111 (MQ=255)
                                                                                                                                                                                                                cGAAACCTACCTGCAGCAGCTCAATGAACGCCTGGAAGGTTTCCAGAGCCACCTGCACGAGGCCAGCGCCGGCCACGCGGACAATAGCACCGCTGCCCGCGAGCTGGATAg                                                                                                                     <  1:398951/111‑1 (MQ=255)
                                                                                                                                                                                                                     cctacctGCAGCAGCGCAATGAACGCCTGGAAGGTTTCCAGAGCCCCCCGCACGAGGCCAGCGCCGGCCACGCGGACAATAGCACCGCTGCCCGCGAGCTGGATAGCCAACTGCGTGAGCACGTCGACGGCCTGCAGAGCAGCGTGCAAGGTg                                                                      <  2:280097/153‑1 (MQ=255)
                                                                                                                                                                                                                                 |                                                                                                                                                                                                                
AGCCGCCTTACAGCCAGGTTGCCGCGCACATCGAGCAGACCCTGATCGGTCTGCTCGATGACCTCAGCCTGCCCGAGCGACACAAGGCTCAGGCGCAGGAAATGCGCGAGCGGGTTGCTCGCGGCTTGAACTGGTATGAACTGATCCCGGTGCTGGATGACTTGGCAGTACTCATGCTGGCGATCACCGACAGCGGGCAGCACGAATTCGAAACCTACCTGCAGCAGCTCAATGAACGCCTGGAAGGTTTCCAGAGCCACCTGCACGAGGCCAGCGCCGGCCACGCGGACAATAGCACCGCTGCCCGCGAGCTGGATAGCCAACTGCGTGAGCACGTCGACGGCCTGCAGAGCAGCGTGCAAGGTGCCGCCGATGTAGACAGCCTCAAGCACATTCTGGAAAACCGCCTGGAAGGCTTGCTGGTGACGATGGAC  >  AE015451/134991‑135424

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