Predicted mutation
evidence seq id position mutation freq annotation gene description
RA AE015451 1,617,614 G→A 7.4% G44G (GGC→GGT PP_1417 ← putative Tricarboxylate transport protein TctB

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*AE0154511,617,6140GA 7.4% 165.6 / 2.2 52G44G (GGC→GGTPP_1417putative Tricarboxylate transport protein TctB
Reads supporting (aligned to +/- strand):  ref base G (24/24);  new base A (2/2);  total (26/26)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.30e-01

CGTACAGCCGAGCCATGGGCAGGCCGACCAGGATGGCGCTGAGGATGAAGCCGAGCGCTTCGAAAGTGGCGGCGAAGACGATCAGCAGGCCGACACAGGCGGTGATCTTGATCAGGGTCTCGCGGTCCAGTTCCGGTTCGTCATCCTTGCGCACGATGGGCGTTGGGCGGAACACCAGGTACAGCAGGCCCAGCCCCATCAGGCCGAGCACCAGCAGCGGGTAAGCGCGCGGGCCCACCGGTTCATAGGAGAATGCAGCCTGATAGGGCCAGGCCATCACGGCCAGGGCGGCACACATCGCCAGCAGGGCCAAGGCGAAAATGCGTTGAAGGATCATTGGGGAATCCTCGTTGGGATTGCTCAAGTGAAATGCAAATCCTGTGGGAGCGGCCTTGCGTCGCGATCGGGCGCGAAGCTGCCGCTAAACCTGCCGATGCGATCTGTACGAAGCATCGTATCGGCTGTGCC  >  AE015451/1617382‑1617849
                                                                                                                                                                                                                                         |                                                                                                                                                                                                                                           
cGTACAGCCGAGCCATGGGCAGGCCGACCAGGATGGCGCTGAGGATGAAGCCGAGCGCTTCGAAAGTGGCGGCGAAGACGATCAGCAGGCCGACACAGGCGGTGATCTTGATCAGGGTCTCGCGGTCCAGTTCCGGTTCGTCATCCTTGCGCACGATGGGCGTTGGGCGGAACACCAGGTACAGCAGGCCCAGCCCCATCAGGCCGAGCACCAGCAGCGGGTAAGCGCGCGGGcc                                                                                                                                                                                                                                           <  1:733680/235‑1 (MQ=255)
      gCCGAGCCATGGGCAGGCCGACCAGGATGGCGCTGAGGATGAAGCCGAGCGCTTCGAAAGTGGCGGCGAAGACGATCAGCAGGCCGACACAGGCGGTGATCTTGATCAGGGTCTCGCGGTCCAGTTCCGGTTCGTCATCCTTGCGCACGATGGGCGTTGGGCGGAACACCAGGTACAGCAGGCCCAGCCCCATCAGGCCGAGCACCAGCAGCGGGTAAGCGCGCGGGcc                                                                                                                                                                                                                                           <  1:155839/229‑1 (MQ=255)
                         gACCAGGATGGCGCTGAGGATGAAGCCGAGCGCTTCGAAAGTGGCGGCGAAGACGATCAGCAGGCCGACACAGGCGGTGATCTTGATCAGGGTCTCGCGGTCCAGTTCCGGTTCGTCATCCTTGCGCACGATGGGCGTTGGGCGGAACACCAGGTACAGCAGGCCCAGCCCCATCAGGCCGAGCACCAGCAGCGGGTAAGCGCGCGGGCCCACCGGTTCATAGGAGAATGCAGCCt                                                                                                                                                                                                                 >  2:446232/1‑236 (MQ=255)
                                       tGAGGATGAAGCCGAGCGCTTCGAAAGTGGCGGCGAAGACGATCAGCAGGCCGACACAGGCGGTGATCTTGATCAGGGTCTCGCGGTCCAGTTCCGGTTCGTCATCCTTGCGCACGATGGGCGTTGGGCGGAACACCAGGTACAGCAGGCCCAGCCCCATCAGGCCGAGCACCAGCAGCGGGTAAGCGCGCGGGCCCACCGGTTCATAGGAGCATGCAGCCTGAt                                                                                                                                                                                                              >  1:728825/1‑225 (MQ=255)
                                           gATGAAGCCGAGCGCTTCGAAAGTGGCGGCGAAGACGATCAGCAGGCCGACACAGGCGGTGATCTTGATCAGGGTCTCGCGGTCCAGTTCCGGTTCGTCATCCTTGCGCACGATGGGCGTTGGGCGGAACACCAGGTACAGCAGGCCCAGCCCCATCAGGCCGAGCACCAGCAGCGGGTAAGCGCGCGGGCCCACCGGTTCAt                                                                                                                                                                                                                                <  2:312325/203‑1 (MQ=255)
                                                  ccGAGCGCTTCGAAAGTGGCGGCGAAGACGATCAGCAGGCCGACACAGGCGGTGATCTTGATCAGGGTCTCGCGGTCCAGTTCCGGTTCGTCATCCTTGCGCACGATGGGCGTTGGGCGGAACACCAGGTACAGCAGGCCCAGCCCCATCAGGCCGAGCACCAGCAGCGGGTAAGCGCGCGGGCCCa                                                                                                                                                                                                                                         <  2:109838/187‑1 (MQ=255)
                                                        gcTTCGAAAGTGGCGGCGAAGACGATCAGCAGGCCGACACAGGCGGTGATCTTGATCAGGGTCTCGCGGTCCAGTTCCGGTTCGTCATCCTTGCGCACGATGGGCGTTGGGCGGAACACCAGGTACAGCAGGCCCAGCCCCATCAGGCCGAGCACCAGCAGCGGGTAAGCGCGCGGGCCCACCGGTTCATAGGAGAATGCAGCCTGATAGGGCCAGGCCATCACgg                                                                                                                                                                                            >  2:522861/1‑226 (MQ=255)
                                                        gcTTCGAAAGTGGCGGCGAAGACGATCAGCAGGCCGACACAGGCGGTGATCTTGATCAGGGTCTCGCGGTCCAGTTCCGGTTCGTCATCCTTGCGCACGATGGGCGTTGGGCGGAACACCAGGTACAGCAGGCCCAGCCCCATCAGGCCGAGCACCAGCAGCGGGTAAGCGCGCGGGCCCACCGGTTCATAGGAGAATGCAGCCTGATAGGGCCAGGCCATCACgg                                                                                                                                                                                            <  1:522862/226‑1 (MQ=255)
                                                        gcTGCGAAAGGGGCGGCGAAGACGATCAGCAGGCCGACACAGGCGGTGATCTTGATCAGGGTCTCGCGGTCCAGTTCCGGTTCGTCATCCTTGCGCACGATGGGCGTTGGGCGGAACACCAGGTACAGCAGGCCCAGCCCCATCAGGCCGAGCACCAGCAGCGGGTAAGCGCGCGGGCCCACCGGTTCATAGGAGAATGCAGCCTGATAGGGCCAGGCCATCACGGCCAGggcgg                                                                                                                                                                                   <  1:184515/235‑1 (MQ=255)
                                                                         gAAGACGATCAGCAGGCCGACACAGGCGGTGATCTTGATCAGGGTCTCGCCGGTCCAGTTCCGGTTCGTCATCCTTGCGCACGATGGGCGTTGGGCGGAACACCAGGTACAGCAGGCCCAGCCCCATCAGGCCGAGCACCAGCAGCGGGTAAGCGCGCGGGCCCACCGGTTCATa                                                                                                                                                                                                                               <  1:247714/175‑1 (MQ=255)
                                                                                     cagGCCGACACAGGCGGTGATCTTGATCAGGGTCTCGCGGTCCAGTTCCGGTTCGTCATCCTTGCGCACGATGGGCGTTGGGCGGAACACCAGGTACAGCAGGCCCAGCCCCATCAGGCCGAGCACCAGCAGCGGGTAAGCGCGCGGGCCCACCGGTTCATAGGAGAATGCAGCCTGATAGGGCCAGGCCATCACGGCCAGGGCGGCACACATCGCCAGCAGGGCCAAGGCGaaaa                                                                                                                                                     <  1:9702/236‑1 (MQ=255)
                                                                                             cacaGGCGGTGATCTTGATCAGGGTCTCGCGGTCCAGTTCCGGTTCGTCATCCTTGCGCACGATGGGCGTTGGGCGGAACACCAGGTACAGCAGGCCCAGCCCCATCAGGCCGAGCACCAGCAGCGGGTAAGCGCGCGGGCCCACCGGTTCATAGGAGAATGCAGCCTg                                                                                                                                                                                                                >  1:529272/1‑169 (MQ=255)
                                                                                                     gTGATCTTGATCAGGGTCTCGCGGTCCAGTTCCGGTTCGTCATCCTTGCGCACGATGGGCGTTGGGCGGAACACCAGGTACAGCAGGCCCAGCCCCATCAGGCCGAGCACCAGCAGCGGGTAAGCGCGCGGGCCCACCGGTTCATAGGAGAATGc                                                                                                                                                                                                                      <  1:577843/155‑1 (MQ=255)
                                                                                                     gTGATCTTGATCAGGGTCTCGCGGTCCAGTTCCGGTTCGTCATCCTTGCGCACGATGGGCGTTGGGCGGAACACCAGGTACAGCAGGCCCAGCCCCATCAGGCCGAGCACCAGCAGCGGGTAAGCGCGCGGGCCCACCGGTTCATAGGAGAATGc                                                                                                                                                                                                                      >  2:577842/1‑155 (MQ=255)
                                                                                                                     tctcGCGGTCCAGTTCCGGTTCGTCATCCTTGCGCACGATGGGCGTTGGGCGGAACACCAGGTACAGCAGGCCCAGCCCCATCAGGCCGAGCACCAGCAGCGGGTAAGCGCGCGGGCCCACCGGTTCATAGGAGAATGCAGCCTGATAGGGCCAGGCCATCACGGCCAGGGCGGCACACATCGCCAGCAGGGCCAAGGCGAAAATGCGTTGAAGGATCATTGGGGAATCCTCGTTg                                                                                                                     >  1:471656/1‑236 (MQ=255)
                                                                                                                         gcgGTCCAGTTCCGGTTCGTCAGCCTTGCGCACGATGGGCGTTGGGCGGAACACCAGGTACAGCAGGCCCAGCCCCATCAGGCCGAGCACCAGCAGCGGGTAAGCGCGCGGGCCCACCGGTTCATAGgaga                                                                                                                                                                                                                          <  2:43221/131‑1 (MQ=255)
                                                                                                                           cgGTCCAGTTCCGGTTCGTCATCCTTGCGCACGATGGGCGTTGGGCGGAACACCAGGTACAGCAGGCCCAGCCCCATCAGGCCGAGCACCAGCAGCGGGTAAGCGCGCGGGCCCACCGGTTCATAGGAGAATGCAGCCTGATAGGGCCAGGCCATCACGGCCAGGGCGGCACACATCGCCAGCAGGGCCAAGGCGAAAATGCGTTGAAGGGTCATTgggg                                                                                                                                >  2:641281/1‑220 (MQ=255)
                                                                                                                                 aGTTCCGGTTCGTCATCCTTGCGCACGATGGGCGTTGGGCGGAACACCAGGTACAGCAGGCCCAGCCCCATCAGGCCGAGCACCAGCAGCGGGTAAGCGCGCGGGCCCACCGGt                                                                                                                                                                                                                                    <  2:611954/114‑1 (MQ=255)
                                                                                                                                 aGTTCCGGTTCGTCATCCTTGCGCACGATGGGCGTTGGGCGGAACACCAGGTACAGCAGGCCCAGCCCCATCAGGCCGAGCACCAGCAGCGGGTAAGCGCGCGGGCCCACCGGt                                                                                                                                                                                                                                    >  1:611955/1‑114 (MQ=255)
                                                                                                                                    tCCGGTTCGTCATCCTTGCGCACGATGGGCGTTGGGCGGAACACCAGGTACAGCAGGCCCAGCCCCATCAGGCCGAGCACCAGCAGCGGGTAAGCGCGCGGGCCCACCGGTTCATAGGAGAATGCAGCCTGATAGGGCCAGGCCATCACGGCCAGGGCGGCACACATCGCCAGCAGGGCCAAGGCGAAAATGCGTTGAAGGATCATTGGGGAATCCTCGTTg                                                                                                                     <  1:20477/222‑1 (MQ=255)
                                                                                                                                      cGGTTCGTCATCCTTGCGCACGATGGGCGTTGGGCGGAACACCAGGTACAGCAGGCCCAGCCCCATCAGGCCGAGCACCAGCAGCGGGTAAGCGCGCGGGCCCACCGGTTCATAGGAGAATGCAGCCTGATAGGGCCAGGCCATCACGGCCAGGGCGGCACACATCGCCAGCAGGGCCAAGGCGAAAATGCGTTGAAGGATCATTGGGGa                                                                                                                               >  1:379719/1‑210 (MQ=255)
                                                                                                                                           cGTCATCCTTGCGCACGATGGGCGTTGGGCGGAACACCAGGTACAGCAGGCCCAGCCCCATCAGGCCGAGCACCAGCAGCGGGTAAGCGCGCGGGCCCACCGGTTCATCGGAGAATGCAGCCTGATAGGGCCGGGCCATCACGGCCAGggcggcc                                                                                                                                                                                 >  2:83881/1‑154 (MQ=255)
                                                                                                                                           cGTCATCCTTGCGCACGATGGGCGTTGGGCGGAACACCAGGTACAGCAGGCCCAGCCCCATCAGGCCGAGCACCAGCAGCGGGTAAGCGCGCGGGCCCACCGGTTCATAGGAGAATGCAGCCTGATAGGGCCAGGCCATCACGGCCAGGGCGGCACACATCGCCAGCAGGGCCAAGGCGAAAATGCGTTGAAGGATCATTGGGGAATCCTCGTTGGGATTg                                                                                                               >  2:298469/1‑221 (MQ=255)
                                                                                                                                           cGTCATCCTTGCGCACGATGGGCGTTGGGCGGAACACCAGGTACAGCAGGCCCAGCCCCATCAGGCCGAGCACCAGCAGCGGGTAAGCGCGCGGGCCCACCGGTTCATAGGAGAATGCAGCCTGATAGGGCCAGGCCATCACGGCCAGGGCGGCACACATCGCCAGCAGGGCCAAGGCGAAAATGCGTTGAAGGATCATTGGGGAATCCTCGTTGGGATTGCTCAAGTGAAAt                                                                                                   <  2:323347/233‑1 (MQ=255)
                                                                                                                                                tCCTTGCGCACGATGGGCGTTGGGCGGAACACCCGGTACAGCAGGCCCAGCCCCATCAGGCCGAGCACCAGCAGCGGGTAAGCGCGCGGGCCCACCGGTTCATAGGAGAATGCAGCCTGATAGGGCCAGGCCATCACGGCCAGGGCGGCACACATCGCCAGCAGGGCCAAGGCGAAAATGCGTTGAAGGATCCTTGGGGAATCCGCGTTGGGATTGCTCAAGTGa                                                                                                      >  2:629998/1‑225 (MQ=255)
                                                                                                                                                tCCTTGCGCACGATGGGCGTTGGGCGGAACACCAGGTACAGCAGGCCCAGCCCCATCAGGCCGAGCACCAGCAGCGGGTAAGCGCGCGGGCCCACCGGTTCATAGGAGAATGCAGCCTGATAGGGCCAGGCCATCACGGCCAGGGCGGCACACATCGCCAGCAGGGCCAAGGCGAAAATGCGTTGAAGGATCATTGGGGAATCCTCGTTGGGATTGCTCAAGTGAAATGCAAATcc                                                                                           <  1:701730/236‑1 (MQ=255)
                                                                                                                                                 ccTTGCGCACGATGGGCGTTGGGCGGAACACCAGGTACAGCGGGCCCAGCCCCATCAGGCCGAGCACCAGCAGCGGGTAAGCGCGCGGGcc                                                                                                                                                                                                                                           <  2:221532/91‑1 (MQ=255)
                                                                                                                                                   ttGCGCACGATGGGCGTTGGGCGGAACACCAGGTACAGCAGGCCCAGCCCCATCAGGCCGAGCACCAGCAGCGGGTAAGCGCGCGGGCCCACCGGTTCATAGGAGAATGCAGCCTGATAGGGCCAGGCCATCACgg                                                                                                                                                                                            <  2:654152/136‑1 (MQ=255)
                                                                                                                                                        cACGATGGGCGTTGGGCGGAACACCAGGTACAGCAGGCCCAGCCCCATCAGGCCGAGCACCAGCAGCGGGTAAGCGCGCGGGCCCACCGGTTCATAGGAGAATGCAGCCTGATAGGGCCAGGCCATCACGGCCAGGGCGGCACACATCGCCAGCAGGGCCAAGGc                                                                                                                                                          >  2:248736/1‑165 (MQ=255)
                                                                                                                                                         aCGATGGGCGTTGGGCGGAACACCAGGTACAGCAGGCCCAGCCCCATCAGGCCGAGCACCAGCAGCGGGTAAGCGCGCGGGCCCACCGGTTCATAGGAGAATg                                                                                                                                                                                                                       <  1:404162/103‑1 (MQ=255)
                                                                                                                                                         aCGATGGGCGTTGGGCGGAACACCAGGTACAGCAGGCCCAGCCCCATCAGGCCGAGCACCAGCAGCGGGTAAGCGCGCGGGCCCACCGGTTCATAGGAGAATg                                                                                                                                                                                                                       >  2:404161/1‑103 (MQ=255)
                                                                                                                                                              gggCGTTGGGCGGAACACCAGGTACAGCAGGCCCAGCCCCATCAGGCCGAGCACCAGCAGCGGGTAAGCGCGCGGGCCCACCGGTTCATAGGAGAATGCAGCCTGATAGGGCCAGGCCGTCACGGCCAGGGCGGCACACATCGCCAGCAGGGCCAAGGCGAAAATGCGTTGAAGGATCATTGGGGAATCCTCGTTGGGATTGCTCAAGTGAAATGCAAATCCTGTGGGAGct                                                                                 >  1:374002/1‑231 (MQ=255)
                                                                                                                                                                             cacCAGGTACAGCAGGCCCAGCCCCATCAGGCCGAGCACCAGCAGCGGGTAAGCGCGCGGGCCCACCGGTTCATAGGAGAATGCAGCCTGATAGggccaggcc                                                                                                                                                                                                   <  2:154198/103‑1 (MQ=255)
                                                                                                                                                                                   gTACAGCAGGCCCAGCCCCATCAGGCCGAGCACCAGCAGCGGGTAAGCGCGCGGGCCCACCGGTTCATAGGAGAATGCAGCCTGATAGGGCCAGGCCATCACGGCCAGGGCGGCACACATCGCCAGCAGGGCCAAGGCGAAAATGCGTTGAAGGATCATTGGGGAATCCTCGTTGGGATTGCTCAAGTGAAATGCAAATCCTGTGGGAGCGGCCTTGCGTCGCGATCGGGCGCGaa                                                        <  1:735228/236‑1 (MQ=255)
                                                                                                                                                                                        gcagGCCCAGCCCCATCAGGCCGAGCACCAGCAGCGGGTAAGCGCGCGGGCCCACCGGTTCATAGGAGAATGCAGCCTGATAGGGCCAGGCCATCACGGCCAGGGCGGCACACATCGCCAGCAGGGCCAAGGCGAAAATGCGTTGAAGGATCATTGGGGAATCCTCGTTGGGATTGCTCAAGTGAAATGCAAATCCTGTGGGAGCGGCCTTGCGTCGCGATCGGGCGCGAAGCTga                                                   >  2:581412/1‑235 (MQ=255)
                                                                                                                                                                                                    ccATCAGGCCGAGCACCAGCAGCGGGTAAGCGCGCGGGCCCACCGGTTCATAGGAGAATGCAGCCTGATAGGGCCAGGCCATCACGGCCAGGGCGGCACACATCGCCAGCAGGGCCAAGGCGAAAATGCGTTGAAGGATCATTGGGGAATCCTCGTTGGGATTGCTCAAGTGAAATGCAAATCCTGTGGGa                                                                                    >  2:683627/1‑191 (MQ=255)
                                                                                                                                                                                                    ccATCAGGCCGAGCACCAGCAGCGGGTAAGCGCGCGGGCCCACCGGTTCATAGGAGAATGCAGCCTGATAGGGCCAGGCCATCACGGCCAGGGCGGCACACATCGCCAGCAGGGCCAAGGCGAAAATGCGTTGAAGGATCATTGGGGAATCCTCGTTGGGATTGCTCAAGTGAAATGCAAATCCTGTGGGAGCGGCCTTGCGTCGCGATCGGGCGCGAAGCTGCCGCTAAACCTGc                                       >  1:469363/1‑236 (MQ=255)
                                                                                                                                                                                                      aTCAGGCCGAGCACCAGCAGCGGGTAAGCGCGCGGGCCCACCGGTTCATAGGAGAATGCAGCCTGATAGGGCCAGGCCATCACGGCCAGGGCGGCACACATCGCCAGCAGGGCCAAGGCGAAAATGCGTTGAAGGATCATTGGGGAATCCTCGTTGGGATTGCTCAAGTGAAATGCAAATCCTGTGGGAGCGGCCTTGCGTCGCGATCGGGCGCGAAGCTGCCGCTAAACCTg                                        >  2:598418/1‑233 (MQ=255)
                                                                                                                                                                                                         aGGCCGAGCACCAGCAGCGGGTAAGCGCGCGGGCCCACCGGTTCATAGGAGAATGCAGCCTGATAGGGCCAGGCCATCACGGc                                                                                                                                                                                           >  2:64616/1‑83 (MQ=255)
                                                                                                                                                                                                          ggCCGAGCACCAGCAGCGGGTAAGCGCGCGGGCCCACCGGTTCATAGGAGAATGCAGCCt                                                                                                                                                                                                                 <  2:495021/60‑1 (MQ=255)
                                                                                                                                                                                                          ggCCGAGCACCAGCAGCGGGTAAGCGCGCGGGCCCACCGGTTCATAGGAGAATGCAGCCt                                                                                                                                                                                                                 >  1:495022/1‑60 (MQ=255)
                                                                                                                                                                                                                  aCCAGCAGCGGGTAAGCGCGCGGGCCCACCGGTTCATAGGAGAATGCAGCCTGATAGGGCCAGGCCATCACGGCCAGGGCGGCACACATCGCCAGCAGGGCCAAGGCGAAAATGCGTTGAAGGATCATTGGGGAATCCTCGTTGGGCTTGCTCAAGTGAAAGGCAAATCCTGTGGGAGCGGCCGTGCGGCGCGATCGGgcggg                                                          >  2:219141/1‑201 (MQ=255)
                                                                                                                                                                                                                    cagcagCGGGTAAGCGCGCGGACCCACCGGTTCATAGGAGAATGCAGCCTGATAGGGCCAGGCCATCACGGCCAGGGCGGCACACATCGCCAGCAGGGCCAAGGCGAAAATGCGTTGAAGGATCATTGGGGAATCCTCGTTGGGATTGCTCAAGTGAAAt                                                                                                   >  1:685410/1‑160 (MQ=255)
                                                                                                                                                                                                                    cagcagCGGGTAAGCGCGCGGACCCACCGGTTCATAGGAGAATGCAGCCTGATAGGGCCAGGCCATCACGGCCAGGGCGGCACACATCGCCAGCAGGGCCAAGGCGAAAATGCGTTGAAGGATCATTGGGGAATCCTCGTTGGGATTGCTCAAGTGAAAt                                                                                                   <  2:685409/160‑1 (MQ=255)
                                                                                                                                                                                                                    cagcagCGGGTAAGCGCGCGGACCCACCGGTTCATAGGAGAATGCAGCCTGATAGGGCCAGGCCACCACGGCCAGGGCGGCACACATCGCCAGCAGGGCCAAGGCGAAAATGCGTTGAAGGATCATTGGGGAATCCTCGTTGGGATTGCTCAAGTGAAATGCAAATCCTGTGGGAGCGGCCTTGCGTCGCGATCGGGCGCGAAGCTGCCGCTAAACCTGCCGATGc                                 <  1:306922/226‑1 (MQ=255)
                                                                                                                                                                                                                    cagcagCGGGTAAGCGCGCGGACCCACCGGTTCATAGGAGAATGCAGCCTGATAGGGCCAGGCCACCACGGCCAGGGCGGCACACATCGCCAGCAGGGCCAAGGCGAAAATGCGTTGAAGGATCATTGGGGAATCCTCGTTGGGATTGCTCAAGTGAAATGCAAATCCTGTGGGAGCGGCCTTGCGTCGCGATCGGGCGCGAAGCTGCCGCTAAACCTGCCGATGc                                 >  2:306921/1‑226 (MQ=255)
                                                                                                                                                                                                                               aaGCGCGCGGGCCCACCGGTTCATAGGAGAATGCAGCCTGATAGGGCCAGGCCATCACGGCCAGGGCGGCACACATCGCCAGCAGGGCCAAGGCGAAAATGCGTTGAAGGATCATTGGGGAATCCTCGTTGGGATTGCTCAAGTGAAATGCAAATCCTGTGGGAGCGGCCTTGCGTCGCGATCGGGCGCGAAGCTGCCGCTAAACCTGCCg                                     >  1:334386/1‑211 (MQ=255)
                                                                                                                                                                                                                                 gcgcgcGGGCCCACCGGTTCATAGGAGAATGCAGCCTGATAGGGCCAGGCCATCACGGCCAGGGCGGCACACATCGCCAGCAGGGCCAAGGCGAAAATGCGTTGAAGGATCATTGGGGAATCCTCGTTGGGATTGCTCAAGTGAAATGCAAATCCTGTGGGAGCGGCCTTGCGTCGCGATCGGGCGCGAAGCTGCCGCTAAACCTGCCGATGCGATCTGTACGAAGCATCGTATCg          <  1:396191/236‑1 (MQ=255)
                                                                                                                                                                                                                                    cgcgGGCCCACCGGTTCATAGGAGAATGCAGCCTGATAGGGCCAGGCCATCACGGCCAGGGCGGCACACATCGCCAGCAGGGCCAAGGCGAAAATGCGTTGAAGGATCATTGGGGAATCCTCGTTGGGATTGCTCAAGTGAAATGCAAATCCTGTGGGAGCGGCCTTGc                                                                          >  2:199889/1‑169 (MQ=255)
                                                                                                                                                                                                                                      cgGGCCCACCGGTTCATAGGAGAATGCAGCCTGATAGGGCCAGGCCATCACGGCCAGGGCGGCACACATCGCCAGCAGGGCCAAGGCGAAAATGCGTTGAAGGATCATTGGGGAATCCTCGTTGGGATTGCTCAAGTGAAATGCAAATCCTGTGGGAGCGGCCTTGCGTCGCGATCGGGCGCGAAGCTGCCGCTAAACCTGCCGATGCGATCTGTACGAAGCATCGTATCGGCtgt     <  1:581413/236‑1 (MQ=255)
                                                                                                                                                                                                                                      cgGGCCCACCGGTTCATAGGAGAATGCAGCCTGATAGGGCCAGGCCATCACGGCCAGGGCGGCACACATCGCCAGCAGGGCCAAGGCGAAAATGCGTTGAAGGATCATTGGGGAATCCTCGTTGGGATTGCTCAAGTGAAATGCAAATCCTGTGGGAGCGGCCTAGCGTCGCGATCGGGCGCGAAGCTGCCGCTAAACCTGCCGATGCGATCTGTACGAAGCATCGTATCGGCtgt     <  1:446233/236‑1 (MQ=255)
                                                                                                                                                                                                                                      cgGGACCACCGGTTCAGAGGAGAATGAAGACGGAGAGGGCCAGGCAATCACGGCCAGGGCGGCACACATCGACAGCAGGGCCAAGGCGAAAATGCGTTGAAGGATCATTGGGGAATCCTCGTTGGGATTGCTCAAg                                                                                                         <  2:72405/136‑1 (MQ=255)
                                                                                                                                                                                                                                         gCCCACCGGTTCATAGGAGAATGCAGCCTGATAGGGCCAGGCCATCACGGCCAGGGCGGCACACATCGCCAGCAGGGCCAAGGCGAAAATGCGTTGAAGGATCATTGGGGAATCCTCGTTGGGATTGCTCAAGTGAAATGCAAATCCTGTGGGAGCGGCCTTGCGTCGCGATCGGGCGCGAAGCTGCCGCTAAACCTGCCGATGCGATCTGTACGAAGCATCGTATCGGCTGTGcc  >  1:693633/1‑236 (MQ=255)
                                                                                                                                                                                                                                         |                                                                                                                                                                                                                                           
CGTACAGCCGAGCCATGGGCAGGCCGACCAGGATGGCGCTGAGGATGAAGCCGAGCGCTTCGAAAGTGGCGGCGAAGACGATCAGCAGGCCGACACAGGCGGTGATCTTGATCAGGGTCTCGCGGTCCAGTTCCGGTTCGTCATCCTTGCGCACGATGGGCGTTGGGCGGAACACCAGGTACAGCAGGCCCAGCCCCATCAGGCCGAGCACCAGCAGCGGGTAAGCGCGCGGGCCCACCGGTTCATAGGAGAATGCAGCCTGATAGGGCCAGGCCATCACGGCCAGGGCGGCACACATCGCCAGCAGGGCCAAGGCGAAAATGCGTTGAAGGATCATTGGGGAATCCTCGTTGGGATTGCTCAAGTGAAATGCAAATCCTGTGGGAGCGGCCTTGCGTCGCGATCGGGCGCGAAGCTGCCGCTAAACCTGCCGATGCGATCTGTACGAAGCATCGTATCGGCTGTGCC  >  AE015451/1617382‑1617849

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