Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 3,751,414 T→G 13.2% H603P (CAT→CCT)  selB ← selenocysteinyl‑tRNA‑specific translation factor

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092733,751,4140TG13.2% 941.6 / 139.8 334H603P (CAT→CCT) selBselenocysteinyl‑tRNA‑specific translation factor
Reads supporting (aligned to +/- strand):  ref base T (163/127);  new base G (25/19);  total (188/146)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.29e-01

Only 100 of 334 total aligned reads displayed.

TATCAAAATTAATCTAATAAAAAGCACATTGTATTCATTAAACAATGTGCTTTTTTTAAATTCATAAGAATTTTGTATTAATTTATTTTAAAACTTAATTAATCATTTCCTTATTTTTCCGGAAATAATAATGCGTCGCGCAATAAATGATCATTTCCACGACGACGCGTAAAGCCAATGCGGTCAAAATATTCCAGAATCTGAATTGCCAGCTTTCGGCCTACGCCTAAGCGATCGCGGAAATCCGCCGCGCAGGTTGAACCACACTCCTGATCGAGATCGCGGATCATATTGG  >  CP009273/3751267‑3751561
                                                                                                                                                   |                                                                                                                                                   
tATCAAAATTAAACTAAAAAAAAGCACATTGTATTCATTAAACAATGTGCTTTTTTTAAATTCATAAGAATTTTGTATTAATTTATTTTAAAACTTAATTAATCATTTCCTTATTTTTCCGGAAATAATAATGCGTCGCGCAATAAATg                                                                                                                                                    <  1:3351286/149‑1 (MQ=255)
    aaaaTTAATCTAATAAAAAGCACATTGTATTCATTAAACAATGTGCTTTTTTTAAATTCATAAGAATTTTGTATTAATTTATTTTAAAACTTAATTAATCATTTCCTTATTTTTCCGGAAATAATAATGCGTCGCGCAATAAATGATCa                                                                                                                                                <  2:1072203/149‑1 (MQ=255)
     aaaaTAAACTAATAAAAAGCACATTGTATTCATTAAACAATGTGCTTTTTTTAAATTCATAAGAATTTTGTATTAATTTATTTTAAAACTTAATTAATCATTTCCTTATTTTTCCGGAAATAATAATGCGTCGCGCAATAAATGATCAt                                                                                                                                               <  1:4950155/149‑1 (MQ=255)
       aTTAATCTAATAAAAAGCACATTGTATTCATTAAACAATGTGCTTTTTTTAAATTCATAAGAATTTTGTATTAATTTATTTTAAAACTTAATTAATCATTTCCTTATTTTTCCGGAAATAATAATGCGTCGCGCAATAAATGATCAttt                                                                                                                                             <  1:3272523/149‑1 (MQ=255)
          aaTCTAATAAAAAGCACATTGTATTCATTAAACAATGTGCTTTTTTTAAATTCATAAGAATTTTGTATTAATTTATTTTAAAACTTAATTAATCATTTCCTTATTTTTCCGGAAATAATAATGCGTCGCGCAATAAATGATCATTTCCa                                                                                                                                          >  1:882973/1‑149 (MQ=255)
                   aaaaGCACATTGTATTCATTAAACAATGTGCTTTTTTTAAATTCATAAGAATTTTGTATTAATTTATTTTAAAACTTAATTAATCATTTCCTTATTTTTCCGGAAATAATAATGCGTCGCGCAATAAATGATCATTTCCACGACGAcgc                                                                                                                                 <  2:1414763/149‑1 (MQ=255)
                   aaaaGCACATTGTATTCATTAAACAATGTGCTTTTTTTAAATTCATAAGAATTTTGTATTAATTTATTTTAAAACTTAATTAATCATTTCCTTATTTTTCCGGAAATAATAATGCGTCGCGCAATAAATGATCATTTCCACGACGAcgc                                                                                                                                 <  2:2954917/149‑1 (MQ=255)
                         acaTTGTATTCATTAAACAATGTGCTTTTTTTAAATTCATAAGAATTTTGTATTAATTTATTTTAAAACTTAATTAATCATTTCCTTATTTTTCCGGAAATAATAATGCGTCGCGCAATAAATGATCATTTCCACGACGACGCGTAAAg                                                                                                                           >  1:1303016/1‑149 (MQ=255)
                           aTTGTATTCATTAAACAATGTGCTTTTTTTAAATTCATAAGAATTTTGTATTAATTTATTTTAAAACTTAATTAATCATTTCCTTATTTTTCCGGAAATAATAATGCGTCGCGCAATAAATGATCATTTCCACGACGACGCGTAAAGcc                                                                                                                         >  2:2608300/1‑149 (MQ=255)
                            ttGTATTCATTAAACAATGTGCTTTTTTTAAATTCATAAGAATTTTGTATTAATTTATTTTAAAACTTAATTAATCATTTCCTTATTTTTCCGGAAATAATAATGCGTCGCGCAATAAATGATCATTTCCACGACGACGCGTAAAGCCa                                                                                                                        >  2:2576332/1‑149 (MQ=255)
                             tGTATTCATTAAACAATGTGCTTTTTTTAAATTCATAAGAATTTTGTATTAATTTATTTTAAAACTTAATTAATCATTTCCTTATTTTTCCGGAAATAATAATGCGTCGCGCAATAAATGATCATTTCCACGACGACGCGTAAAGCCaa                                                                                                                       >  2:1956815/1‑149 (MQ=255)
                             tGTATTCATTAAACAATGTGCTTTTTTTAAATTCATAAGAATTTTGTATTAATTTATTTTAAAACTTAATTAATCATTTCCTTATTTTTCCGGAAATAATAATGCGTCGCGCAATAAATGATCATTTCCACGACGACGCGTAAAGCCaa                                                                                                                       >  2:2849663/1‑149 (MQ=255)
                                    aTTAAACAATGTGCTTTTTTTAAATTCATAAGAATTTTGTATTAATTTATTTTAAAACTTAATTAATCATTTCCTTATTTTTCCGGAAATAATAATGCGTCGCGCAATAAAGGATCATTTCCACGACGACGCGTAAAGCCAATGCGGTc                                                                                                                <  1:420195/149‑1 (MQ=255)
                                     ttAAACAATGTGCTTTTTTTAAATTCATAAGAATTTTGTATTAATTTATTTTAAAACTTAATTAATCATTTCCTTATTTTTCCGGAAATAATAATGCGTCGCGCAATAAATGATCATTTCCACGACGACGCGTAAAGCCAATGCGGTCa                                                                                                               >  2:59299/1‑149 (MQ=255)
                                     ttAAACAATGTGCTTTTTTTAAATTCATAAGAATTTTGTATTAATTTATTTTAAAACTTAATTAATCATTTCCTTATTTTTCCGGAAATAATAATGCGTCGCGCAATAAAGGATCATTTCCACGACGACGCGTAAAGCCAATGCGGTCa                                                                                                               >  2:3153164/1‑149 (MQ=255)
                                      tAAACAATGTGCTTTTTTTAAATTCATAAGAATTTTGTATTAATTTATTTTAAAACTTAATTAATCATTTCCTTATTTTTCCGGAAATAATAATGCGTCGCGCAATAAATGATCATTTCCACGACGACGCGTAAAGCCAATGCGGTCaa                                                                                                              >  1:2917324/1‑149 (MQ=255)
                                        aaataTGTGCTTTTTTTAAATTCATAAGAATTTTGTATTAATTTATTTTAAAACTTAATTAATCATTTCCTTATTTTTCCGGAAATAATAATGCGTCGCGCAATAAAGGATCATTTCCACGACGACGCGTAAAGCCAATGCGGTCaaaa                                                                                                            <  2:1944055/145‑1 (MQ=255)
                                        aaCAATGTGCTTTTTTTAAATTCATAAGAATTTTGTATTAATTTATTTTAAAACTTAATTAATCATTTCCTTATTTTTCCGGAAATAATAATGCGTCGCGCAATAAATGATCATTTCCACGACGACGCGTAAAGCCAATGCGGTCaaaa                                                                                                            >  1:3888196/1‑149 (MQ=255)
                                        aaCAATGTGCTTTTTTTAAATTCATAAGAATTTTGTATTAATTTATTTTAAAACTTAATTAATCATTTCCTTATTTTTCCGGAAATAATAATGCGTCGCGCAATAAAGGATCATTTCCACGACGACGCGTAAAGCCCATGCGGTCaaaa                                                                                                            >  1:3102874/1‑149 (MQ=255)
                                         aCAATGTGCTTTTTTTAAATTCATAAGAATTTTGTATTAATTTATTTTAAAACTTAATTAATCATTTCCTTATTTTTCCGGAAATAATAATGCGTCGCGCAATAAATGATCATTTCCACGACGACGCGTAAAGCCAATGCGGTCAAAat                                                                                                           >  2:3385933/1‑149 (MQ=255)
                                         aCAATGTGCTTTTTTTAAATTCATAAGAATTTTGTATTAATTTATTTTAAAACTTAATTAATCATTTCCTTATTTTTCCGGAAATAATAATGCGTCGCGCAATAAATGATCATTTCCACGACGACGCGTAAAGCCAATGCGGTCAAAat                                                                                                           >  1:4398170/1‑149 (MQ=255)
                                         aCAATGTGCATTTTTTAAATTCATAAGAATTTTGTATTAATTTATTTTAAAACTTAATTAATCATTTCCTTATTTTTCCGGAAATAATAATGCGTCGCGCAATAAATGATCATTTCCACGACGACGCGTAAAGCCAATGCGGTCAAAat                                                                                                           <  2:2770514/149‑1 (MQ=255)
                                          cAATGTGCTTTTTTTAAATTCATAAGAATTTTGTATTAATTTATTTTAAAACTTAATTAATCATTTCCTTATTTTTCCGGAAATAATAATGCGTCGCGCAATAAATGATCATTTCCACGACGACGCGTAAAGCCAATGCGGTCAAAata                                                                                                          >  2:170910/1‑149 (MQ=255)
                                          cAATGTGATTTTTTTAAATTCATAAGAATTTTGTATTAATTTATTTTAAAACTTAATTAATCATTTCCTTATTTTTCCGGAAATAATAATGCGTCGCGCAATAAATGATCATTTCCACGACGACGCGTAAAGCCAATGCGGTCAAAata                                                                                                          <  1:3475410/149‑1 (MQ=255)
                                             tgtgCTTTTTTTAAATTCATAAGAATTTTGTATTAATTTATTTTAAAACTTAATTAATCATTTCCTTATTTTTCCGGAAATAATAATGCGTCGCGCAATAAATGATCATTTCCACGACGACGCGTAAAGCCAATGCGGTCAAAATATTc                                                                                                       >  2:2653060/1‑149 (MQ=255)
                                                gCTTTTTTTAAATTCATAAGAATTTTGTATTAATTTATTTTAAAACTTAATTAATCATTTCCTTATTTTTCCGGAAATAATAATGCGTCGCGCAATAAAGGATCATTTCCACGACGACGCGTAAAGCCAATGCGGTCAAAATATTCCAg                                                                                                    >  1:4067813/1‑149 (MQ=255)
                                                   ttttttAAATTCATAAGAATTTTGTATTAATTTATTTTAAAACTTAATTAATCATTTCCTTATTTTTCCGGAAATAATAATGCGTCGCGCAATAAATGATCATTTCCACGACGACGCGTAAAGCCAATGCGGTCAAAATATTCCAGAAt                                                                                                 <  2:4626134/149‑1 (MQ=255)
                                                   ttttttAAATTCATAAGAATTTTGTATTAATTTATTTTAAAACTTAATTAATCATTTCCTTATTTTTCCGGAAATAATAATGCGTCGCGCAATAAATGATCATTTCCACGACGACGCGTAAAGCCAATGCGGTCAAAATATTCCAGAAt                                                                                                 <  2:2324521/149‑1 (MQ=255)
                                                    tttttAAATTCATAAGAATTTTGTATTAATTTATTTTAAAACTTAATTAATCATTTCCTTATTTTTCCGGAAATAATAATGCGTCGCGCAATAAATGATCATTTCCACGACGACGCGTAAAGCCAATGCGGTCAAAATATTCCAGAATc                                                                                                <  1:686288/149‑1 (MQ=255)
                                                     ttttAAATTCATAAGAATTTTGTATTAATTTATTTTAAAACTTAATTAATCATTTCCTTATTTTTCCGGAAATAATAATGCGTCGCGCAATAAATGATCATTTCCACGACGACGCGTAAAGCCAATGCGGTCAAAATATTCCAGAATCt                                                                                               <  2:3503555/149‑1 (MQ=255)
                                                        tAAATTCATAAGAATTTTGTATTAATTTATTTTAAAACTTAATTAATCATTTCCTTATTTTTCCGGAAATAATAATGCGTCGCGCAATAAATGATCATTTCCACGACGACGCGTAAAGCCAATGCGGTCAAAATATTCCAGAATCTGaa                                                                                            <  1:2272766/149‑1 (MQ=255)
                                                           aTTCATAAGAATTTTGTATTAATTTATTTTAAAACTTAATTAATCATTTCCTTATTTTTCCGGAAATAATAATGCGTCGCGCAATAAAGGATCATTTCCACGACGACGCGTAAAGCCAATGCGGTCAAAATATTCCAGACTCTGAATTg                                                                                         >  2:4232253/1‑149 (MQ=255)
                                                              cATAAGAATTTTGTATTAATTTATTTTAAAACTTAATTAATCATTTCCTTATTTTTCCGGAAATAATAATGCGTCGCGCAATAAATGATCATTTCCACGACGACGCGTAAAGCCAATGCGGTCAAAATATTCCAGAATCTGAATTGCCa                                                                                      >  1:1673281/1‑149 (MQ=255)
                                                               aTAAGAATTTTGTATTAATTTATTTTAAAACTTAATTAATCATTTCCTTATTTTTCCGGAAATAATAATGCGTCGCGCAATAAATGATCATTTCCACGACGACGCGTAAAGCCAATGCGGTCAAAATATTCCAGAATCTGAATTGCCAg                                                                                     <  1:1966859/149‑1 (MQ=255)
                                                                tAAGAATTTTGTATTAATTTATTTTAAAACTTAATTAATCATTTCCTTATTTTTCCGGAAATAATAATGCGTCGCGCAATAAATGATCATTTCCACGACGACGCGTAAAGCCAATGCGGTCAAAATATTCCAGAATCTGAATTGCCAGc                                                                                    >  1:3314383/1‑149 (MQ=255)
                                                                 aaGAATTTTGTATTAATTTATTTTAAAACTTAATTAATCATTTCCTTATTTTTCCGGAAATAATAATGCGTCGCGCAATAAATGATCATTTCCACGACGACGCGTAAAGCCAATGCGGTCAAAATATTCCAGAATCTGAATTGCCAGCt                                                                                   >  2:2579982/1‑149 (MQ=255)
                                                                 aaGAATTTTGTATTAATTTATTTTAAAACTTAATTAATCATTTCCTTATTTTTCCGGAAATAATAATGCGTCGCGCAATAAATGATCATTTCCACGACGACGCGTAAAGCCAATGCGGTCAAAATATTCCAGAATCTGAATTGCCAGCt                                                                                   <  1:4128012/149‑1 (MQ=255)
                                                                 aaGAATTTTGTATTAATTTATTTTAAAACTTAATTAATCATTTCCTTATTTTTCCGGAAATAATAATGCGTCGCGCAATAAAGGATCATTTCCACGACGACGCGTAAAGCCAATGCGGTCAAAATATTCCAGAATCTGAATTGCCAGCt                                                                                   <  2:2728820/149‑1 (MQ=255)
                                                                  aGAATTTTGTATTAATTTATTTTAAAACTTAATTAATCATTTCCTTATTTTTCCGGAAATAATAATGCGTCGCGCAATAAAGGATCATTTCCACGACGACGCGTAAAGCCAATGCGGTCAAAATATTCCAGAATCTGAATTGCCAGCtt                                                                                  <  2:4524173/149‑1 (MQ=255)
                                                                   gAATTTTGTATTAATTTATTTTAAAACTTAATTAATCATTTCCTTATTTTTCCGGAAATAATAATGCGTCGCGCAATAAATGATCATTTCCACGACGACGCGTAAAGCCAATGCGGTCAAAATATTCCAGAATCTGAATTGCCAGCttt                                                                                 <  2:4105347/149‑1 (MQ=255)
                                                                   gAATTTTGTATTAATTTATTTTAAAACTTAATTAATCATTTCCTTATTTTTCCGGAAATAATAATGCGTCGCGCAATAAATGATCATTTCCACGACGACGCGTAAAGCCAATGCGGTCAAAATATTCCAGAATCTGAATTGCCAGCttt                                                                                 >  2:361574/1‑149 (MQ=255)
                                                                    aaTTTTGTATTAATTTATTTTAAAACTTAATTAATCATTTCCTTATTTTTCCGGAAATAATAATGCGTCGCGCAATAAATGATCATTTCCACGACGACGCGTAAAGCCAATGCGGTCAAAATATTCCAGAATCTGAATTGCCAGCTTTc                                                                                <  1:3597279/149‑1 (MQ=255)
                                                                     aTTTTGTATTAATTTATTTTAAAACTTAATTAATCATTTCCTTATTTTTCCGGAAATAATAATGCGTCGCGCAATAAATGATCATTTCCACGACGACGCGTAAAGCCAATGCGGTCAAAATATTCCAGAATCTGAATTGCCAGCTTTCg                                                                               >  1:1006345/1‑149 (MQ=255)
                                                                         tGTATTAATTTATTTTAAAACTTAATTAATCATTTCCTTATTTTTCCGGAAATAATAATGCGTCGCGCAATAAATGATCATTTCCACGACGACGCGTAAAGCCAATGCGGTCAAAATATTCCAGAATCTGAATTGCCAGCTTTCGGCCt                                                                           <  2:1536916/149‑1 (MQ=255)
                                                                           tATTAATTTATTTTAAAACTTAATTAATCATTTCCTTATTTTTCCGGAAATAATAATGCGTCGCGCAATAAAGGATCATTTCCACGACGACGCGTAAAGCCAATGCGGTCAAAATATTCCAGAATCTGAATTGCCAGCTTTCGGCCTAc                                                                         <  2:220321/149‑1 (MQ=255)
                                                                              tAATTTATTTTAAAACTTAATTAATCATTTCCTTATTTTTCCGGAAATAATAATGCGTCGCGCAATAAATGATCATTTCCACGACGACGCGTAAAGCCAATGCGGTCAAAATATTCCAGAATCTGAATTGCCAGCTTTCGGCCTACGcc                                                                      <  2:411702/149‑1 (MQ=255)
                                                                               aaTTTATTTTAAAACTTAATTAATCATTTCCTTATTTTTCCGGAAATAATAATGCGTCGCGCAATAAATGATCATTTCCACGACGACGCGTAAAGCCAATGCGGTCAAAATATTCCAGAATCTGAATTGCCAGCTTTCGGCCTACGCCt                                                                     >  2:1019/1‑149 (MQ=255)
                                                                               aaTTTATTTTAAAACATAATTAATCATTTCCTTATTTTTCCGGAAATAATAATGCGTCGCGCAATAAATGATCATTTCCACGACGACGCGTAAAGCCAATGCGGTCAAAATATTCCAGAATCTGAATTGCCAGCTTTCGGCCTACGCCt                                                                     >  2:2091/1‑149 (MQ=255)
                                                                                atttatttTAAAACTTAATTAATCATTTCCTTATTTTTCCGGAAATAATAATGCGTCGCGCAATAAATGATCATTTCCACGACGACGCGTAAAGCCAATGCGGTCAAAATATTCCAGAATCTGAATTGCCAGCTTTCGGCCTACGCCTa                                                                    <  1:840561/149‑1 (MQ=255)
                                                                                atttatttTAAAACTTAATTAATCATTTCCTTATTTTTCCGGAAATAATAATGCGTCGCGCAATAAATGATCATTTCCACGACGACGCGTAAAGCCAATGCGGTCAAAATATTCCAGAATCTGAATTGCCAGCTTTCGGCCTACGCCTa                                                                    <  2:374830/149‑1 (MQ=255)
                                                                                 tttatttTAAAACTTAATTAATCATTTCCTTATTTTTCCGGAAATAATAATGCGTCGCGCAATAAATGATCATTTCCACGACGACGCGTAAAGCCAATGCGGTCAAAATATTCCAGAATCTGAATTGCCAGCTTTCGGCCTACGCCTaa                                                                   <  2:3050093/149‑1 (MQ=255)
                                                                                           aaCTTAATTAATCATTTCCTTATTTTTCCGGAAATAATAATGCGTCGCGCAATAAATGATCATTTCCACGACGACGCGTAAAGCCAATGCGGTCAAAATATTCCAGAATCTGAATTGCCAGCTTTCGGCCTACGCCTAAGCGATCGCgg                                                         >  1:2007668/1‑149 (MQ=255)
                                                                                           aaCTTAATTAATCATTTCCTTATTTTTCCGGAAATAATAATGCGTCGCGCAATAAATGATCATTTCCACGACGACGCGTAAAGCCAATGCGGTCAAAATATTCCAGAATCTGAATTGCCAGCTTTCGGCCTACGCCTAAGCGATCGCgg                                                         >  2:2971622/1‑149 (MQ=255)
                                                                                           aaCTTAATTAATCATTTCCTTATTTTTCCGGAAATAATAATGCGTCGCGCAATAAATGATCATTTCCACGACGACGCGTAAAGCCAATGCGGTCAAAATATTCCAGAATCTGAATTGCCAGCTTTCGGCCTACGCCTAAGCGATCGCgg                                                         >  2:738102/1‑149 (MQ=255)
                                                                                            acataattaatCATTTCCTTATTTTTCCGGAAATAATAATGCGTCGCGCAATAAATGATCATTTCCACGACGACGCGTAAAGCCAATGCGGTCAAAATATTCCAGAATCTGAATTGCCAGCTTTCGGCCTACGCCTAAGCGATCGCGGa                                                        <  2:205755/146‑1 (MQ=255)
                                                                                            aCTTAATTAATCATTTCCTTATTTTTCCGGAAATAATAATGCGTCGCGCAATAAATGATCATTTCCACGACGACGCGTAAAGCCAATGCGGTCAAAATATTCCAGAATCTGAATTGCCAGCTTTCGGCCTACGCCTAAGCGATCGCGGa                                                        >  1:4850683/1‑149 (MQ=255)
                                                                                                aattaatCATTTCCTTATTTTTCCGGAAATAATAATGCGTCGCGCAATAAATGATCATTTCCACGACGACGCGTAAAGCCAATGCGGTCAAAATATTCCAGAATCTGAATTGCCAGCTTTCGGCCTACGCCTAAGCGATCGCGGAAATc                                                    >  2:3286847/1‑149 (MQ=255)
                                                                                                  ttaatCATTTCCTTATTTTTCCGGAAATAATAATGCGTCGCGCAATAAATGATCATTTCCACGACGACGCGTAAAGCCAATGCGGTCAAAATATTCCAGAATCTGAATTGCCAGCTTTCGGCCTACGCCTAAGCGATCGCGGAAATccg                                                  >  2:4335858/1‑149 (MQ=255)
                                                                                                   taatCATTTCCTTATTTTTCCGGAAATAATAATGCGTCGCGCAATAAATGATCATTTCCACGACGACGCGTAAAGCCAATGCGGTCAAAATATTCCAGAATCTGAATTGCCAGCTTTCGGCCTACGCCTAAGCGATCGCGGAAATccgc                                                 >  2:345422/1‑149 (MQ=255)
                                                                                                    aatCATTTCCTTATTTTTCCGGAAATAATAATGCGTCGCGCAATAAATGATCATTTCCACGACGACGCGTAAAGCCAATGCGGTCAAAATATTCCAGAATCTGAATTGCCAGCTTTCGGCCTACGCCTAAGCGATCGCGGAAATccgcc                                                >  2:508125/1‑149 (MQ=255)
                                                                                                    aatCATTTCCTTATTTTTCCGGAAATAATAATGCGTCGCGCAATAAATGATCATTTCCACGACGACGCGTAAAGCCAATGCGGTCAAAATATTCCAGAATCTGAATTGCCAGCTTTCGGCCTACGCCTAAGCGATCGCGGAAATccgcc                                                >  1:595652/1‑149 (MQ=255)
                                                                                                        aTTTCCTTATTTTTCCGGAAATAATAATGCGTCGCGCAATAAATGATCATTTCCACGACGACGCGTAAAGCCAATGCGGTCAAAATATTCCAGAATCTGAATTGCCAGCTTTCGGCCTACGCCTAAGCGATCGCGGAAATCCGCCgcgc                                            >  1:3422838/1‑149 (MQ=255)
                                                                                                         tttCCTTATTTTTCCGGAAATAATAATGCGTCGCGCAATAAATGATCATTTCCACGACGACGCGTAAAGCCAATGCGGTCAAAATATTCCAGAATCTGAATTGCCAGCTTTCGGCCTACGCCTAAGCGATCGCGGAAATCCGCCGCGCa                                           <  1:3604683/149‑1 (MQ=255)
                                                                                                          ttCCTTATTTTTCCGGAAATAATAATGCGTCGCGCAATAAAGGATCATTTCCACGACGACGCGTAAAGCCAATGCGGTCAAAATATTCCAGAATCTGAATTGCCAGCTTTCGGCCTACGCCTAAGCGATCGCGGAAATCCGCCGCGCAg                                          >  1:4674860/1‑149 (MQ=255)
                                                                                                             cTTATTTTTCCGGAAATAATAATGCGTCGCGCAATAAATGATCATTTCCACGACGACGCGTAAAGCCAATGCGGTCAAAATATTCCAGAATCTGAATTGCCAGCTTTCGGCCTACGCCTAAGCGATCGCGGAAATCCGCCGCGCAGGtt                                       >  1:1849858/1‑149 (MQ=255)
                                                                                                             cTTATTTTTCCGGAAATAATAATGCGTCGCGCAATAAATGATCATTTCCACGACGACGCGTAAAGCCAATGCGGTCAAAATATTCCAGAATCTGAATTGCCAGCTTTCGGCCTACGCCTAAGCGATCGCGGAAATCCGCCGCGCAGGtt                                       >  2:1786645/1‑149 (MQ=255)
                                                                                                              ttATTTTTCCGGAAATAATAATGCGTCGCGCAATAAATGATCATTTCCACGACGACGCGTAAAGCCAATGCGGTCAAAATATTCCAGAATCTGAATTGCCAGCTTTCGGCCTACGCCTAAGCGATCGCGGAAATCCGCCGCGCAGGTTg                                      >  2:2429100/1‑149 (MQ=255)
                                                                                                              ttATTTTTCCGGAAATAATAATGCGTCGCGCAATAAAGGATCATTTCCACGACGACGCGTAAAGCCAATGCGGTCAAAATATTCCAGAATCTGAATTGCCAGCTTTCGGCCTACGCCTAAGCGATCGCGGAAATCCGCCGCGCAGGTTg                                      >  2:4034771/1‑149 (MQ=255)
                                                                                                              ttATTTTTACGGAAATAATAATGCGTCGCGCAATATAGTATCATTTCCACGACGACGCGTAAAGCCAATGCGGTCAAAATATTCCAGAATCTGAATTGCCAGCTTTCGGCCTACGCCTAAGCGATCGCGGAAATCCGCCGCGCAGGTTg                                      <  1:4792852/149‑1 (MQ=255)
                                                                                                                aTTTTTCCGGAAATAATAATGCGTCGCGCAATAAATGATCATTTCCACGACGACGCGTAAAGCCAATGCGGTCAAAATATTCCAGAATCTGAATTGCCAGCTTTCGGCCTACGCCTAAGCGATCGCGGAAATCCGCCGCGCAGGTTGaa                                    >  2:1608743/1‑149 (MQ=255)
                                                                                                                  ttttCCGGAAATAATAATGCGTCGCGCAATAAATGATCATTTCCACGACGACGCGTAAAGCCAATGCGGTCAAAATATTCCAGAATCTGAATTGCCAGCTTTCGGCCTACGCCTAAGCGATCGCGGAAATCCGCCGCGCAGGTTGAAcc                                  <  2:3314383/149‑1 (MQ=255)
                                                                                                                    ttCCGGAAATAATAATGCGTCGCGCAATAAAGGATCATTTCCACGACGACGCGTAAAGCCAATGCGGTCAAAATATTCCAGAATCTGAATTGCCAGCTTTCGGCCTACGCCTAAGCGATCGCGGAAATCCGCCGCGCAGGTTGAACcac                                >  2:3327895/1‑149 (MQ=255)
                                                                                                                      ccGGAAATAATAATGCGTCGCGCAATAAATGATCATTTCCACGACGACGCGTAAAGCCAATGCGGTCAAAATATTCCAGAATCTGAATTGCCAGCTTTCGGCCTACGCCTAAGCGATCGCGGAAATCCGCCGCGCAGGTTGAACCacac                              <  1:4184552/149‑1 (MQ=255)
                                                                                                                      ccGGAAATAATAATGCGTCGCGCAATAAATGATCATTTCCACGACGACGCGTAAAGCCAATGCGGTCAAAATATTCCAGAATCTGAATTGCCAGCTTTCGGCCTACGCCTAAGCGATCGCGGAAATCCGCCGCGCAGGTTGAACCacac                              >  1:3705041/1‑149 (MQ=255)
                                                                                                                        ggAAATAATAATGCGTCGCGCAATAAATGATCATTTCCACGACGACGCGTAAAGCCAATGCGGTCAAAATATTCCAGAATCTGAATTGCCAGCTTTCGGCCTACGCCTAAGCGATCGCGGAAATCCGCCGCGCAGGTTGAACCACACTc                            >  1:3811503/1‑149 (MQ=255)
                                                                                                                         gAAATAATAATGCGTCGCGCAATAAATGATCATTTCCACGACGACGCGTAAAGCCAATGCGGTCAAAATATTCCAGAATCTGAATTGCCAGCTTTCGGCCTACGCCTAAGCGATCGCGGAAATCCGCCGCGCAGGTTGAACCACACTcc                           <  1:2700872/149‑1 (MQ=255)
                                                                                                                         gAAATAATAATGCGTCGCGCAATAAATGATCATTTCCACGACGACGCGTAAAGCCAATGCGGTCAAAATATTCCAGAATCTGAATTGCCAGCTTTCGGCCTACGCCTAAGCGATCGCGGAAATCCGCCGCGCAGGTTGAACCACACTcc                           >  1:1214397/1‑149 (MQ=255)
                                                                                                                          aaaTAATAATGCGTCGCGCAATAAATGATCATTTCCACGACGACGCGTAAAGCCAATGCGGTCAAAATATTCCAGAATCTGAATTGCCAGCTTTCGGCCTACGCCTAAGCGATCGCGGAAATCCGCCGCGCAGGTTGAACCACACTCCt                          >  1:3043746/1‑149 (MQ=255)
                                                                                                                           aataataatGCGTCGCGCAATAAAGGATCATTTCCACGACGACGCGTAAAGCCAATGCGGTCAAAATATTCCAGAATCTGAATTGCCAGCTTTCGGCCTACGCCTAAGCGATCGCGGAAATCCGCCGCGCAGGTTGAACCACACTCCTg                         <  2:4147448/149‑1 (MQ=255)
                                                                                                                             taataatGCGTCGCGCAATAAATGATCATTTCCACGACGACGCGTAAAGCCAATGCGGTCAAAATATTCCAGAATCTGAATTGCCAGCTTTCGGCCTACGCCTAAGCGATCGCGGAAATCCGCCGCGCAGGTTGAACCACACTCCTGAt                       <  2:480046/149‑1 (MQ=255)
                                                                                                                             taataatGCGTCGCGCAATAAATGATCATTTCCACGACGACGCGTAAAGCCAATGCGGTCAAAATATTCCAGAATCTGAATTGCCAGCTTTCGGCCTACGCCTAAGCGATCGCGGAAATCCGCCGCGCAGGTTGAACCACACTCCTGAt                       >  1:2378018/1‑149 (MQ=255)
                                                                                                                             taataatGCGTCGCGCAATAAATGATCATTTCCACGACGACGCGTAAAGCCAATGCGGTCAAAATATTCCAGAATCTGAATTGCCAGCTTTCGGCCTACGCCTAAGCGATCGCGGAAATCCGCCGCGCAGGTTGAACCACACTCCTGAt                       >  1:4417910/1‑149 (MQ=255)
                                                                                                                              aataatGCGTCGCGCAATAAATGATCATTTCCACGACGACGCGTAAAGCCAATGCGGTCAAAATATTCCAGAATCTGAATTGCCAGCTTTCGGCCTACGCCTAAGCGATCGCGGAAATCCGCCGCGCAGGTTGAACCACACTCCTGATc                      <  2:4431563/149‑1 (MQ=255)
                                                                                                                               ataatGCGTCGCGCAATAAATGATCATTTCCACGACGACGCGTAAAGCCAATGCGGTCAAAATATTCCAGAATCTGAATTGCCAGCTTTCGGCCTACGCCTAAGCGATCGCGGAAATCCGCCGCGCAGGTTGAACCACACTCCTGATCg                     >  1:2092127/1‑149 (MQ=255)
                                                                                                                                taatGCGTCGCGCAATAAATGATCATTTCCACGACGACGCGTAAAGCCAATGCGGTCAAAATATTCCAGAATCTGAATTGCCAGCTTTCGGCCTACGCCTAAGCGATCGCGGAAATCCGCCGCGCAGGTTGAACCACACTCCTGATCga                    <  2:1183108/149‑1 (MQ=255)
                                                                                                                                 aatGCGTCGCGCAATAAATGATCATTTCCACGACGACGCGTAAAGCCAATGCGGTCAAAATATTCCAGAATCTGAATTGCCAGCTTTCGGCCTACGCCTAAGCGATCGCGGAAATCCGCCGCGCAGGTTGAACCACACTCCTGATCgag                   >  1:2188264/1‑149 (MQ=255)
                                                                                                                                 aatGCGTCGCGCAATAAAGGATCATTTCCACGACGACGCGTAAAGCCAATGCGGTCAAAATATTCCAGAATCTGAATTGCCAGCTTTCGGCCTACGCCTAAGCGATCGCGGAAATCCGCCGCGCAGGTTGAACCACACTCCTGATCgag                   <  1:4816715/149‑1 (MQ=255)
                                                                                                                                  atGCGTCGCGCAATAAATGATCATTTCCACGACGACGCGTAAAGCCAATGCGGTCAAAATATTCCAGAATCTGAATTGCCAGCTTTCGGCCTACGCCTAAGCGATCGCGGAAATCCGCCGCGCAGGTTGAACCACACTCCTGATCgaga                  >  1:5169443/1‑149 (MQ=255)
                                                                                                                                  atGCGTCGCGCAATAAATGATCATTTCCACGACGACGCGTAAAGCCAATGCGGTCAAAATATTCCAGAATCTGAATTGCCAGCTTTCGGCCTACGCCTAAGCGATCGCGGAAATCCGCCGCGCAGGTTGAACCACACTCCTGATCgaga                  >  1:2518219/1‑149 (MQ=255)
                                                                                                                                   tGCGTCGCGCAATAAATGATCATTTCCACGACGACGCGTAAAGCCAATGCGGTCAAAATATTCCAGAATCTGAATTGCCAGCTTTCGGCCTACGCCTAAGCGATCGCGGAAATCCGCCGCGCAGGTTGAACCACACTCCTGATCGAGAt                 <  2:5187951/149‑1 (MQ=255)
                                                                                                                                   tGCGTCGCGCAATAAAGGATCATTTCCACGACGACGCGTAAAGCCAATGCGGTCAAAATATTCCAGAATCTGAATTGCCAGCTTTCGGCCTACGCCTAAGCGATCGCGGAAATCCGCCGCGCAGGTTGAACCACACTCCTGATCGAGAt                 <  1:4589693/149‑1 (MQ=255)
                                                                                                                                     cGTCGCGCAATAAATGATCATTTCCACGACGACGCGTAAAGCCAATGCGGTCAAAATATTCCAGAATCTGAATTGCCAGCTTTCGGCCTACGCCTAAGCGATCGCGGAAATCCGCCGCGCAGGTTGAACCACACTCCTGATCGAGATcg               <  2:2466275/149‑1 (MQ=255)
                                                                                                                                      gTCGCGCAATAAATGATCATTTCCACGACGACGCGTAAAGCCAATGCGGTCAAAATATTCCAGAATCTGAATTGCCAGCTTTCGGCCTACGCCTAAGCGATCGCGGAAATCCGCCGCGCAGGTTGAACCACACTCCTGATCGAGATcgc              <  1:177738/149‑1 (MQ=255)
                                                                                                                                      gTCGCGCAATAAATGATCATTTCCACGACGACGCGTAAAGCCAATGCGGTCAAAATATTCCAGAATCTGAATTGCCAGCTTTCGGCCTACGCCTAAGCGATCGCGGAAATCCGCCGCGCAGGTTGAACCACACTCCTGATCGAGATcgc              >  1:548827/1‑149 (MQ=255)
                                                                                                                                      gTCGAGAAATAAATGATCCTTTACACGACGACGCGTAAAGCCAATGCGGTCAAAATATTCCAGAATCTGAATTGCCAGCTTTCGGCCTACGCATAATCGATCGCGGAAATCCGCCGCGCAGGTTGAACAACACTACTTATCGAGATcgc              <  2:1370719/149‑1 (MQ=255)
                                                                                                                                           gcAATAAATGATCATTTCCACGACGACGCGTAAAGCCAATGCGGTCAAAATATTCCAGAATCTGAATTGCCAGCTTTCGGCCTACGCCTAAGCGATCGCGGAAATCCGCCGCGCAGGTTGAACCACACTCCTGATCGAGATCGCGGATc         <  2:2370260/149‑1 (MQ=255)
                                                                                                                                              aTAAATGATCATTTCCACGACGACGCGTAAAGCCAATGCGGTCAAAATATTCCAGAATCTGAATTGCCAGCTTTCGGCCTACGCCTAAGCGATCGCGGAAATCCGCCGCGCAGGTTGAACCACACTCCTGATCGAGATCGCGGATCata      >  1:115381/1‑149 (MQ=255)
                                                                                                                                               tAAATGATCATTTCCACGACGACGCGTAAAGCCAATGCGGTCAAAATATTCCAGAATCTGAATTGCCAGCTTTCGGCCTACGCCTAAGCGATCGCGGAAATCCGCCGCGCAGGTTGAACCACACTCCTGATCGAGATCGCGGATCatat     <  1:2604344/149‑1 (MQ=255)
                                                                                                                                                aaaTGATCATTTCCACGACGACGCGTAAAGCCAATGCGGTCAAAATATTCCAGAATCTGAATTGCCAGCTTTCGGCCTACGCCTAAGCGATCGCGGAAATCCGCCGCGCAGGTTGAACCACACTCCTGATCGAGATCGCGGATCATAtt    >  2:88203/1‑149 (MQ=255)
                                                                                                                                                  aTGATCATTTCCACGACGACGCGTAAAGCCAATGCGGTCAAAATATTCCAGAATCTGAATTGCCAGCTTTCGGCCTACGCCTAAGCGATCGCGGAAATCCGCCGCGCAGGTTGAACCACACTCCTGATCGAGATCGCGGATCATATTgg  <  1:163750/149‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                   
TATCAAAATTAATCTAATAAAAAGCACATTGTATTCATTAAACAATGTGCTTTTTTTAAATTCATAAGAATTTTGTATTAATTTATTTTAAAACTTAATTAATCATTTCCTTATTTTTCCGGAAATAATAATGCGTCGCGCAATAAATGATCATTTCCACGACGACGCGTAAAGCCAATGCGGTCAAAATATTCCAGAATCTGAATTGCCAGCTTTCGGCCTACGCCTAAGCGATCGCGGAAATCCGCCGCGCAGGTTGAACCACACTCCTGATCGAGATCGCGGATCATATTGG  >  CP009273/3751267‑3751561

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: