Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 4,127,022 T→G 100% E188D (GAA→GAC)  yijF ← DUF1287 family protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092734,127,0220TG99.0% 875.3 / ‑3.9 216E188D (GAA→GAC) yijFDUF1287 family protein
Reads supporting (aligned to +/- strand):  ref base T (3/0);  new base G (131/82);  total (134/82)
Fisher's exact test for biased strand distribution p-value = 2.90e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
Rejected as polymorphism: E-value score below prediction cutoff.
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 225 total aligned reads displayed.

TTGTGCTGGCACTGGCGCTGGTGAGCCGTAAGAAAAAAGAAGCCGTCAGTGTAAAAAGGATCTGAATTTTTTCTTCATGTGGGGCGATCTCTTATTTAACAAAATAACGATAATGCCCCACCATCCGCCAGTTAAACAGCACATCTTCTTCCTGCGCGCCTGCGCCAATGTTATGTATCACCAGCGGCGTACCGTCGCGGGCGAAGCCATCTGAAACCACCCCAATATGTGCCAGCCCGTTATCCAGTCGCCAGGAGACAATATCGCCCGCTTGATAGTCACTGGGGT  >  CP009273/4126877‑4127164
                                                                                                                                                 |                                                                                                                                              
ttgtgctgGCACTGGCGCTGGTGAGCCGTAAGAAAAAAGAAGCCGTCAGTGTAAAAAGGATCTGAATTTTTTCTTCATGTGGGGCGATCTCTTATTTAACAAAATAACGATAATGCCCCACCATCCGCCAGTTAAACAGCACAtcgtc                                                                                                                                              >  1:2032838/1‑148 (MQ=255)
ttgtgctgGCACTGGCGCTGGTGAGCCGTAAGAAAAAAGAAGCCGTCAGTGTAAAAAGGATCTGAATTTTTTCTTCATGTGGGGCGATCTCTTATTTAACAAAATAACGATAATGCCCCACCATCCGCCAGTTAAACAGCACAtcgtc                                                                                                                                              >  2:1994198/1‑148 (MQ=255)
 tgtgctgGCACTGGCGCTGGTGAGCCGTAAGAAAAAAGAAGCCGTCAGTGTAAAAAGGATCTGAATTTTTTCTTCATGTGGGGCGATCTCTTATTTAACAAAATAACGATAATGCCCCACCATCCGCCAGTTAAACAGCACAtcgtct                                                                                                                                             >  1:2957231/1‑148 (MQ=255)
 tgtgctgGCACTGGCGCTGGTGAGCCGTAAGAAAAAAGAAGCCGTCAGTGTAAAAAGGATCTGAATTTTTTCTTCATGTGGGGCGATCTCTTATTTAACAAAATAACGATAATGCCCCACCATCCGCCAGTTAAACAGCACAtcgtct                                                                                                                                             >  2:414793/1‑148 (MQ=255)
   tgctgGCACTGGCGCTGGTGAGCCGTAAGAAAAAAGAAGCCGTCAGTGTAAAAAGGATCTGAATTTTTTCTTCATGTGGGGCGATCTCTTATTTAACAAAATAACGATAATGCCCCACCATCCGCCAGTTAAACAGCACATCgtcttc                                                                                                                                           <  1:984372/148‑1 (MQ=255)
       ggCACTGGCGCTGGTGAGCCGTAAGAAAAAAGAAGCCGTCAGTGTAAAAAGGATCTGAATTTTTTCTTCATGTGGGGCGATCTCTTATTTAACAAAATAACGATAATGCCCCACCATCCGCCAGTTAAACAGCACATCGTCTTCCTgc                                                                                                                                       >  2:2778409/1‑148 (MQ=255)
        gCACTGGCGCTGGTGAGCCGTAAGAAAAAAGAAGCCGTCAGTGTAAAAAGGATCTGAATTTTTTCTTCATGTGGGGCGATCTCTTATTTAACAAAATAACGATAATGCCCCACCATCCGCCAGTTAAACAGCACATCGTCTTCCTgcg                                                                                                                                      >  2:315886/1‑148 (MQ=255)
         cACTGGCGCTGGTGAGCCGTAAGAAAAAAGAAGCCGTCAGTGTAAAAAGGATCTGAATTTTTTCTTCATGTGGGGCGATCTCTTATTTAACAAAATAACGATAATGCCCCACCATCCGCCAGTTAAACAGCACATCGTCTTCCTgcgc                                                                                                                                     >  1:2352951/1‑148 (MQ=255)
           cTGGCGCTGGTGAGCCGTAAGAAAAAAGAAGCCGTCAGTGTAAAAAGGATCTGAATTTTTTCTTCATGTGGGGCGATCTCTTATTTAACAAAATAACGATAATGCCCCACCATCCGCCAGTTAAACAGCACATCGTCTTCCTgcgcgc                                                                                                                                   <  2:1639516/148‑1 (MQ=255)
                 cTGGTGAGCCGTAAGAAAAAAGAAGCCGTCAGTGTAAAAAGGATCTGAATTTTTTCTTCATGTGGGGCGATCTCTTATTTAACAAAATAACGATAATGCCCCACCATCCGCCAGTTAAACAGCACATCGTCTTCCTGCGCGCCTgcgc                                                                                                                             <  1:1811308/148‑1 (MQ=255)
                   ggTGAGCCGTAAGAAAAAAGAAGCCGTCAGTGTAAAAAGGATCTGAATTTTTTCTTCATGTGGGGCGATCTCTTATTTAACAAAATAACGATAATGCCCCACCATCCGCCAGTTAAACAGCACATCGTCTTCCTGCGCGCCTGCGCCa                                                                                                                           >  1:1745701/1‑148 (MQ=255)
                   ggTGAGCCGTAAGAAAAAAGAAGCCGTCAGTGTAAAAAGGATCTGAATTTTTTCTTCATGTGGGGCGATCTCTTATTTAACAAAATAACGATAATGCCCCACCATCCGCCAGTTAAACAGCACATCGTCTTCCTGCGCGCCTGCGCCa                                                                                                                           <  2:490263/148‑1 (MQ=255)
                    gTGAGCCGTAAGAAAAAAGAAGCCGTCAGTGTAAAAAGGATCTGAATTTTTTCTTCATGTGGGGCGATCTCTTATTTAACAAAATAACGATAATGCCCCACCATCCGCCAGTTAAACAGCACATCGTCTTCCTGCGCGCCTGCGCCaa                                                                                                                          >  1:1207784/1‑148 (MQ=255)
                      gAGCCGTAAGAAAAAAGAAGCCGTCAGTGTAAAAAGGATCTGAATTTTTTCTTCATGTGGGGCGATCTCTTATTTAACAAAATAACGATAATGCCCCACCATCCGCCAGTTAAACAGCACATCGTCTTCCTGCGCGCCTGCGCCAATg                                                                                                                        <  1:624546/148‑1 (MQ=255)
                       aGCCGTAAGAAAAAAGAAGCCGTCAGTGTAAAAAGGATCTGAATTTTTTCTTCATGTGGGGCGATCTCTTATTTAACAAAATACCGATAATGCCCCACCATCCGCCAGTTAAACAGCACATCGTCTTCCTGCGCGCCTGCGCCACTGt                                                                                                                       >  2:1525765/1‑148 (MQ=255)
                       aGCCGTAAGAAAAAAGAAGCCGTCAGTGTAAAAAGGATCTGAATTTTTTCTTCATGTGGGGCGATCTCTTATTTAACAAAATAACGATAATGCCCCACCATCCGCCAGTTAAACAGCACATCGTCTTCCTGCGCGCCTGCGCCAATGt                                                                                                                       >  1:3292628/1‑148 (MQ=255)
                        gCCGTAAGAAAAAAGAAGCCGTCAGTGTAAAAAGGATCTGAATTTTTTCTTCATGTGGGGCGATCTCTTATTTAACAAAATAACGATAATGCCCCACCATCCGCCAGTTAAACAGCACATCGTCTTCCTGCGCGCCTGCGCCAATGtt                                                                                                                      <  1:2831368/148‑1 (MQ=255)
                        gCCGTAAGAAAAAAGAAGCCGTCAGTGTAAAAAGGATCTGAATTTTTTCTTCATGTGGGGCGATCTCTTATTTAACAAAATAACGATAATGCCCCACCATCCGCCAGTTAAACAGCACATCGTCTTCCTGCGCGCCTGCGCCAATGtt                                                                                                                      <  2:823940/148‑1 (MQ=255)
                          cGTAAGAAAAAAGAAGCCGTCAGTGTAAAAAGGATCTGAATTTTTTCTTCATGTGGGGCGATCTCTTATTTAACAAAATAACGATAATGCCCCACCATCCGCCAGTTAAACAGCACATCGTCTTCCTGCGCGCCTGCGCCAATGTTAt                                                                                                                    <  2:216279/148‑1 (MQ=255)
                          cGTAAGAAAAAAGAAGCCGTCAGTGTAAAAAGGATCTGAATTTTTTCTTCATGTGGGGCGATCTCTTATTTAACAAAATAACGATAATGCCCCACCATCCGCCAGTTAAACAGCACATCGTCTTCCTGCGCGCCTGCGCCAATGTTAt                                                                                                                    <  2:2219820/148‑1 (MQ=255)
                            tAAGAAAAAAGAAGCCGTCAGTGTAAAAAGGATCTGAATTTTTTCTTCATGTGGGGCGATCTCTTATTTAACAAAATAACGATAATGCCCCACCATCCGCCAGTTAAACAGCACATCGTCTTCCTGCGCGCCTGCGCCAATGTTATGt                                                                                                                  <  2:1008118/148‑1 (MQ=255)
                            tAAGAAAAAAGAAGCCGTCAGTGTAAAAAGGATCTGAATTTTTTCTTCATGTGGGGCGATCTCTTATTTAACAAAATAACGATAATGCCCCACCATCCGCCAGTTAAACAGCACATCGTCTTCCTGCGCGCCTGCGCCAATGTTATGt                                                                                                                  <  2:2178847/148‑1 (MQ=255)
                            tAAGAAAAAAGAAGCCGTCAGTGTAAAAAGGATCTGAATTTTTTCTTCATGTGGGGCGATCTCTTATTTAACAAAATAACGATAATGCCCCACCATCCGCCAGTTAAACAGCACATCGTCTTCCTGCGCGCCTGCGCCAATGTTATGt                                                                                                                  >  1:536258/1‑148 (MQ=255)
                            tAAGAAAAAAGAAGCCGTCAGTGTAAAAAGGATCTGAATTTTTTCTTCATGTGGGGCGATCTCTTATTTAACAAAATAACGATAATGCCCCACCATCCGCCAGTTAAACAGCACATCGTCTTCCTGCGCGCCTGCGCCAATGTTATGt                                                                                                                  <  2:3241344/148‑1 (MQ=255)
                            tAAGAAAAAAGAAGCCGTCAGTGTAAAAAGGATCTGAATTTTTTCTTCATGTGGGGCGATCTCTTATTTAACAAAATAACGATAATGCCCCACCATCCGCCAGTTAAACAGCACATCGTCTTCCTGCGCGCCTGCGCCAATGTTATGt                                                                                                                  >  1:1110785/1‑148 (MQ=255)
                              aGAAAAAAGAAGCCGTCAGTGTAAAAAGGATCTGAATTTTTTCTTCATGTGGGGCGATCTCTTATTTAACAAAATAACGATAATGCCCCACCATCCGCCAGTTAAACAGCACATCGTCTTCCTGCGCGCCTGCGCCAATGTTATGTAt                                                                                                                <  1:276952/148‑1 (MQ=255)
                                aaaaaaGAAGCCGTCAGTGTAAAAAGGATCTGAATTTTTTCTTCATGTGGGGCGATCTCTTATTTAACAAAATAACGATAATGCCCCACCATCCGCCAGTTAAACAGCACATCGTCTTCCTGCGCGCCTGCGCCAATGTTATGTATCa                                                                                                              >  1:1344173/1‑148 (MQ=255)
                                   aaaGAAGCCGTCAGTGTAAAAAGGATCTGAATTTTTTCTTCATGTGGGGCGATCTCTTATTTAACAAAATAACGATAATGCCCCACCATCCGCCAGTTAAACAGCACATCGTCTTCCTGCGCGCCTGCGCCAATGTTATGTATCACCa                                                                                                           >  2:668244/1‑148 (MQ=255)
                                      gaagCCGTCAGTGTAAAAAGGATCTGAATTTTTTCTTCATGTGGGGCGATCTCTTATTTAACAAAATAACGATAATGCCCCACCATCCGCCAGTTAAACAGCACATCGTCTTCCTGCGCGCCTGCGCCAATGTTATGTATCACCAgcg                                                                                                        <  2:1670321/148‑1 (MQ=255)
                                       aagCCGTCAGTGTAAAAAGGATCTGAATTTTTTCTTCATGTGGGGCGATCTCTTATTTAACAAAATAACGATAATGCCCCACCATCCGCCAGTTAAACAGCACATCGTCTTCCTGCGCGCCTGCGCCAATGTTATGTATCACCAgcgg                                                                                                       >  1:1883898/1‑148 (MQ=255)
                                        agCCGTCAGTGTAAAAAGGATCTGAATTTTTTCTTCATGTGGGGCGATCTCTTATTTAACAAAATAACGATAATGCCCCACCATCCGCCAGTTAAACAGCACATCGTCTTCCTGCGCGCCTGCGCCAATGTTATGTATCACCAgcggc                                                                                                      >  2:1736997/1‑148 (MQ=255)
                                          ccGTCAGTGTAAAAAGGATCTGAATTTTTTCTTCATGTGGGGCGATCTCTTATTTAACAAAATAACGATAATGCCCCACCATCCGCCAGTTAAACAGCACATCGTCTTCCTGCGCGCCTGCGCCAATGTTATGTATCACCAGCGGCGt                                                                                                    >  2:2350843/1‑148 (MQ=255)
                                              cAGTGTAAAAAGGATCTGAATTTTTTCTTCATGTGGGGCGATCTCTTATTTAACAAAATAACGATAATGCCCCACCATCCGCCAGTTAAACAGCACATCGTCTTCCTGCGCGCCTGCGCCAATGTTATGTATCACCAGCGGCGTACCg                                                                                                >  2:3261779/1‑148 (MQ=255)
                                                 tgtAAAAAGGATCTGAATTTTTTCTTCATGTGGGGCGATCTCTTATTTAACAAAATAACGATAATGCCCCACCATCCGCCAGTTAAACAGCACATCGTCTTCCTGCGCGCCTGCGCCAATGTTATGTATCACCAGCGGCGTACCGTcg                                                                                             <  1:1507864/148‑1 (MQ=255)
                                                       aaGGATCTGAATTTTTTCTTCATGTGGGGCGATCTCTTATTTAACAAAATAACGATAATGCCCCACCATCCGCCAGTTAAACAGCACATCGTCTTCCTGCGCGCCTGCGCCAATGTTATGTATCACCAGCGGCGTACCGTCGCGGGCg                                                                                       >  2:2632448/1‑148 (MQ=255)
                                                         ggATCTGAATTTTTTCTTCATGTGGGGCGATCTCTTATTTAACAAAATAACGATAATGCCCCACCATCCGCCAGTTAAACAGCACATCGTCTTCCTGCGCGCCTGCGCCAATGTTATGTATCACCAGCGGCGTACCGTCGCGGGCGaa                                                                                     >  1:686417/1‑148 (MQ=255)
                                                          gATCTGAATTTTTTCTTCATGTGGGGCGATCTCTTATTTAACAAAATAACGATAATGCCCCACCATCCGCCAGTTAAACAGCACATCGTCTTCCTGCGCGCCTGCGCCAATGTTATGTATCACCAGCGGCGTACCGTCGCGGGCGAAg                                                                                    >  1:347380/1‑148 (MQ=255)
                                                           aTCTGAATTTTTTCTTCATGTGGGGCGATCTCTTATTTAACAAAATAACGATAATGCCCCACCATCCGCCAGTTAAACAGCACATCGTCTTCCTGCGCGCCTGCGCCAATGTTATGTATCACCAGCGGCGTACCGTCGCGGGCGAAGc                                                                                   >  1:198489/1‑148 (MQ=255)
                                                           aTCTGAATTTTTTCTTCATGTGGGGCGATCTCTTATTTAACAAAATAACGATAATGCCCCACCATCCGCCAGTTAAACAGCACATCGTCTTCCTGCGCGCCTGCGCCAATGTTATGTATCACCAGCGGCGTACCGTCGCGGGCGAAGc                                                                                   >  1:2262954/1‑148 (MQ=255)
                                                              tGAATTTTTTCTTCATGTGGGGCGATCTCTTATTTAACAAAATAACGATAATGCCCCACCATCCGCCAGTTAAACAGCACATCGTCTTCCTGCGCGCCTGCGCCAATGTTATGTATCACCAGCGGCGTACCGTCGCGGGCGAAGCCAt                                                                                <  1:3173395/148‑1 (MQ=255)
                                                              tGAATTTTTTCTTCATGTGGGGCGATCTCTTATTTAACAAAATAACGATAATGCCCCACCATCCGCCAGTTAAACAGCACATCGTCTTCCTGCGCGCCTGCGCCAATGTTATGTATCACCAGCGGCGTACCGTCGCGGGCGAAGCCAt                                                                                <  2:350333/148‑1 (MQ=255)
                                                                  ttttttCTTCATGTGGGGCGATCTCTTATTTAACAAAATAACGATAATGCCCCACCATCCGCCAGTTAAACAGCACATCGTCTTCCTGCGCGCCTGCGCCAATGTTATGTATCACCAGCGGCGTACCGTCGCGGGCGAAGCCATCTGa                                                                            >  2:25349/1‑148 (MQ=255)
                                                                   tttttCTTCATGTGGGGCGATCTCTTATTTAACAAAATAACGATAATGCCCCACCATCCGCCAGTTAAACAGCACATCGTCTTCCTGCGCGCCTGCGCCAATGTTATGTATCACCAGCGGCGTACCGTCGCGGGCGAAGCCATCTGaa                                                                           >  1:2364686/1‑148 (MQ=255)
                                                                   tttttCTTCATGTGGGGCGATCTCTTATTTAACAAAATAACGATAATGCCCCACCATCCGCCAGTTAAACAGCACATCGTCTTCCTGCGCGCCTGCGCCAATGTTATGTATCACCAGCGGCGTACCGTCGCGGGCGAAGCCATCTGaa                                                                           >  1:2364637/1‑148 (MQ=255)
                                                                     tttCTTCATGTGGGGCGATCTCTTATTTAACAAAATAACGATAATGCCCCACCATCCGCCAGTTAAACAGCACATCGTCTTCCTGCGCGCCTGCGCCAATGTTATGTATCACCAGCGGCGTACCGTCGCGGGCGAAGCCATCTGAAac                                                                         <  2:2727931/148‑1 (MQ=255)
                                                                       tcttcATGTGGGGCGATCTCTTATTTAACAAAATAACGATAATGCCCCACCATCCGCCAGTTAAACAGCACATCGTCTTCCTGCGCGCCTGCGCCAATGTTATGTATCACCAGCGGCGTACCGTCGCGGGCGAAGCCATCTGAAacca                                                                       >  1:1617389/1‑148 (MQ=255)
                                                                        cttcATGTGGGGCGATCTCTTATTTAACAAAATAACGATAATGCCCCACCATCCGCCAGTTAAACAGCACATCGTCTTCCTGCGCGCCTGCGCCAATGTTATGTATCACCAGCGGCGTACCGTCGCGGGCGAAGCCATCTGAAaccac                                                                      >  2:2485557/1‑148 (MQ=255)
                                                                         ttcATGTGGGGCGATCTCTTATTTAACAAAATAACGATAATGCCCCACCATCCGCCAGTTAAACAGCACATCGTCTTCCTGCGCGCCTGCGCCAATGTTATGTATCACCAGCGGCGTACCGTCGCGGGCGAAGCCATCTGAAaccacc                                                                     <  1:2793535/148‑1 (MQ=255)
                                                                             tgtgGGGCGATCTCTTATTTAACAAAATAACGATAATGCCCCACCATCCGCCAGTTAAACAGCACATCGTCTTCCTGCGCGCCTGCGCCAATGTTATGTATCACCAGCGGCGTACCGTCGCGGGCGAAGCCATCTGAAACCACCCCaa                                                                 >  2:847004/1‑148 (MQ=255)
                                                                              gtgGGGCGATCTCTTATTTAACAAAATAACGATAATGCCCCACCATCCGCCAGTTAAAAAGCACATCGTCTTCCTGCGCGACTGCGCCAATGTTATGTATCACCAGCGGCGTACCGTCGCGGGCGAAGCCATCTGAAACCACCCCAat                                                                >  1:704489/1‑148 (MQ=255)
                                                                                 gggCGATCTCTTATTTAACAAAATAACGATAATGCCCCACCATCCGCCAGTTAAACAGCACATCGTCTTCCTGCGCGCCTGCGCCAATGTTATGTATCACCAGCGGCGTACCGTCGCGGGCGAAGCCATCTGAAACCACCCCAATAtg                                                             >  2:3012613/1‑148 (MQ=255)
                                                                                       tctcTTATTTAACAAAATAACGATAATGCCCCACCATCCGCCAGTTAAACAGCACATCTTCTTCCTGCGCGCCTGCGCCACTGGTATGTATCACCAGCGGCGTACCGTCGCGGGCGAAGCCAGCTGAAACCACCCCACTCTGTGCCAg                                                       >  1:1925979/1‑148 (MQ=255)
                                                                                       tctcTTATTTAACAAAATAACGATAATGCCCCACCATCCGCCAGTTAAACAGCACATCGTCTTCCTGCGCGCCTGCGCCAATGTTATGTATCACCAGCGGCGTACCGTCGCGGGCGAAGCCATCTGAAACCACCCCAATATGTGCCAg                                                       >  1:1035081/1‑148 (MQ=255)
                                                                                       tctcTTATTTAACAAAATAACGATAATGCCCCACCATCCGCCAGTTAAACAGCACATCGTCTTCCTGCGCGCCTGCGCCAATGTTATGTATCACCAGCGGCGTACCGTCGCGGGCGAAGCCATCTGAAACCACCCCAATATGTGCCAg                                                       >  2:328471/1‑148 (MQ=255)
                                                                                        ctctTATTTAACAAAATAACGATAATGCCCCACCATCCGCCAGTTAAACAGCACATCGTCTTCCTGCGCGCCTGCGCCAATGTTATGTATCACCAGCGGCGTACCGTCGCGGGCGATGCCATCTGAAACCACCCCAATATGTGCCAGc                                                      <  2:2566612/148‑1 (MQ=255)
                                                                                        ctctTATTTAACAAAATAACGATAATGCCCCACCATCCGCCAGTTAAACAGCACATCGTCTTCCTGCGCGCCTGCGCCAATGTTATGTATCACCAGCGGCGTACCGTCGCGGGCGAAGCCATCTGAAACCACCCCAATATGTGCCAGc                                                      <  2:3258626/148‑1 (MQ=255)
                                                                                        ctctTATTTAACAAAATAACGATAATGCCCCACCATCCGCCAGTTAAACAGCACATCGTCTTCCTGCGCGCCTGCGCCAATGTTATGTATCACCAGCGGCGTACCGTCGCGGGCGAAGCCATCTGAAACCACCCCAATATGTGCCAGc                                                      >  2:1513966/1‑148 (MQ=255)
                                                                                         tctTATTTAACAAAATAACGATAATGCCCCACCATCCGCCAGTTAAACAGCACATCGTCTTCCTGCGCGCCTGCGCCAATGTTATGTATCACCAGCGGCGTACCGTCGCGGGCGAAGCCATCTGAAACCACCCCAATATGTGCCAGcc                                                     >  1:1147883/1‑148 (MQ=255)
                                                                                         tctTATTTAACAAAATAACGATAATGCCCCACCATCCGCCAGTTAAACAGCACATCGTCTTCCTGCGCGCCTGCGCCAATGTTATGTATCACCAGCGGCGTACCGTCGCGGGCGAAGCCATCTGAAACCACCCCAATATGTGCCAGcc                                                     >  1:2202765/1‑148 (MQ=255)
                                                                                          ctTATTTAACAAAATAACGATAATGCCCCACCATCCGCCAGTTAAACAGCACATCGTCTTCCTGCGCGCCTGCGCCAATGTTATGTATCACCAGCGGCGTACCGTCGCGGGCGAAGCCATCTGAAACCACCCCAATATGTGCCAGccc                                                    >  1:2725664/1‑148 (MQ=255)
                                                                                            tATTTAACAAAATAACGATAATGCCCCACCATCCGCCAGTTAAACAGCACATCGTCTTCCTGCGCGCCTGCGCCAATGTTATGTATCACCAGCGGCGTACCGTCGCGGGCGAAGCCATCTGAAACCACCCCAATATGTGCCAGCCCGt                                                  >  2:1939596/1‑148 (MQ=255)
                                                                                              tttAACAAAATAACGATAATGCCCCACCATCCGCCAGTTAAACAGCACATCGTCTTCCTGCGCGCCTGCGCCAATGTTATGTATCACCAGCGGCGTACCGTCGCGGGCGAAGCCATCTGAAACCACCCCAATATGTGCCAGCCCGTTa                                                >  2:1196558/1‑148 (MQ=255)
                                                                                              tttAACAAAATAACGATAATGCCCCACCATCCGCCAGTTAAACAGCACATCGTCTTCCTGCGCGCCTGCGCCAATGTTATGTATCACCAGCGGCGTACCGTCGCGGGCGAAGCCATCTGAAACCACCCCAATATGTGCCAGCCCGTTa                                                <  1:1137120/148‑1 (MQ=255)
                                                                                              tttAACAAAATAACGATAATGCCCCACCATCCGCCAGTTAAACAGCACATCGTCTTCCTGCGCGCCTGCGCCAATGTTATGTATCACCAGCGGCGTACCGTCGCGGGCGAAGCCATCTGAAACCACCCCAATATGTGCCAGCCCGTTa                                                >  2:2606334/1‑148 (MQ=255)
                                                                                              tttAACAAAATAACGATAATGCCCCACCATCCGCCAGTTAAACAGCACATCGTCTTCCTGCGCGCCTGCGCCAATGTTATGTATCACCAGCGGCGTACCGTCGCGGGCGAAGCCATCTGAAACCACCCCAATATGTGCCAGCCCGTTa                                                >  2:711008/1‑148 (MQ=255)
                                                                                                tAACAAAATAACGATAATGCCCCACCATCCGCCAGTTAAACAGCACATCGTCTTCCTGCGCGCCTGCGCCAATGTTATGTATCACCAGCGGCGTACCGTCGCGGGCGAAGCCATCTGAAACCACCCCAATATGTGCCAGCCCGTTATc                                              <  2:536258/148‑1 (MQ=255)
                                                                                                  aCAAAATAACGATAATGCCCCACCATCCGCCAGTTAAACAGCACATCGTCTTCCTGCGCGCCTGCGCCAATGTTATGTATCACCAGCGGCGTACCGTCGCGGGCGAAGCCATCTGAAACCACCCCAATATGTGCCAGCCCGTTATCCa                                            >  2:3155842/1‑148 (MQ=255)
                                                                                                    aaaaTAACGATAATGCCCCACCATCCGCCAGTTAAACAGCACATCGTCTTCCTGCGCGCCTGCGCCAATGTTATGTATCACCAGCGGCGTACCGTCGCGGGCGAAGCCATCTGAAACCACCCCAATATGTGCCAGCCCGTTATCCAGt                                          <  1:2875393/148‑1 (MQ=255)
                                                                                                     aaaTAACGATAATGCCCCACCATCCGCCAGTTAAACAGCACATCGTCTTCCTGCGCGCCTGCGCCAATGTTATGTATCACCAGCGGCGTACCGTCGCGGGCGAAGCCATCTGAAACCACCCCAATATGTGCCAGCCCGTTATCCAGTc                                         >  2:2118821/1‑148 (MQ=255)
                                                                                                         aaCGATAATGCCCCACCATCCGCCAGTTAAACAGCACATCTTCTTCCTGCGCGCCTGCGCCAATGTTATGTATCACCAGCGGCGTACCGTCGCGGGCGAAGCCATCTGAAACCACCCCAATATGTGCCAGCCCGTTATCCAGTCGCCa                                     >  1:3206425/1‑148 (MQ=255)
                                                                                                            gATAATGCCCCACCATCCGCCAGTTAAACAGCACATCGTCTTCCTGCGCGCCTGCGCCAATGTTATGTATCACCAGCGGCGTACCGTCGCGGGCGAAGCCATCTGAAACCACCCCAATATGTGCCAGCCCGTTATCCAGTCGCCAGga                                  >  2:112916/1‑148 (MQ=255)
                                                                                                             aTAATGCCCCACCATCCGCCAGTTAAACAGCACATCGTCTTCCTGCGCGCCTGCGCCAATGTTATGTATCACCAGCGGCGTACCGTCGCGGGCGAAGCCATCTGAAACCACCCCAATATGTGCCAGCCCGTTATCCAGTCGCCAGgag                                 >  1:2674020/1‑148 (MQ=255)
                                                                                                             aTAATGCCCCACCATCCGCCAGTTAAACAGCACATCGTCTTCCTGCGCGCCTGCGCCAATGTTATGTATCACCAGCGGCGTACCGTCGCGGGCGAAGCCATCTGAAACCACCCCAATATGTGCCAGCCCGTTATCCAGTCGCCAGgag                                 <  1:2973768/148‑1 (MQ=255)
                                                                                                             aTAATGCCCCACCATCCGCCAGTTAAACAGCACATCGTCTTCCTGCGCGCCTGCGCCAATGTTATGTATCACCAGCGGCGTACCGTCGCGGGCGAAGCCATCTGAAACCACCCCAATATGTGCCAGCCCGTTATCCAGTCGCCAGgag                                 >  1:1739609/1‑148 (MQ=255)
                                                                                                             aTAATGCCCCACCATCCGCCAGTTAAACAGCACATCGTCTTCCTGCGCGCCTGCGCCAATGTTATGTATCACCAGCGGCGTACCGTCGCGGGCGAAGCCATCTGAAACCACCCCAATATGTGCCAGCCCGTTATCCAGTCGCCAGgag                                 <  1:2973824/148‑1 (MQ=255)
                                                                                                             aTAATGCCCCACCATCCGCCAGTTAAACAGCACATCGTCTTCCTGCGCGCCTGCGCCAATGTTATGTATCACCAGCGGCGTACCGTCGCGGGCGAAGCCATCTGAAACCACCCCAATATGTGCCAGCCCGTTATCCAGTCGCCAGgag                                 >  1:410710/1‑148 (MQ=255)
                                                                                                              tAATGCCCCACCATCCGCCAGTTAAACAGCACATCGTCTTCCTGCGCGCCTGCGCCAATGTTATGTATCACCAGCGGCGTACCGTCGCGGGCGAAGCCATCTGAAACCACCCCAATATGTGCCAGCCCGTTATCCAGTCGCCAGgaga                                <  1:2528552/148‑1 (MQ=255)
                                                                                                              tAATGCCCCACCATCCGCCAGTTAAACAGCACATCGTCTTCCTGCGCGCCTGCGCCAATGTTATGTATCACCAGCGGCGTACCGTCGCGGGCGAAGCCATCTGAAACCACCCCAATATGTGCCAGCCCGTTATCCAGTCGCCAGgaga                                >  1:633373/1‑148 (MQ=255)
                                                                                                               aaTGCCCCACCATCCGCCAGTTAAACAGCACATCGTCTTCCTGCGCGCCTGCGCCAATGTTATGTATCACCAGCGGCGTACCGTCGCGGGCGAAGCCATCTGAAACCACCCCAATATGTGCCAGCCCGTTATCCAGTCGCCAGGAGAc                               >  2:1076212/1‑148 (MQ=255)
                                                                                                                aTGCCCCACCATCCGCCAGTTAAACAGCACATCGTCTTCCTGCGCGCCTGCGCCAATGTTATGTATCACCAGCGGCGTACCGTCGCGGGCGAAGCCATCTGAAACCACCCCAATATGTGCCAGCCCGTTATCCAGTCGCCAGGAGACa                              >  1:831818/1‑148 (MQ=255)
                                                                                                                 tGCCCCACCATCCGCCAGTTAAACAGCACATCGTCTTCCTGCGCGCCTGCGCCAATGTTATGTATCACCAGCGGCGTACCGTCGCGGGCGAAGCCATCTGAAACCACCCCAATATGTGCCAGCCCGTTATCCAGTCGCCAGGAGACaa                             <  2:633015/148‑1 (MQ=255)
                                                                                                                 tGCCCCACCATCCGCCAGTTAAACAGCACATCGTCTTCCTGCGCGCCTGCGCCAATGTTATGTATCACCAGCGGCGTACCGTCGCGGGCGAAGCCATCTGAAACCACCCCAATATGTGCCAGCCCGTTATCCAGTCGCCAGGAGACaa                             <  1:2975654/148‑1 (MQ=255)
                                                                                                                 tGCCCCACCATCCGCCAGTTAAACAGCACATCGTCTTCCTGCGCGCCTGCGCCAATGTTATGTATCACCAGCGGCGTACCGTCGCGGGCGAAGCCATCTGAAACCACCCCAATATGTGCCAGCCCGTTATCCAGTCGCCAGGAGACaa                             >  1:258669/1‑148 (MQ=255)
                                                                                                                  gCCCCACCATCCGCCAGTTAAACAGCACATCGTCTTCCTGCGCGCCTGCGCCAATGTTATGTATCACCAGCGGCGTACCGTCGCGGGCGAAGCCATCTGAAACCACCCCAATATGTGCCAGCCCGTTATCCAGTCGCCAGGAGACAat                            >  2:2664036/1‑148 (MQ=255)
                                                                                                                  gCCCCACCATCCGCCAGTTAAACAGCACATCGTCTTCCTGCGCGCCTGCGCCAATGTTATGTATCACCAGCGGCGTACCGTCGCGGGCGAAGCCATCTGAAACCACCCCAATATGTGCCAGCCCGTTATCCAGTCGCCAGGAGACAat                            <  1:1268362/148‑1 (MQ=255)
                                                                                                                    cccACCATCCGCCAGTTAAACAGCACATCGTCTTCCTGCGCGCCTGCGCCAATGTTATGTATCACCAGCGGCGTACCGTCGCGGGCGAAGCCATCTGAAACCACCCCAATATGTGCCAGCCCGTTATCCAGTCGCCAGGAGACAatat                          <  2:2484377/148‑1 (MQ=255)
                                                                                                                     ccaccaTCCGCCAGTTAAACAGCACATCGTCTTCCTGCGCGCCTGCGCCAATGTTATGTATCACCAGCGGCGTACCGTCGCGGGCGAAGCCATCTGAAACCACCCCAATATGTGCCAGCCCGTTATCCAGTCGCCAGGAGACAATATc                         >  2:914338/1‑148 (MQ=255)
                                                                                                                     ccaccaTCCGCCAGTTAAACAGCACATCGTCTTCCTGCGCGCCTGCGCCAATGTTATGTATCACCAGCGGCGTACCGTCGCGGGCGAAGCCATCTGAAACCACCCCAATATGTGCCAGCCCGTTATCCAGTCGCCAGGAGACAATATc                         >  2:2688290/1‑148 (MQ=255)
                                                                                                                      caccaTCCGCCAGTTAAACAGCACATCGTCTTCCTGCGCGCCTGCGCCAATGTTATGTATCACCAGCGGCGTACCGTCGCGGGCGAAGCCATCTGAAACCACCCCAATATGTGCCAGCCCGTTATCCAGTCGCCAGGAGACAATATCg                        >  2:1256937/1‑148 (MQ=255)
                                                                                                                       accaTCCGCCAGTTAAACAGCACATCGTCTTCCTGCGCGCCTGCGCCAATGTTATGTATCACCAGCGGCGTACCGTCGCGGGCGAAGCCATCTGAAACCACCCCAATATGTGCCAGCCCGTTATCCAGTCGCCAGGAGACAATATCGc                       >  1:886436/1‑148 (MQ=255)
                                                                                                                          aTCCGCCAGTTAAACAGCACATCGTCTTCCTGCGCGCCTGCGCCAATGTTATGTATCACCAGCGGCGTACCGTCGCGGGCGAAGCCATCTGAAACCACCCCAATATGTGCCAGCCCGTTATCCAGTCGCCAGGAGACAATATCGCCCg                    >  1:2615857/1‑148 (MQ=255)
                                                                                                                               ccAGTTAAACAGCACATCGTCTTCCTGCGCGCCTGCGCCAATGTTATGTATCACCAGCGGCGTACCGTCGCGGGCGAAGCCATCTGAAACCACCCCAATATGTGCCAGCCCGTTATCCAGTCGCCAGGAGACAATATCGCCCGCTTGa               <  2:2957231/148‑1 (MQ=255)
                                                                                                                               ccAGTTAAACAGCACATCGTCTTCCTGCGCGCCTGCGCCAATGTTATGTATCACCAGCGGCGTACCGTCGCGGGCGAAGCCATCTGAAACCACCCCAATATGTGCCAGCCCGTTATCCAGTCGCCAGGAGACAATATCGCCCGCTTGa               >  2:1769806/1‑148 (MQ=255)
                                                                                                                                cAGTTAAACAGCACATCGTCTTCCTGCGCGCCTGCGCCAATGTTATGTATCACCAGCGGCGTACCGTCGCGGGCGAAGCCATCTGAAACCACCCCAATATGTGCCAGCCCGTTATCCAGTCGCCAGGAGACAATATCGCCCGCTTGAt              >  1:3248490/1‑148 (MQ=255)
                                                                                                                                cAGTTAAACAGCACATCGTCTTCCTGCGCGCCTGCGCCAATGTTATGTATCACCAGCGGCGTACCGTCGCGGGCGAAGCCATCTGAAACCACCCCAATATGTGCCAGCCCGTTATCCAGTCGCCAGGAGACAATATCGCCCGCTTGAt              <  2:2352953/148‑1 (MQ=255)
                                                                                                                                  gTTAAACAGCACATCGTCTTCCTGCGCGCCTGCGCCAATGTTATGTATCACCAGCGGCGTACCGTCGCGGGCGAAGCCATCTGAAACCACCCCAATATGTGCCAGCCCGTTATCCAGTCGCCAGGAGACAATATCGCCCGCTTGATAg            >  1:115377/1‑148 (MQ=255)
                                                                                                                                   ttAAACAGCACATCGTCTTCCTGCGCGCCTGCGCCAATGTTATGTATCACCAGCGGCGTACCGTCGCGGGCGAAGCCATCTGAAACCACCCCAATATGTGCCAGCCCGTTATCCAGTCGCCAGGAGACAATATCGCCCGCTTGATAGt           >  2:2044759/1‑148 (MQ=255)
                                                                                                                                      aaCAGCACATCGTCTTCCTGCGCGCCTGCGCCAATGTTATGTATCACCAGCGGCGTACCGTCGCGGGCGAAGCCATCTGAAACCACCCCAATATGTGCCAGCCCGTTATCCAGTCGCCAGGAGACAATATCGCCCGCTTGATAGTCAc        >  1:2137737/1‑148 (MQ=255)
                                                                                                                                        cAGCACATCGTCTTCCTGCGCGCCTGCGCCAATGTTATGTATCACCAGCGGCGTACCGTCGCGGGCGAAGCCATCTGAAACCACCCCAATATGTGCCAGCCCGTTATCCAGTCGCCAGGAGACAATATCGCCCGCTTGATAGTCACTg      <  2:243603/148‑1 (MQ=255)
                                                                                                                                         aGCACATCGTCTTCCTGCGCGCCTGCGCCAATGTTATGTATCACCAGCGGCGTACCGTCGCGGGCGAAGCCATCTGAAACCACCCCAATATGTGCCAGCCCGTTATCCAGTCGCCAGGAGACAATATCGCCCGCTTGATAGTCACTgg     >  1:404764/1‑148 (MQ=255)
                                                                                                                                                 |                                                                                                                                              
TTGTGCTGGCACTGGCGCTGGTGAGCCGTAAGAAAAAAGAAGCCGTCAGTGTAAAAAGGATCTGAATTTTTTCTTCATGTGGGGCGATCTCTTATTTAACAAAATAACGATAATGCCCCACCATCCGCCAGTTAAACAGCACATCTTCTTCCTGCGCGCCTGCGCCAATGTTATGTATCACCAGCGGCGTACCGTCGCGGGCGAAGCCATCTGAAACCACCCCAATATGTGCCAGCCCGTTATCCAGTCGCCAGGAGACAATATCGCCCGCTTGATAGTCACTGGGGT  >  CP009273/4126877‑4127164

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