Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 2,920,050 C→A 87.5% I128I (ATC→ATA)  ygdH → UPF0717 family protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092732,920,0500CA87.5% 480.5 / 71.4 206I128I (ATC→ATA) ygdHUPF0717 family protein
Reads supporting (aligned to +/- strand):  ref base C (17/9);  new base A (94/86);  total (111/95)
Fisher's exact test for biased strand distribution p-value = 2.93e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.32e-01

Only 100 of 207 total aligned reads displayed.

GGAAGAGGCTTTTGTCGATGGGCGAATTATTCGCGCTTTGCAGGCCAACTTGTTCGCGGTCCTGCGTGACATTCTCTTCGTTTACGGGCAAATCCATAACACCGTTCGTTTTCCCAACCTGAATCTCGACAACTCCGTCCACATCACTAACCTGGTCTTTTCCATCTTGCGTAACGCTCGCGCGCTGCATGTGGGTGAAGCGCCAAATATGGTGGTCTGCTGGGGCGGTCACTCAATTAACGAAAACGAGTATTTGTATGCCCGTCGCGTCGGAAACCAGCTGGGCCTG  >  CP009273/2919906‑2920194
                                                                                                                                                |                                                                                                                                                
ggAAGAGGCTTTTGTCTATGGGCGAATTATTCGCGCTTTGCAGGCCAACTTGTTCGCGGTCCTGCGTGACATTCTCTTCGTTTACGGGCAAATCCATAACACCGTTCGTTTTCCCAACCTGAATCTCGACAACTCCGTCCACATAACt                                                                                                                                               <  2:231151/148‑1 (MQ=255)
 gAAGAGGCTTTTGTCTATGGGCGAATTATTCGCGCTTTGCAGGCCAACTTGTTCGCGGTCCTGCGTGACATTCTCTTCGTTTACGGGCAAATCCATAACACCGTTCGTTTTCCCAACCTGAATCTCGACAACTCCGTCCACATAACTa                                                                                                                                              <  1:1063958/148‑1 (MQ=255)
  aaGAGGCTTTTGTCTATGGGCGAATTATTCGCGCTTTGCAGGCCAACTTGTTCGCGGTCCTGCGTGACATTCTCTTCGTTTACGGGCAAATCCATAACACCGTTCGTTTTCCCAACCTGAATCTCGACAACTCCGTCCACATAACTaa                                                                                                                                             >  2:2330578/1‑148 (MQ=255)
  aaGAGGCTTTTGTCTATGGGCGAATTATTCGCGCTTTGCAGGCCAACTTGTTCGCGGTCCTGCGTGACATTCTCTTCGTTTACGGGAAAATCCATAACACCGTTCGTTTTCCCAACCTGAATCTCGACAACTCCGTCCACATAACTaa                                                                                                                                             >  2:1093048/1‑148 (MQ=255)
  aaGAGGCTTTTGTCTATGGGCGAATTATTCGCGCTTTGCAGGCCAACTTGTTCGCGGTCCTGCGTGACATTCTCTTCGTTTACGGGAAAATCCATAACACCGTTCGTTTTCCCAACCTGAATCTCGACAACTCCGTCCACATAACTaa                                                                                                                                             >  2:1093138/1‑148 (MQ=255)
      ggCTTTTGTCTATGGGCGAATTATTCGCGCTTTGCAGGCCAACTTGTTCGCGGTCCTGCGTGACATTCTCTTCGTTTACGGGCAAATCCATAACACCGTTCGTTTTCCCAACCTGAATCTCGACAACTCCGTCCACATAACTAACCTg                                                                                                                                         <  2:1992562/148‑1 (MQ=255)
      ggCTTTTGTCGATGGGCGAATTATTCGCGCTTTGCAGGCCAACTTGTTCGCGGTCCTGCGTGACATTCTCTTCGTTTACGGGCAAATCCATAACACCGTTCGTTTTCCCAACCTGAATCTCGACAACTCCGTCCACATCACTAACCTg                                                                                                                                         >  2:2134812/1‑148 (MQ=255)
       gCTTTTGTCTATGGGCGAATTATTCGCGCTTTGCAGGCCAACTTGTTCGCGGTCCTGCGTGACATTCTCTTCGTTTACGGGCAAATCCATAACACCGTTCGTTTTCCCAACCTGAATCTCGACAACTCCGTCCACATAACTAACCTgg                                                                                                                                        >  1:672433/1‑148 (MQ=255)
       gCTTTTGTCTATGGGCGAATTATTCGCGCTTTGCAGGCCAACTTGTTCGCGGTCCTGCGTGACATTCTCTTCGTTTACGGGCAAATCCATAACACCGTTCGTTTTCCCAACCTGAATCTCGACAACTCCGTCCACATAACTAACCTgg                                                                                                                                        >  1:1996702/1‑148 (MQ=255)
       gCTTTTGTCGATGGGCGAATTATTCGCGCTTTGCAGGCCAACTTGTTCGCGGTCCTGCGTGACATTCTCTTCGTTTACGGGCAAATCCATAACACCGTTCGTTTTCCCAACCTGAATCTCGACAACTCCGTCCACATCACTAACCTgg                                                                                                                                        >  1:1872834/1‑148 (MQ=255)
        cTTTTGTCTATGGGCGAATTATTCGCGCTTTGCAGGCCAACTTGTTCGCGGTCCTGCGTGACATTCTCTTCGTTTACGGGCAAATCCATAACACCGTTCGTTTTCCCAACCTGAATCTCGACAACTCCGTCCACATAACTAACCTGGt                                                                                                                                       <  1:2690896/148‑1 (MQ=255)
          tttGTCTATGGGCGAATTATTCGCGCTTTGCAGGCCAACTTGTTCGCGGTCCTGCGTGACATTCTCTTCGTTTACGGGCAAATCCATAACACCGTTCGTTTTCCCAACCTGAATCTCGACAACTCCGTCCACATAACTAACCTGGTCt                                                                                                                                     >  1:713906/1‑148 (MQ=255)
          tttGTCTATGGGCGAATTATTCGCGCTTTGCAGGCCAACTTGTTCGCGGTCCTGCGTGACATTCTCTTCGTTTACGGGCAAATCCATAACACCGTTCGTTTTCCCAACCTGAATCTCGACAACTCCGTCCACATAACTAACCTGGTCt                                                                                                                                     >  2:593018/1‑148 (MQ=255)
          tttGTCGATGGGCGAATTATTCGCGCTTTGCAGGCCAACTTGTTCGCGGTCCTGCGTGACATTCTCTTCGTTTACGGGCAAATCCATAACACCGTTCGTTTTCCCAACCTGAATCTCGACAACTCCGTCCACATCACTAACCTGGTCt                                                                                                                                     <  2:1539920/148‑1 (MQ=255)
              tctaTGGGCGAATTATTCGCGCTTTGCAGGCCTACTTGTTCGCGGTCCTGCGTGACATTCTCTTCGTTTACGGGCAAATCCATAACACCGTTCGTTTTCCCAACCTGAATCTCGACAACTCCGTCCACATAACTAACCTGGTCTTTTc                                                                                                                                 <  2:1006816/145‑1 (MQ=255)
              tctaTGGGCGAATTATTCGCGCTTTGCAGGCCAACTTGTTCGCGGTCCTGCGTGACATTCTCTTCGTTTACGGGCAAATCCATAACACCGTTCGTTTTCCCAACCTGAATCTCGACAACTCCGTCCACATAACTAACCTGGTCTTTTc                                                                                                                                 <  2:532868/145‑1 (MQ=255)
               ctaTGGGCGAATTATTCGCGCTTTGCAGGCCAACTTGTTCGCGGTCCTGCGTGACATTCTCTTCGTTTACGGGCAAATCCATAACACCGTTCGTTTTCCCAACCTGAATCTCGACAACTCCGTCCACCTAACTAACCTGGTCTTTTcc                                                                                                                                >  1:3142189/3‑148 (MQ=255)
               ctaTGGGCGAATTATTCGCGCTTTGCAGGCCAACTTGTTCGCGGTCCTGCGTGACATTCTCTTCGTTTACGGGCAAATCCATAACACCGTTCGTTTTCCCAACCTGAATCTCGACAACTCCGTCCACATAACTAACCTGGTCTTTTcc                                                                                                                                <  2:2712786/146‑1 (MQ=255)
                 aTGGGCGAATTATTCGCGCTTTGCAGGCCAACTTGTTCGCGGTCCTGCGTGACATTCTCTTCGTTTACGGGCAAATCCATAACACCGTTCTTTTTCCCAACCTGAATCTCGACAACTCCGTCCACATAACTAACCTGGTCTTTTCCAt                                                                                                                              >  2:1150561/1‑148 (MQ=255)
                   gggCGAATTATTCGCGCTTTGCAGGCCAACTTGTTCGCGGTCCTGCGTGACATTCTCTTCGTTTACGGGCAAATCCATAACACCGTTCGTTTTCCCAACCTGAATCTCGACAACTCCGTCCACATCACTAACCTGGTCTTTTCCATCt                                                                                                                            >  2:2435977/1‑148 (MQ=255)
                   gggCGAATTATTCGCGCTTTGCAGGCCAACTTGTTCGCGGTCCTGCGTGACATTCTCTTCGTTTACGGGCAAATCCATAACACCGTTCGTTTTCCCAACCTGAATCTCGACAACTCCGTCCACATAACTAACCTGGTCTTTTCCATCt                                                                                                                            >  1:2787963/1‑148 (MQ=255)
                     gCGAATTATTCGCGCTTTGCAGGCCAACTTGTTCGCGGTCCTGCGTGACATTCTCTTCGTTTACGGGCAAATCCATAACACCGTTCGTTTTCCCAACCTGAATCTCGACAACTCCGTCCACATAACTAACCTGGTCTTTTCCATCTTg                                                                                                                          <  2:681937/148‑1 (MQ=255)
                         attattCGCGCTTTGCAGGCCAACTTGTTCGCGGTCCTGCGTGACATTCTCTTCGTTTACGGGCAAATCCATAACACCGTTCGTTTTCCCAACCTGAATCTCGACAACTCCGTCCACATAACTAACCTGGTCTTTTCCATCTTGCGTa                                                                                                                      <  1:2572158/148‑1 (MQ=255)
                              tCGCGCTTTGCAGGCCAACTTGTTCGCGGTCCTGCGTGACATTCTCTTCGTTTACGGGCAAATCCATAACACCGTTCGTTTTCCCAACCTGAATCTCGACAACTCCGTCCACATAACTAACCTGGTCTTTTCCATCTTGCGTAACGCt                                                                                                                 >  2:2642249/1‑148 (MQ=255)
                                    tttGCAGGCCAACTTGTTCGCGGTCCTGCGTGACATTCTCTTCGTTTACGGGCAAATCCATAACACCGTTCGTTTTCCCAACCTGAATCTCGACAACTCCGTCCACATCACTAACCTGGTCTTTTCCATCTTGCGTAACGCTcgcgcg                                                                                                           <  1:1607718/148‑1 (MQ=255)
                                    tttGCAGGCCAACTTGTTCGCGGTCCTGCGTGACATTCTCTTCGTTTACGGGCAAATCCATAACACCGTTCGTTTTCCCAACCTGAATCTCGACAACTCCGTCCACATAACTAACCTGGTCTTTTCCATCTTGCGTAACGCTcgcgcg                                                                                                           <  1:2565338/148‑1 (MQ=255)
                                    tttGCAGGCCAACTTGTTCGCGGTCCTGCGTGACATTCTCTTCGTTTACGGGCAAATCCATAACACCGTTCGTTTTCCCAACCTGAATCTCGACAACTCCGTCCACATAACTAACCTGGTCTTTTCCATCTTGCGTAACGCTcgcgcg                                                                                                           <  2:3012242/148‑1 (MQ=255)
                                    tttGCAGGCCAACTTGTTCGCGGTCCTGCGTGACATTCTCTTCGTTTACGGGCAAATCCATAACACCGTTCGTTTTACCAACCTGAATCTCGACAACTCCGTCCACATCACTAACCTGGTCTTTTCCATCTTGCGTAACGCTcgcgcg                                                                                                           >  1:2334043/1‑148 (MQ=255)
                                     ttGCAGGCCAACTTGTTCGCGGTCCTGCGTGACATTCTCTTCGTTTACGGGCAAATCCATAACACCGTTCGTTTTCCCAACCTGAATCTCGACAACTCCGTCCACATAACTAACCTGGTCTTTTCCATCTTGCGTAACGCTCgcgcgc                                                                                                          <  1:2270353/148‑1 (MQ=255)
                                     ttGCAGGCCAACTTGTTCGCGGTCCTGCGTGACATTCTCTTCGTTTACGGGCAAATCCATAACACCGTTCGTTTTCCCAACCTGAATCTCGACAACTCCGTCCACATAACTAACCTGGTCTTTTCCATCTTGCGTAACGCTCgcgcgc                                                                                                          >  1:1142970/1‑148 (MQ=255)
                                          ggCCAACTTGTTCGCGGTCCTGCGTGACATTCTCTTCGTTTACGGGCAAATCCATAACACCGTTCGTTTTCCCAACCTGAATCTCGACAACTCCGTCCACATAACTAACCTGGTCTTTTCCATCTTGCGTAACGCTCGCGCGCTGCAt                                                                                                     <  1:1022426/148‑1 (MQ=255)
                                          ggCCAACTTGTTCGCGGTCCTGCGTGACATTCTCTTCGTTTACGGGCAAATCCATAACACCGTTCGTTTTCCCAACCTGAATCTCGACAACTCCGTCCACATAACTAACCTGGTCTTTTCCATCTTGCGTAACGCTCGCGCGCTGCAt                                                                                                     <  2:2749479/148‑1 (MQ=255)
                                          ggCCAACTTGTTCGCGGTCCTGCGTGACACTCTCTTCGTTTACGGGCAAATCCATAACACCGTTCGTTTTCCCAACCTGAATCTCGACAACTCCGTCCACATAACTAACCTGGTCTTTTCCATCTTGCGTAACGCTCGCGCGCTGCAt                                                                                                     >  1:1170706/1‑148 (MQ=255)
                                             cAACTTGTTCGCGGTCCTGCGTGACATTCTCTTCGTTTACGGGCAAATCCATAACACCGTTCGTTTTCCCAACCTGAATCTCGACAACTCCGTCCACATAACTAACCTGGTCTTTTCCATCTTGCGTAACGCTCGCGCGCTGCAtgtg                                                                                                  >  2:441899/1‑148 (MQ=255)
                                             cAACTTGTTCGCGGTCCTGCGTGACATTCTCTTCGTTTACGGGCAAATCCATAACACCGTTCGTTTTCCCAACCTGAATCTCGACAACTCCGTCCACATAACTAACCTGGTCTTTTCCATCTTGCGTAACGCTCGCGCGCTGCAtgtg                                                                                                  >  2:807256/1‑148 (MQ=255)
                                              aaCTTGTTCGCGGTCCTGCGTGACATTCTCTTCGTTTACGGGCAAATCCATAACACCGTTCGTTTTCCCAACCTGAATCTCGACAACTCCGTCCACATAACTAACCTGGTCTTTTCCATCTTGCGTAACGCTCGCGCGCTGCATGTgg                                                                                                 >  1:785868/1‑148 (MQ=255)
                                                cTTGTTCGCGGTCCTGCGTGACATTCTCTTCGTTTACGGGCAAATCCATAACACCGTTCGTTTTCCCAACCTGAATCTCGACAACTCCGTCCACATAACTAACCTGGTCTTTTCCATCTTGCGTAACGCTCGCGCGCTGCATGTGGGt                                                                                               <  2:2406154/148‑1 (MQ=255)
                                                  tGTTCGCGGTCCTGCGTGACATTCTCTTCGTTTACGGGCAAATCCATAACACCGTTCGTTTTCCCAACCTGAATCTCGACAACTCCGTCCACATAACTAACCTGGTCTTTTCCATCTTGCGTAACGCTCGCGCGCTGCATGTGGGTGa                                                                                             <  1:64688/148‑1 (MQ=255)
                                                    ttCGCGGTCCTGCGTGACATTCTCTTCGTTTACGGGCAAATCCATAACACCGTTCGTTTTCCCAACCTGAATCTCGACAACTCCGTCCACATAACTAACCTGGTCTTTTCCATCTTGCGTAACGCTCGCGCGCTGCATGTGGGTGAAg                                                                                           >  2:2726815/1‑148 (MQ=255)
                                                       gcgGTCCTGCGTGACATTCTCTTCGTTTACGGGCAAATCCATAACACCGTTCGTTTTCCCAACCTGAATCTCGACAACTCCGTCCACATCACTAACCTGGTCTTTTCCATCTTGCGTAACGCTCGCGCGCTGCATGTGGGTGAAgcgc                                                                                        <  2:2774257/148‑1 (MQ=255)
                                                       gcgGTCCTGCGTGACATTCTCTTCGTTTACGGGCAAATCCATAACACCGTTCGTTTTCCCAACCTGAATCTCGACAACTCCGTCCACATCACTAACCTGGTCTTTTCCATCTTGCGTAACGCTCGCGCGCTGCATGTGGGTGAAgcgc                                                                                        <  2:2774228/148‑1 (MQ=255)
                                                        cgGTCCTGCGTGACATTCTCTTCGTTTACGGGCAAATCCATAACACCGTTCGTTTTCCCAACCTGAATCTCGACAACTCCGTCCACATAACTAACCTGGTCTTTTCCATCTTGCGTAACGCTCGCGCGCTGCATGTGGGTGAAGCGcc                                                                                       >  2:1913338/1‑148 (MQ=255)
                                                           tCCTGCGTGACATTCTCTTCGTTTACGGGCAAATCCATAACACCGTTCGTTTTCCCAACCTGAATCTCGACAACTCCGTCCACATAACTAACCTGGTCTTTTCCATCTTGCGTAACGCTCGCGCGCTGCATGTGGGTGAAGCGCCaaa                                                                                    <  2:1287035/148‑1 (MQ=255)
                                                            ccTGCGTGACATTCTCTTCGTTTACGGGCAAATCCATAACACCGTTCGTTTTCCCAACCTGAATCTCGACAACTCCGTCCACATAACTAACCTGGTCTTTTCCATCTTGCGTAACGCTCGCGCGCTGCATGTGGGTGAAGCGCCAAat                                                                                   <  2:3020953/148‑1 (MQ=255)
                                                               gCGTGACATTCTCTTCGTTTACGGGCAAATCCATAACACCGTTCGTTTTCCCAACCTGAATCTCGACAACTCCGTCCACATAACTAACCTGGTCTTTTCCATCTTGCGTAACGCTCGCGCGCTGCATGTGGGTGAAGCGCCAAATAtg                                                                                <  2:1269760/148‑1 (MQ=255)
                                                               gCGTGACATTCTCTTCGTTTACGGGCAAATCCATAACACCGTTCGTTTTCCCAACCTGAATCTCGACAACTCCGTCCACATAACTAACCTGGTCTTTTCCATCTTGCGTAACGCTCGCGCGCTGCATGTGGGTGAAGCGCCAAATAtg                                                                                >  2:459922/1‑148 (MQ=255)
                                                                  tGACATTCTCTTCGTTTACGGGCAAATCCATAACACCGTTCGTTTTCCCAACCTGAATCTCGACAACTCCGTCCACATCACTAACCTGGTCTTTTCCATCTTGCGTAACGCTCGCGCGCTGCATGTGGGTGAAGCGCCAAATAtggtg                                                                             <  2:1919312/148‑1 (MQ=255)
                                                                  tGACATTCTCTTCGTTTACGGGCAAATCCATAACACCGTTCGTTTTCCCAACCTGAATCTCGACAACTCCGTCCACATACCTAACCTGGTCTTTTCCATCTTGCGTAACGCTCGCGCGCTGCATGTGGGTGAAGCGCCAAATAtggtg                                                                             <  2:373818/148‑1 (MQ=255)
                                                                  tGACATTCTCTTCGTTTACGGGCAAATCCATAACACCGTTCGTTTTCCCAACCTGAATCTCGACAACTCCGTCCACATAACTAACCTGGTCTTTTCCATCTTGCGTAACGCTCGCGCGCTGCATGTGGGTGAAGCGCCAAATAtggtg                                                                             <  2:3254107/148‑1 (MQ=255)
                                                                   gACATTCTCTTCGTTTACGGGCAAATCCATAACACCGTTCGTTTTCCCAACCTGAATCTCGACAACTCCGTCCACATCACTAACCTGGTCTTTTCCATCTTGCGTAACGCTCGCGCGCTGCATGTGGGTGAAGCGCCAAATAtggtgg                                                                            >  1:3022259/1‑148 (MQ=255)
                                                                   gACATTCTCTTCGTTTACGGGCAAATCCATAACACCGTTCGTTTTCCCAACCTGAATCTCGACAACTCCGTCCACATAACTAACCTGGTCTTTTCCATCTTGCGTAACGCTCGCGCGCTGCATGTGGGTGAAGCGCCAAATAtggtgg                                                                            >  1:688505/1‑148 (MQ=255)
                                                                    aCATTCTCTTCGTTTACGGGCAAATCCATAACACCGTTCGTTTTCCCAACCTGAATCTCGACAACTCCGTCCACATAACTAACCTGGTCTTTTCCATCTTGCGTAACGCTCGCGCGCTGCATGTGGGTGAAGCGCCAAATATggtggt                                                                           <  1:244811/148‑1 (MQ=255)
                                                                       ttCTCTTCGTTTACGGGCAAATCCATAACACCGTTCGTTTTCCCAACCTGAATCTCGACAACTCCGTCCACATAACTAACCTGGTCTTTTCCATCTTGCGTAACGCTCGCGCGCTGCATGTGGGTGAAGCGCCAAATATGGTGGTctg                                                                        >  1:2722724/1‑148 (MQ=255)
                                                                             tCGTTTACGGGCAAATCCATAACACCGTTCGTTTTCCCAACCTGAATCTCGACAACTCCGTCCACATAACTAACCTGGTCTTTTCCATCTTGCGTAACGCTCGCGCGCTGCATGTGGGTGAAGCGCCAAATATGGTGGTCTGCTgggg                                                                  <  2:2984677/148‑1 (MQ=255)
                                                                                tttACGGGCAAATCCATAACACCGTTCGTTTTCCCAACCTGAATCTCGACAACTCCGTCCACATAACTAACCTGGTCTTTTCCATCTTGCGTAACGCTCGCGCGCTGCATGTGGGTGAAGCGCCAAATATGGTGGTCTGCTGGGGCgg                                                               <  2:1657125/148‑1 (MQ=255)
                                                                                 ttACGGGCAAATCCATAACACCGTTCGTTTTCCCAACCTGAATCTCGACAACTCCGTCCACATAACTAACCTGGTCTTTTCCATCTTGCGTAACGCTCGCGCGCTGCATGTGGGTGAAGCGCCAAATATGGTGGTCTGCTGGGGCGGt                                                              <  2:1853419/148‑1 (MQ=255)
                                                                                      ggaaaaTCCATAACACCGTTCGTTTTCCCAACCTGAATCTCGACAACTCCGTCCACATAACTAACCTGGTCTTTTCCATCTTGCGTAACGCTCGCGCGCTGCATGTGGGTGAAGCGCCAAATATGGTGGTCTGCTGGGGCGGTCACTc                                                         <  1:1093115/145‑1 (MQ=255)
                                                                                        cAAATCCATAACACCGTTCGTTTTCCCAACCTGAATCTCGACAACTCCGTCCACATAACTAACCTGGTCTTTTCCATCTTGCGTAACGCTCGCGCGCTGCATGTGGGTGAAGCGCCAAATATGGTGGTCTGCTGGGGCGGTCACTCaa                                                       >  1:1195506/1‑148 (MQ=255)
                                                                                           aTCCATAACACCGTTCGTTTTCCCAACCTGAATCTCGACAACTCCGTCCACATAACTAACCTGGTCTTTTCCATCTTGCGTAACGCTCGCGCGCTGCATGTGGGTGAAGCGCCAAATATGGTGGTCTGCTGGGGCGGTCACTCAATTa                                                    <  2:292702/148‑1 (MQ=255)
                                                                                            tCCATAACACCGTTCGTTTTCCCAACCTGAATCTCGACAACTCCGTCCACATAACTAACCTGGTCTTTTCCATCTTGCGTAACGCTCGCGCGCTGCATGTGGGTGAAGCGCCAAATATGGTGGTCTGCTGGGGCGGTCACTCAATTaa                                                   >  1:3105008/1‑148 (MQ=255)
                                                                                            tCCATAACACCGTTCGTTTTCCCAACCTGAATCTCGACAACTCCGTCCACATAACTAACCTGGTCTTTTCCATCTTGCGTAACGCTCGCGCGCTGCATGTGGGTGAAGCGCCAAATATGGTGGTCTGCTGGGGCGGTCACTCAATTaa                                                   >  2:284609/1‑148 (MQ=255)
                                                                                            tCCATAACACCGTTCGTTTTCCCAACCTGAATCTCGACAACTCCGTCCACATAACTAACCTGGTCTTTTCCATCTTGCGTAACGCTCGCGCGCTGCATGTGGGTGAAGCGCCAAATATGGTGGTCTGCTGGGGCGGTCACTCAATTaa                                                   >  1:894878/1‑148 (MQ=255)
                                                                                              cATAACACCGTTCGTTTTCCCAACCTGAATCTCGACAACTCCGTCCACATAACTAACCTGGTCTTTTCCATCTTGCGTAACGCTCGCGCGCTGCATGTGGGTGAAGCGCCAAATATGGTGGTCTGCTGGGGCGGTCACTCAATTAACg                                                 >  1:1834235/1‑148 (MQ=255)
                                                                                                tAACACCGTTCGTTTTCCCAACCTGAATCTCGACAACTCCGTCCACATCACTAACCTGGTCTTTTCCATCTTGCGTAACGCTCGCGCGCTGCATGTGGGTGAAGCGCCAAATATGGTGGTCTGCTGGGGCGGTCACTCAATTAACGaa                                               >  2:2888331/1‑148 (MQ=255)
                                                                                                 aaCACCGTTCGTTTTCCCAACCTGAATCTCGACAACTCCGTCCACATAACTAACCTGGTCTTTTCCATCTTGCGTAACGCTCGCGCGCTGCATGTGGGTGAAGCGCCAAATATGGTGGTCTGCTGGGGCGGTCACTCAATTAACGaaa                                              <  1:441010/148‑1 (MQ=255)
                                                                                                  acacCGTTCGTTTTCCCAACCTGAATCTCGACAACTCCGTCCACATAACTAACCTGGTCTTTTCCATCTTGCGTAACGCTCGCGCGCTGCATGTGGGTGAAGCGCCAAATATGGTGGTCTGCTGGGGCGGTCACTCAATTAACGaaaa                                             <  2:3218709/148‑1 (MQ=255)
                                                                                                  acacCGTTCGTTTTCCCAACCTGAATCTCGACAACTCCGTCCACATAACTAACCTGGTCTTTTCCATCTTGCGTAACGCTCGCGCGCTGCATGTGGGTGAAGCGCCAAATATGGTGGTCTGCTGGGGCGGTCACTCAATTAACGaaaa                                             <  2:3281058/148‑1 (MQ=255)
                                                                                                  acacCGTTCGTTTTCCCAACCTGAATCTCGACAACTCCGTCCACATAACTAACCTGGTCTTTTCCATCTTGCGTAACGCTCGCGCGCTGCATGTGGGTGAAGCGCCAAATATGGTGGTCTGCTGGGGCGGTCACTCAATTAACGaaaa                                             >  2:1332154/1‑148 (MQ=255)
                                                                                                  acacCGTTCGTTTTCCCAACCTGAATCTCGACAACTCCGTCCACATAACTAACCTGGTCTTTTCCATCTTGCGTAACGCTCGCGCGCTGCATGTGGGTGAAGCGCCAAATATGGTGGTCTGCTGGGGCGGTCACTCAATTAACGaaaa                                             <  2:2104158/148‑1 (MQ=255)
                                                                                                   cacCGTTCGTTTTCCCAACCTGAATCTCGACAACTCCGTCCACATAACTAACCTGGTCTTTTCCATCTTGCGTAACGCTCGCGCGCTGCATGTGGGTGAAGCGCCAAATATGGTGGTCTGCTGGGGCGGTCACTCAATTAACGAAAAc                                            >  1:1722956/1‑148 (MQ=255)
                                                                                                     ccGTTCGTTTTCCCAACCTGAATCTCGACAACTCCGTCCACATAACTAACCTGGTCTTTTCCATCTTGCGTAACGCTCGCGCGCTGCATGTGGGTGAAGCGCCAAATATGGTGGTCTGCTGGGGCGGTCACTCAATTAACGAAAACGa                                          >  1:1286456/1‑148 (MQ=255)
                                                                                                     ccGTTCGTTTTCCCAACCTGAATCTCGACAACTCCGTCCACATAACTAACCTGGTCTTTTCCATCTTGCGTAACGCTCGCGCGCTGCATGTGGGTGAAGCGCCAAATATGGTGGTCTGCTGGGGCGGTCACTCAATTAACGAAAACGa                                          >  2:1067904/1‑148 (MQ=255)
                                                                                                         tcgttTTCCCAACCTGAATCTCGACAACTCCGTCCACATCACTAACCTGGTCTTTTCCATCTTGCGTAACGCTCGCGCGCTGCATGTGGGTGAAGCGCCAAATATGGTGGTCTGCTGGGGCGGTCACTCAATTAACGAAAACGAGTAt                                      >  1:2098002/1‑148 (MQ=255)
                                                                                                         tcgttTTCCCAACCTGAATCTCGACAACTCCGTCCACATAACTAACCTGGTCTTTTCCATCTTGCGTAACGCTCGCGCGCTGCATGTGGGTGAAGCGCCAAATATGGTGGTCTGCTGGGGCGGTCACTCAATTAACGAAAACGAGTAt                                      <  1:2654044/148‑1 (MQ=255)
                                                                                                         tcgttTTCCCAACCTGAATCTCGACAACTCCGTCCACATAACTAACCTGGTCTTTTCCATCTTGCGTAACGCTCGCGCGCTGCATGTGGGTGAAGCGCCAAATATGGTGGTCTGCTGGGGCGGTCACTCAATTAACGAAAACGAGTAt                                      <  2:1192108/148‑1 (MQ=255)
                                                                                                          cgatttCCCAACCTGAATCTCGACAACTCCGTCCACATAACTAACCTGGTCTTTTCCATCTTGCGTAACGCTCGCGCGCTGCATGTGGGTGAAGCGCCAAATATGGTGGTCTGCTGGGGCGGTCACTCAATTAACGAAAACGAGTAtt                                     <  2:1370651/145‑1 (MQ=255)
                                                                                                           gttTTCCCAACCTGAATCTCGACAACTCCGTCCACATAACTAACCTGGTCTTTTCCATCTTGCGTAACGCTCGCGCGCTGCATGTGGGTGAAGCGCCAAATATGGTGGTCTGCTGGGGCGGTCACTCAATTAACGAAAACGAGTAttt                                    <  1:766871/148‑1 (MQ=255)
                                                                                                           gttTTCCCAACCTGAATCTCGACAACTCCGTCCACATAACTAACCTGGTCTTTTCCATCTTGCGTAACGCTCGCGCGCTGCATGTGGGTGAAGCGCCAAATATGGTGGTCTGCTGGGGCGGTCACTCAATTAACGAAAACGAGTAttt                                    <  2:1062471/148‑1 (MQ=255)
                                                                                                            ttttCCCAACCTGAATCTCGACAACTCCGTCCACATAACTAACCTGGTCTTTTCCATCTTGCGTAACGCTCGCGCGCTGCATGTGGGTGAAGCGCCAAATATGGTGGTCTGCTGGGGCGGTCACTCAATTAACGAAAACGAGTATTTg                                   <  1:1368292/148‑1 (MQ=255)
                                                                                                            ttttCCCAACCTGAATCTCGACAACTCCGTCCACATAACTAACCTGGTCTTTTCCATCTTGCGTAACGCTCGCGCGCTGCATGTGGGTGAAGCGCCAAATATGGTGGTCTGCTGGGGCGGTCACTCAATTAACGAAAACGAGTATTTg                                   <  1:252211/148‑1 (MQ=255)
                                                                                                             tttCCCAACCTGAATCTCGACAACTCCGTCCACATAACTAACCTGGTCTTTTCCATCTTGCGTAACGCTCGCGCGCTGCATGTGGGTGAAGCGCCAAATATGGTGGTCTGCTGGGGCGGTCACTCAATTAACGAAAACGAGTATTTGt                                  <  2:713906/148‑1 (MQ=255)
                                                                                                             tttCCCAACCTGAATCTCGACAACTCCGTCCACATAACTAACCTGGTCTTTTCCATCTTGCGTAACGCTCGCGCGCTGCATGTGGGTGAAGCGCCAAATATGGTGGTCTGCTGGGGCGGTCACTCAATTAACGAAAACGAGTATTTGt                                  >  2:1690042/1‑148 (MQ=255)
                                                                                                             tttCCCAACCTGAATCTCGACAACTCCGTCCACATAACTAACCTGGTCTTTTCCATCTTGCGTAACGCTCGCGCGCTGCATGTGGGTGAAGCGCCAAATATGGTGGTCTGCTGGGGCGGTCACTCAATTAACGAAAACGAGTATTTGt                                  >  2:1690325/1‑148 (MQ=255)
                                                                                                                cccAACCTGAATCTCGACAACTCCGTCCACATAACTAACCTGGTCTTTTCCATCTTGCGTAACGCTCGCGCGCTGCATGTGGGTGAAGCGCCAAATATGGTGGTCTGCTGGGGCGGTCACTCAATTAACGAAAACGAGTATTTGTATg                               >  1:2644497/1‑148 (MQ=255)
                                                                                                                cccAACCTGAATCTCGACAACTCCGTCCACATAACTAACCTGGTCTTTTCCATCTTGCGTAACGCTCGCGCGCTGCATGTGGGTGAAGCGCCAAATATGGTGGTCTGCTGGGGCGGTCACTCAATTAACGAAAACGAGTATTTGTATg                               <  1:2994841/148‑1 (MQ=255)
                                                                                                                   aaCCTGAATCTCGACAACTCCGTCCACATAACTAACCTGGTCTTTTCCATCTTGCGTAACGCTCGCGCGCTGCATGTGGGTGAAGCGCCAAATATGGTGGTCTGCTGGGGCGGTCACTCAATTAACGAAAACGAGTATTTGTATGccc                            >  1:2908038/1‑148 (MQ=255)
                                                                                                                    aCCTGAATCTCGACAACTCCGTCCACATAACTAACCTGGTCTTTTCCATCTTGCGTAACGCTCGCGCGCTGCATGTGGGTGAAGCGCCAAATATGGTGGTCTGCTGGGGCGGTCACTCAATTAACGAAAACGAGTATTTGTATGCCcg                           <  1:2615389/148‑1 (MQ=255)
                                                                                                                              cGACAACTCCGTCCACATAACTAACCTGGTCTTTTCCATCTTGCGTAACGCTCGCGCGCTGCATGTGGGTGAAGCGCCAAATATGGTGGTCTGCTGGGGCGGTCACTCAATTAACGAAAACGAGTATTTGTATGCCCGTCGCGTCGGa                 >  2:1575180/1‑148 (MQ=255)
                                                                                                                              cGACAACTCCGTCCACATAACTAACCTGGTCTTTTCCATCTTGCGTAACGCTCGCGCGCTGCATGTGGGTGAAGCGCCAAATATGGTGGTCTGCTGGGGCGGTCACTCAATTAACGAAAACGAGTATTTGTATGCCCGTCGCGTCGGa                 >  2:326206/1‑148 (MQ=255)
                                                                                                                              cGACAACTCCGTCCACATAACTAACCTGGTCTTTTCCATCTTGCGTAACGCTCGCGCGCTGCATGTGGGTGAAGCGCCAAATATGGTGGTCTGCTGGGGCGGTCACTCAATTAACGAAAACGAGTATTTGTATGCCCGTCGCGTCGGa                 >  2:578775/1‑148 (MQ=255)
                                                                                                                                aCAACTCCGTCCACATAACTAACCTGGTCTTTTCCATCTTGCGTAACGCTCGCGCGCTGCATGTGGGTGAAGCGCCAAATATGGTGGTCTGCTGGGGCGGTCACTCAATTAACGAAAACGAGTATTTGTATGCCCGTCGCGTCGGaaa               >  1:2405894/1‑148 (MQ=255)
                                                                                                                                aCAACTCCGTCCACATAACTAACCTGGTCTTTTCCATCTTGCGTAACGCTCGCGCGCTGCATGTGGGTGAAGCGCCAAATATGGTGGTCTGCTGGGGCGGTCACTCAATTAACGAAAACGAGTATTTGTATGCCCGTCGCGTCGGaaa               <  2:488642/148‑1 (MQ=255)
                                                                                                                                 cAACTCCGTCCACATCACTAACCTGGTCTTTTCCATCTTGCGTAACGCTCGCGCGCTGCATGTGGGTGAAGCGCCAAATATGGTGGTCTGCTGGGGCGGTCACTCAATTAACGAAAACGAGTATTTGTATGCCCGTCGCGTCGGAAAc              >  2:888018/1‑148 (MQ=255)
                                                                                                                                 cAACTCCGTCCACATAACTAACCTGGTCTTTTCCATCTTGCGTAACGCTCGCGCGCTGCATGTGGGTGAAGCGCCAAATATGGTGGTCTGCTGGGGCGGTCACTCAATTAACGAAAACGAGTATTTGTATGCCCGTCGCGTCGGAAAc              >  2:218678/1‑148 (MQ=255)
                                                                                                                                 cAACTCCGTCCACATAACTAACCTGGTCTTTTCCATCTTGCGTAACGCTCGCGCGCTGCATGTGGGTGAAGCGCCAAATATGGTGGTCTGCTGGGGCGGTCACTCAATTAACGAAAACGAGTATTTGTATGCCCGTCGCGTCGGAAAc              >  1:1068936/1‑148 (MQ=255)
                                                                                                                                   aCTCCGTCCACATCACTAACCTGGTCTTTTCCATCTTGCGTAACGCTCGCGCGCTGCATGTGGGTGAAGCGCCAAATATGGTGGTCTGCTGGGGCGGTCACTCAATTAACGACAACGAGTATTTGTATGCCCGTCGCGTCGGAAACCa            >  2:2043344/1‑148 (MQ=255)
                                                                                                                                     tCCGTCCACATAACTAACCTGGTCTTTTCCATCTTGCGTAACGCTCGCGCGCTGCATGTGGGTGAAGCGCCAAATATGGTGGTCTGCTGGGGCGGTCACTCAATTAACGAAAACGAGTATTTGTATGCCCGTCGCGTCGGAAACCAGc          >  2:2501172/1‑148 (MQ=255)
                                                                                                                                         tCCACATCACTAACCTGGTCTTTTCCATCTTGCGTAACGCTCGCGCGCTGCATGTGGGTGAAGCGCCAAATATGGTGGTCTGCTGGGGCGGTCACTCAATTAACGAAAACGAGTATTTGTATGCCCGTCGCGTCGGAAACCAGCTggg      <  2:631031/148‑1 (MQ=255)
                                                                                                                                             caTCACTAACCTGGTCTTTTCCATCTTGCGTAACGCTCGCGCGCTGCATGTGGGTGAAGCGCCAAATATGGTGGTCTGCTGGGGCGGTCACTCAATTAACGAAAACGAGTATTTGTATGCCCGTCGCGTCGGAAACCAGCTGGGCCTg  >  1:1850411/1‑148 (MQ=255)
                                                                                                                                             caTAACTAACCTGGTCTTTTCCATCTTGCGTAACGCTCGCGCGCTGCATGTGGGTGAAGCGCCAAATATGGTGGTCTGCTGGGGCGGTCACTCAATTAACGAAAACGAGTATTTGTATGCCCGTCGCGTCGGAAACCAGCTGGGCCTg  >  2:429486/1‑148 (MQ=255)
                                                                                                                                                |                                                                                                                                                
GGAAGAGGCTTTTGTCGATGGGCGAATTATTCGCGCTTTGCAGGCCAACTTGTTCGCGGTCCTGCGTGACATTCTCTTCGTTTACGGGCAAATCCATAACACCGTTCGTTTTCCCAACCTGAATCTCGACAACTCCGTCCACATCACTAACCTGGTCTTTTCCATCTTGCGTAACGCTCGCGCGCTGCATGTGGGTGAAGCGCCAAATATGGTGGTCTGCTGGGGCGGTCACTCAATTAACGAAAACGAGTATTTGTATGCCCGTCGCGTCGGAAACCAGCTGGGCCTG  >  CP009273/2919906‑2920194

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