Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 2,785,712 C→A 8.7% 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,7120CA 8.7% 1242.7 / 77.9 331H361N (CAC→AAC) gabDsuccinate‑semialdehyde dehydrogenase I, NADP‑dependent
Reads supporting (aligned to +/- strand):  ref base C (162/140);  new base A (19/10);  total (181/150)
Fisher's exact test for biased strand distribution p-value = 2.46e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 335 total aligned reads displayed.

GGCAGTGAGCAAACTGCACATCGGCGACGGGCTGGATAACGGCGTCACCATCGGGCCGCTGATCGATGAAAAAGCGGTAGCAAAAGTGGAAGAGCATATTGCCGATGCGCTGGAGAAAGGCGCGCGCGTGGTTTGCGGCGGTAAAGCGCACGAACGCGGCGGCAACTTCTTCCAGCCGACCATTCTGGTGGACGTTCCGGCCAACGCCAAAGTGTCGAAAGAAGAGACGTTCGGCCCCCTCGCCCCGCTGTTCCGCTTTAAAGATGAAGCTGATGTGATTGCGCAAGCCAATGA  >  CP009273/2785564‑2785857
                                                                                                                                                    |                                                                                                                                                 
ggcagTGAGCAAACTGCACATCGGCGACGGGCTGGATAACGGCGTCACCATCGGGCCGCTGATCGATGAAAAAGCGGTAGCAAAAGTGGAAGAGCATATTGCCGATGCGCTGGAGAAAGGCGCGCGCGTGGTTTGCGGCGGTAAAgcgc                                                                                                                                                   >  1:3594548/1‑149 (MQ=255)
     tGAGCAAACTGCACATCGGCGACGGGCTGGATAACGGCGTCACCATCGGGCCGCTGATCGATGAAAAAGCGGTAGCAAAAGTGGAAGAGCATATTGCCGATGCGCTGGAGAAAGGCGCGCGCGTGGTTTGCGGCGGTAAAGcgaacgaa                                                                                                                                              <  2:1429883/149‑1 (MQ=255)
     tGAGCAAACTGCACATCGGCGACGGGCTGGATAACGGCGTCACCATCGGGCCGCTGATCGATGAAAAAGCGGTAGCAAAAGTGGAAGAGCATATTGCCGATGCGCTGGAGAAAGGCGCGCGCGTGGTTTGCGGCGGTAAAGCGCACGaa                                                                                                                                              <  1:4781115/149‑1 (MQ=255)
     tGAGCAAACTGCACATCGGCGACGGGCTGGATAACGGCGTCACCATCGGGCCGCTGATCGATGAAAAAGCGGTAGCAAAAGTGGAAGAGCATATTGCCGATGCGCTGGAGAAAGGCGCGCGCGTGGTTTGCGGCGGTAAAGCGCACGaa                                                                                                                                              <  1:2437939/149‑1 (MQ=255)
         cAAACTGCACATCGGCGACGGGCTGGATAACGGCGTCACCATCGGGCCGCTGATCGATGAAAAAGCGGTAGCAAAAGTGGAAGAGCATATTGCCGATGCGCTGGAGAAAGGCGCGCGCGTGGTTTGCGGCGGTAAAGCGCACGAAcgcg                                                                                                                                          >  1:2368944/1‑149 (MQ=255)
          aaaCTGCACATCGGCGACGGGCTGGATAACGGCGTCACCATCGGGCCGCTGATCGATGAAAAAGCGGTAGCAAAAGTGGAAGAGCATATTGCCGATGCGCTGGAGAAAGGCGCGCGCGTGGTTTGCGGCGGTAAAGCGCACGAACgcgg                                                                                                                                         >  1:2243868/1‑149 (MQ=255)
            aCTGCACATCGGCGACGGGCTGGATAACGGCGTCACCATCGGGCCGCTGATCGATGAAAAAGCGGTAGCAAAAGTGGAAGAGCATATTGCCGATGCGCTGGAGAAAGGCGCGCGCGTGGTTTGCGGCGGTAAAGCGCACGAACgcggcg                                                                                                                                       >  2:1295138/1‑149 (MQ=255)
            aCTGCACATCGGCGACGGGCTGGATAACGGCGTCACCATCGGGCCGCTGATCGATGAAAAAGCGGTAGCAAAAGTGGAAGAGCATATTGCCGATGCGCTGGAGAAAGGCGCGCGCGTGGTTTGCGGCGGTAAAGCGCACGAACgcggcg                                                                                                                                       >  1:3330456/1‑149 (MQ=255)
              tGCACATCGGCGACGGGCTGGATAACGGCGTCACCATCGGGCCGCTGATCGATGAAAAAGCGGTAGCAAAAGTGGAAGAGCATATTGCCGATGCGCTGGAGAAAGGCGCGCGCGTGGTTTGCGGCGGTAAAGCGCACGAACGCggcggc                                                                                                                                     >  2:1747222/1‑149 (MQ=255)
               gCACATCGGCGACGGGCTGGATAACGGCGTCACCATCGGGCCGCTGATCGATGAAAAAGCGGTAGCAAAAGTGGAAGAGCATATTGCCGATGCGCTGGAGAAAGGCGCGCGCGTGGTTTGCGGCGGTAAAGCGCACGAACGCGGCGGCa                                                                                                                                    <  2:1505459/149‑1 (MQ=255)
               gCACATCGGCGACGGGCTGGATAACGGCGTCACCATCGGGCCGCTGATCGATGAAAAAGCGGTAGCAAAAGTGGAAGAGCATATTGCCGATGCGCTGGAGAAAGGCGCGCGCGTGGTTTGCGGCGGTAAAGCGCACGAACGCGGCGGCa                                                                                                                                    >  1:1428897/1‑149 (MQ=255)
               gCACATCGGCGACGGGCTGGATAACGGCGTCACCATCGGGCCGCTGATCGATGAAAAAGCGGTAGCAAAAGTGGAAGAGCATATTGCCGATGCGCTGGAGAAAGGCGCGCGCGTGGTTTGCGGCGGTAAAGCGCACGAACGCGGCGGCa                                                                                                                                    <  2:2492299/149‑1 (MQ=255)
               gCACATCGGCGACGGGCTGGATAACGGCGTCACCATCGGGCCGCTGATCGATGAAAAAGCGGTAGCAAAAGTGGAAGAGCATATTGCCGATGCGCTGGAGAAAGGCGCGCGCGTGGTTTGCGGCGGTAAAGCGAACGAACGCGGCGGCa                                                                                                                                    >  2:3754561/1‑149 (MQ=255)
                 acaTCGGCGACGGGCTGGATAACGGCGTCACCATCGGGCCGCTGATCGATGAAAAAGCGGTAGCAAAAGTGGAAGAGCATATTGCCGATGCGCTGGAGAAAGGCGCGCGCGTGGTTTGCGGCGGTAAAGCGCACGAACGCGGCGGCAAc                                                                                                                                  >  2:719019/1‑149 (MQ=255)
                  caTCGGCGACGGGCTGGATAACGGCGTCACCATCGGGCCGCTGATCGATGAAAAAGCGGTAGCAAAAGTGGAAGAGCATATTGCCGATGCGCTGGAGAAAGGCGCGCGCGTGGTTTGCGGCGGTAAAGCGCACGAACGCGGCGGCAAct                                                                                                                                 >  2:190595/1‑149 (MQ=255)
                   aTCGGCGACGGGCTGGATAACGGCGTCACCATCGGGCCGCTGATCGATGAAAAAGCGGTAGCAAAAGTGGAAGAGCATATTGCCGATGCGCTGGAGAAAGGCGCGCGCGTGGTTTGCGGCGGTAAAGCGCACGAACGCGGCGGCAActt                                                                                                                                <  2:1968492/149‑1 (MQ=255)
                    tCGGCGACGGGCTGGATAACGGCGTCACCATCGGGCCGCTGATCGATGAAAAAGCGGTAGCAAAAGTGGAAGAGCATATTGCCGATGCGCTGGAGAAAGGCGCGCGCGTGGTTTGCGGCGGTAAAGCGCACGAACGCGGCGGCAActtc                                                                                                                               >  1:4372106/1‑149 (MQ=255)
                       gCGACGGGCTGGATAACGGCGTCACCATCGGGCCGCTGATCGATGAAAAAGCGGTAGCAAAAGTGGACGAGCATATTGCCGATGCGCTGGAGAAAGGCGCGCGCGTGGTTTGCGGCGGTAAAGCGCACGAACGCGGCGGCAACttcttc                                                                                                                            <  1:3908112/149‑1 (MQ=255)
                        cGACGGGCTGGATAACGGCGTCACCATCGGGCCGCTGATCGATGAAAAAGCGGTAGCAAAAGTGGAAGAGCATATTGCCGATGCGCTGGAGAAAGGCGCGCGCGTGGTTTGCGGCGGTAAAGCGCACGAACGCGGCGGCAACTTCTTcc                                                                                                                           <  1:3419186/149‑1 (MQ=255)
                         gACGGGCTGGATAACGGCGTCACCATCGGGCCGCTGATCGATGAAAAAGCGGTAGCAAAAGTGGAAGAGCATATTGCCGATGCGCTGGAGAAAGGCGCGCGCGTGGTTTGCGGCGGTAAAGCGCACGAACGCGGCGGCAACTTCTTCCa                                                                                                                          <  2:3874307/149‑1 (MQ=255)
                            gggCTGGATAACGGCGTCACCATCGGGCCGCTGATCGATGAAAAAGCGGTAGCAAAAGTGGAAGAGCATATTGCCGATGCGCTGGAGAAAGGCGCGCGCGTGGTTTGCGGCGGTAAAGCGCACGAACGCGGCGGCAACTTCTTCCAGcc                                                                                                                       >  2:936334/1‑149 (MQ=255)
                                    tAACGGCGTCACCATCGGGCCGCTGATCGATGAAAAAGCGGTAGCAAAAGTGGAAGAGCATATTGCCGATGCGCTGGAGAAAGGCGCGCGCGTGGTTTGCGGCGGTAAAGCGCACGAACGCGGCGGCAACTTCTTCCAGCCGACCATTc                                                                                                               >  2:4150154/1‑149 (MQ=255)
                                    gcccGGCGTCACCATCGGGCCGCTGATCGATGACAAAGCGGTAGCAAAAGGCGACGAGCATATTGCCGCTGCGCTGGAGAAAGGCGCGCGCGTGGATTGCGGCGGTAAAGCGCACGAACGCGGCGGCAACTTCTTCCAGCCGACCATTc                                                                                                               <  2:4956228/146‑1 (MQ=255)
                                      aCGGCGTCACCATCGGGCCGCTGATCGATGAAAAAGCGGTAGCAAAAGTGGAAGAGCATATTGCCGATGCGCTGGAGAAAGGCGCGCGCGTGGTTTGCGGCGGTAAAGCGCACGAACGCGGCGGCAACTTCTTCCAGCCGACCATTCtg                                                                                                             <  1:4707839/149‑1 (MQ=255)
                                       cGGCGTCACCATCGGGCCGCTGATCGATGAAAAAGCGGTAGCAAAAGTGGAAGAGCATATTGCCGATGCGCTGGAGAAAGGCGCGCGCGTGGTTTGCGGCGGTAAAGCGCACGAACGCGGCGGCAACTTCTTCCAGCCGACCATTCtgg                                                                                                            <  1:4078014/149‑1 (MQ=255)
                                       cGGCGTCACCATCGGGCCGCTGATCGATGAAAAAGCGGTAGCAAAAGTGGAAGAGCATAGTGCCGATGCGCTGGAGAAAGGCGCGCGCGTGGTTTGCGGCGGTAAAGCGCACGAACGCGGCGGCAACTTCTTCCAGCCGACCATTCtgg                                                                                                            <  2:1869744/149‑1 (MQ=255)
                                           gTCACCATCGGGCCGCTGATCGATGAAAAAGCGGTAGCAAAAGTGGAAGAGCATATTGCCGATGCGCTGGAGAAAGGCGCGCGCGTGGTTTGCGGCGGTAAAGCGCACGAACGCGGCGGCAACTTCTTCCAGCCGACCATTCTGGTGGa                                                                                                        <  2:1160042/149‑1 (MQ=255)
                                           gTCACCATCGGGCCGCTGATCGATGAAAAAGCGGTAGCAAAAGTGGAAGAGCATATTGCCGATGCGCTGGAGAAAGGCGCGCGCGTGGTTTGCGGCGGTAAAGCGCACGAACGCGGCGGCAACTTCTTCCAGCCGACCATTCTGGTGGa                                                                                                        >  1:1913379/1‑149 (MQ=255)
                                              aCCATCGGGCCGCTGATCGATGAAAAAGCGGTAGCAAAAGTGGAAGAGCATATTGCCGATGCGCTGGAGAAAGGCGCGCGCGTGGTTTGCGGCGGTAAAGCGCACGAACGCGGCGGCAACTTCTTCCAGCCGACCATTCTGGTGGACGt                                                                                                     <  2:852087/149‑1 (MQ=255)
                                              aCCATCGGGCCGCTGATCGATGAAAAAGCGGTAGCAAAAGTGGAAGAGCATATTGCCGATGCGCTGGAGAAAGGCGCGCGCGTGGTTTGCGGCGGTAAAGCGCACGAACGCGGCGGCAACTTCTTCCAGCCGACCATTCTGGTGGACGt                                                                                                     >  1:4520319/1‑149 (MQ=255)
                                              aCCATCGGGCCGCTGATCGATGAAAAAGCGGTAGCAAAAGTGGAAGAGCATATTGCCGATGCGCTGGAGAAAGGCGCGCGCGTGGTTTGCGGCGGTAAAGCGAACGAACGCGGCGGCAACTTCTTCCAGCCGACCAGTCTGGTGGACGt                                                                                                     >  2:4195475/1‑149 (MQ=255)
                                              aCCATCGGGCCGCTGATCGATGAAAAAGCGGGAGCAAAAGTGGAAGAGCATATTGCCGATGCGCTGGAGAAAGGCGCGCGCGTGGTTTGCGGCGGTAAAGCGCACGAACGCGGCGGCAACTTCTTCCAGCCGACCATTCTGGTGGACGt                                                                                                     <  1:1549241/149‑1 (MQ=255)
                                                cATCGGGCCGCTGATCGATGAAAAAGCGGTAGCAAAAGTGGAAGAGCATATTGCCGATGCGCTGGAGAAAGGCGCGCGCGTGGTTTGCGGCGGTAAAGCGCACGAACGCGGCGGCAACTTCTTCCAGCCGACCATTCTGGTGGACGTTc                                                                                                   >  2:1099008/1‑149 (MQ=255)
                                                 aTCGGGCCGCTGATCGATGAAAAAGCGGTAGCAAAAGTGGAAGAGCATATTGCCGATGCGCTGGAGAAAGGCGCGCGCGTGGTTTGCGGCGGTAAAGCGCACGAACGCGGCGGCAACTTCTTCCAGCCGACCATTCTGGTGGACGTTcc                                                                                                  <  2:3823356/149‑1 (MQ=255)
                                                    gggCCGCTGATCGATGAAAAAGCGGTAGCAAAAGTGGAAGAGCATATTGCCGATGCGCTGGAGAAAGGCGCGCGCGTGGTTTGCGGCGGTAAAGCGCACGAACGCGGCGGCAACTTCTTCCAGCCGACCATTCTGGTGGACGTTCCGGc                                                                                               >  2:5033209/1‑149 (MQ=255)
                                                     ggCCGCTGATCGATGAAAAAGCGGTAGCAAAAGTGGAAGAGCATATTGCCGATGCGCTGGAGAAAGGCGCGCGCGTGGTTTGCGGCGGTAAAGCGAACGAACGCGGCGGCAACTTCTTCCAGCCGACCATTCTGGTGGACGTTCCGGcc                                                                                              <  2:526565/149‑1 (MQ=255)
                                                          cTGATCGATGAAAAAGCGGTAGCAAAAGTGGAAGAGCATATTGCCGATGCGCTGGAGAAAGGCGCGCGCGTGGTTTGCGGCGGTAAAGCGCACGAACGCGGCGGCAACTTCTTCCAGCCGACCATTCTGGTGGACGTTCCGGCCAACGc                                                                                         <  1:2099342/149‑1 (MQ=255)
                                                            gATCGATGAAAAAGCGGTAGCAAAAGTGGAAGAGCATATTGCCGATGCGCTGGAGAAAGGCGCGCGCGTGGTTTGCGGCGGTAAAGCGCACGAACGCGGCGGCAACTTCTTCCAGCCGACCATTCTGGTGGACGTTCCGGCCAACGCCa                                                                                       >  1:744952/1‑149 (MQ=255)
                                                               cGATGAAAAAGCGGTAGCAAAAGTGGAAGAGCATATTGCCGATGCGCTGGAGAAAGGCGCGCGCGTGGTTTGCGGCGGTAAAGCGCACGAACGCGGCGGCAACTTCTTCCAGCCGACCATTCTGGTGGACGTTCCGGCCAACGCCAAAg                                                                                    >  2:3656797/1‑149 (MQ=255)
                                                                gATGAAAAAGCGGTAGCAAAAGTGGAAGAGCATATTGCCGATGCGCTGGAGAAAGGCGCGCGCGTGGTTTGCGGCGGTAAAGCGCACGAACGCGGCGGCAACTTCTTCCAGCCGACCATTCTGGTGGACGTTCCGGCCAACGCCAAAgt                                                                                   >  2:3329265/1‑149 (MQ=255)
                                                                gATGAAAAAGCGGTAGCAAAAGTGGAAGAGCATATTGCCGATGCGCTGGAGAAAGGCGCGCGCGTGGTTTGCGGCGGTAAAGCGCACGAACGCGGCGGCAACTTCTTCCAGCCGACCATTCTGGTGGACGTTCCGGCCAACGCCAAAgt                                                                                   <  2:1104860/149‑1 (MQ=255)
                                                                 aTGAAAAAGCGGTAGCAAAAGTGGAAGAGCATATTGCCGATGCGCTGGAGAAAGGCGCGCGCGTGGTTTGCGGCGGTAAAGCGCACGAACGCGGCGGCAACTTCTTCCAGCCGACCATTCTGGTGGACGTTCCGGCCAACGCCAAAgtg                                                                                  >  1:4600766/1‑149 (MQ=255)
                                                                 aTGAAAAAGCGGTAGCAAAAGTGGAAGAGCATATTGCCGATGCGCTGGAGAAAGGCGCGCGCGTGGTTTGCGGCGGTAAAGCGCACGAACGCGGCGGCAACTTCTTCCAGCCGACCATTCTGGTGGACGTTCCGGCCAACGCCAAAgtg                                                                                  >  1:1614296/1‑149 (MQ=255)
                                                                 aTGAAAAAGCGGTAGCAAAAGTGGAAGAGCATATTGCCGATGCGCTGGAGAAAGGCGCGCGCGTGGTTTGCGGCGGTAAAGCGCACGAACGCGGCGGCAACTTCTTCCAGCCGACCATTCTGGTGGACGTTCCGGCCAACGCCAAAgtg                                                                                  >  1:1358313/1‑149 (MQ=255)
                                                                     aaaaGCGGTAGCAAAAGTGGAAGAGCATATTGCCGATGCGCTGGAGAAAGGCGCGCGCGTGGTTTGCGGCGGTAAAGCGCACGAACGCGGCGGCAACTTCTTCCAGCCGACCATTCTGGTGGACGTTCCGGCCAACGCCAAAGTGTCGa                                                                              >  1:4827933/1‑149 (MQ=255)
                                                                     aaaaGCGGTAGCAAAAGTGGAAGAGCATATTGCCGATGCGCTGGAGAAAGGCGCGCGCGTGGTTTGCGGCGGTAAAGCGCACGAACGCGGCGGCAACTTCTTCCAGCCGACCATTCTGGTGGACGTTCCGGCCAACGCCAAAGTGTCGa                                                                              >  1:551083/1‑149 (MQ=255)
                                                                       aaGCGGTAGCAAAAGTGGAAGAGCATATTGCCGATGCGCTGGAGAAAGGCGCGCGCGTGGTTTGCGGCGGTAAAGCGCACGAACGCGGCGGCAACTTCTTCCAGCCGACCATTCTGGTGGACGTTCCGGCCAACGCCAAAGTGTCGaaa                                                                            >  2:1382894/1‑149 (MQ=255)
                                                                       aaGCGGTAGCAAAAGTGGAAGAGCATATTGCCGATGCGCTGGAGAAAGGCGCGCGCGTGGTTTGCGGCGGTAAAGCGCACGAACGCGGCGGCAACTTCTTCCAGCCGACCATTCTGGTGGACGTTCCGGCCAACGCCAAAGTGTCGaaa                                                                            >  1:2127033/1‑149 (MQ=255)
                                                                       aaGCGGTAGCAAAAGTGGAAGAGCATATTGCCGATGCGCTGGAGAAAGGCGCGCGCGTGGTTTGCGGCGGTAAAGCGCACGAACGCGGCGGCAACTTCTTCCAGCCGACCATTCTGGTGGACGTTCCGGCCAACGCCAAAGTGTCGaaa                                                                            <  1:4178361/149‑1 (MQ=255)
                                                                        aGCGGTAGCAAAAGTGGAAGAGCATATTGCCGATGCGCTGGAGAAAGGCGCGCGCGTGGTTTGCGGCGGTAAAGCGCACGAACGCGGCGGCAACTTCTTCCAGCCGACCATTCTGGTGGACGTTCCGGCCAACGCCAAAGTGTCGAaag                                                                           >  1:2919135/1‑149 (MQ=255)
                                                                        aGCGGTAGCAAAAGTGGAAGAGCATATTGCCGATGCGCTGGAGAAAGGCGCGCGCGTGGTTTGCGGCGGTAAAGCGCACGAACGCGGCGGCAACTTCTTCCAGCCGACCATTCTGGTGGACGTTCCGGCCAACGCCAAAGTGTCGAaag                                                                           >  1:1868123/1‑149 (MQ=255)
                                                                         gCGGTAGCAAAAGTGGAAGAGCATATTGCCGATGCGCTGGAGAAAGGCGCGCGCGTGGTTTGCGGCGGTAAAGCGCACGAACGCGGCGGCAACTTCTTCCAGCCGACCATTCTGGTGGACGTTCCGGCCAACGCCAAAGTGTCGAaaga                                                                          >  1:1012207/1‑149 (MQ=255)
                                                                           ggTAGCAAAAGTGGAAGAGCATATTGCCGATGCGCTGGAGAAAGGCGCGCGCGTGGTTTGCGGCGGTAAAGCGCACGAACGCGGCGGCAACTTCTTCCAGCCGACCATTCTGGTGGACGTTCCGGCCAACGCCAAAGTGTCGAaagaag                                                                        <  2:813429/149‑1 (MQ=255)
                                                                           ggTAGCAAAAGTGGAAGAGCATATTGCCGATGCGCTGGAGAAAGGCGCGCGCGTGGTTTGCGGCGGTAAAGCGCACGAACGCGGCGGCAACTTCTTCCAGCCGACCATTCTGGTGGACGTTCCGGCCAACGCCAAAGTGTCGAaagaag                                                                        >  1:351138/1‑149 (MQ=255)
                                                                           ggTAGCAAAAGTGGAAGAGCATATTGCCGATGCGCTGGAGAAAGGCGCGCGCGTGGTTTGCGGCGGTAAAGCGCACGAACGCGGCGGCAACTTCTTCCAGCCGACCATTCTGGTGGACGTTCCGGCCAACGCCAAAGTGTCGAaagaag                                                                        <  2:4850794/149‑1 (MQ=255)
                                                                            gTAGCAAAAGTGGAAGAGCATATTGCCGATGCGCTGGAGAAAGGCGCGCGCGTGGTTTGCGGCGGTAAAGCGAACGAACGCGGCGGCAACTTCTTCCAGCCGACCATTCTGGTGGACGTTCCGGCCAACGCCAAAGTGTCGAAagaaga                                                                       >  2:26976/1‑149 (MQ=255)
                                                                             tAGCAAAAGTGGAAGAGCATATTGCCGATGCGCTGGAGAAAGGCGCGCGCGTGGTTTGCGGCGGTAAAGCGAACGAACGCGGCGGCAACTTCTTCCAGCCGACCATTCTGGTGGACGTTCCGGCCAACGCCAAAGTGTCGAAAGAagag                                                                      <  2:3587585/149‑1 (MQ=255)
                                                                               gCAAAAGTGGAAGAGCATATTGCCGATGCGCTGGAGAAAGGCGCGCGCGTGGTTTGCGGCGGTAAAGCGCACGAACGCGGCGGCAACTTCTTCCAGCCGACCATTCTGGTGGACGTTCCGGCCAACTCCAAAGTGTCGAAAGAAGAGAc                                                                    >  2:3395536/1‑149 (MQ=255)
                                                                               gCAAAAGTGGAAGAGCATATTGCCGATGCGCTGGAGAAAGGCGCGCGCGTGGTTTGCGGCGGTAAAGCGAACGAACGCGGCGGCAACTTCTTCCAGCCGACCATTCTGGTGGACGTTCCGGCCAACGCCAAAGTGTCGAAAGAAGAGAc                                                                    >  2:4681497/1‑149 (MQ=255)
                                                                                cAAAAGTGGAAGAGCATATTGCCGATGCGCTGGAGAAAGGCGCGCGCGTGGTTTGCGGCGGTAAAGCGCACGAACGCGGCGGCAACTTCTTCCAGCCGACCATTCTGGTGGACGTTCCGGCCAACGCCAAAGTGTCGAAAGAAGAGACg                                                                   <  2:2608235/149‑1 (MQ=255)
                                                                                 aaaaGTGGAAGAGCATATTGCCGATGCGCTGGAGAAAGGCGCGCGCGTGGTTTGCGGCGGTAAAGCGCACGAACGCGGCGGCAACTTCTTCCAGCCGACCATTCTGGTGGACGTTCCGGCCAACGCCAAAGTGTCGAAAGAAGAGACGt                                                                  >  2:613750/1‑149 (MQ=255)
                                                                                  aaaGTGGAAGAGCATATTGCCGATGCGCTGGAGAAAGGCGCGCGCGTGGTTTGCGGCGGTAAAGCGCACGAACGCGGCGGCAACTTCTTCCAGCCGACCATTCTGGTGGACGTTCCGGCCAACGCCAAAGTGTCGAAAGAAGAGACGtt                                                                 >  2:1613746/1‑149 (MQ=255)
                                                                                  aaaGTGGAAGAGCATATTGCCGATGCGCTGGAGAAAGGCGCGCGCGTGGTTTGCGGCGGTAAAGCGCACGAACGCGGCGGCAACTTCTTCCAGCCGACCATTCTGGTGGACGTTCCGGCCAACGCCAAAGTGTCGAAAGAAGAGACGtt                                                                 >  2:3039499/1‑149 (MQ=255)
                                                                                   aaGTGGAAGAGCATATTGCCGATGCGCGGGAGAAAGGCGCGCGCGTGGTTTGCGGCGGTAAAGCGCACGAACGCGGCGGCAACTTCTTCCAGCCGACCATTCTGGTGGACGTTCCGGCCCACGCCAAAGTGTCGAAAGAAGAGACGtgc                                                                >  2:1690665/1‑147 (MQ=255)
                                                                                    aGTGGAAGAGCATATTGCCGATGCGCTGGAGAAAGGCGCGCGCGTGGTTTGCGGCGGTAAAGCGCACGAACGCGGCGGCAACTTCTTCCAGCCGACCATTCTGGTGGACGTTCCGGCCAACGCCAAAGTGTCGAAAGAAGAGACGTTCg                                                               >  2:2262756/1‑149 (MQ=255)
                                                                                      tGGAAGAGCATATTGCCGATGCGCTGGAGAAAGGCGCGCGCGTGGTTTGCGGCGGTAAAGCGCACGAACGCGGCGGCAACTTCTTCCAGCCGACCATTCTGGTGGACGTTCCGGCCAACGCCAAAGTGTCGAAAGAAGAGACGTTCGGc                                                             >  1:3654157/1‑149 (MQ=255)
                                                                                       ggAAGAGCATATTGCCGATGCGCTGGAGAAAGGCGCGCGCGTGGTTTGCGGCGGTAAAGCGCACGAACGCGGCGGCAACTTCTTCCAGCCGACCATTCTGGTGGACGTTCCGGCCAACGCCAAAGTGTCGAAAGAAGAGACGTTCGGcc                                                            >  2:4211445/1‑149 (MQ=255)
                                                                                         aaGAGCATATTGCCGATGCGCTGGAGAAAGGCGCGCGCGTGGTTTGCGGCGGTAAAGCGCACGAACGCGGCGGCAACTTCTTCCAGCCGACCATTCTGGTGGACGTTCCGGCCAACGCCAAAGTGTCGAAAGAAGAGACGTTCGGcccc                                                          <  2:4372106/149‑1 (MQ=255)
                                                                                          agagCATATTGCCGATGCGCTGGAGAAAGGCGCGCGCGTGGTTTGCGGCGGTAAAGCGCACGAACGCGGCGGCAACTTCTTCCAGCCGACCATTCTGGTGGACGTTCCGGCCAACGCCACAGTGTCGAAAGACGAGACGTTCGGccccc                                                         >  2:460715/1‑149 (MQ=255)
                                                                                           gagCATATTGCCGATGCGCTGGAGAAAGGCGCGCGCGTGGTTTGCGGCGGTAAAGCGCACGAACGCGGCGGCAACTTCTTCCAGCCGACCATTCTGGTGGACGTTCCGGCCAACGCCAAAGTGTCGAAAGAAGAGACGTTCGGCCCCCt                                                        <  2:2404927/149‑1 (MQ=255)
                                                                                           gagCATATTGCCGATGCGCTGGAGAAAGGCGCGCGCGTGGTTTGCGGCGGTAAAGCGCACGAACGCGGCGGCAACTTCTTCCAGCCGACCATTCTGGTGGACGTTCCGGCCAACGCCAAAGTGTCGAAAGAAGAGACGTTCGGCCCCCt                                                        <  1:2687807/149‑1 (MQ=255)
                                                                                           gagCATATTGCCGATGCGCTGGAGAAAGGCGCGCGCGTGGTTTGCGGCGGTAAAGCGCACGAACGCGGCGGCAACTTCTTCCAGCCGACCATTCTGGTGGACGTTCCGGCCAACGCCAAAGTGTCGAAAGAAGAGACGTTCGGCCCCCt                                                        <  1:3018584/149‑1 (MQ=255)
                                                                                           gagCATATTGCCGATGCGCTGGAGAAAGGCGCGCGCGTGGTTTGCGGCGGTAAAGCGCACGAACGCGGCGGCAACTTCTTCCAGCCGACCATTCTGGTGGACGTTCCGGCCAACGCCAAAGTGTCGAAAGAAGAGACGTTCGGCCCCCt                                                        <  1:926634/149‑1 (MQ=255)
                                                                                              cATATTGCCGATGCGCTGGAGAAAGGCGCGCGCGTGGTTTGCGGCGGTAAAGCGCACGAACGCGGCGGCAACTTCTTCCAGCCGACCATTCTGGTGGACGTTCCGGCCAACGCCAAAGTGTCGAAAGAAGAGACGTTCGGCCCCCTCGc                                                     >  1:31881/1‑149 (MQ=255)
                                                                                              cATATTGCCGATGCGCTGGAGAAAGGCGCGCGCGTGGTTTGCGGCGGTAAAGCGCACGAACGCGGCGGCAACTTCTTCCAGCCGACCATTCTGGTGGACGTTCCGGCCAACGCCAAAGTGTCGAAAGAAGAGACGTTCGGCCCCCTCGc                                                     >  1:4284107/1‑149 (MQ=255)
                                                                                                tatTGCCGATGCGCTGGAGAAAGGCGCGCGCGTGGTTTGCGGCGGTAAAGCGCACGAACGCGGCGGCAACTTCTTCCAGCCGACCATTCTGGTGGACGTTCCGGCCAACGCCAAAGTGTCGAAAGAAGAGACGTTCGGCCCCCTCGccc                                                   <  2:5080344/149‑1 (MQ=255)
                                                                                                tagtGCCGATGCGCTGGAGAAAGGCGCGCGCGTGGTTTGCGGCGGTAAAGCGCACGAACGCGGCGGCAACTTCTTCCAGCCGACCATTCTGGTGGACGTTCCGGCCAACGCCAAAGTGTCGAAAGAAGAGACGTTCGGCCCCCTCGccc                                                   <  2:44430/146‑1 (MQ=255)
                                                                                                   tGCCGATGCGCTGGAGAAAGGCGCGCGCGTGGTTTGCGGCGGTAAAGCGCACGAACGCGGCGGCAACTTCTTCCAGCCGACCATTCTGGTGGACGTTCCGGCCAACGCCAAAGTGTCGAAAGAAGAGACGTTCGGCCCCCTCGCCCCGc                                                <  1:5117159/149‑1 (MQ=255)
                                                                                                    gCCGATGCGCTGGAGAAAGGCGCGCGCGTGGTTTGCGGCGGTAAAGCGCACGAACGCGGCGGCAACTTCTTCCAGCCGACCATTCTGGTGGACGTTCCGGCCAACGCCAAAGTGTCGAAAGAAGAGACGTTCGGCCCCCTCGCCCCGCt                                               >  1:3155122/1‑149 (MQ=255)
                                                                                                      cGATGCGCTGGAGAAAGGCGCGCGCGTGGTTTGCGGCGGTAAAGCGCACGAACGCGGCGGCAACTTCTTCCAGCCGACCATTCTGGTGGACGTTCCGGCCAACGCCAAAGTGTCGAAAGAAGAGACGTTCGGCCCCCTCGCCCCGCTGt                                             <  2:1847733/149‑1 (MQ=255)
                                                                                                        aTGCGCTGGAGAAAGGCGCGCGCGTGGTTTGCGGCGGTAAAGCGCACGAACGCGGCGGCACCTTCTTCCAGCCGACCATTCTGGTGGACGTTCCGGCCAACGCCAAAGTGTCGAAAGCAGAGACGTTCGGCCCCCTCGCCCCGCTGTTc                                           >  2:3707030/1‑149 (MQ=255)
                                                                                                         tGCGCTGGAGAAAGGCGCGCGCGTGGTTTGCGGCGGTAAAGCGCACGAACGCGGCGGCAACTTCTTCCAGCCGACCATTCTGGTGGACGTTCCGGCCAACGCCAAAGTGTCGAAAGAAGAGACGTTCGGCCCCCTCGCCCCGCTGTTcc                                          >  2:4050471/1‑149 (MQ=255)
                                                                                                           cgcTGGAGAAAGGCGCGCGCGTGGTTTGCGGCGGTAAAGCGCACGAACGCGGCGGCAACTTCTTCCAGCCGACCATTCTGGTGGACGTTCCGGCCAACGCCAAAGTGTCGAAAGAAGAGACGTTCGGCCCCCTCGCCCCGCTGTTCCGc                                        <  2:953971/149‑1 (MQ=255)
                                                                                                            gcTGGAGAAAGGCGCGCGCGTGGTTTGCGGCGGTAAAGCGCACGAACGCGGCGGCAACTTCTTCCAGCCGACCATTCTGGTGGACGTTCCGGCCAACGCCAAAGTGTCGAAAGAAGAGACGTTCGGCCCCCTCGCCCCGCTGTTCCGCt                                       <  2:3783766/149‑1 (MQ=255)
                                                                                                             cTGGAGAAAGGCGCGCGCGTGGTTTGCGGCGGTAAAGCGCACGAACGCGGCGGCAACTTCTTCCAGCCGACCATTCTGGTGGACGTTCCGGCCAACGCCAAAGTGTCGAAAGAAGAGACGTTCGGCCCCCTCGCCCCGCTGTTCCGCtt                                      <  1:3030482/149‑1 (MQ=255)
                                                                                                                gagaAAGGCGCGCGCGTGGTTTGCGGCGGTAAAGCGCACGAACGCGGCGGCAACTTCTTCCAGCCGACCATTCTGGTGGACGTTCCGGCCAACGCCAAAGTGTCGAAAGAAGAGACGTTCGGCCCCCTCGCCCCGCTGTTCCGCTTTaa                                   >  2:4476127/1‑149 (MQ=255)
                                                                                                                gagaAAGGCGCGCGCGTGGTTTGCGGCGGTAAAGCGCACGAACGCGGCGGCAACTTCTTCCAGCCGACCATTCTGGTGGACGTTCCGGCCAACGCCAAAGTGTCGAAAGAAGAGACGTTCGGCCCCCTCGCCCCGCTGTTCCGCTTTaa                                   >  1:3245581/1‑149 (MQ=255)
                                                                                                                 agaAAGGCGCGCGCGTGGTTTGCGGCGGTAAAGCGAACGAACGCGGCGGCAACTTCTTCCAGCCGACCATTCTGGTGGACGTTCCGGCCAACGCCAAAGTGTCGAAAGAAGAGACGTTCGGCCCCCTCGCCCCGCTGTTCCGCTTTaaa                                  <  2:3102504/149‑1 (MQ=255)
                                                                                                                  gaAAGGCGCGCGCGTGGTTTGCGGCGGTAAAGCGCACGAACGCGGCGGCAACTTCTTCCAGCCGACCATTCTGGTGGACGTTCCGGCCAACGCCAAAGTGTCGAAAGAAGAGACGTTCGGCCCCCTCGCCCCGCTGTTCCGCTTTAAAg                                 >  1:4425910/1‑149 (MQ=255)
                                                                                                                   acaGGCGCGCGCGTGGTTTGCGGCGGTAAAGCGCACGAACGCGGCGGCAACTTCTTCCAGCCGACCATTCTGGTGGACGTTCCGGCCAACGCCAAAGTGTCGAAAGAAGAGACGTTCGGCCCCCTCGCCCCGCTGTTCCGCTTTAAAGa                                <  2:2190674/147‑1 (MQ=255)
                                                                                                                   aaaGGCGCGCGCGTGGTTTGCGGCGGTAAAGCGCACGAACGCGGCGGCAACTTCTTCCAGCCGACCATTCTGGTGGACGTTCCGGCCAACGCCAAAGTGTCGAAAGAAGAGACGTTCGGCCCCCTCGCCCCGCTGTTCCGCTTTAAAGa                                <  1:1885589/149‑1 (MQ=255)
                                                                                                                    aaGGCGCGCGCGTGGTTTGCGGCGGTAAAGCGCACGAACGCGGCGGCAACTTCTTCCAGCCGACCATTCTGGTGGACGTTCCGGCCAACGCCAAAGTGTCGAAAGAAGAGACGTTCGGCCCCCTCGCCCCGCTGTTCCGCTTTAAAGAt                               >  1:413887/1‑149 (MQ=255)
                                                                                                                     aGGCGCGCGCGTGGTTTGCGGCGGTAAAGCGCACGAACGCGGCGGCAACTTCTTCCAGCCGACCATTCTGGTGGACGTTCCGGCCAACGCCAAAGTGTCGAAAGAAGAGACGTTCGGCCCCCTCGCCCCGCTGTTCCGCTTTAAAGATg                              >  1:4927113/1‑149 (MQ=255)
                                                                                                                         gcgcgcgTGGTTTGCGGCGGTAAAGCGCACGAACGCGGCGGCAACTTCTTCCAGCCGACCATTCTGGTGGACGTTCCGGCCAACGCCAAAGTGTCGAAAGAAGAGACGTTCGGCCCCCTCGCCCCGCTGTTCCGCTTTAAAGATGAAGc                          >  2:3474422/1‑149 (MQ=255)
                                                                                                                          cgcgcgTGGTTTGCGGCGGTAAAGCGAACGAACGCGGCGGCAACTTCTTCCAGCCGACCATTCTGGTGGACGTTCCGGCCAACGCCAAAGTGTCGAAAGAAGAGACGTTCGGCCCCCTCGCCCCGCTGTTCCGCTTTAAAGATGAAGCt                         >  2:2788413/1‑149 (MQ=255)
                                                                                                                                     tGCGGCGGTAAAGCGCACGAACGCGGCGGCAACTTCTTCCAGCCGACCATTCTGGTGGACGTTCCGGCCAACGCCAAAGTGTCGAAAGAAGAGACGTTCGGCCCCCTCGCCCCGCTGTTCCGCTTTAAAGATGAAGCTGATGTGATTgc              <  2:3929437/149‑1 (MQ=255)
                                                                                                                                         gcggTAAAGCGCACGAACGCGGCGGCAACTTCTTCCAGCCGACCATTCTGGTGGACGTTCCGGCCAACGCCAAAGTGTCGAAAGAAGAGACGTTCGGCCCCCTCGCCCCGCTGTTCCGCTTTAAAGATGAAGCTGATGTGATTGCGCaa          >  2:2146408/1‑149 (MQ=255)
                                                                                                                                           ggTAAAGCGCACGAACGCGGCGGCAACTTCTTCCAGCCGACCATTCTGGTGGACGTTCCGGCCAACGCCAAAGTGTCGAAAGAAGAGTCGTTCGGCCCCCTCGCCCCGCTGTTCCGCTTTAAAGATGAAGCTGATGTGATTGCGCAAGc        <  1:96972/149‑1 (MQ=255)
                                                                                                                                               aaGCGAACGAACGCGGCGGCAACTTCTTCCAGCCGACCATTCTGGTGGACGTTCCGGCCAACGCCAAAGTGTCGAAAGAAGAGACGTTCGGCCCCCTCGCCCCGCTGTTCCGCTTTAAAGATGAACCTGATGTGATAGCGCAAGCCAAt    >  2:4325833/1‑149 (MQ=255)
                                                                                                                                                aGCGCACGAACGCGGCGGCAACTTCTTCCAGCCGACCATTCTGGTGGACGTTCCGGCCAACGCCAAAGTGTCGAAAGAAGAGACGTTCGGCCCCCTCGCCCCGCTGTTCCGCTTTAAAGATGAAGCTGATGTGATTGCGCAAGCCAATg   >  1:1207356/1‑149 (MQ=255)
                                                                                                                                                    |                                                                                                                                                 
GGCAGTGAGCAAACTGCACATCGGCGACGGGCTGGATAACGGCGTCACCATCGGGCCGCTGATCGATGAAAAAGCGGTAGCAAAAGTGGAAGAGCATATTGCCGATGCGCTGGAGAAAGGCGCGCGCGTGGTTTGCGGCGGTAAAGCGCACGAACGCGGCGGCAACTTCTTCCAGCCGACCATTCTGGTGGACGTTCCGGCCAACGCCAAAGTGTCGAAAGAAGAGACGTTCGGCCCCCTCGCCCCGCTGTTCCGCTTTAAAGATGAAGCTGATGTGATTGCGCAAGCCAATGA  >  CP009273/2785564‑2785857

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: ‑