Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 4,405,730 G→T 100% intergenic (+13/+104) aidB → / ← yjfN DNA alkylation damage repair protein; flavin‑containing DNA binding protein, weak isovaleryl CoA dehydrogenase/DUF1471 family periplasmic protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092734,405,7300GT97.2% 603.5 / 6.3 196intergenic (+13/+104)aidB/yjfNDNA alkylation damage repair protein; flavin‑containing DNA binding protein, weak isovaleryl CoA dehydrogenase/DUF1471 family periplasmic protein
Reads supporting (aligned to +/- strand):  ref base G (1/5);  new base T (78/112);  total (79/117)
Fisher's exact test for biased strand distribution p-value = 4.04e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
Rejected as polymorphism: Frequency below/above cutoff threshold.
Rejected as polymorphism: Variant not supported by required number of reads on each strand.

Only 100 of 196 total aligned reads displayed.

TTGAAATATGCTTCTCCGCCAATGGCGCAGGCGTGGTGTCAGGTGATGTTAGATACGCGCGGCGGCGTACGGTTGTCAGAGCAGATCCAGAATGATTTATTGCTGCGGGCGACGGGGGGAGTGTGTGTGTAAGCGTATACGACTGATGCGACGCTGGTTTCGATTAACTAAATGAAATATGTGAAAATTGTAGGCCGGACAAGGCGCTCGCGCCGCATCCGGCATTGTTCATTTGCTGTAGATGACACTCATGCATACAGTATCGCCTGTACGCGCCAGTTATCAGG  >  CP009273/4405586‑4405872
                                                                                                                                                |                                                                                                                                              
