Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 2,785,712 C→A 29.3% H361N (CAC→AAC)  gabD → succinate‑semialdehyde dehydrogenase I, NADP‑dependent

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092732,785,7120CA29.3% 396.7 / 195.3 201H361N (CAC→AAC) gabDsuccinate‑semialdehyde dehydrogenase I, NADP‑dependent
Reads supporting (aligned to +/- strand):  ref base C (89/53);  new base A (38/21);  total (127/74)
Fisher's exact test for biased strand distribution p-value = 8.73e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 207 total aligned reads displayed.

GCAGTGAGCAAACTGCACATCGGCGACGGGCTGGATAACGGCGTCACCATCGGGCCGCTGATCGATGAAAAAGCGGTAGCAAAAGTGGAAGAGCATATTGCCGATGCGCTGGAGAAAGGCGCGCGCGTGGTTTGCGGCGGTAAAGCGCACGAACGCGGCGGCAACTTCTTCCAGCCGACCATTCTGGTGGACGTTCCGGCCAACGCCAAAGTGTCGAAAGAAGAGACGTTCGGCCCCCTCGCCCCGCTGTTCCGCTTTAAAGATGAAGCTGATGTGATTGCGCAAGCCAATGACAC  >  CP009273/2785565‑2785860
                                                                                                                                                   |                                                                                                                                                    
gcagTGAGCAAACTGCACATCGGCGACGGGCTGGATAACGGCGTCACCATCGGGCCGCTGATCGATGAAAAAGCGGTAGCAAAAGTGGAAGAGCATATTGCCGATGCGCTGGAGAAAGGCGCGCGCGTGGTTTGCGGCGGTAAAGCGCa                                                                                                                                                     >  2:447416/1‑149 (MQ=255)
gcagTGAGCAAACTGCACATCGGCGACGGGCTGGATAACGGCGTCACCATCGGGCCGCTGATCGATGAAAAAGCGGTAGCAAAAGTGGAAGAGCATATTGCCGATGCGCTGGAGAAAGGCGCGCGCGTGGTTTGCGGCGGTAAAGCGCa                                                                                                                                                     <  1:3094416/149‑1 (MQ=255)
 cagTGAGCAAACTGCACATCGGCGACGGGCTGGATAACGGCGTCACCATCGGGCCGCTGATCGATGAAAAAGCGGTAGCAAAAGTGGAAGAGCATATTGCCGATGCGCTGGAGAAAGGCGCGCGCGTGGTTTGCGGCGGTAAAGCGCAc                                                                                                                                                    <  2:1385166/149‑1 (MQ=255)
 cagTGAGCAAACTGCACATCGGCGACGGGCTGGATAACGGCGTCACCATCGGGCCGCTGATCGATGAAAAAGCGGTAGCAAAAGTGGAAGAGCATATTGCCGATGCGCTGGAGAAAGGCGCGCGCGTGGTTTGCGGCGGTAAAGCGCAc                                                                                                                                                    <  2:3564338/149‑1 (MQ=255)
 cagTGAGCAAACTGCACATCGGCGACGGGCTGGATAACGGCGTCACCATCGGGCCGCTGATCGATGAAAAAGCGGTAGCAAAAGTGGAAGAGCATATTGCCGATGCGCTGGAGAAAGGCGCGCGCGTGGTTTGCGGCGGTAAAGCGAAc                                                                                                                                                    <  1:138524/149‑1 (MQ=255)
    tGAGCAAACTGCACATCGGCGACGGGCTGGATAACGGCGTCACCATCGGGCCGCTGATCGATGAAAAAGCGGTAGCAAAAGTGGAAGAGCATATTGCCGATGCGCTGGAGAAAGGCGCGCGCGTGGTTTGCGGCGGTAAAGCGCACGaa                                                                                                                                                 <  1:4005132/149‑1 (MQ=255)
     gAGCAAACTGCACATCGGCGACGGGCTGGATAACGGCGTCACCATCGGGCCGCTGATCGATGAAAAAGCGGTAGCAAAAGTGGAAGAGCATATTGCCGATGCGCTGGAGAAAGGCGCGCGCGTGGTTTGCGGCGGTAAAGCGCACGAAc                                                                                                                                                <  1:996595/149‑1 (MQ=255)
        cAAACTGCACATCGGCGACGGGCTGGATAACGGCGTCACCATCGGGCCGCTGATCGATGAAAAAGCGGTAGCAAAAGTGGAAGAGCATATTGCCGATGCGCTGGAGAAAGGCGCGCGCGTGGTTTGCGGCGGTAAAGCGCACGAAcgcg                                                                                                                                             >  2:3514492/1‑149 (MQ=255)
           aCTGCACATCGGCGACGGGCTGGATAACGGCGTCACCATCGGGCCGCTGATCGATGAAAAAGCGGTAGCAAAAGTGGAAGAGCATATTGCCGATGCGCTGGAGAAAGGCGCGCGCGTGGTTTGCGGCGGTAAAGCGCACGAACgcggcg                                                                                                                                          >  1:3153314/1‑149 (MQ=255)
              gCACATCGGCGACGGGCTGGATAACGGCGTCACCATCGGGCCGCTGATCGATGAAAAAGCGGTAGCAAAAGTGGAAGAGCATATTGCCGATGCGCTGGAGAAAGGCGCGCGCGTGGTTTGCGGCGGTAAAGCGCACGAACGCGGCGGCa                                                                                                                                       <  1:3886324/149‑1 (MQ=255)
                acaTCGGCGACGGGCTGGATAACGGCGTCACCATCGGGCCGCTGATCGATGAAAAAGCGGTAGCAAAAGTGGAAGAGCATATTGCCGATGCGCTGGAGAAAGGCGCGCGCGTGGTTTGCGGCGGTAAAGCGCACGAACGCGGCGGCAAc                                                                                                                                     <  1:4087053/149‑1 (MQ=255)
                acaTCGGCGACGGGCTGGATAACGGCGTCACCATCGGGCCGCTGATCGATGAAAAAGCGGTAGCAAAAGTGGAAGAGCATATTGCCGATGCGCTGGAGAAAGGCGCGCGCGTGGTTTGCGGCGGTAAAGCGAACGAACGCGGCGGCAAc                                                                                                                                     >  2:4008493/1‑149 (MQ=255)
                  aTCGGCGACGGGCTGGATAACGGCGTCACCATCGGGCCGCTGATCGATGAAAAAGCGGTAGCAAAAGTGGAAGAGCATATTGCCGATGCGCTGGAGAAAGGCGCGCGCGTGGTTTGCGGCGGTAAAGCGCACGAACGCGGCGGCAActt                                                                                                                                   >  2:3560647/1‑149 (MQ=255)
                    cGGCGACGGGCTGGATAACGGCGTCACCATCGGGCCGCTGATCGATGAAAAAGCGGTAGCAAAAGTGGAAGAGCATATTGCCGATGCGCTGGAGAAAGGCGCGCGCGTGGTTTGCGGCGGTAAAGCGAACGAACGCGGCGGCAActtct                                                                                                                                 >  1:953001/1‑149 (MQ=255)
                    cGGCGACGGGCTGGATAACGGCGTCACCATCGGGCCGCTGATCGATGAAAAAGCGGTAGCAAAAGTGGAAGAGCATATTGCCGATGCGCTGGAGAAAGGCGCGCGCGTGGTTTGCGGCGGTAAAGCGAACGAACGCGGCGGCAActtct                                                                                                                                 >  1:2769661/1‑149 (MQ=255)
                       cGACGGGCTGGATAACGGCGTCACCATCGGGCCGCTGATCGATGAAAAAGCGGTAGCAAAAGTGGAAGAGCATATTGCCGATGCGCTGGAGAAAGGCGCGCGCGTGGTTTGCGGCGGTAAAGCGAACGAACGCGGCGGCAACTTCTTcc                                                                                                                              <  1:3348641/149‑1 (MQ=255)
                        gACGGGCTGGATAACGGCGTCACCATCGGGCCGCTGATCGATGAAAAAGCGGTAGCAAAAGTGGAAGAGCATATTGCCGATGCGCTGGAGAAAGGCGCGCGCGTGGTTTGCGGCGGTAAAGCGCACGAACGCGGCGGCAACTTCTTCCa                                                                                                                             >  1:3572762/1‑149 (MQ=255)
                        gACGGGCTGGATAACGGCGTCACCATCGGGCCGCTGATCGATGAAAAAGCGGTAGCAAAAGTGGAAGAGCATATTGCCGATGCGCTGGAGAAAGGCGCGCGCGTGGTTTGCGGCGGTAAAGCGCACGAACGCGGCGGCAACTTCTTCCa                                                                                                                             >  1:3264608/1‑149 (MQ=255)
                                   tAACGGCGTCACCATCGGGCCGCTGATCGATGAAAAAGCGGTAGCAAAAGTGGAAGAGCATATTGCCGATGCGCTGGAGAAAGGCGCGCGCGTGGTTTGCGGCGGTAAAGCGCACGAACGCGGCGGCAACTTCTTCCAGCCGACCATTc                                                                                                                  >  1:3768796/1‑149 (MQ=255)
                                    aaCGGCGTCACCATCGGGCCGCTGATCGATGAAAAAGCGGTAGCAAAAGTGGAAGAGCATATTGCCGATGCGCTGGAGAAAGGCGCGCGCGTGGTTTGCGGCGGTAAAGCGCACGAACGCGGCGGCAACTTCTTCCAGCCGACCATTCt                                                                                                                 <  2:2852208/149‑1 (MQ=255)
                                    aaCGGCGTCACCATCGGGCCGCTGATCGATGAAAAAGCGGTAGCAAAAGTGGAAGAGCATATTGCCGATGCGCTGGAGAAAGGCGCGCGCGTGGTTTGCGGCGGTAAAGCGAACGAACGCGGCGGCAACTTCTTCCAGCCGACCATTCt                                                                                                                 >  2:1527861/1‑149 (MQ=255)
                                     aCGGCGTCACCATCGGGCCGCTGATCGATGAAAAAGCGGTAGCAAAAGTGGAAGAGCATATTGCCGATGCGCTGGAGAAAGGCGCGCGCGTGGTTTGCGGCGGTAAAGCGCACGAACGCGGCGGCAACTTCTTCCAGCCGACCAGTCtg                                                                                                                >  2:3018532/1‑149 (MQ=255)
                                       ggCGTCACCATCGGGCCGCTGATCGATGAAAAAGCGGTAGCAAAAGTGGAAGAGCATATTGCCGATGCGCTGGAGAAAGGCGCGCGCGTGGTTTGCGGCGGTAAAGCGCACGAACGCGGCGGCAACTTCTTCCAGCCGACCATTCtggt                                                                                                              <  1:379502/149‑1 (MQ=255)
                                           tCACCATCGGGCCGCTGATCGATGAAAAAGCGGTAGCAAAAGTGGAAGAGCATATTGCCGATGCGCTGGAGAAAGGCGCGCGCGTGGTTTGCGGCGGTAAAGCGCACGAACGCGGCGGCAACTTCTTCCAGCCGACCATTCTGGTGGAc                                                                                                          <  2:2946644/149‑1 (MQ=255)
                                            cACCATCGGGCCGCTGATCGATGAAAAAGCGGTAGCAAAAGTGGAAGAGCATATTGCCGATGCGCTGGAGAAAGGCGCGCGCGTGGTTTGCGGCGGTAAAGCGCACGAACGCGGCGGCAACTTCTTCCAGCCGACCATTCTGGTGGACg                                                                                                         >  1:2386855/1‑149 (MQ=255)
                                              ccATCGGGCCGCTGATCGATGAAAAAGCGGTAGCAAAAGTGGAAGAGCATATTGCCGATGCGCTGGAGAAAGGCGCGCGCGTGGTTTGCGGCGGTAAAGCGAACGAACGCGGCGGCAACTTCTTCCAGCCGACCATTCTGGTGGACGtt                                                                                                       >  2:1856392/1‑149 (MQ=255)
                                              ccATCGGGCCGCTGATCGATGAAAAAGCGGTAGCAAAAGTGGAAGAGCATATTGCCGATGCGCTGGAGAAAGGCGCGCGCGTGGTTTGCGGCGGTAAAGCGAACGAACGCGGCGGCAACTTCTTCCAGCCGACCATTCTGGTGGACGtt                                                                                                       >  2:1195370/1‑149 (MQ=255)
                                                aTCGGGCCGCTGATCGATGAAAAAGCGGTAGCAAAAGTGGAAGAGCATATTGCCGATGCGCTGGAGAAAGGCGCGCGCGTGGTTTGCGGCGGTAAAGCGCACGAACGCGGCGGCAACTTCTTCCAGCCGACCATTCTGGTGGACGTTcc                                                                                                     >  1:852006/1‑149 (MQ=255)
                                                    ggCCGCTGATCGATGAAAAAGCGGTAGCAAAAGTGGAAGAGCATATTGCCGATGCGCTGGAGAAAGGCGCGCGCGTGGTTTGCGGCGGTAAAGCGCACGAACGCGGCGGCAACTTCTTCCAGCCGACCATTCTGGTGGACGTTCCGGcc                                                                                                 <  1:1854576/149‑1 (MQ=255)
                                                      ccGCTGATCGATGAAAAAGCGGTAGCAAAAGTGGAAGAGCATATTGCCGATGCGCTGGAGAAAGGCGCGCGCGTGGTTTGCGGCGGTAAAGCGCACGAACGCGGCGGCAACTTCTTCCAGCCGACCATTCTGGTGGACGTTCCGGCCaa                                                                                               >  2:50017/1‑149 (MQ=255)
                                                        gCTGATCGATGAAAAAGCGGTAGCAAAAGTGGAAGAGCATATTGCCGATGCGCTGGAGAAAGGCGCGCGCGTGGTTTGCGGCGGTAAAGCGAACGAACGCGGCGGCAACTTCTTCCAGCCGACCATTCTGGTGGACGTTCCGGCCAACg                                                                                             >  2:1595147/1‑149 (MQ=255)
                                                         cTGATCGATGAAAAAGCGGTAGCAAAAGTGGAAGAGCATATTGCCGATGCGCTGGAGAAAGGCGCGCGCGTGGTTTGCGGCGGTAAAGCGCACGAACGCGGCGGCAACTTCTTCCAGCCGACCATTCTGGTGGACGTTCCGGCCAACGc                                                                                            >  1:1582899/1‑149 (MQ=255)
                                                          tGATCGATGAAAAAGCGGTAGCAAAAGTGGAAGAGCATATTGCCGATGCGCTGGAGAAAGGCGCGCGCGTGGTTTGCGGCGGTAAAGCGAACGAACGCGGCGGCAACTTCTTCCAGCCGACCATTCTGGTGGACGTTCCGGCCAACGcc                                                                                           >  2:1388060/1‑149 (MQ=255)
                                                           taTCGATGAAAAATCGGTATCAAAAGTGGAAGATCAAATTTCCGATGCGATGGATAAAGGCGCGCGCGTGGTTAGAGTCGGTAAATAGAACGAACGCGTCGGCAACTTCTTCCAGCCGACCAAACTGGTGTACGTTCCGGCCAACGCCa                                                                                          >  2:200410/2‑149 (MQ=255)
                                                                aTGAAAAAGCGGTAGCAAAAGTGGAAGAGCATATTGCCGATGCGCTGGAGAAAGGCGCGCGCGTGGTTTGCGGCGGTAAAGCGCACGAACGCGGCGGCAACTTCTTCCAGCCGACCATTCTGGTGGACGTTCCGGCCAACGCCAAAgtg                                                                                     >  2:1519266/1‑149 (MQ=255)
                                                                 tGAAAAAGCGGTAGCAAAAGTGGAAGAGCATATTGCCGATGCGCTGGAGAAAGGCGCGCGCGTGGTTTGCGGCGGTAAAGCGCACGAACGCGGCGGCAACTTCTTCCAGCCGACCATTCTGGTGGACGTTCCGGCCAACGCCAAAgtgt                                                                                    >  2:4249505/1‑149 (MQ=255)
                                                                  gAAAAAGCGGTAGCAAAAGTGGAAGAGCATATTGCCGATGCGCTGGAGAAAGGCGCGCGCGTGGTTTGCGGCGGTAAAGCGCACGAACGCGGCGGCAACTTCTTCCAGCCGACCATTCTGGTGGACGTTCCGGCCAACGCCAAAGTGTc                                                                                   >  1:1165274/1‑149 (MQ=255)
                                                                  gAAAAAGCGGTAGCAAAAGTGGAAGAGCATATTGCCGATGCGCTGGAGAAAGGCGCGCGCGTGGTTTGCGGCGGTAAAGCGAACGAACGCGGCGGCAACTTCTTCCAGCCGACCATTCTGGTGGACGTTCCGGCCAACGCCAAAGTGTc                                                                                   <  2:1551472/149‑1 (MQ=255)
                                                                   aaacaGCGGTAGCAAACGTGGAAGAGCATATTGCCGCTGCGCTGGAGAAAGGCGCGCGCGTGGTTTGCGGCGGTAAAGCGCACGCACGCGGCGGCAACTTCTTCCAGCCGACCATTCTGGTGGACGTTCCGGCCACCGCCAAAGTGTCg                                                                                  <  1:3891359/149‑1 (MQ=255)
                                                                   aaaaaGCGGTAGCAAAAGTGGAAGAGCATATTGCCGATGCGCTGGAGAAAGGCGCGCGCGTGGTTTGCGGCGGTAAAGCGCACGAACGCGGCGGCAACTTCTTCCAGCCGACCATTCTGGTGGACGTTCCGGCCAACGCCAAAGTGTCg                                                                                  >  2:1011488/1‑149 (MQ=255)
                                                                     aaaGCGGTAGCAAAAGTGGAAGAGCATATTGCCGATGCGCTGGAGAAAGGCGCGCGCGTGGTTTGCGGCGGTAAAGCGCACGAACGCGGCGGCAACTTCTTCCAGCCGACCATTCTGGTGGACGTTCCGGCCAACGCCAAAGTGTCGaa                                                                                >  1:3910520/1‑149 (MQ=255)
                                                                      aaGCGGTAGCAAAAGTGGAAGAGCATATTGCCGATGCGCTGGAGAAAGGCGCGCGCGTGGTTTGCGGCGGTAAAGCGCACGAACGCGGCGGCAACTTCTTCCAGCCGACCATTCTGGTGGACGTTCCGGCCAACGCCAAAGTGTCGaaa                                                                               >  2:2770192/1‑149 (MQ=255)
                                                                      aaGCGGTAGCAAAAGTGGAAGAGCATATTGCCGATGCGCTGGAGAAAGGCGCGCGCGTGGTTTGCGGCGGTAAAGCGCACGAACGCGGCGGCAACTTCTTCCAGCCGACCATTCTGGTGGACGTTCCGGCCAACGCCAAAGTGTCGaaa                                                                               >  2:613859/1‑149 (MQ=255)
                                                                        gCGGTAGCAAAAGTGGAAGAGCATATTGCCGATGCGCTGGAGAAAGGCGCGCGCGTGGTTTGCGGCGGTAAAGCGCACGAACGCGGCGGCAACTTCTTCCAGCCGACCATTCTGGTGGACGTTCCGGCCAACGCCAAAGTGTCGAaaga                                                                             >  1:2844701/1‑149 (MQ=255)
                                                                        gCGGTAGCAAAAGTGGAAGAGCATATTGCCGATGCGCTGGAGAAAGGCGCGCGCGTGGTTTGCGGCGGTAAAGCGCACGAACGCGGCGGCAACTTCTTCCAGCCGACCATTCTGGTGGACGTTCCGGCCAACGCCAAAGTGTCGAaaga                                                                             <  2:1777652/149‑1 (MQ=255)
                                                                        gCGGTAGCAAAAGTGGAAGAGCATATTGCCGATGCGCTGGAGAAAGGCGCGCGCGTGGTTTGCGGCGGTAAAGCGCACGAACGCGGCGGCAACTTCTTCCAGCCGACCATTCTGGTGGACGTTCCGGCCAACGCCAAAGTGTCGAaaga                                                                             <  2:3951543/149‑1 (MQ=255)
                                                                            tAGCAAAAGTGGAAGAGCATATTGCCGATGCGCTGGAGAAAGGCGCGCGCGTGGTTTGCGGCGGTAAAGCGCACGAACGCGGCGGCAACTTCTTCCAGCCGACCATTCTGGTGGACGTTCCGGCCAACGCCAAAGTGTCGAAAGAagag                                                                         <  1:2222814/149‑1 (MQ=255)
                                                                            tAGCAAAAGTGGAAGAGCATATTGCCGATGCGCTGGAGAAAGGCGCGCGCGTGGTTTGCGGCGGTAAAGCGCACGAACGCGGCGGCAACTTCTTCCAGCCGACCATTCTGGTGGACGTTCCGGCCAACGCCAAAGTGTCGAAAGAagag                                                                         <  1:834588/149‑1 (MQ=255)
                                                                             aGCAAAAGTGGAAGAGCATATTGCCGATGCGCTGGAGAAAGGCGCGCGCGTGGTTTGCGGCGGTAAAGCGCACGAACGCGGCGGCAACTTCTTCCAGCCGACCATTCTGGTGGACGTTCCGGCCAACGCCAAAGTGTCGAAAGAAgaga                                                                        >  1:762036/1‑149 (MQ=255)
                                                                              gCAAAAGTGGAAGAGCATATTGCCGATGCGCTGGAGAAAGGCGCGCGCGTGGTTTGCGGCGGTAAAGCGCACGAACGCGGCGGCAACTTCTTCCAGCCGACCATTCTGGTGGACGTTCCGGCCAACGCCAAAGTGTCGAAAGAAGAGAc                                                                       <  1:3514492/149‑1 (MQ=255)
                                                                                 aaaGTGGAAGAGCATATTGCCGATGCGCTGGAGAAAGGCGCGCGCGTGGTTTGCGGCGGTAAAGCGCACGAACGCGGCGGCAACTTCTTCCAGCCGACCATTCTGGTGGACGTTCCGGCCAACGCCAAAGTGTCGAAAGAAGAGACGtt                                                                    <  1:1675327/149‑1 (MQ=255)
                                                                                 aaaGTGGAAGAGCATATTGCCGATGCGCTGGAGAAAGGCGCGCGCGTGGTTTGCGGCGGTAAAGCGAACGAACGCGGCGGCAACTTCTTCCAGCCGACCATTCTGGTGGACGTTCCGGCCAACGCCAAAGTGTCGAAAGAAGAGACGtt                                                                    <  1:640939/149‑1 (MQ=255)
                                                                                  aaGTGGAAGAGCATATTGCCGATGCGCTGGAGAAAGGCGCGCGCGTGGTTTGCGGCGGTAAAGCGCACGAACGCGGCGGCAACTTCTTCCAGCCGACCATTCTGGTGGACGTTCCGGCCAACGCCAAAGTGTCGAAAGAAGAGACGTTc                                                                   <  1:1464339/149‑1 (MQ=255)
                                                                                  aaGTGGAAGAGCATATTGCCGATGCGCTGGAGAAAGGCGCGCGCGTGGTTTGCGGCGGTAAAGCGCACGAACGCGGCGGCAACTTCTTCCAGCCGACCATTCTGGTGGACGTTCCGGCCAACGCCAAAGTGTCGAAAGAAGAGACGTTc                                                                   <  2:1440107/149‑1 (MQ=255)
                                                                                  aaGTGGAAGAGCATATTGCCGATGCGCTGGAGAAAGGCGCGCGCGTGGTTTGCGGCGGTAAAGCGCACGAACGCGGCGGCAACTTCTTCCAGCCGACCATTCTGGTGGACGTTCCGGCCAACGCCAAAGTGTCGAAAGAAGAGACGTTc                                                                   <  1:4197161/149‑1 (MQ=255)
                                                                                  aaGTGGAAGAGCATATTGCCGATGCGCTGGAGAAAGGCGCGCGCGTGGTTTGCGGCGGTAAAGCGCACGAACGCGGCGGCAACTTCTTCCAGCCGACCATTCTGGTGGACGTTCCGGCCAACGCCAAAGTGTCGAAAGAAGAGACGTTc                                                                   <  1:1765984/149‑1 (MQ=255)
                                                                                  aaGTGGAAGAGCATATTGCCGATGCGCTGGAGAAAGGCGCGCGCGAGGTTTGCGGCGGTAAAGCGCACGAACGCGGCGGCAACTTCTTCCAGCCGACCATTCTGGTGGACGTTCCGGCCAACGCCCAAGTGTCGAAAGAAGAGACGtta                                                                   >  2:3512918/1‑148 (MQ=255)
                                                                                   aGTGGAAGAGCATATTGCCGATGCGCTGGAGAAAGGCGCGCGCGTGGTTTGCGGCGGTAAAGCGCACGAACGCGGCGGCAACTTCTTCCAGCCGACCATTCTGGTGGACGTTCCGGCCAACGCCAAAGTGTCGAAAGAAGAGACGTTCg                                                                  <  2:793375/149‑1 (MQ=255)
                                                                                     tGGAAGAGCATATTGCCGATGCGCTGGAGAAAGGCGCGCGCGTGGTTTGCGGCGGTAAAGCGAACGAACGCGGCGGCAACTTCTTCCAGCCGACCATTCTGGTGGACGTTCCGGCCAACGCCAAAGTGTCGAAAGAAGAGACGTTCGGc                                                                >  1:485841/1‑149 (MQ=255)
                                                                                      ggAAGAGCATATTGCCGATGCGCTGGAGAAAGGCGCGCGCGTGGTTTGCGGCGGTAAAGCGCACGAACGCGGCGGCAACTTCTTCCAGCCGACCATTCTGGTGGACGTTCCGGCCAACGCCAAAGTGTCGAAAGAAGAGACGTTCGGcc                                                               >  1:4304885/1‑149 (MQ=255)
                                                                                        aaGAGCATAGTGCCGATGCGCTGGAGAAAGGCGCGCGCGTGGTTTGCGGCGGTAAAGCGCACGAACGCGGCGGCAACTTCTTCCAGCCGACCATTCTGGTGGACGTTCCGGCCAACGCCAAAGTGTCGAAAGAAGAGACGTTCGGcccc                                                             <  2:318792/149‑1 (MQ=255)
                                                                                         agagCATATTGCCGATGCGCTGGAGAAAGGCGCGCGCGTGGTTTGCGGCGGTAAAGCGCACGAACGCGGCGGCAACTTCTTCCAGCCGACCATTCTGGTGGACGTTCCGGCCAACGCCAAAGTGTCGAAAGAAGAGACGTTCGGccccc                                                            >  2:753482/1‑149 (MQ=255)
                                                                                          gagCATATTGCCGATGCGCTGGAGAAAGGCGCGCGCGTGGTTTGCGGCGGTAAAGCGAACGAACGCGGCGGCAACTTCTTCCAGCCGACCATTCTGGTGGACGTTCCGGCCAACGCCAAAGTGTCGAAAGAAGAGACGTTCGGCCCCCt                                                           >  1:1732806/1‑149 (MQ=255)
                                                                                            gCATATTGCCGATGCGCTGGAGAAAGGCGCGCGCGTGGTTTGCGGCGGTAAAGCGCACGAACGCGGCGGCAACTTCTTCCAGCCGACCATTCTGGTGGACGTTCCGGCCAACGCCAAAGTGTCGAAAGAAGAGACGTTCGGCCCCCTCg                                                         >  2:549175/1‑149 (MQ=255)
                                                                                            gCATATTGCCGATGCGCTGGAGAAAGGCGCGCGCGTGGTTTGCGGCGGTAAAGCGCACGAACGCGGCGGCAACTTCTTCCAGCCGACCATTCTGGTGGACGTTCCGGCCAACGCCAAAGTGTCGAAAGAAGAGACGTTCGGCCCCCTCg                                                         >  2:3673745/1‑149 (MQ=255)
                                                                                              atatTGCCGATGCGCTGGAGAAAGGCGCGCGCGTGGTTTGCGGCGGTAAAGCGCACGAACGCGGCGGCAACTTCTTCCAGCCGACCATTCTGGTGGACGTTCCGGCCAACGCCAAAGTGTCGAAAGAAGAGACGTTCGGCCCCCTCGcc                                                       >  1:3965205/1‑149 (MQ=255)
                                                                                              atatTGCCGATGCGCTGGAGAAAGGCGCGCGCGTGGTTTGCGGCGGTAAAGCGCACGAACGCGGCGGCAACTTCTTCCAGCCGACCATTCTGGTGGACGTTCCGGCCAACGCCAAAGTGTCGAAAGAAGAGACGTTCGGCCCCCTCGcc                                                       >  1:3965227/1‑149 (MQ=255)
                                                                                              atatTGCCGATGCGCTGGAGAAAGGCGCGCGCGTGGTTTGCGGCGGTAAAGCGCACGAACGCGGCGGCAACTTCTTCCAGCCGACCATTCTGGTGGACGTTCCGGCCAACGCCAAAGTGTCGAAAGAAGAGACGTTCGGCCCCCTCGcc                                                       >  2:1965534/1‑149 (MQ=255)
                                                                                              atatTGCCGATGCGCTGGAGAAAGGCGCGCGCGTGGTTTGCGGCGGTAAAGCGAACGAACGCGGCGGCAACTTCTTCCAGCCGACCATTCTGGTGGACGTTCCGGCCAACGCCAAAGTGTCGAAAGAAGAGACGTTCGGCCCCCTCGcc                                                       >  1:1030855/1‑149 (MQ=255)
                                                                                               tagtGCCGATGCGCTGGAGAAAGGCGCGCGCGTGGTTTGCGGCGGTAAAGCGCACGAACGCGGCGGCAACTTCTTCCAGCCGACCATTCTGGTGGACGTTCCGGCCAACGCCAAAGTGTCGAAAGAAGAGACGTTCGGCCCCCTCGccc                                                      <  1:4140452/146‑1 (MQ=255)
                                                                                                atTGCCGATGCGCTGGAGAAAGGCGCGCGCGTGGTTTGCGGCGGTAAAGCGCACGAACGCGGCGGCAACTTCTTCCAGCCGACCATTCTGGTGGACGTTCCGGCCAACGCCAAAGTGTCGAAAGAAGAGACGTTCGGCCCCCTCGcccc                                                     >  2:3775765/1‑149 (MQ=255)
                                                                                                atTGCCGATGCGCTGGAGAAAGGCGCGCGCGTGGTTTGCGGCGGTAAAGCGCACGAACGCGGCGGCAACTTCTTCCAGCCGACCATTCTGGTGGACGTTCCGGCCAACGCCAAAGTGTCGAAAGAAGAGACGTTCGGCCCCCTCGcccc                                                     >  2:215996/1‑149 (MQ=255)
                                                                                                atTGCCGATGCGCTGGAGAAAGGCGCGCGCGTGGTTTGCGGCGGTAAAGCGAACGAACGCGGCGGCAACTTCTTCCAGCCGACCATTCTGGTGGACGTTCCGGCCAACGCCAAAGTGTCGAAAGAAGAGACGTTCGGCCCCCTCGcccc                                                     >  2:1647882/1‑149 (MQ=255)
                                                                                                      gATGCGCTGGAGAAAGGCGCGCGCGTGGTTTGCGGCGGTAAAGCGCACGAACGCGGCGGCAACTTCTTCCAGCCGACCATTCTGGTGGACGTTCCGGCCAACGCCAAAGTGTCGAAAGAAGAGACGTTCGGCCCCCTCGCCCCGCTGtt                                               <  1:583960/149‑1 (MQ=255)
                                                                                                      gATGCGCTGGAGAAAGGCGCGCGCGTGGTTTGCGGCGGTAAAGCGAACGAACGCGGCGGCAACTTCTTCCAGCCGACCATTCTGGTGGACGTTCCGGCCAACGCCAAAGTGTCGAAAGAAGAGACGTTCGGCCCCCTCGCCCCGCTGtt                                               >  1:3602934/1‑149 (MQ=255)
                                                                                                       aggcgcTGGAGAACGGCGCGCGCGTGGTTTGCGGCGGTAAAGCGAACGACCGCGGCGGCAACTTCGTCCAGCCGACCATTCTGGTGGACGTTCCGGCCAACGCCAAAGTGTCGAAAGAAGAGACGTTCGGCCCCCTCGCCCCGCTGTTc                                              <  2:4148819/147‑1 (MQ=255)
                                                                                                       aTGCGCTGGAGAAAGGCGCGCGCGTGGTTTGCGGCGGTAAAGCGAACGAACGCGGCGGCAACTTCTTCCAGCCGACCATTCTGGTGGACGTTCCGGCCAACGCCAAAGTGTCGAAAGAAGAGACGTTCGGCCCCCTCGCCCCGCTGTTc                                              <  1:3886167/149‑1 (MQ=255)
                                                                                                        tGCGCTGGAGAAAGGCGCGCGCGTGGTTTGCGGCGGTAAAGCGCACGAACGCGGCGGCAACTTCTTCCAGCCGACCATTCTGGTGGACGTTCCGGCCAACGCCAAAGTGTCGAAAGAAGAGACGTTCGGCCCCCTCGCCCCGCTGTTcc                                             >  2:1877111/1‑149 (MQ=255)
                                                                                                          cgcTGGAGAAAGGCGCGCGCGTGGTTTGCGGCGGTAAAGCGCACGAACGCGGCGGCAACTTCTTCCAGCCGACCATTCTGGTGGACGTTCCGGCCAACGCCAAAGTGTCGAAAGAAGAGACGTTCGGCCCCCTCGCCCCGCTGTTCCGc                                           <  2:2364832/149‑1 (MQ=255)
                                                                                                           gcTGGAGAAAGGCGCGCGCGTGGTTTGCGGCGGTAAAGCGCACGAACGCGGCGGCAACTTCTTCCAGCCGACCATTCTGGTGGACGTTCCGGCCAACGCCAAAGTGTCGAAAGAAGAGACGTTCGGCCCCCTCGCCCCGCTGTTCCGCt                                          <  1:3863875/149‑1 (MQ=255)
                                                                                                            cTGGAGAAAGGCGCGCGCGTGGTTTGCGGCGGTAAAGCGCACGAACGCGGCGGCAACTTCTTCCAGCCGACCATTCTGGTGGACGTTCCGGCCAACGCCAAAGTGTCGAAAGAAGAGACGTTCGGCCCCCTCGCCCCGCTGTTCCGCtt                                         >  2:1074119/1‑149 (MQ=255)
                                                                                                            cTGGAGAAAGGCGCGCGCGTGGTTTGCGGCGGTAAAGCGAACGAACGCGGCGGCAACTTCTTCCAGCCGACCATTCTGGTGGACGTTCCGGCCAACGCCAAAGTGTCGAAAGAAGAGACGTTCGGCCCCCTCGCCCCGCTGTTCCGCtt                                         >  2:3919837/1‑149 (MQ=255)
                                                                                                                agaAAGGCGCGCGCGTGGTTTGCGGCGGTAAAGCGCACGAACGCGGCGGCAACTTCTTCCAGCCGACCATTCTGGTGGACGTTCCGGCCAACGCCAAAGTGTCGAAAGAAGAGACGTTCGGCCCCCTCGCCCCGCTGTTCCGCTTTaaa                                     <  1:622182/149‑1 (MQ=255)
                                                                                                                agaAAGGCGCGCGCGTGGTTTGCGGCGGTAAAGCGAACGAACGCGGCGGCAACTTCTTCCAGCCGACCATTCTGGTGGACGTTCCGGCCAACGCCAAAGTGTCGAAAGAAGAGACGTTCGGCCCCCTCGCCCCGCTGTTCCGCTTTaaa                                     <  1:3160789/149‑1 (MQ=255)
                                                                                                                    aGGCGCGCGCGTGGTTTGCGGCGGTAAAGCGCACGAACGCGGCGGCAACTTCTTCCAGCCGACCATTCTGGTGGACGTTCCGGCCAACGCCAAAGTGTCGAAAGAAGAGACGTTCGGCCCCCTCGCCCCGCTGTTCCGCTTTAAAGATg                                 <  2:2708286/149‑1 (MQ=255)
                                                                                                                     ggCGCGCGCGTGGTTTGCGGCGGTAAAGCGCACGAACGCGGCGGCAACTTCTTCCAGCCGACCAGTCTGGTGGACGTTCCGGCCAACGCCAAAGTGTCGAAAGAAGTGACGTTCGGCCCCCTCGCCCCGCTGTTCCGCTTTAACGATGa                                >  2:4314213/1‑149 (MQ=255)
                                                                                                                      gcgcgcgcGTGGTTTGCGGCGGTAAAGCGCACGAACGCGGCGGCAACTTCTTCCAGCCGACCATTCTGGTGGACGTTCCGGCCAACGCCAAAGTGTCGAAAGAAGAGACGTTCGGCCCCCTCGCCCCGCTGTTCCGCTTTAAAGATGaa                               <  2:3575498/149‑1 (MQ=255)
                                                                                                                          gcgcgTGGTTTGCGGCGGTAAAGCGCACGAACGCGGCGGCAACTTCTTCCAGCCGACCATTCTGGTGGACGTTCCGGCCAACGCCAAAGTGTCGAAAGAAGAGACGTTCGGCCCCCTCGCCCCGCTGTTCCGCTTTAAAGATGAAGCTg                           >  1:775290/1‑149 (MQ=255)
                                                                                                                          gcgcgTGGTTTGCGGCGGTAAAGCGCACGAACGCGGCGGCAACTTCTTCCAGCCGACCATTCTGGTGGACGTTCCGGCCAACGCCAAAGTGTCGAAAGAAGAGACGTTCGGCCCCCTCGCCCCGCTGTTCCGCTTTAAAGATGAAGCTg                           >  1:261560/1‑149 (MQ=255)
                                                                                                                          gcgcgTGGTTTGCGGCGGTAAAGCGAACGAACGCGGCGGCAACTTCTTCCAGCCGACCATTCTGGTGGACGTTCCGGCCAACGCCAAAGTGTCGAAAGAAGAGACGTTCGGCCCCCTCGCCCCGCTGTTCCGCTTTAAAGATGAAGCTg                           >  2:2721262/1‑149 (MQ=255)
                                                                                                                            gcgTGGTTTGCGGCGGTAAAGCGAACGAACGCGGCGGCAACTTCTTCCAGCCGACCATTCTGGTGGACGTTCCGGCCAACGCCAAAGTGTCGAAAGAAGAGACGTTCGGCCCCCTCGCCCCGCTGTTCCGCTTTAAAGATGAAGCTGAt                         <  1:4008932/149‑1 (MQ=255)
                                                                                                                                ggTTTGCGGCGGTAAAGCGCACGAACGCGGCGGCAACTTCTTCCAGCCGACCATTCTGGTGGACGTTCCGGCCAACGCCAAAGTGTCGAAAGAAGAGACGTTCGGCCCCCTCGCCCCGCTGTTCCGCTTTAAAGATGAAGCTGATGTGa                     >  2:2517427/1‑149 (MQ=255)
                                                                                                                                 gTTTGCGGCGGTAAAGCGCACGAACGCGGCGGCAACTTCTTCCAGCCGACCATTCTGGTGGACGTTCCGGCCAACGCCAAAGTGTCGAAAGAAGAGACGTTCGGCCCCCTCGCCCCGCTGTTCCGCTTTAAAGATGAAGCTGATGTGAt                    >  1:1598404/1‑149 (MQ=255)
                                                                                                                                 gTTTGCGGCGGTAAAGCGCACGAACGCGGCGGCAACTTCTTCCAGCCGACCATTCTGGTGGACGTTCCGGCCAACGCCAAAGTGTCGAAAGAAGAGACGTTCGGCCCCCTCGCCCCGCTGTTCCGCTTTAAAGATGAAGCTGATGTGAt                    >  1:960355/1‑149 (MQ=255)
                                                                                                                                    tGCGGCGGTAAAGCGAACGAACGCGGCGGCAACTTCTTCCAGCCGACCATTCTGGTGGACGTTCCGGCCAACGCCAAAGTGTCGAAAGAAGAGACGTTCGGCCCCCTCGCCCCGCTGTTCCGCCTTAAAGATGAAGCTGATGTGATTgc                 <  2:3263761/149‑1 (MQ=255)
                                                                                                                                       ggcggTAAAGCGAACGAACGCGGCGGCAACTTCTTCCAGCCGACCATTCTGGTGGACGTTCCGGCCAACGCCAAAGTGTCGAAAGAAGAGACGTTCGGCCCCCTCGCCCCGCTGTTCCGCTTTAAAGATGAAGCTGATGTGATTGCGCa              >  1:2689613/1‑149 (MQ=255)
                                                                                                                                       ggcggTAAAGCGAACGAACGCGGCGGCAACTTCTTCCAGCCGACCATTCTGGTGGACGTTCCGGCCAACGCCAAAGTGTCGAAAGAAGAGACGTTCGGCCCCCTCGCCCCGCTGTTCCGCTTTAAAGATGAAGCTGATGTGATTGCGCa              >  1:533116/1‑149 (MQ=255)
                                                                                                                                             aaaGCGCACGAACGCGGCGGCAACTTCTTCCAGCCGACCATTCTGGTGGACGTTCCGGCCAACGCCAAAGTGTCGAAAGAAGAGACGTTCGGCCCCCTCGCCCCGCTGTTCCGCTTTAAAGATGAAGCTGATGTGATTGCGCAAGCCaa        <  2:2870405/149‑1 (MQ=255)
                                                                                                                                               aGCGAACGAACGCGGCGGCAACTTCTTCCAGCCGACCATTCTGGTGGACGTTCCGGCCAACGCCAAAGTGTCGAAAGAAGAGACGTTCGGCCCCCTCGCCCCGCTGTTCCGCTTTAAAGATGAAGCTGATGTGATTGCGCAAGCCAATg      >  1:1426093/1‑149 (MQ=255)
                                                                                                                                                   cACGAACGCGGCGGCAACTTCTTCCAGCCGACCATTCTGGTGGACGTTCCGGCCAACGCCAAAGTGTCGAAAGAAGAGACGTTCGGCCCCCTCGCCCCGCTGTTCCGCTTTAAAGATGAAGCTGATGTGATTGCGCAAGCCAATGacac  >  1:2020163/1‑149 (MQ=255)
                                                                                                                                                   |                                                                                                                                                    
GCAGTGAGCAAACTGCACATCGGCGACGGGCTGGATAACGGCGTCACCATCGGGCCGCTGATCGATGAAAAAGCGGTAGCAAAAGTGGAAGAGCATATTGCCGATGCGCTGGAGAAAGGCGCGCGCGTGGTTTGCGGCGGTAAAGCGCACGAACGCGGCGGCAACTTCTTCCAGCCGACCATTCTGGTGGACGTTCCGGCCAACGCCAAAGTGTCGAAAGAAGAGACGTTCGGCCCCCTCGCCCCGCTGTTCCGCTTTAAAGATGAAGCTGATGTGATTGCGCAAGCCAATGACAC  >  CP009273/2785565‑2785860

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