ttGAAATATGCTTCTCCGCCAATGGCGCAGGCGTGGTGTCAGGTGATGTTAGATACGCGCGGCGGCGTACGGTTGTCAGAGCAGATCCAGAATGATTTATTGCTGCGGGCGACGGGGGGAGTGTGTGTGTAAGCGTATACGACTTATg                                                                                                                                             >  2:393916/1‑148 (MQ=255)
 tGAAATATGCTTCTCCGCCAATGGCGCAGGCGTGGTGTCAGGTGATGTTAGATACGCGCGGCGGCGTACGGTTGTCAGAGCAGATCCAGAATGATTTATTGCTGCGGGCGACGGGGGGAGTGTGTGTGTAAGCGTATACGACTTATGc                                                                                                                                            >  2:2374163/1‑148 (MQ=255)
  gAAATATGCTTCTCCGCCAATGGCGCAGGCGTGGTGTCAGGTGATGTTAGATACGCGCGGCGGCGTACGGTTGTCAGAGCAGATCCAGAATGATTTATTGCTGCGGGCGACGGGGGGAGTGTGTGTGTAAGCGTATACGACTTATGCg                                                                                                                                           <  1:1758818/148‑1 (MQ=255)
   aaaTATGCTTCTCCGCCAATGGCGCAGGCGTGGTGTCAGGTGATGTTAGATACGCGCGGCGGCGTACGGTTGTCAGAGCAGATCCAGAATGATTTATTGCTGCGGGCGACGGGGGGAGTGTGTGTGTAAGCGTATACGACTTATGCGa                                                                                                                                          >  2:3527362/1‑148 (MQ=255)
       atGCTTCTCCGCCAATGGCGCAGGCGTGGTGTCAGGTGATGTTAGATACGCGCGGCGGCGTACGGTTGTCAGAGCAGATCCAGAATGATTTATTGCTGCGGGCGACGGGGGGAGTGTGTGTGTAAGCGTATACGACTTATGCGACGCt                                                                                                                                      >  1:3280322/1‑148 (MQ=255)
          cTTCTCCGCCAATGGCGCAGGCGTGGTGTCAGGTGATGTTAGATACGCGCGGCGGCGTACGGTTGTCAGAGCAGATCCAGAATGATTTATTGCTGCGGGCGACGGTGGGAGTGTGTGTGTAAGCGTATACGACTTATGCGACGCTGGt                                                                                                                                   <  2:3370271/148‑1 (MQ=255)
          cTTCTCCGCCAATGGCGCAGGCGTGGTGTCAGGTGATGTTAGATACGCGCGGCGGCGTACGGTTGTCAGAGCAGATCCAGAATGATTTATTGCTGCGGGCGACGGGGGGAGTGTGTGTGTAAGCGTATACGACTTATGCGACGCTGGt                                                                                                                                   <  1:2460874/148‑1 (MQ=255)
          cTTCTCCGCCAATGGCGCAGGCGTGGTGTCAGGTGATGTTAGATACGCGCGGCGGCGTACGGTTGTCAGAGCAGATCCAGAATGATTTATTGCTGCGGGCGACGGGGGGAGTGTGTGTGTAAGCGTATACGACTTATGCGACGCTGGt                                                                                                                                   >  1:1860885/1‑148 (MQ=255)
           gtctcCGCCAATGGCGCAGGCGTGGTGTCAGGTGATGTTAGATACGCGCGGCGGCGTACGGTTGTCAGAGCAGATCCAGAATGATTTATTGCTGCGGGCGACGGGGGGAGTGTGTGTGTAAGCGTATACGACTTATGCGACGCTGGtt                                                                                                                                  <  2:3252845/147‑1 (MQ=255)
            tctcCGCCAATGGCGCAGGCGTGGTGTCAGGTGATGTTAGATACGCGCGGCGGCGTACGGTTGTCAGAGCAGATCCAGAATGATTTATTGCTGCGGGCGACGGGGGGAGTGTGTGTGTACGCGTATACGACTTATGCGCCGCTGGttt                                                                                                                                 >  2:2412221/1‑148 (MQ=255)
              tcCGCCAATGGCGCAGGCGTGGTGTCAGGTGATGTTAGATACGCGCGGCGGCGTACGGTTGTCAGAGCAGATCCAGAATGATTTATTGCTGCGGGCGACGGGGGGAGTGTGTGTGTAAGCGTATACGACTTATGCGACGCTGGTTTCg                                                                                                                               >  1:1707818/1‑148 (MQ=255)
               ccGCCAATGGCGCAGGCGTGGTGTCAGGTGATGTTAGATACGCGCGGCGGCGTACGGTTGTCAGAGCAGATCCAGAATGATTTATTGCTGCGGGCGACGGGGGGAGTGTGTGTGTAAGCGTATACGACTTATGCGACGCTGGTTTCGa                                                                                                                              <  2:2765670/148‑1 (MQ=255)
                cGCCAATGGCGCAGGCGTGGTGTCAGGTGATGTTAGATACGCGCGGCGGCGTACGGTTGTCAGAGCAGATCCAGAATGATTTATTGCTGCGGGCGACGGGGGGAGTGTGTGTGTAAGCGTATACGACTTATGCGACGCTGGTTTCGAt                                                                                                                             <  2:2186622/148‑1 (MQ=255)
                   cAATGGCGCAGGCGTGGTGTCAGGTGATGTTAGATACGCGCGGCGGCGTACGGTTGTCAGAGCAGATCCAGAATGATTTATTGCTGCGGGCGACGGGGGGAGTGTGTGTGTAAGCGTATACGACTTATGCGACGCTGGTTTCGATTaa                                                                                                                          <  2:365815/148‑1 (MQ=255)
                   cAATGGCGCAGGCGTGGTGTCAGGTGATGTTAGATACGCGCGGCGGCGTACGGTTGTCAGAGCAGATCCAGAATGATTTATTGCTGCGGGCGACGGGGGGAGTGTGTGTGTAAGCGTATACGACTTATGCGACGCTGGTTTCGATTaa                                                                                                                          <  1:710424/148‑1 (MQ=255)
                   cAATGGCGCAGGCGTGGTGTCAGGTGATGTTAGATACGCGCGGCGGCGTACGGTTGTCAGAGCAGATCCAGAATGATTTATTGCTGCGGGCGACGGGGGGAGTGTGTGTGTAAGCGTATACGACTTATGCGACGCTGGTTTCGATTaa                                                                                                                          <  1:710223/148‑1 (MQ=255)
                    aaTGGCGCAGGCGTGGTGTCAGGTGATGTTAGATACGCGCGGCGGCGTACGGTTGTCAGAGCAGATCCAGAATGATTTATTGCTGCGGGCGACGGGGGGAGTGTGTGTGTAAGCGTATACGACTTATGCGACGCTGGTTTCGATTAAc                                                                                                                         >  1:3692341/1‑148 (MQ=255)
                    aaTGGCGCAGGCGTGGTGTCAGGTGATGTTAGATACGCGCGGCGGCGTACGGTTGTCAGAGCAGATCCAGAATGATTTATTGCTGCGGGCGACGGGGGGAGTGTGTGTGTAAGCGTATACGACTTATGCGACGCTGGTTTCGATTAAc                                                                                                                         <  1:3680885/148‑1 (MQ=255)
                    aaTGGCGCAGGCGTGGTGTCAGGTGATGTTAGATACGCGCGGCGGCGTACGGTTGTCAGAGCAGATCCAGAATGATTTATTGCTGCGGGCGACGGGGGGAGTGTGTGTGTAAGCGTATACGACTTATGCGACGCTGGTTTCGATTAAc                                                                                                                         <  1:228433/148‑1 (MQ=255)
                    aaTGGCGCAGGCGTGGTGTCAGGTGATGTTAGATACGCGCGGCGGCGTACGGTTGTCAGAGCAGATCCAGAATGATTTATTGCTGCGGGCGACGGGGGGAGTGTGTGTGTAAGCGTATACGACTTATGCGACGCTGGTTTCGATTAAc                                                                                                                         >  2:1630604/1‑148 (MQ=255)
                     aTGGCGCAGGCGTGGTGTCAGGTGATGTTAGATACGCGCGGCGGCGTACGGTTGTCAGAGCAGATCCAGAATGATTTATTGCTGCGGGCGACGGGGGGAGTGTGTGTGTAAGCGTATACGACTTATGCGACGCTGGTTTCGATTAACt                                                                                                                        >  1:2941569/1‑148 (MQ=255)
                      tGGCGCAGGCGTGGTGTCAGGTGATGTTAGATACGCGCGGCGGCGTACGGTTGTCAGAGCAGATCCAGAATGATTTATTGCTGCGGGCGACGGGGGGAGTGTGTGTGTAAGCGTATACGACTTATGCGACGCTGGTTTCGATTAACTa                                                                                                                       <  2:1543906/148‑1 (MQ=255)
                      tGGCGCAGGCGTGGTGTCAGGTGATGTTAGATACGCGCGGCGGCGTACGGTTGTCAGAGCAGATCCAGAATGATTTATTGCTGCGGGCGACGGGGGGAGTGTGTGTGTAAGCGTATACGACTTATGCGACGCTGGTTTCGATTAACTa                                                                                                                       >  2:2926676/1‑148 (MQ=255)
                             ggCGTGGTGTCAGGTGATGTTAGATACGCGCGGCGGCGTACGGTTGTCAGAGCAGATCCAGAATGATTTATTGCTGCGGGCGACGGGGGGAGTGTGTGTGTAAGCGTATACGACTTATGCGACGCTGGTTTCGATTAACTAAATGaaa                                                                                                                >  1:1099649/1‑148 (MQ=255)
                              gCGTGGTGTCAGGTGATGTTAGATACGCGCGGCGGCGTACGGTTGTCAGAGCAGATCCAGAATGATTTATTGCTGCGGGCGACGGGGGGAGTGTGTGTGTAAGCGTATACGACTTATGCGACGCTGGTTTCGATTAACTAAATGAAat                                                                                                               <  1:3771384/148‑1 (MQ=255)
                                 tggtgTCAGGTGATGTTAGATACGCGCGGCGGCGTACGGTTGTCAGAGCAGATCCAGAATGATTTATTGCTGCGGGCGACGGGGGGAGTGTGTGTGTAAGCGTATACGACTTATGCGACGCTGGTTTCGATTAACTAAATGAAATAtg                                                                                                            <  1:3379955/148‑1 (MQ=255)
                                      tCAGGTGATGTTAGATACGCGCGGCGGCGTACGGTTGTCAGAGCAGATCCAGAATGATTTATTGCTGCGGGCGACGGGGGGAGTGTGTGTGTAAGCGTATACGACTTATGCGACGCTGGTTTCGATTAACTAAATGAAATATGTGaaa                                                                                                       <  1:950507/148‑1 (MQ=255)
                                      tCAGGTGATGTTAGATACGCGCGGCGGCGTACGGTTGTCAGAGCAGATCCAGAATGATTTATTGCTGCGGGCGACGGGGGGAGTGTGTGTGTAAGCGTATACGACTTATGCGACGCTGGTTTCGATTAACTAAATGAAATATGTGaaa                                                                                                       >  1:3334456/1‑148 (MQ=255)
                                      tCAGGTGATGTTAGATACGCGCGGCGGCGTACGGTTGTCAGAGCAGATCCAGAATGATTTATTGCTGCGGGCGACGGGGGGAGTGTGTGTGTAAGCGTATACGACTTATGCGACGCTGGTTTCGATTAACTAAATGAAATATGTGaaa                                                                                                       >  1:3334508/1‑148 (MQ=255)
                                      tCAGGTGATGTTAGATACGCGCGGCGGCGTACGGTTGTCAGAGCAGATCCAGAATGATTTATTGCTGCGGGCGACGGGGGGAGTGTGTGTGTAAGCGTATACGACTTATGCGACGCTGGTTTCGATTAACTAAATGAAATATGTGaaa                                                                                                       <  1:3584943/148‑1 (MQ=255)
                                       cAGGTGATGTTAGATACGCGCGGCGGCGTACGGTTGTCAGAGCAGATCCAGAATGATTTATTGCTGCGGGCGACGGGGGGAGTGTGTGTGTAAGCGTATACGACTTATGCGACGCTGGTTTCGATTAACTAAATGAAATATGTGaaaa                                                                                                      <  1:3166968/148‑1 (MQ=255)
                                       cAGGTGATGTTAGATACGCGCGGCGGCGTACGGTTGTCAGAGCAGATCCAGAATGATTTATTGCTGCGGGCGACGGGGGGAGTGTGTGTGTAAGCGTATACGACTTATGCGACGCTGGTTTCGATTAACTAAATGAAATATGTGaaaa                                                                                                      <  1:613327/148‑1 (MQ=255)
                                        aGGTGATGTTAGATACGCGCGGCGGCGTACGGTTGTCAGAGCAGATCCAGAATGATTTATTGCTGCGGGCGACGGGGGGAGTGTGTGTGTAAGCGTATACGACTTATGCGACGCTGGTTTCGATTAACTAAATGAAATATGTGAAAAt                                                                                                     <  2:267492/148‑1 (MQ=255)
                                        aGGTGATGTTAGATACGCGCGGCGGCGTACGGTTGTCAGAGCAGATCCAGAATGATTTATTGCTGCGGGCGACGGGGGGAGTGTGTGTGTAAGCGTATACGACTTATGCGACGCTGGTTTCGATTAACTAAATGAAATATGTGAAAAt                                                                                                     <  1:1474645/148‑1 (MQ=255)
                                          gTGATGTTAGATACGCGCGGCGGCGTACGGTTGTCAGAGCAGATCCAGAATGATTTATTGCTGCGGGCGACGGGGGGAGTGTGTGTGTAAGCGTATACGACTTATGCGACGCTGGTTTCGATTAACTAAATGAAATATGTGAAAATTg                                                                                                   >  1:143158/1‑148 (MQ=255)
                                           tGATGTTAGATACGCGCGGCGGCGTACGGTTGTCAGAGCAGATCCAGAATGATTTATTGCTGCGGGCGACGGGGGGAGTGTGTGTGTAAGCGTATACGACTTATGCGACGCTGGTTTCGATTAACTAAATGAAATATGTGAAAATTGt                                                                                                  <  1:2444010/148‑1 (MQ=255)
                                           tGATGTTAGATACGCGCGGCGGCGTACGGTTGTCAGAGCAGATCCAGAATGATTTATTGCTGCGGGCGACGGGGGGAGTGTGTGTGTAAGCGTATACGACTTATGCGACGCTGGTTTCGATTAACTAAATGAAATATGTGAAAATTGt                                                                                                  <  2:1223095/148‑1 (MQ=255)
                                           tGATGTTAGATACGCGCGGCGGCGTACGGTTGTCAGAGCAGATCCAGAATGATTTATTGCTGCGGGCGACGGGGGGAGTGTGTGTGTAAGCGTATACGACTTATGCGACGCTGGTTTCGATTAACTAAATGAAATATGTGAAAATTGt                                                                                                  <  1:2806517/148‑1 (MQ=255)
                                           tGATGTTAGATACGCGCGGCGGCGTACGGTTGTCAGAGCAGATCCAGAATGATTTATTGCTGCGGGCGACGGGGGGAGGGTGTGTGTAAGCGTATACGACTTATGCGACGCTGGTTTCGATTAACTAAATGAAATATGTGAAAATTGt                                                                                                  >  2:361320/1‑148 (MQ=255)
                                              tGTTAGATACGCGCGGCGGCGTACGGTTGTCAGAGCAGATCCAGAATGATTTATTGCTGCGGGCGACGGGGGGAGTGTGTGTGTAAGCGTATACGACTTATGCGACGCTGGTTTCGATTAACTAAATGAAATATGTGAAAATTGTAgg                                                                                               <  2:1605955/148‑1 (MQ=255)
                                               gTTAGATACGCGCGGCGGCGTACGGTTGTCAGAGCAGATCCAGAATGATTTATTGCTGCGGGCGACGGGGGGAGTGTGTGTGTAAGCGTATACGACTTATGCGACGCTGGTTTCGATTAACTAAATGAAATATGTGAAAATTGTAGGc                                                                                              >  1:2397841/1‑148 (MQ=255)
                                               gTTAGATACGCGCGGCGGCGTACGGTTGTCAGAGCAGATCCAGAATGATTTATTGCTGCGGGCGACGGGGGGAGTGTGTGTGTAAGCGTATACGACTTATGCGACGCTGGTTTCGATTAACTAAATGAAATATGTGAAAATTGTAGGc                                                                                              >  2:2863895/1‑148 (MQ=255)
                                               gTTAGATACGCGCGGCGGCGTACGGTTGTCAGAGCAGATCCAGAATGATTTATTGCTGCGGGCGACGGGGGGAGTGTGTGTGTAAGCGTATACGACTGATGCGACGCTGGTTTCGATTAACTAAATGAAATATGTGAAAATTGTAGGc                                                                                              >  2:2863924/1‑148 (MQ=255)
                                               gTTAGATACGCGCGGCGGCGTACGGTTGTCAGAGCAGATCCAGAAGGATTTATTGCTGCGGGCGACGGGGGGAGTGTGTGTGTAAGCGTATACTACTTATGCGACGCTGGTTTCGATTAACTAAATGAAGTATGTGAAAAAGGTAGGc                                                                                              >  2:3442028/1‑148 (MQ=255)
                                                ttAGATACGCGCGGCGGCGTACGGTTGTCAGAGCAGATCCAGAATGATTTATTGCTGCGGGCGACGGGGGGAGTGTGTGTGTAAGCGTATACGACTTATGCGACGCTGGTTTCGATTAACTAAATGAAATATGTGAAAATTGTAGGcc                                                                                             >  1:3793052/1‑148 (MQ=255)
                                                 tAGATACGCGCGGCGGCGTACGGTTGTCAGAGCAGATCCAGAATGATTTATTGCTGCGGGCGACGGGGGGAGTGTGTGTGTAAGCGTATACGACTTATGCGACGCTGGTTTCGATTAACTAAATGAAATATGTGAAAATTGTAGGCCg                                                                                            <  1:51908/148‑1 (MQ=255)
                                                    aTACGCGCGGCGGCGTACGGTTGTCAGAGCAGATCCAGAATGATTTATTGCTGCGGGCGACGGGGGGAGTGTGTGTGTAAGCGTATACGACTTATGCGACGCTGGTTTCGATTAACTAAATGAAATATGTGAAAATTGTAGGCCGGAc                                                                                         >  2:872861/1‑148 (MQ=255)
                                                      aCGCGCGGCGGCGTACGGTTGTCAGAGCAGATCCAGAATGATTTATTGCTGCGGGCGACGGGGGGAGTGTGTGTGTAAGCGTATACGACTTATGCGACGCTGGTTTCGATTAACTAAATGAAATATGTGAAAATTGTAGGCCGGACaa                                                                                       >  1:2283542/1‑148 (MQ=255)
                                                       cgcgcgGCGGCGTACGGTTGTCAGAGCAGATCCAGAATGATTTATTGCTGCGGGCGACGGGGGGAGTGTGTGTGTAAGCGTATACGACTTATGCGACGCTGGTTTCGATTAACTAAATGAAATATGTGAAAATTGTAGGCCGGACAAg                                                                                      <  1:1966723/148‑1 (MQ=255)
                                                        gcgagGCGGCGTACGGTTGTCAGAGCAGATCCAGAATGATTTATTGCTGCGGGCGACGGGGGGAGTGTGTGTGTAAGCGTATACGACTTATGCGACGCTGGTTTCGATTAACTAAATGAAATATGTGAAAATTGTAGGCCGGACAAgg                                                                                     >  2:898715/1‑148 (MQ=255)
                                                          gcggcggcgTACGGTTGTCAGAGCAGATCCAGAATGATTTATTGCTGCGGGCGACGGGGGGAGTGTGTGTGTAAGCGTATACGACTTATGCGACGCTGGTTTCGATTAACTAAATGAAATATGTGAAAATTGTAGGCCGGACAAGgcg                                                                                   <  1:962547/148‑1 (MQ=255)
                                                           cggcggcgTACGGTTGTCAGAGCAGATCCAGAATGATTTATTGCTGCGGGCGACGGGGGGAGTGTGTGTGTAAGCGTATACGACTGATGCTACGCTGGTTTCGATTAACTAAATGAAATATGTGAAAATTGTAGGCCGGACAAGgcgc                                                                                  <  2:3457374/148‑1 (MQ=255)
                                                                   tACGGTTGTCAGAGCAGATCCAGAATGATTTATTGCTGCGGGCGACGGGGGGAGTGTGTGTGTAAGCGTATACGACTTATGCGACGCTGGTTTCGATTAACTAAATGAAATATGTGAAAATTGTAGGCCGGACAAGGCGCTCGcgccg                                                                          >  2:1731419/1‑148 (MQ=255)
                                                                      ggTTGTCAGAGCAGATCCAGAATGATTTATTGCTGCGGGCGACGGGGGGAGTGTGTGTGTAAGCGTATACGACTTATGCGACGCTGGTTTCGATTAACTAAATGAAATATGTGAAAATTGTAGGCCGGACAAGGCGCTCGCGCCGCAt                                                                       <  1:1949183/148‑1 (MQ=255)
                                                                            cAGAGCAGATCCAGAATGATTTATTGCTGCGGGCGACGGGGGGAGTGTGTGTGTAAGCGTATACGACTTATGCGACGCTGGTTTCGATTAACTAAATGAAATATGTGAAAATTGTAGGCCGGACAAGGCGCTCGCGCCGCATCCGGCa                                                                 <  1:217289/148‑1 (MQ=255)
                                                                             agagCAGATCCAGAATGATTTATTGCTGCGGGCGACGGGGGGAGTGTGTGTGTAAGCGTATACGACTTATGCGACGCTGGTTTCGATTAACTAAATGAAATATGTGAAAATTGTAGGCCGGACAAGGCGCTCGCGCCGCATCCGGCAt                                                                >  1:3155333/1‑148 (MQ=255)
                                                                               agCAGATCCAGAATGATTTATTGCTGCGGGCGACGGGGGGAGTGTGTGTGTAAGCGTATACGACTTATGCGACGCTGGTTTCGATTAACTAAATGAAATATGTGAAAATTGTAGGCCGGACAAGGCGCTCGCGCCGCATCCGGCATTg                                                              >  1:2113911/1‑148 (MQ=255)
                                                                               agCAGATCCAGAATGATTTATTGCTGCGGGCGACGGGGGGAGTGTGTGTGTAAGCGTATACGACTTATGCGACGCTGGTTTCGATTAACTAAATGAAATATGTGAAAATTGTAGGCCGGACAAGGCGCTCGCGCCGCATCCGGCATTg                                                              <  1:2794293/148‑1 (MQ=255)
                                                                                gCAGATCCAGAATGATTTATTGCTGCGGGCGACGGGGGGAGTGTGTGTGTAAGCGTATACGACTTATGCGACGCTGGTTTCGATTAACTAAATGAAATATGTGAAAATTGTAGGCCGGACAAGGCGCTCGCGCCGCATCCGGCATTGt                                                             >  1:3085195/1‑148 (MQ=255)
                                                                                  aGATCCAGAATGATTTATTGCTGCGGGCGACGGGGGGAGTGTGTGTGTAAGCGTATACGACTTATGCGACGCTGGTTTCGATTAACTAAATGAAATATGTGAAAATTGTAGGCCGGACAAGGCGCTCGCGCCGCATCCGGCATTGTTc                                                           <  2:145832/148‑1 (MQ=255)
                                                                                   gATCCAGAATGATTTATTGCTGCGGGCGACGGGGGGAGTGTGTGTGTAAGCGTATACGACTTATGCGACGCTGGTTTCGATTAACTAAATGAAATATGTGAAAATTGTAGGCCGGACAAGGCGCTCGCGCCGCATCCGGCATTGTTCa                                                          <  1:1063587/148‑1 (MQ=255)
                                                                                   gATCCAGAATGATTTATTGCTGCGGGCGACGGGGGGAGTGTGTGTGTAAGCGTATACGACTTATGCGACGCTGGTTTCGATTAACTAAATGAAATATGTGAAAATTGTAGGCCGGACAAGGCGCTCGCGCCGCATCCGGCATTGTTCa                                                          <  1:1497177/148‑1 (MQ=255)
                                                                                     tCCAGAATGATTTATTGCTGCGGGCGACGGGGGGAGTGTGTGTGTAAGCGTATACGACTTATGCGACGCTGGTTTCGATTAACTAAATGAAATATGTGAAAATTGTAGGCCGGACAAGGCGCTCGCGCCGCATCCGGCATTGTTCAtt                                                        <  2:1363857/148‑1 (MQ=255)
                                                                                     tCCAGAATGATTTATTGCTGCGGGCGACGGGGGGAGTGTGTGTGTAAGCGTATACGACTTATGCGACGCTGGTTTCGATTAACTAAATGAAATATGTGAAAATTGTAGGCCGGACAAGGCGCTCGCGCCGCATCCGGCATTGTTCAtt                                                        <  2:2079382/148‑1 (MQ=255)
                                                                                      ccAGAATGATTTATTGCTGCGGGCGACGGGGGGAGTGTGTGTGTAAGCGTATACGACTTATGCGACGCTGGTTTCGATTAACTAAATGAAATATGTGAAAATTGTAGGCCGGACAAGGCGCTCGCGCCGCATCCGGCATTGTTCAttt                                                       >  2:2541722/1‑148 (MQ=255)
                                                                                      ccAGAATGATTTATTGCTGCGGGCGACGGGGGGAGGGTGTGTGTAAGCGTATACGACTTATGCGACGCTGGTTTCGATTAACTAAATGAAATATGTGAAACTTGTAGGCCGGACAAGGCGCTCGCGCCGCATCCGGCAGTGTTCAgtt                                                       >  2:2541728/1‑148 (MQ=39)
                                                                                              aTTTATTGCTGCGGGCGACGGGGGGAGTGTGTGTGTAAGCGTATACGACTTATGCGACGCTGGTTTCGATTAACTAAATGAAATATGTGAAAATTGTAGGCCGGACAAGGCGCTCGCGCCGCATCCGGCATTGTTCATTTGCTGTAGa                                               <  2:2481324/148‑1 (MQ=255)
                                                                                               tttATTGCTGCGGGCGACGGGGGGAGTGTGTGTGTAAGCGTATACGACTTATGCGACGCTGGTTTCGATTAACTAAATGAAATATGTGAAAATTGTAGGCCGGACAAGGCGCTCGCGCCGCATCCGGCATTGTTCATTTGCTGTAGAt                                              <  2:3334508/148‑1 (MQ=255)
                                                                                               tttATTGCTGCGGGCGACGGGGGGAGTGTGTGTGTAAGCGTATACGACTTATGCGACGCTGGTTTCGATTAACTAAATGAAATATGTGAAAATTGTAGGCCGGACAAGGCGCTCGCGCCGCATCCGGCATTGTTCATTTGCTGTAGAt                                              <  2:3334456/148‑1 (MQ=255)
                                                                                                       tgcGGGCGACGGGGGGAGTGTGTGTGTAAGCGTATACGACTTATGCGACGCTGGTTTCGATTAACTAAATGAAATATGTGAAAATTGTAGGCCGGACAAGGCGCTCGCGCCGCATCCGGCATTGTTCATTTGCTGTAGATGACACTCa                                      <  2:1666385/148‑1 (MQ=255)
                                                                                                        gcGGGCGACGGGGGGAGTGTGTGTGTAAGCGTATACGACTTATGCGACGCTGGTTTCGATTAACTAAATGAAATATGTGAAAATTGTAGGCCGGACAAGGCGCTCGCGCCGCATCCGGCATTGTTCATTTGCTGTAGATGACACTCAt                                     >  2:850578/1‑148 (MQ=255)
                                                                                                          gggCGACGTGGGGAGTGTGTGTGTAAGCGTATACGACTTATGCGAAGCTTGTTTCGATTAACTAATTGAAATATGTGAAAATTGTAGGCCGGACAAGGCGCTCGCGCCGCATCCGGCATTGTTCATTTGCTGTAGATGACAATCATGc                                   <  2:3036884/148‑1 (MQ=39)
                                                                                                          gggCGACGGGGGGAGTGTGTGTGTAAGCGTATACGACTTATGCGACGCTGGTTTCGATTAACTAAATGAAATATGTGAAAATTGTAGGCCGGACAAGGCGCTCGCGCCGCATCCGGCATTGTTCATTTGCTGTAGATGACACTCATGc                                   <  1:2072290/148‑1 (MQ=255)
                                                                                                           ggCGACGGGGGGAGTGTGTGTGTAAGCGTATACGACTTATGCGACGCTGGTTTCGATTAACTAAATGAAATATGTGAAAATTGTAGGCCGGACAAGGCGCTCGCGCCGCATCCGGCATTGTTCATTTGCTGTAGATGACACTCATGCa                                  >  1:2575193/1‑148 (MQ=255)
                                                                                                             cGACGGGGGGAGTGTGTGTGTAAGCGTATACGACTTATGCGACGCTGGTTTCGATTAACTAAATGAAATATGTGAAAATTGTAGGCCGGACAAGGCGCTCGCGCCGCATCCGGCATTGTTCATTTGCTGTAGATGACACTCATGCATa                                <  1:2374163/148‑1 (MQ=255)
                                                                                                              gACGGGGGGAGTGTGTGTGTAAGCGTATACGACTTATGCGACGCTGGTTTCGATTAACTAAATGAAATATGTGAAAATTGTAGGCCGGACAAGGCGCTCGCGCCGCATCCGGCATTGTTCATTTGCTGTAGATGACACTCATGCATAc                               <  1:3112915/148‑1 (MQ=255)
                                                                                                              gACGGGGGGAGTGTGTGTGTAAGCGTATACGACTTATGCGACGCTGGTTTCGATTAACTAAATGAAATATGTGAAAATTGTAGGCCGGACAAGGCGCTCGCGCCGCATCCGGCATTGTTCATTTGCTGTAGATGACACTCATGCATAc                               <  2:1766482/148‑1 (MQ=255)
                                                                                                                cGGGGGGAGTGTGTGTGTAAGCGTATACGACTTATGCGACGCTGGTTTCGATTAACTAAATGAAATATGTGAAAATTGTAGGCCGGACAAGGCGCTCGCGCCGCATCCGGCATTGTTCATTTGCTGTAGATGACACTCATGCATACAg                             >  1:2893248/1‑148 (MQ=255)
                                                                                                                 ggggggAGTGTGTGTGTAAGCGTATACGACTGATGCTACGCTGGTTTCGATTAACTAAATGAAATATGTGAAAATTGTAGGCCGGACAAGGCGCTCGCGCCGCATCCGGCATTGTTCATTTGCTGTAGATGACACTCATGCATACAGt                            <  2:3698178/148‑1 (MQ=255)
                                                                                                                      gAGTGTGTGTGTAAGCGTATACGACTTATGCGACGCTGGTTTCGATTAACTAAATGAAATATGTGAAAATTGTAGGCCGGACAAGGCGCTCGCGCCGCATCCGGCATTGTTCATTTGCTGTAGATGACACTCATGCATACAGTATCGc                       <  1:1048619/148‑1 (MQ=255)
                                                                                                                      gAGTGTGTGTGTAAGCGTATACGACTTATGCGACGCTGGTTTCGATTAACTAAATGAAATATGTGAAAATTGTAGGCCGGACAAGGCGCTCGCGCCGCATCCGGCATTGTTCATTTGCTGTAGATGACACTCATGCATACAGTATCGc                       <  2:2154935/148‑1 (MQ=255)
                                                                                                                      gAGTGTGTGTGTAAGCGTATACGACTTATGCGACGCTGGTTTCGATTAACTAAATGAAATATGTGAAAATTGTAGGCCGGACAAGGCGCTCGCGCCGCATCCGGCATTGTTCATTTGCTGTAGATGACACTCATGCATACAGTATCGc                       <  2:3772280/148‑1 (MQ=255)
                                                                                                                      gAGTGTGTGTGTAAGCGTATACGACTTATGCGACGCTGGTTTCGATTAACTAAATGAAATATGTGAAAATTGTAGGCCGGACAAGGCGCTCGCGCCGCATCCGGCATTGTTCATTTGCTGTAGATGACACTCATGCATACAGTATCGc                       <  2:2155423/148‑1 (MQ=255)
                                                                                                                       aGTGTGTGTGTAAGCGTATACGACTTATGCGACGCTGGTTTCGATTAACTAAATGAAATATGTGAAAATTGTAGGCCGGACAAGGCGCTCGCGCCGCATCCGGCATTGTTCATTTGCTGTAGATGACACTCATGCATACAGTATCGcc                      <  1:2857024/148‑1 (MQ=255)
                                                                                                                             tgtgtAAGCGTATACGACTTATGCGACGCTGGTTTCGATTAACTAAATGAAATATGTGAAAATTGTAGGCCGGACAAGGCGCTCGCGCCGCATCCGGCATTGTTCATTTGCTGTAGATGACACTCATGCATACAGTATCGCCTGTAcg                >  1:2711664/1‑148 (MQ=255)
                                                                                                                             tgtgtAAGCGTATACGACTTATGCGACGCTGGTTTCGATTAACTAAATGAAATATGTGAAAATTGTAGGCCGGACAAGGCGCTCGCGCCGCATCCGGCATTGTTCATTTGCTGTAGATGACACTCATGCATACAGTATCGCCTGTAcg                <  1:2076433/148‑1 (MQ=255)
                                                                                                                              gtgtAAGCGTATACGACTTATGCGACGCTGGTTTCGATTAACTAAATGAAATATGTGAAAATTGTAGGCCGGACAAGGCGCTCGCGCCGCATCCGGCATTGTTCATTTGCTGTAGATGACACTCATGCATACAGTATCGCCTGTAcgc               >  1:1504815/1‑148 (MQ=255)
                                                                                                                               tgtAAGCGTATACGACTTATGCGACGCTGGTTTCGATTAACTAAATGAAATATGTGAAAATTGTAGGCCGGACAAGGCGCTCGCGCCGCATCCGGCATTGTTCATTTGCTGTAGATGACACTCATGCATACAGTATCGCCTGTAcgcg              <  2:2264633/148‑1 (MQ=255)
                                                                                                                                  aaGCGTATACGACTTATGCGACGCTGGTTTCGATTAACTAAATGAAATATGTGAAAATTGTAGGCCGGACAAGGCGCTCGCGCCGCATCCGGCATTGTTCATTTGCTGTAGATGACACTCATGCATACAGTATCGCCTGTACGCGCCa           >  2:412022/1‑148 (MQ=255)
                                                                                                                                   aGCGTATACGACTTATGCGACGCTGGTTTCGATTAACTAAATGAAATATGTGAAAATTGTAGGCCGGACAAGGCGCTCGCGCCGCATCCGGCATTGTTCATTTGCTGTAGATGACACTCATGCATACAGTATCGCCTGTACGCGCCAg          <  2:2524195/148‑1 (MQ=255)
                                                                                                                                   aGCGTATACGACTTATGCGACGCTGGTTTCGATTAACTAAATGAAATATGTGAAAATTGTAGGCCGGACAAGGCGCTCGCGCCGCATCCGGCATTGTTCATTTGCTGTAGATGACACTCATGCATACAGTATCGCCTGTACGCGCCAg          <  1:1199473/148‑1 (MQ=255)
                                                                                                                                   aGCGTATACGACTTATGCGACGCTGGTTTCGATTAACTAAATGAAATATGTGAAAATTGTAGGCCGGACAAGGCGCTCGCGCCGCATCCGGCATTGTTCATTTGCTGTAGATGACACTCATGCATACAGTATCGCCTGTACGCGCCAg          <  1:1630604/148‑1 (MQ=255)
                                                                                                                                   aGCGTATACGACTTATGCGACGCTGGTTTCGATTAACTAAATGAAATATGTGAAAATTGTAGGCCGGACAAGGCGCTCGCGCCGCATCCGGCATTGTTCATTTGCTGTAGATGACACTCATGCATACAGTATCGCCTGTACGCGCCAg          <  2:1189058/148‑1 (MQ=255)
                                                                                                                                      gTATACGACTTATGCGACGCTGGTTTCGATTAACTAAATGAAATATGTGAAAATTGTAGGCCGGACAAGGCGCTCGCGCCGCATCCGGCATTGTTCATTTGCTGTAGATGACACTCATGCATACAGTATCGCCTGTACGCGCCAGTTa       <  2:1099649/148‑1 (MQ=255)
                                                                                                                                      gTATACGAATTATGCGACGCTGGTTTCGATTAACTAAATGAATTATGTTAAAATTGTAGGCCGGACAAGGCGCTCGCGCCGCATCCGGCATTGTTCATTTGCTGTAGATGACACTCATGCATACAGTATCGCCTGTACGCGCCAGTTa       <  1:3502388/148‑1 (MQ=255)
                                                                                                                                        ataCGACTTATGCGACGCTGGTTTCGATTAACTAAATGAAATATGTGAAAATTGTAGGCCGGACAAGGCGCTCGCGCCGCATCCGGCATTGTTCATTTGCTGTAGATGACACTCATGCATACAGTATCGCCTGTACGCGCCAGTTATc     >  2:65277/1‑148 (MQ=255)
                                                                                                                                        ataCGACTTATGCGACGCTGGTTTCGATTAACTAAATGAAATATGTGAAAATTGTAGGCCGGACAAGGCGCTCGCGCCGCATCCGGCATTGTTCATTTGCTGTAGATGACACTCATGCATACAGTATCGCCTGTACGCGCCAGTTATc     <  1:872861/148‑1 (MQ=255)
                                                                                                                                         taCGACTTATGCGACGCTGGTTTCGATTAACTAAATGAAATATGTGAAAATTGTAGGCCGGACAAGGCGCTCGCGCCGCATCCGGCATTGTTCATTTGCTGTAGATGACACTCATGCATACAGTATCGCCTGTACGCGCCAGTTATCa    <  1:1780995/148‑1 (MQ=255)
                                                                                                                                          aCGACTTATGCGACGCTGGTTTCGATTAACTAAATGAAATATGTGAAAATTGTAGGCCGGACAAGGCGCTCGCGCCGCATCCGGCATTGTTCATTTGCTGTAGATGACACTCATGCATACAGTATCGCCTGTACGCGCCAGTTATCAg   <  1:2294699/148‑1 (MQ=255)
                                                                                                                                           cGACTTATGCGACGCTGGTTTCGATTAACTAAATGAAATATGTGAAAATTGTAGGCCGGACAAGGCGCTCGCGCCGCATCCGGCATTGTTCATTTGCTGTAGATGACACTCATGCATACAGTATCGCCTGTACGCGCCAGTTATCAgg  >  1:2193445/1‑148 (MQ=255)
                                                                                                                                                |                                                                                                                                              
TTGAAATATGCTTCTCCGCCAATGGCGCAGGCGTGGTGTCAGGTGATGTTAGATACGCGCGGCGGCGTACGGTTGTCAGAGCAGATCCAGAATGATTTATTGCTGCGGGCGACGGGGGGAGTGTGTGTGTAAGCGTATACGACTGATGCGACGCTGGTTTCGATTAACTAAATGAAATATGTGAAAATTGTAGGCCGGACAAGGCGCTCGCGCCGCATCCGGCATTGTTCATTTGCTGTAGATGACACTCATGCATACAGTATCGCCTGTACGCGCCAGTTATCAGG  >  CP009273/4405586‑4405872

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: